ОВР (IDW) (Geostatistical Analyst)
Краткая информация
Использует измеренные значения вокруг интерполируемого местоположения для прогнозирования значения для любого неисследованного места, исходя из предположения, что близко расположенные объекты похожи больше, чем более удаленные объекты.
Использование
-
Проинтерполированные значения ограничиваются диапазоном значений, используемых для интерполяции. Поскольку метод ОВР находит среднее значение, основываясь на взвешенном расстоянии до опорных точек, среднее не может превышать самое высокое входное значение и не может быть меньше самого низкого значения. Следовательно, этот метод не может создавать хребты или долины, если эти орографические элементы не были определены с помощью опорных точек.
-
Метод ОВР может создать эффект «бычьих глаз» вокруг точек с измеренными данными.
-
Обязательных допущений для входных данных не существует.
-
Этот метод хорошо подходит для очень больших наборов входных данных.
Синтаксис
| Параметр | Объяснение | Тип данных |
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 (дополнительно) |
Расстояние, контролирующее влияние соседних точек на интерполируемое значение. Более высокое значение степени приводит к меньшему влиянию удаленных точек. | Double |
search_neighborhood (дополнительно) |
Определяет, какие точки, находящиеся в окрестности, будут использованы для вычисления результата. По умолчанию – стандартные. Класс функции Поиска соседних точек (Search Neighborhood) SearchNeighborhoodStandard, SearchNeighborhoodSmooth, SearchNeighborhoodStandardCircular и SearchNeighborhoodSmoothCircular. Standard
Сглаживание
Стандартная окружность (StandardCircular)
Круговое сглаживание (SmoothCircular)
| Geostatistical Search Neighborhood |
weight_field (дополнительно) |
Используется для усиления влияния наблюдаемого значения. Чем больше вес, тем больше точка влияет на интерполируемое значение. При совпадении точек наблюдения присвойте наибольший вес более достоверному измерению. | Field |
Пример кода
Интерполировать последовательность точечных объектов в растр.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.IDW_ga("ca_ozone_pts", "OZONE", "outIDW", "C:/gapyexamples/output/idwout", "2000", "2",
arcpy.SearchNeighborhoodStandard(300000, 300000, 0, 15, 10, "ONE_SECTOR"), "")
Интерполировать последовательность точечных объектов в растр.
# Name: InverseDistanceWeighting_Example_02.py
# Description: Interpolate a series of point features onto a rectangular raster
# using Inverse Distance Weighting (IDW).
# 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 = "outIDW"
outRaster = "C:/gapyexamples/output/idwout"
cellSize = 2000.0
power = 2
# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
maxNeighbors = 15
minNeighbors = 10
sectorType = "ONE_SECTOR"
searchNeighbourhood = arcpy.SearchNeighborhoodStandard(majSemiaxis, minSemiaxis,
angle, maxNeighbors,
minNeighbors, sectorType)
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute IDW
arcpy.IDW_ga(inPointFeatures, zField, outLayer, outRaster, cellSize,
power, searchNeighbourhood)