Доступно с лицензией Spatial Analyst.
Доступно с лицензией 3D Analyst.
Краткая информация
Интерполирует поверхность растра на основании значений в точках с использованием метода обратно взвешенных расстояний (ОВР).
Использование
Выходное значение для ячейки, вычисляемое с использованием метода ОВР, ограничено диапазоном применяемых для интерполяции значений. Поскольку метод ОВР находит среднее значение, основываясь на взвешенном расстоянии до опорных точек, среднее не может превышать самое высокое входное значение и не может быть меньше самого низкого значения. Следовательно, этот метод не может создавать хребты или долины, если эти орографические элементы не были определены с помощью опорных точек (Watson and Philip, 1985).
Лучшие результаты при использовании метода ОВР могут быть получены в том случае, если сеть опорных точек достаточно плотная для того, чтобы отразить локальную вариацию, которую вы хотите смоделировать. Если сеть входных опорных точек редкая или нерегулярная, результирующая поверхность будет недостаточно адекватно представлять требуемую поверхность (Watson and Philip, 1985).
Влияние входной точки на интерполируемое значение изотропно. Поскольку влияние входной точки на интерполируемое значение зависит от расстояния, метод ОВР не сохраняет хребты (Philip and Watson, 1982).
Некоторые входные данные могут содержать несколько точек с одинаковыми x,y координатами. Если значения точек в одних и тех же местоположениях совпадают, такие точки считаются идентичными, и предполагается, что они не влияют на выходные данные. Если значения различаются, такие точки рассматриваются как совпадающие.
Различные инструменты интерполяции могут использовать это условие для данных по-разному. Например, в некоторых случаях для вычисления используется первая встреченная совпадающая точка; в других случаях используется последняя точка. Такая ситуация может приводить к тому, что в некоторых точках выходного растра значения будут отличаться от ожидаемых. Решением этой проблемы может стать предварительная подготовка данных, при которой совпадающие точки будут удалены. Инструмент Собрать события в наборе инструментов Пространственная статистика полезен для идентификации совпадающих точек в ваших данных.
Опция, задающая применение барьеров, используется для определения местоположения линейных объектов, о которых известно, что они нарушают непрерывность поверхности. Эти объекты не имеют z-значений. Типичными примерами барьеров являются скалы, разломы и дамбы. Барьеры ограничивают выбранный поднабор входных опорных точек, используемых для интерполяции выходных z-значений, точками, расположенными по ту же сторону от барьера, что и текущая обрабатываемая ячейка. Разделение барьером определяется анализом линии прямой видимости между каждой парой точек. Это означает, что для исключения из зоны влияния каждой точки соседней точки для пары точек не нужно проводить топологическое разделение. Входные опорные точки, которые лежат непосредственно на линии барьера, будут включены в поднаборы точек, расположенные по обе стороны от барьера.
Пространственные объекты барьеров вводятся в виде полилиний. ОВР использует только координаты x,y линейных объектов, поэтому вводить z-значения по левую и правую сторону от барьера необязательно. Любые введенные z-значения будут проигнорированы.
Использование барьеров значительно увеличивает время, затрачиваемое на обработку.
Инструмент имеет ограничения, так как не может обработать более 45 миллионов входных точек. Если ваш входной класс пространственных объектов содержит более 45 миллионов точек, инструмент может не выдать результат. Вы можете обойти это ограничение путем интерполяции поверхности для области изучения по кускам, с сохранением области перекрытия по краям фрагментов, а затем создании мозаики из результатов для получения единого большого набора растровых данных. Также можно использовать набор данных terrain для хранения и отображения точек и поверхностей, состоящих из миллиардов точек измерения.
Если у вас имеется дополнительный модуль Geostatistical Analyst, с помощью входящей в него версии инструмента ОВР вы можете обрабатывать более крупные наборы данных.
Входные векторные данные должны включать по крайней мере одно корректное поле таблицы атрибутов.
См. раздел Параметры среды анализа и Spatial Analyst для получения дополнительной информации о среде геообработки данного инструмента.
Литература:
Дж. М. Филип (Philip, G. M.) и Д. Ф. Уотсон (D. F. Watson), "A Precise Method for Determining Contoured Surfaces" (Точный метод определения поверхностей с изолиниями). Australian Petroleum Exploration Association Journal (Журнал австралийской ассоциации нефтеразведки) 22:205-212.1982.
Д. Ф. Уотсон (D. F. Watson) и Дж. М. Филип (Philip, G. M.), "A Refinement of Inverse Distance Weighted Interpolation" (Уточнение интерполяции по методу обратно взвешенных расстояний). Geoprocessing 2:315-327.1985.
Синтаксис
Idw (in_point_features, z_field, {cell_size}, {power}, {search_radius}, {in_barrier_polyline_features})| Параметр | Объяснение | Тип данных |
in_point_features | Входные точечные объекты, содержащие z-значения, с помощью которых будет построен растр поверхности. | Feature Layer |
z_field | Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле формы Shape, если входные точечные объекты содержат z-значения. | Field |
cell_size (Дополнительный) | Размер ячейки, который будет использован при создании выходного растра. Будет использоваться значение в параметрах среды, если оно было явно задано; в противном случае, будет использован меньший из размеров ширины или высоты экстента входной пространственной привязки, деленный на 250. | Analysis Cell Size |
power (Дополнительный) | Экспонента, используемая для возведения расстояния в степень. Контролирует значимость точек, находящихся в окрестности интерполируемого значения. Более высокое значение степени приводит к меньшему влиянию удаленных точек. Это может быть любое действительное число больше 0, но наиболее корректные результаты будут получены с использованием значений в диапазоне от 0.5 до 3. Значение по умолчанию равно 2. Значение по умолчанию равно 2. | Double |
search_radius (Дополнительный) | Класс Radius определяет, какие из входных точек будут использоваться для интерполяции значения для каждой ячейки в выходном растре. Есть два типа классов радиуса: RadiusVariable и RadiusFixed. Переменный радиус поиска используется для нахождения определенного количества входных опорных точек для интерполяции. Фиксированный тип радиуса использует заданное фиксированное расстояние, в пределах которого все входные точки будут использоваться для интерполяции. По умолчанию используется тип Переменный.
| Radius |
in_barrier_polyline_features (Дополнительный) | Полилинейные объекты, которые будут использоваться в качестве разрыва или границы при поиске входных опорных точек. | Feature Layer |
Возвращаемое значение
| Название | Объяснение | Тип данных |
| out_raster | Выходной растр интерполированной поверхности. Это всегда растр с плавающей точкой. | Raster |
Пример кода
IDW, пример 1 (окно Python)
В этом примере входными данными является точечный шейп-файл, интерполируется выходная поверхность - растр TIFF.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outIDW = Idw("ozone_pts.shp", "ozone", 2000, 2, RadiusVariable(10, 150000))
outIDW.save("C:/sapyexamples/output/idwout.tif")
IDW, пример 2 (автономный скрипт)
В этом примере входными данными является точечный шейп-файл, интерполируется выходная поверхность - растр Grid.
# Name: IDW_Ex_02.py
# Description: Interpolate a series of point features onto a rectangular
# raster using Inverse Distance Weighting (IDW).
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
cellSize = 2000.0
power = 2
searchRadius = RadiusVariable(10, 150000)
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute IDW
outIDW = Idw(inPointFeatures, zField, cellSize, power, searchRadius)
# Save the output
outIDW.save("C:/sapyexamples/output/idwout02")
Параметры среды
Информация о лицензиях
- ArcGIS Desktop Basic: Требует Spatial Analyst или 3D Analyst
- ArcGIS Desktop Standard: Требует Spatial Analyst или 3D Analyst
- ArcGIS Desktop Advanced: Требует Spatial Analyst или 3D Analyst