Интерполяция по методу локальных полиномов (Local Polynomial Interpolation) (Geostatistical Analyst)
Резюме
Подбирает полином указанного порядка (нулевого, первого, второго, третьего и т. д.), каждый в пределах определенных перекрывающихся окрестностей, чтобы получить результирующую поверхность.
Использование
-
Следует использовать Интерполяцию по методу локальных полиномов (Local Polynomial Interpolation), если в наборе данных есть слабая вариация.
-
Интерполяция по методу глобального полинома (Global Polynomial Interpolation) полезна для создания сглаженных поверхностей и выявления трендов в данных. Однако в науках о Земле исследуемая переменная обычно имеет вариацию с малым диапазоном наряду с трендом в данных. Если в наборе данных есть вариация с малым диапазоном, ее можно обнаружить с помощью инструмента Интерполяция по методу локальных полиномов (Local Polynomial Interpolation).
Синтаксис
| Параметр | Объяснение | Тип данных |
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 |
power (дополнительно) |
Степень полинома. | Long |
search_neighborhood (дополнительно) |
Определяет, какие точки, находящиеся в окрестности, будут использованы для вычисления результата. По умолчанию – стандартные. Это функции класса Поиска соседних точек (Search Neighborhood) SearchNeighborhoodStandard, SearchNeighborhoodSmooth), SearchNeighborhoodStandardCircular и SearchNeighborhoodSmoothCircular. Standard
Сглаживание
Стандартная окружность (StandardCircular)
Круговое сглаживание (SmoothCircular)
| Geostatistical Search Neighborhood |
kernel_function (дополнительно) |
Функция ядра, которая используется в моделировании.
| String |
use_condition_number (дополнительно) |
Опция для управления вычислением интерполируемых значений и стандартными ошибками интерполяции в случае нестабильности интерполяции. Эта опция доступна только для полиномиальных функций порядка 1, 2 и 3.
| Boolean |
bandwidth (дополнительно) |
Применяется для задания максимального расстояния, при котором точки данных используются для интерполяции. С увеличением ширины полосы сдвиг проинтерполированных значений увеличивается, а дисперсия интерполяции уменьшается. | Double |
condition_number (дополнительно) |
У каждой инвертируемой квадратной матрицы имеется число обусловленности, которое указывает, насколько неточным может быть решение линейных уравнений с небольшим изменением в коэффициентах матрицы (это может быть вызвано неточными данными). Если число обусловленности является большим, небольшое изменение в коэффициентах матрицы приведет к большому изменению в векторе решения. | Double |
weight_field (дополнительно) |
Используется для усиления влияния наблюдаемого значения. Чем больше вес, тем больше точка влияет на интерполируемое значение. При совпадении точек наблюдения присвойте наибольший вес более достоверному измерению. | Field |
output_type (дополнительно) |
Тип поверхности для хранения результатов интерполяции.
| String |
Пример кода
Интерполяция точечных объектов в прямоугольный растр.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.LocalPolynomialInterpolation_ga(
"ca_ozone_pts", "OZONE", "outLPI", "C:/gapyexamples/output/lpiout", "2000",
"2", arcpy.SearchNeighborhoodSmooth(300000, 300000, 0, 0.5), "QUARTIC",
"", "", "", "", "PREDICTION")
Интерполяция точечных объектов в прямоугольный растр.
# Name: LocalPolynomialInterpolation_Example_02.py
# Description: Local Polynomial interpolation fits many polynomials, each
# within specified overlapping neighborhoods.
# 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 = "outLPI"
outRaster = "C:/gapyexamples/output/lpiout"
cellSize = 2000.0
power = 2
kernelFunction = "QUARTIC"
bandwidth = ""
useConNumber = ""
conNumber = ""
weightField = ""
outSurface = "PREDICTION"
# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
smoothFactor = 0.5
searchNeighbourhood = arcpy.SearchNeighborhoodSmooth(majSemiaxis, minSemiaxis,
angle, smoothFactor)
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute LocalPolynomialInterpolation
arcpy.LocalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer, outRaster,
cellSize, power, searchNeighbourhood,
kernelFunction, bandwidth, useConNumber,
conNumber, weightField, outSurface)