需要 Geostatistical Analyst 许可。
使用方法
使用此工具可获得表示感兴趣区域表面渐进趋势的平滑表面。
当数据中存在短程变化时，应使用局部多项式插值法工具。
语法
GlobalPolynomialInterpolation(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 |
代码示例
将点要素插值成矩形栅格。
import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.GlobalPolynomialInterpolation_ga("ca_ozone_pts", "OZONE", "outGPI",
"C:/gapyexamples/output/gpiout", "2000", "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
# Execute GlobalPolynomialInterpolation
arcpy.GlobalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer,
outRaster, cellSize, power)许可信息
- Basic: 需要 地统计分析
- Standard: 需要 地统计分析
- Advanced: 需要 地统计分析