全局多项式插值法 (地统计分析)
用法
-
使用此工具可获得表示感兴趣区域表面渐进趋势的平滑表面。
-
而在数据中存在短程变化时,应使用局部多项式插值法工具。
语法
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)
相关主题
7/10/2012