趋势 (空间分析)

摘要

使用趋势面法将点插值成栅格曲面。

了解有关趋势面法工作原理的详细信息

用法

语法

Trend (in_point_features, z_field, {cell_size}, {order}, {regression_type}, {out_rms_file})
参数说明数据类型
in_point_features

包含要插值到表面栅格中的 z 值的输入点要素。

Feature Layer
z_field

存放每个点的高度值或量级值的字段。

如果输入点要素包含 z 值,则该字段可以是数值型字段或者 Shape 字段。

如果回归类型为 Logistic,则该字段的值只能为 0 或 1。

Field
cell_size
(可选)

要创建的输出栅格的像元大小。

如果明确设置该值,则它将是环境中的值,否则,它是输入空间参考中输入点要素范围的宽度或高度除以 250 之后得到的较小值。

Analysis Cell Size
order
(可选)

多项式的阶。

该值必须是介于 1 和 12 之间的整数。值为 1 会对点进行平面拟合,而较高的值则会拟合出更为复杂的曲面。默认值为 1。

Long
regression_type
(可选)

要执行的回归类型。

  • LINEAR 执行多项式回归,对输入点进行最小二乘曲面拟合。这种类型适用于连续型数据。
  • LOGISTIC 执行逻辑趋势面分析。为二元数据生成连续的概率曲面。
String
out_rms_file
(可选)

包含插值的 RMS 误差和卡方相关信息的输出文本文件的文件名。

扩展名必须为“.txt”。

File

返回值

名称说明数据类型
out_raster

输出插值后的表面栅格。

Raster

代码示例

趋势面法示例 1(Python 窗口)

该示例输入一个点 shapefile,然后通过对表面插值,输出得到 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")
趋势面法示例 2(独立脚本)

该示例输入一个点 shapefile,然后通过对表面插值,输出得到 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")

环境

相关主题

许可信息

ArcView: 需要 Spatial Analyst 或 3D Analyst
ArcEditor: 需要 Spatial Analyst 或 3D Analyst
ArcInfo: 需要 Spatial Analyst 或 3D Analyst

7/10/2012