含障碍的核插值法 (地统计分析)

摘要

一个移动窗口预测器,它使用两点之间的最短距离,这样可以将线障碍任意一侧的点都连接起来。

用法

语法

KernelInterpolationWithBarriers_ga (in_features, z_field, out_ga_layer, {out_raster}, {cell_size}, {in_barrier_features}, {kernel_function}, {bandwidth}, {power}, {ridge}, {output_type})
参数说明数据类型
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
in_barrier_features
(可选)

使用“非欧氏”距离而非“通视分析”距离的绝对障碍要素。

Feature Layer
kernel_function
(可选)

模拟中所使用的核函数。

  • QUARTIC 四阶多项式函数。
  • EXPONENTIAL 函数按比例增长或衰减。
  • GAUSSIAN 朝正/负无穷方向快速跌落的钟形函数。
  • EPANECHNIKOV 不连续的抛物线函数。
  • POLYNOMIAL5 五阶多项式函数。
  • CONSTANT指示函数。
核函数
String
bandwidth
(可选)

用于指定预测所用数据点之间的最大距离。随着带宽的增加,预测偏差将增加,而预测方差会减少。

Double
power
(可选)

设置多项式的阶数。

Long
ridge
(可选)

用于线性方程组解的数值稳定性。对于不含障碍的规则分布数据,它不会影响预测。对于数据位于要素障碍附近或被障碍隔离的区域,预测可能不稳定,而往往需要相对较大的岭参数值。

Double
output_type
(可选)

用于存储插值结果的表面类型。

  • PREDICTION预测表面可通过内插值生成。
  • PREDICTION_STANDARD_ERROR 标准误差表面可通过内插值的标准误差生成。
String

代码示例

KernelInterpolationWithBarriers(Python 窗口)

使用障碍要素类将点要素插入矩形栅格。

import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.KernelInterpolationWithBarriers_ga("ca_ozone_pts", "OZONE", "outKIWB",
                                         "C:/gapyexamples/output/kiwbout", "2000",
                                         "ca_outline", "QUARTIC", "", "", "50", "PREDICTION")
KernelInterpolationWithBarriers(独立脚本)

使用障碍要素类将点要素插入矩形栅格。

# Name: KernelInterpolationWithBarriers_Example_02.py
# Description: Kernel Interpolation with Barriers is a moving window predictor
#   that uses non-Euclidean distances.
# 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 = "outKIWB"
outRaster = "C:/gapyexamples/output/kiwbout"
cellSize = 2000.0
inBarrier = "ca_outline.shp"
kernelFunction = "QUARTIC"
bandwidth = ""
power = ""
ridgeParam = "50"
outputType = "PREDICTION"

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

# Execute KernelInterpolationWithBarriers
arcpy.KernelInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer, outRaster,
                                         cellSize, inBarrier, kernelFunction, bandwidth,
                                         power, ridgeParam, outputType)



环境

相关主题


7/10/2012