Интерполяция по методу глобального полинома (Global Polynomial Interpolation) (Geostatistical Analyst)
Резюме
Подгоняет сглаженную поверхность, определенную математической функцией (полиномом) к входным опорным точкам.
Использование
-
В результате выполнения этого инструмента строится сглаженная поверхность, которая представляет постепенные тренды поверхности в исследуемой области.
-
Инструмент Интерполяция по методу локальных полиномов (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 |
weight_field (дополнительно) |
Используется для усиления влияния наблюдаемого значения. Чем больше вес, тем больше точка влияет на интерполируемое значение. При совпадении точек наблюдения присвойте наибольший вес более достоверному измерению. | Field |
Пример кода
Интерполяция точечных объектов в прямоугольный растр.
import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.GlobalPolynomialInterpolation_ga("ca_ozone_pts", "OZONE", "outGPI",
"C:/gapyexamples/output/gpiout", "2000", "2", "")
Интерполяция точечных объектов в прямоугольный растр.
# Name: GlobalPolynomialInterpolation_Example_02.py
# Description: Global Polynomial interpolation fits a smooth surface that is
# defined by a mathematical function (a polynomial) to the input
# sample points. The Global Polynomial surface changes gradually
# and captures coarse-scale pattern in the data. Global Polynomial
# interpolation is like taking a piece of paper and fitting it
# between the raised points (raised to the height of value).
# 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 = "outGPI"
outRaster = "C:/gapyexamples/output/gpiout"
cellSize = 2000.0
power = 2
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute GlobalPolynomialInterpolation
arcpy.GlobalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer,
outRaster, cellSize, power)