Интерполяция ядра с барьерами (Geostatistical Analyst)
Краткая информация
Предиктор со скользящим окном, который использует самое короткое расстояние между точками, в результате чего связанными оказываются точки с любой из сторон барьеров.
Использование
-
Абсолютный векторный барьер рассматривается с позиции неевклидовой геометрии, а не с позиции кратчайшего расстояния (линии взгляда). При рассмотрении расстояния с позиции линии взгляда подразумевается, что прямая линия, находящаяся между измеряемым положением и положением, для которого выполняется интерполяция, не пересекается с объектом барьера. Если расстояние вокруг барьера находится в пределах характеристик радиуса поиска, оно будет рассматриваться с позиции неевклидовой геометрии.
-
Время обработки зависит от сложности геометрии классов объектов барьеров. Для создания нового класса объектов путем сглаживая или удаления некоторых объектов, можно использовать инструменты из группы Генерализация (Generalization).
Для функций ядра EXPONENTIAL, GAUSSIAN и CONSTANT коэффициент сглаживания применяется так, чтобы ядра имели конечный радиус, равный значению bandwidth.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_features |
Входные точечные объекты, содержащие z-значения для интерполяции. | Feature Layer |
z_field |
Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле Shape, если входные объекты содержат z-значения или m-значения. | Field |
out_ga_layer (дополнительно) |
Создаваемый геостатический слой. Этот слой является обязательными выходными данными, только если не запрошен выходной растр. | Geostatistical Layer |
out_raster (дополнительно) |
Выходной растр. Этот растр является обязательными выходными данными, только если не запрошен выходной геостатистический слой. | Raster Dataset |
cell_size (дополнительно) |
Размер ячейки, который будет использован при создании выходного растра. Это значение можно явно задать на закладке Анализ растра (Raster Analysis) окна Настройки параметров среды (Environment Settings). Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250. | Analysis Cell Size |
in_barrier_features (дополнительно) |
Пространственные объекты абсолютных барьеров с использованием неевклидовых расстояний вместо расстояний по линии прямой видимости. | Feature Layer |
kernel_function (дополнительно) |
Функция ядра, которая используется в моделировании.
| String |
bandwidth (дополнительно) |
Применяется для задания максимального расстояния, при котором точки данных используются для интерполяции. С увеличением ширины полосы сдвиг проинтерполированных значений увеличивается, а дисперсия интерполяции уменьшается. | Double |
power (дополнительно) |
Задает порядок полиномиальной функции. | Long |
ridge (дополнительно) |
Используется для числовой стабилизации решения системы линейных уравнений. Это не влияет на интерполируемые значения в случае регулярно распределенных данных без барьеров. Интерполяция для областей, в которых данные расположены возле барьера, представленного пространственным объектом, или изолированы барьерами, могут быть неустойчивыми и зачастую требуют относительно больших значений параметров хребта. | Double |
output_type (дополнительно) |
Тип поверхности для хранения результатов интерполяции.
| String |
Пример кода
Интерполяция точечных объектов в прямоугольный растр с использованием класса объекта барьеров.
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")
Интерполяция точечных объектов в прямоугольный растр с использованием класса объекта барьеров.
# 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"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute KernelInterpolationWithBarriers
arcpy.KernelInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer, outRaster,
cellSize, inBarrier, kernelFunction, bandwidth,
power, ridgeParam, outputType)