全局多项式插值法 (地统计分析)

摘要

将使用数学函数（多项式）定义的平滑表面与输入采样点拟合。

全局多项式插值法的工作原理

用法

• 使用此工具可获得表示感兴趣区域表面渐进趋势的平滑表面。

• 当数据中存在短程变化时，应使用局部多项式插值法工具。

语法

GlobalPolynomialInterpolation_ga (in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {power}, {weight_field})

参数	说明	数据类型
in_features	包含要插入的 z 值的输入点要素。	Feature Layer
z_field	表示每个点的高度或量级值的字段。如果输入要素包含 z 值或 m 值，则该字段可以是数值字段或 Shape 字段。	Field
out_ga_layer (可选)	生成的地统计图层。只有未请求任何输出栅格时才需要输出该图层。	Geostatistical Layer
out_raster (可选)	输出栅格。只有未请求任何输出地统计图层时才需要输出该栅格。	Raster Dataset
cell_size (可选)	要创建的输出栅格的像元大小。 可在“环境设置”的“栅格分析”下显式设置此值。如果未设置，则该值为输入空间参考中输入点要素范围的宽度与高度中的较小值除以 250。	Analysis Cell Size
power (可选)	多项式的阶。	Long
weight_field (可选)	用于强调某个观测。权重越大，对预测的影响就越大。对于重合的观测，为最可靠的测量值分配最大权重。	Field

代码实例

GlobalPolynomialInterpolation 示例 1（Python 窗口）

将点要素插值成矩形栅格。

import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.GlobalPolynomialInterpolation_ga("ca_ozone_pts", "OZONE", "outGPI",
                                       "C:/gapyexamples/output/gpiout", "2000", "2", "")

GlobalPolynomialInterpolation 示例 2（独立脚本）

将点要素插值成矩形栅格。

# Name: GlobalPolynomialInterpolation_Example_02.py
# Description: Global Polynomial interpolation fits a smooth surface that is
#              defined by a mathematical function (a polynomial) to the input 
#              sample points. The Global Polynomial surface changes gradually 
#              and captures coarse-scale pattern in the data.  Global Polynomial
#              interpolation is like taking a piece of paper and fitting it 
#              between the raised points (raised to the height of value).
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outLayer = "outGPI"
outRaster = "C:/gapyexamples/output/gpiout"
cellSize = 2000.0
power = 2

# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")

# Execute GlobalPolynomialInterpolation
arcpy.GlobalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer, 
                                       outRaster, cellSize, power)

环境

像元大小, 当前工作空间, 掩膜, 输出坐标系, 范围, 临时工作空间, 捕捉栅格, 重合点

许可信息

ArcGIS for Desktop Basic:需要 Geostatistical Analyst

ArcGIS for Desktop Standard:需要 Geostatistical Analyst

ArcGIS for Desktop Advanced:需要 Geostatistical Analyst