Доступно с лицензией Geostatistical Analyst.
Сводка
Повторно вычисляет параметры вариограммы Диапазон, Самородок и Частичный порог на основании данных, полученных в пределах меньшей окрестности, перемещаясь по всем местоположениям.
Более подробно о работе инструмента Движущееся окно кригинга
Использование
Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML), представляющая модель кригинга, отличающуюся от эмпирического байесовского кригинга.
Для работы инструмента необходимо чтобы набор входных данных содержал более 10 точек. Однако инструмент наиболее эффективен с большими наборами данных с нестационарными трендами.
В скриптах Python класс GeostatisticalDatasets из ArcPy будет удобен для внесения значений параметра Входные наборы данных (Input dataset(s)).
Для форматов данных, которые поддерживают пустые значения (Null), например для классов пространственных объектов файловых баз геоданных, значение Null используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают пустые значения, например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
arcpy.ga.GAMovingWindowKriging(in_ga_model_source, in_datasets, in_locations, neighbors_max, out_featureclass, {cell_size}, {out_surface_grid})| Parameter | Объяснение | Тип данных |
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 |
Пример кода
Прогнозирует значения в выбранных местоположениях точек.
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"
# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
outPoints)Environments
Информация о лицензиях
- Basic: Требуется Geostatistical Analyst
- Standard: Требуется Geostatistical Analyst
- Advanced: Требуется Geostatistical Analyst