Скользящее окно кригинга (Geostatistical Analyst)
Краткая информация
Повторно вычисляет параметры вариограммы Диапазон (Range), Наггет (Nugget) и Частичный порог (Partial Sill) на основании данных, полученных в пределах меньшей окрестности, перемещаясь по всем местоположениям.
Использование
-
Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML).
-
Для работы инструмента необходимо чтобы набор входных данных содержал более 10 точек. Однако инструмент наиболее эффективен с большими наборами данных с нестационарными трендами.
В скриптах Python класс GeostatisticalDatasets из ArcPy будет удобен для внесения значений параметра Входные наборы данных (Input dataset(s)).
-
Для форматов данных, которые поддерживают нулевые значения (Null), например для классов пространственных объектов файловых и персональных баз геоданных, нулевое значение (Null) используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают нулевые значения (Null), например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_ga_model_source |
Анализируемый источник геостатистической модели. | File; Geostatistical Layer |
in_datasets | Объект GeostatisticalDatasets. Или, это может быть строка элементов, разделенных точкой с запятой. Каждый элемент состоит из следующих компонентов:
| Geostatistical Value Table |
in_locations |
Точки, для которых будут проинтерполированны значения. | Feature Layer |
neighbors_max |
Число соседей для использования в скользящем окне. | Long |
out_featureclass |
Класс пространственных объектов, хранящий результаты. | Feature Class |
cell_size (дополнительно) |
Размер ячейки, который будет использован при создании выходного растра. Это значение можно явно задать на закладке Анализ растра (Raster Analysis) окна Настройки параметров среды (Environment Settings). Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250. | Analysis Cell Size |
out_surface_grid (дополнительно) |
Проинтерполированные значения в выходном классе пространственных объектов вычисляются в растр на основе метода Интерполяции локальных полиномов. | Raster Dataset |
Пример кода
Прогнозирует значения в выбранных местоположениях точек.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GAMovingWindowKriging_ga(
"C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
"C:/gapyexamples/data/obs_pts.shp", "10", "C:/gapyexamples/output/outMWK", "", "")
Прогнозирует значения в выбранных местоположениях точек.
# Name: MovingWindowKriging_Example_02.py
# Description: The kriging model is automatically estimated for each neighborhood
# as the kriging interpolation moves through all the location points.
# Requirements: Geostatistical Analyst Extension
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"
# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
inPoints = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
obsPoints = "C:/gapyexamples/data/obs_pts.shp"
maxNeighbors = 10
outPoints = "C:/gapyexamples/output/outMWK"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
outPoints)