Тренд (Trend) (Spatial Analyst)
Резюме
Интерполирует поверхность растра на основании значений в точках с использованием метода тренда.
Использование
-
По мере увеличения порядка (степени) полинома подбираемая поверхность становится все более сложной. Полином высокой степени не всегда позволит получить наиболее точную поверхность; точность поверхности будет зависеть от данных.
-
Дополнительный выходной файл RMS содержит информацию о среднеквадратической ошибке интерполяции. Эта информация может быть использована для определения лучшего значения для степени полинома, через изменение значения степени до тех пор, пока вы не получите наименьшее значение среднеквадратической ошибки. Более подробно о том, как работает Тренд (Trend), смотрите в файле RMS.
-
Для использования опции LOGISTIC Типа регрессии, поле z-значения входных точечных объектов должно иметь коды ноль (0) и единица (1).
-
Некоторые входные данные могут иметь несколько точек с одинаковыми x,y координатами. Если значения точек в одних и тех же местоположениях совпадают, такие точки считаются идентичными и предполагается, что они не влияют на выходные данные. Если значения различаются, такие точки рассматриваются как совпадающие.
Различные инструменты интерполяции могут использовать это условие для данных по-разному. Например, в некоторых случаях для вычисления используется первая совпадающая точка, которую вы встретили; в других случаях используется последняя точка. Такая ситуация может приводить к тому, что в некоторых точках выходного растра значения будут отличаться от ожидаемых. Решением этой проблемы может стать предварительная подготовка данных, при которой совпадающие точки будут удалены. Инструмент Собрать события (Collect Events) в наборе инструментов Пространственная статистика (Spatial Statistics) полезен для идентификации совпадающих точек в ваших данных.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_point_features |
Входные точечные объекты, содержащие z-значения, с помощью которых будет построен растр поверхности. | Feature Layer |
z_field |
Поле, в котором хранится значение высоты или величины для каждой точки. Это может быть числовое поле или поле формы Shape, если входные точечные объекты содержат z-значения. Если тип регрессии — Logistic, значения в поле могут быть только 0 или 1. | Field |
cell_size (дополнительно) |
Размер ячейки, который будет использован при создании выходного растра. Будет использоваться значение в параметрах среды, если оно было явно задано; в противном случае, будет использован меньший из размеров ширины или высоты экстента входной пространственной привязки, деленный на 250. | Analysis Cell Size |
order (дополнительно) |
Степень полинома. Это должно быть целое значение от 1 до 12. Значение, равное 1, будет означать, что для точек будет подобрана плоская поверхность, а более высокое значение приведет к построению поверхности более сложной формы. Значение по умолчанию равно 1. | Long |
regression_type (дополнительно) |
Используемый тип уравнения регрессии.
| String |
out_rms_file (дополнительно) |
Имя выходного текстового файла, который содержит информацию о среднеквадратической ошибке (RMS) и хи-квадрат интерполяции. Расширение файла должно быть .txt. | File |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_raster |
Выходной растр интерполированной поверхности. | Raster |
Пример кода
В этом примере вводится точечный шейп-файл и интерполируется выходная поверхность в качестве растра TIFF.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outTrend = Trend("ozone_pts.shp", "ozone", 2000, 2, "LINEAR")
outTrend.save("C:/sapyexamples/output/trendout.tif")
В этом примере вводится точечный шейп-файл и интерполируется выходная поверхность в качестве растра Grid.
# Name: Trend_Ex_02.py
# Description: Interpolate a series of point features
# onto a rectangular raster using a trend technique.
# 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
PolynomialOrder = 2
regressionType = "LINEAR"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute Trend
outTrend = Trend(inPointFeatures, zField, cellSize,
PolynomialOrder, regressionType)
# Save the output
outTrend.save("C:/sapyexamples/output/trendout02")