Естественная окрестность (Spatial Analyst)
Резюме
Интерполируется поверхность растра на основании значений точек с использованием метода естественной окрестности.
Более подробно об инструменте Естественная окрестность (Natural Neighbor)
Использование
-
Если центр ячейки периметра выходного растра выходит за пределы выпуклой оболочки (определённой входными точками), то этим ячейкам будут присвоены значения NoData. Если входная точка попадает на одну из ячеек из расположенных по периметру, и центр ячейки находится вне выпуклой оболочки, ячейке все равно будет присвоено значение NoData.
-
Некоторые входные данные могут иметь несколько точек с одинаковыми x,y координатами. Если значения точек в одних и тех же местоположениях совпадают, такие точки считаются идентичными и предполагается, что они не влияют на выходные данные. Если значения различаются, такие точки рассматриваются как совпадающие.
Различные инструменты интерполяции могут использовать это условие для данных по-разному. Например, в некоторых случаях для вычисления используется первая совпадающая точка, которую вы встретили; в других случаях используется последняя точка. Такая ситуация может приводить к тому, что в некоторых точках выходного растра значения будут отличаться от ожидаемых. Решением этой проблемы может стать предварительная подготовка данных, при которой совпадающие точки будут удалены. Инструмент Собрать события (Collect Events) в наборе инструментов Пространственная статистика (Spatial Statistics) полезен для идентификации совпадающих точек в ваших данных.
-
Инструмент имеет ограничения, так как не может обработать более 15 миллионов входных точек. Если входной класс объектов содержит очень большое количество точек (более 15 миллионов), инструмент может прекратить обработку и результат не будет получен.
Это ограничение можно обойти путем обработки поверхности и создания мозаики из результатов для получения единого большого набора растровых данных. Убедитесь, что между частями есть наложение. Также можно использовать набор данных Terrain (Terrain dataset) для хранения и визуализации точек и поверхностей, состоящих из миллиардов точек измерения.
Рекоммендуется, чтобы входные данные были в проекционной системе координат, а не в географической.
-
Если имеется модуль Дополнительный модуль ArcGIS 3D Analyst, в качестве альтернативы можно воспользоваться набором данных TIN. Сначала создайте TIN из исходных данных. Затем преобразуйте результирующую сеть TIN в растр с помощью инструмента TIN в растр (TIN To Raster) и параметра Естественные окрестности (Natural Neighbors). Такая методика особенно эффективна при наличии линий перегиба или области данных неправильной формы.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_point_features |
Входные точечные объекты, содержащие z-значения, с помощью которых будет построен растр поверхности. | Feature Layer |
z_field |
Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле формы Shape, если входные точечные объекты содержат z-значения. | Field |
cell_size (дополнительно) |
Размер ячейки, который будет использован при создании выходного растра. Будет использоваться значение в параметрах среды, если оно было явно задано; в противном случае, будет использован меньший из размеров ширины или высоты экстента входной пространственной привязки, деленный на 250. | Analysis Cell Size |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_raster |
Выходной растр интерполированной поверхности. | Raster |
Пример кода
В этом примере вводится точечный шейп-файл и интерполируется выходная поверхность в качестве растра TIFF.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outNaturalNeighbor = NaturalNeighbor("ozone_pts.shp", "ozone", 2000)
outNaturalNeighbor.save("C:/sapyexamples/output/nnout.tif")
В этом примере вводится точечный шейп-файл и интерполируется выходная поверхность в качестве растра Grid.
# Name: NaturalNeighbor_Ex_02.py
# Description: Interpolate a series of point features onto a
# rectangular raster using Natural Neighbor interpolation.
# 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 = 40000
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute NaturalNeighbor
outNatNbr = NaturalNeighbor(inPointFeatures, zField, cellSize)
# Save the output
outNatNbr.save("C:/sapyexamples/output/nnout02")