Статистические исследования природы и структуры пространственных данных, моделирование и построение поверхностей при использовании модуля Geostatistical Analyst программного пакета ArcGIS 9.3
Курсовой проект - Компьютеры, программирование
Другие курсовые по предмету Компьютеры, программирование
сти между измеренными величинами) и ковариация для всех пар опорных точек из набора данных; на облаке вариограммы эти значения показаны как функция расстояния между двумя точками. Облако вариограммы/ковариации может быть использовано для изучения локальных характеристик пространственной автокорреляции в наборе данных и определения выпадающих значений в выборке.
Для варианта 12 не удалось исследовать облако вариограммы поскольку количество точек превышает допустимое значение 300 точек.
.3 Картографирование значений глубин
1.На панели инструментов Geostatistical Analyst выберем опцию Geostatistical Wizard... (Мастер операций геостатистики...).
.В окне Layer (Слой) выберем слой варианта. В окне Attribute (Атрибут) выберем depth. В окне Methods (Методы) выберем Kriging (Кригинг). Нажмем Next (Далее).
3.В диалоговом окне Geostatistical Method Selection, в строке Order of trend removal (Порядок вычитаемого тренда) выберем Second (Второй).
Рисунок 11 - Выбор геостатистического метода.
4.В диалоге Geostatistical Method Selection (Выбор геостатистического метода) нажмем Next (Далее). В диалоге Detrending (Вычитание тренда) нажмем Next (Далее).
Рисунок 12 - Вычитание тренда.
5.Наберем новое значение размера лага. Сначала поставим значение побольше, например 10000, а затем уменьшайте его, пока крайние точки на графике не приблизятся к крайним значениям на оси distance (рис. 13).
Рисунок 13 - Моделирование вариограммы/ковариации.
6.В строке Number of lags (Количество лагов) введем 9.
7.Отметим галочкой опцию Show search direction (Показать направление поиска).
.Включим опцию Anisotropy (Анизотропия). Нажмем Next (Далее).
.Щелкнем мышью на графическом изображении, чтобы выбрать точку, для которой будет выполняться интерполяция (см. рис.16). Объекты, расположенные ближе друг к другу, должны быть более похожими - по мере удаления опорных точек от точки с неизвестным значением, они могут оказаться не слишком полезными для вычисления значения в искомой точке.
Рисунок 16 - Поиск соседства.
10.Нажмем Next (Далее) в диалоге Searching Neighborhood (Поиск соседства).
11.Выберем закладку QQPlot (График КК), чтобы просмотреть график квантиль-квантиль.
.Чтобы выделить цветом местоположение определенной точки, выберем строку в таблице, соответствующую интересующей нас точке. Выбранная точка на диаграмме будет выделена зеленым цветом.
картографирование поверхность визуализация рельеф
Рисунок 17- График КК.
13.Нажмем кнопку Save cross validation... (Сохранить результаты перекрестной проверки...) с тем, чтобы использовать таблицу для дальнейшего анализа.
14.Нажмем OK. Щелкнем по названию слоя, чтобы выделить его цветом, затем еще раз щелкнем по нему и изменим название на "Вычтенный тренд".
.Нажмем правую клавишу мыши на созданном слое "Вычтенный тренд" и выберите в открывающемся меню опцию Create Prediction Standart Error Map (Построить карту значений стандартной ошибки интерполяции). Построенная карта представлена на рис. 18.
Рисунок 18 - Карта значений стандартной ошибки интерполяции.
2.4 Сравнение моделей
1.Правой клавишей мыши щелкним на слое "Вычтенный тренд", выберем опцию Compare... (Сравнить).
Рисунок 19 - Сравнение моделей.
2.Закроем диалог Cross Validation Comparison (Сравнение результатов перекрестной проверки).
3.Нажмем правую клавишу мыши на слое Default ("По умолчанию") и выберем Remove (Удалить).
4.Выберем слой "Вычтенный тренд" и переместим его в нижнюю часть таблицы содержания для того, чтобы можно было видеть точки выборки. Итоговый результат отображен на рисунке 20.
Рисунок 20 - Итоговый результат.
2.5 Картографирование вероятности превышения критического значения глубин
1.На панели Geostatistical Analyst выберем Geostatistical Wizard... (Мастер операций геостатистики...).
2.В окне Layer (Слой) выберем слой варианта. В окне Attribute (Атрибут) выберем depth. В окне Method (Метод) выберем Kriging (Кригинг).
3.В диалоге Choose Input Data and Method (Выберем входные данные и метод) нажмем Next (Далее).
4.Выберем Indicator Kriging (Индикаторный кригинг) (рис.21).
Рисунок 21 - Выбор индикаторного кригинга и первичного порогового значения.
5.Установим первичное пороговое значение (Primary Threshold) равным 12.
6.Отметим опцию Exceed (Превышает).
7.В диалоге Geostatistical Method Selection (Выбор геостатистического метода) нажмем Next (Далее).
8.Зададим размер лага равным 5000 и количество лагов, равным 9.
9.Нажмите Next (Далее) в диалоге Semivariogram/Covariance Modeling (Моделирование вариограммы/ковариации)
Рисунок 22 - Окно выбора вариограммы/ковариации.
10.Нажмем Next (Далее) в диалоге Searching Neighborhood (Поиск соседства).
Голубая линия на диаграмме в следующем диалоге показывает пороговое значение 12m. Точкам слева от этой линии присвоен индикатор 0, а точкам справа от линии индикатор 1 (рис. 23).
11.Выделим цветом строку в таблице, со значением индикатора, равным 0. Эта точка будет выделена зеленым цветом на графике, слева от голубой линии, соответствующей пороговому значению.
Рисунок 23 - Интерполяция значений индикатора.
12.Нажмем кнопку Finish (Готово) в диалоге Cross Validation (Перекрестная проверка).
13.Нажмем OK в диалоге Method Summary (Информация о результирующем слое). В виде данных ArcMap верхним слоем отобразится карта вероятности.
На карте показаны предполагаемые значения индикатора, которые могут быть интерпретированы как вероят?/p>