Интерполяция ядра с барьерами (Geostatistical Analyst)

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

Сводка

Предиктор с движущимся окном, который использует самое короткое расстояние между точками, в результате чего связанными оказываются точки с любой из сторон барьеров.

Как действует интерполяция ядра с барьерами

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

  • Абсолютный векторный барьер рассматривается с позиции неевклидовой геометрии, а не с позиции кратчайшего расстояния (линии взгляда). При рассмотрении расстояния с позиции линии взгляда подразумевается, что прямая линия, находящаяся между измеряемым положением и положением, для которого выполняется интерполяция, не пересекается с объектом барьера. Если расстояние вокруг барьера находится в пределах характеристик радиуса поиска, оно будет рассматриваться с позиции неевклидовой геометрии.

  • Время обработки зависит от сложности геометрии классов объектов барьеров. Для создания нового класса объектов путем сглаживая или удаления некоторых объектов, можно использовать инструменты из группы Генерализация (Generalization).

  • Для функций ядра Экспоненциальная, Гауссова и Константа коэффициент сглаживания применяется так, чтобы ядра имели конечный радиус, равный значению bandwidth.

Синтаксис

arcpy.ga.KernelInterpolationWithBarriers(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {in_barrier_features}, {kernel_function}, {bandwidth}, {power}, {ridge}, {output_type})
ParameterОбъяснениеТип данных
in_features

Входные точечные объекты, содержащие z-значения для интерполяции.

Feature Layer
z_field

Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле Shape, если входные объекты содержат z-значения или m-значения.

Field
out_ga_layer
(Дополнительный)

Создаваемый геостатический слой. Этот слой является обязательными выходными данными, только если не запрошен выходной растр.

Geostatistical Layer
out_raster
(Дополнительный)

Выходной растр. Этот растр является обязательными выходными данными, только если не запрошен выходной геостатистический слой.

Raster Dataset
cell_size
(Дополнительный)

Размер ячейки, который будет использован при создании выходного растра.

Это значение можно явно задать в Параметрах среды через Размер ячейки.

Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250.

Analysis Cell Size
in_barrier_features
(Дополнительный)

Пространственные объекты абсолютных барьеров с использованием неевклидовых расстояний вместо расстояний по линии прямой видимости.

Feature Layer
kernel_function
(Дополнительный)

Функция ядра, которая используется в моделировании.

  • EXPONENTIAL Функция растет или убывает пропорционально.
  • GAUSSIAN Колоколообразная функция, которая быстро падает в направлении плюс/минус бесконечности.
  • QUARTIC Полиномиальная функция четвертого порядка.
  • EPANECHNIKOV Прерывистая параболическая функция.
  • POLYNOMIAL5 Полиномиальная функция пятого порядка.
  • CONSTANTФункция индикатора.
String
bandwidth
(Дополнительный)

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

Double
power
(Дополнительный)

Задает порядок полиномиальной функции.

Long
ridge
(Дополнительный)

Используется для числовой стабилизации решения системы линейных уравнений. Это не влияет на интерполируемые значения в случае регулярно распределенных данных без барьеров. Интерполяция для областей, в которых данные расположены возле барьера, представленного пространственным объектом, или изолированы барьерами, могут быть неустойчивыми и зачастую требуют относительно больших значений параметров хребта.

Double
output_type
(Дополнительный)

Тип поверхности для хранения результатов интерполяции.

Дополнительные сведения о типах выходной поверхности см. в Какие типы выходных поверхностей могут генерировать модели интерполяции?

  • PREDICTIONПоверхности интерполяции создаются из проинтерполированных значений.
  • PREDICTION_STANDARD_ERROR Поверхности типа Стандартная ошибка создаются из стандартных ошибок проинтерполированных значений.
String

Пример кода

KernelInterpolationWithBarriers, пример 1 (окно Python)

Интерполяция точечных объектов в прямоугольный растр с использованием класса объекта барьеров.

import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.KernelInterpolationWithBarriers_ga("ca_ozone_pts", "OZONE", "outKIWB",
                                         "C:/gapyexamples/output/kiwbout", "2000",
                                         "ca_outline", "QUARTIC", "", "", "50", "PREDICTION")
KernelInterpolationWithBarriers, пример 2 (автономный скрипт Python)

Интерполяция точечных объектов в прямоугольный растр с использованием класса объекта барьеров.

# Name: KernelInterpolationWithBarriers_Example_02.py
# Description: Kernel Interpolation with Barriers is a moving window predictor
#   that uses non-Euclidean distances.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outLayer = "outKIWB"
outRaster = "C:/gapyexamples/output/kiwbout"
cellSize = 2000.0
inBarrier = "ca_outline.shp"
kernelFunction = "QUARTIC"
bandwidth = ""
power = ""
ridgeParam = "50"
outputType = "PREDICTION"

# Execute KernelInterpolationWithBarriers
arcpy.KernelInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer, outRaster,
                                         cellSize, inBarrier, kernelFunction, bandwidth,
                                         power, ridgeParam, outputType)

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

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

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