Доступно с лицензией Geostatistical Analyst.
Краткая информация
Повторно вычисляет параметры вариограммы Диапазон (Range), Наггет (Nugget) и Частичный порог (Partial Sill) на основании данных, полученных в пределах меньшей окрестности, перемещаясь по всем местоположениям.
Более подробно о работе инструмента Движущееся окно кригинга
Использование
Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML), представляющая модель кригинга, отличающуюся от эмпирического байесовского кригинга.
Для работы инструмента необходимо чтобы набор входных данных содержал более 10 точек. Однако инструмент наиболее эффективен с большими наборами данных с нестационарными трендами.
В скриптах Python класс GeostatisticalDatasets из ArcPy будет удобен для внесения значений параметра Входные наборы данных (Input dataset(s)).
Для форматов данных, которые поддерживают пустые значения (Null), например для классов пространственных объектов файловых баз геоданных, значение Null используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают пустые значения, например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
GAMovingWindowKriging_ga (in_ga_model_source, in_datasets, in_locations, neighbors_max, out_featureclass, {cell_size}, {out_surface_grid})
Параметр | Объяснение | Тип данных |
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 (дополнительно) | Размер ячейки, который будет использован при создании выходного растра. Это значение можно явно задать в разделе Анализ растра в Параметрах среды. Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250. | Analysis Cell Size |
out_surface_grid (дополнительно) | Проинтерполированные значения в выходном классе пространственных объектов вычисляются в растр на основе метода локальных полиномов. | Raster Dataset |
Пример кода
MovingWindowKriging Пример 1 (окно Python)
Прогнозирует значения в выбранных местоположениях точек.
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", "", "")
MovingWindowKriging Пример 2 (автономный скрипт)
Прогнозирует значения в выбранных местоположениях точек.
# 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)
Параметры среды
- Размер ячейки (Cell size)
- Совпадающие точки (Coincident Points)
- Текущая рабочая область (Current Workspace)
- Экстент (Extent)
- Географические преобразования (Geographic Transformations)
- Маска (Mask)
- Выходная система координат (Output Coordinate System)
- Временная рабочая область (Scratch Workspace)
- Растр привязки (Snap Raster)
Информация о лицензировании
- ArcGIS for Desktop Basic: Требует Geostatistical Analyst
- ArcGIS for Desktop Standard: Требует Geostatistical Analyst
- ArcGIS for Desktop Advanced: Требует Geostatistical Analyst