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

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

Краткая информация

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

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

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

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

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

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

Параметры

ПодписьОписаниеТип данных
Входные объекты

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

Feature Layer
Поле значений Z

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

Field
Выходной геостатистический слой
(Дополнительный)

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

Geostatistical Layer
Выходной растр
(Дополнительный)

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

Raster Dataset
Размер выходной ячейки
(Дополнительный)

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

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

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

Analysis Cell Size
Входные объекты абсолютных барьеров
(Дополнительный)

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

Feature Layer
Функции ядра
(Дополнительный)

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

  • Экспоненциальная Функция растет или убывает пропорционально.
  • Гауссова Колоколообразная функция, которая быстро падает в направлении плюс/минус бесконечности.
  • Биквадратная Полиномиальная функция четвертого порядка.
  • Функция Епанечникова Прерывистая параболическая функция.
  • Полином пятого порядка Полиномиальная функция пятого порядка.
  • КонстантФункция индикатора.
String
Ширина канала
(Дополнительный)

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

Double
Порядок полинома
(Дополнительный)

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

Long
Параметр гребень
(Дополнительный)

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

Double
Тип выходной поверхности
(Дополнительный)

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

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

arcpy.ga.KernelInterpolationWithBarriers(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {in_barrier_features}, {kernel_function}, {bandwidth}, {power}, {ridge}, {output_type})
ИмяОписаниеТип данных
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

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