Скользящее окно кригинга (Moving Window Kriging) (Geostatistical Analyst)
Резюме
Повторно вычисляет параметры вариограммы Диапазон (Range), Самородок (Nugget) и Частичный порог (Partial Sill) на основании данных, полученных в пределах меньшей окрестности, перемещаясь по всем местоположениям.
Принцип работы инструмента Скользящее окно кригинга (Moving Window Kriging)
Использование
-
Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML).
-
Набор входных данных должен содержать более 10 точек. Однако рекомендуется использовать его с большими наборами данных с нестационарными трендами.
-
Для форматов данных, которые поддерживают нулевые значения (Null), например для классов пространственных объектов файловых и персональных баз геоданных, нулевое значение (Null) значение используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают нулевые значения (Null), например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_ga_model_source |
Анализируемый источник геостатистической модели. | File; Geostatistical Layer |
in_datasets [in_dataset,...] |
Входные наборы данных отображаются в двух частях:
| 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)