Доступно с лицензией Spatial Analyst.
Доступно с лицензией 3D Analyst.
Сводка
Интерполирует поверхность растра на основании значений в точках с использованием метода тренда.
Использование
По мере увеличения порядка (степени) полинома подбираемая поверхность становится все более сложной. Полином высокой степени не всегда позволит получить наиболее точную поверхность; точность поверхности будет зависеть от данных.
Для использования опции Логический в Типе регрессии, поля z-значений входных точечных объектов должны иметь коды ноль (0) и единица (1).
Выходной размер ячейки можно задать как числовым значением, так и получить из существующего набора растровых данных. Если размер ячейки не был явно задан как значение параметра, он извлекается из параметра среды Размер ячейки, если такой задан. Если же и в параметр среды значение размера ячейки не указано, но настроен параметр Растр привязки, то будет использоваться значение величины его ячейки. Если же ни один из параметров не задан, то размер ячейки будет вычислен как наименьшее значение длины либо ширины экстента растра, поделенное на 250, при этом экстент берется из параметра среды Выходная система координат.
Если размер ячейки указан в виде числа, то инструмент будет использовать это значение для выходного растра.
Если для определения размера ячейки указан другой растр, в параметре будет показан путь к набору растровых данных вместо значения размера ячейки. Размер ячейки этого растра будет использоваться непосредственно в анализе указанная пространственная привязка набора данных будет такой же, как и выходная система координат. Если же пространственная привязка набора данных отличается от выходной пространственной привязки, на будет перепроецирована на основе выбранного Метода проецирования размера ячейки.
Дополнительный выходной файл RMS содержит информацию о среднеквадратической ошибке интерполяции. Эта информация может быть использована для определения лучшего значения для степени полинома, через изменение значения степени до тех пор, пока вы не получите наименьшее значение среднеквадратической ошибки. Более подробно о том, как работает Тренд, смотрите в файле RMS.
Некоторые входные данные могут содержать несколько точек с одинаковыми x,y координатами. Если значения точек в одних и тех же местоположениях совпадают, такие точки считаются идентичными, и предполагается, что они не влияют на выходные данные. Если значения различаются, такие точки рассматриваются как совпадающие.
Различные инструменты интерполяции могут использовать это условие для данных по-разному. Например, в некоторых случаях для вычисления используется первая встреченная совпадающая точка; в других случаях используется последняя точка. Такая ситуация может приводить к тому, что в некоторых точках выходного растра значения будут отличаться от ожидаемых. Решением этой проблемы может стать предварительная подготовка данных, при которой совпадающие точки будут удалены. Инструмент Собрать события в наборе инструментов Пространственная статистика полезен для идентификации совпадающих точек в ваших данных.
При использовании в качестве входных данных форматов, которые поддерживают значения Null, таких как классы пространственных объектов файловой базы геоданных, значение Null будет игнорироваться.
Синтаксис
arcpy.3d.Trend(in_point_features, z_field, out_raster, {cell_size}, {order}, {regression_type}, {out_rms_file})| Parameter | Объяснение | Тип данных |
in_point_features | Входные точечные объекты, содержащие z-значения, с помощью которых будет построен растр поверхности. | Feature Layer |
z_field | Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле формы Shape, если входные точечные объекты содержат z-значения. Если тип регрессии – Логический, значения в поле могут быть только 0 или 1. | Field |
out_raster | Выходной растр интерполированной поверхности. Это всегда растр с плавающей точкой. | Raster Dataset |
cell_size (Дополнительный) | Размер ячейки выходного растра, который будет создан. Этот параметр можно задать как числовым значением, так и получить из существующего набора растровых данных. Если размер ячейки не был в явном виде задан как значение параметра, будет использоваться значение размера ячейки, указанное в параметрах среды геообработки, если же и оно не указано, будут использоваться дополнительные правила вычисления размера ячейки на основе входных данных. Более подробно см. раздел об использовании. | Analysis Cell Size |
order (Дополнительный) | Степень полинома. Должно быть задано целым числом от 1 до 12. Значение, равное 1, будет означать, что для точек будет подобрана плоская поверхность, а более высокое значение приведет к построению поверхности более сложной формы. Значение по умолчанию равно 1. | Long |
regression_type (Дополнительный) | Используемый тип уравнения регрессии.
| String |
out_rms_file (Дополнительный) | Имя выходного текстового файла, который содержит информацию о среднеквадратической ошибке (RMS) и хи-квадрат интерполяции. Расширение файла должно быть .txt. | File |
Пример кода
В этом примере входными данными является точечный шейп-файл, интерполируется выходная поверхность – растр TIFF.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Trend_3d("ca_ozone_pts.shp", "ozone",
"C:/output/trendout.tif", 2000, 2, "LINEAR")В этом примере входными данными является точечный шейп-файл, интерполируется выходная поверхность – растр Grid.
# Name: Trend_3d_Ex_02.py
# Description: Interpolate a series of point features onto a
# rectangular raster using a trend technique.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/sapyexamples/output/trendout02"
cellSize = 2000.0
PolynomialOrder = 2
regressionType = "LINEAR"
# Execute Trend
arcpy.Trend_3d(inPointFeatures, zField, outRaster, cellSize,
PolynomialOrder, regressionType)Environments
Информация о лицензиях
- Basic: Требуется 3D Analyst или Spatial Analyst
- Standard: Требуется 3D Analyst или Spatial Analyst
- Advanced: Требуется 3D Analyst или Spatial Analyst