Движущееся окно кригинга (Geostatistical Analyst)

Доступно с лицензией Geostatistical Analyst.

Сводка

Повторно вычисляет параметры вариограммы Диапазон, Самородок и Частичный порог на основании данных, полученных в пределах меньшей окрестности, перемещаясь по всем местоположениям.

Более подробно о работе инструмента Движущееся окно кригинга

Использование

  • Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML), представляющая модель кригинга, отличающуюся от эмпирического байесовского кригинга.

  • Для работы инструмента необходимо чтобы набор входных данных содержал более 10 точек. Однако инструмент наиболее эффективен с большими наборами данных с нестационарными трендами.

  • В скриптах Python класс GeostatisticalDatasets из ArcPy будет удобен для внесения значений параметра Входные наборы данных (Input dataset(s)).

  • Для форматов данных, которые поддерживают пустые значения (Null), например для классов пространственных объектов файловых баз геоданных, значение Null используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают пустые значения, например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.

Синтаксис

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

Пример кода

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"

# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
                               outPoints)

Информация о лицензиях

  • Basic: Требуется Geostatistical Analyst
  • Standard: Требуется Geostatistical Analyst
  • Advanced: Требуется Geostatistical Analyst

Связанные разделы