含障碍的扩散插值 (地统计分析)

需要 Geostatistical Analyst 许可。

摘要

使用基于热方程的核插值表面,并且允许使用栅格和要素障碍重新定义输入点间的距离。

了解有关“含障碍的扩散插值”工作原理的详细信息

使用情况

  • 绝对要素障碍采用“非欧氏”距离方法而非“通视分析”方法。“通视分析”方法要求测量位置与需要预测的位置间的直线与障碍要素不相交。如果障碍周围的距离位于搜索邻域规范之内,则将对其使用“非欧氏”距离方法。

  • 处理时间取决于障碍要素类几何的复杂性。通过使用制图综合工具集中的工具对要素进行平滑处理或删除某些要素可创建新要素类。

  • 对于输入附加障碍栅格参数,其值必须大于或等于 1。值 1 表示没有障碍。

  • 输入累积障碍栅格参数中值的单位应该与 x,y 坐标单位相同。如果相邻单元的值相同,则表示该位置没有障碍。

  • 输入流动障碍栅格中值的单位应该与 x,y 坐标单位相同。如果相邻像元的值相同,则表示该位置没有障碍。另外,如果值由高变低,则意味着没有障碍。

  • 任何可选栅格障碍输入中存在的 NoData 都等同于包含输入绝对障碍要素。

参数

标注说明数据类型
输入要素

包含要插入的 z 值的输入点要素。

Feature Layer
Z 值字段

表示每个点的高度或量级值的字段。如果输入要素包含 z 值或 m 值,则该字段可以是数值字段或 Shape 字段。

Field
输出地统计图层
(可选)

生成的地统计图层。只有未请求任何输出栅格时才需要输出该图层。

Geostatistical Layer
输出栅格
(可选)

输出栅格。只有未请求任何输出地统计图层时才需要输出该栅格。

Raster Dataset
输出像元大小
(可选)

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

可以通过像元大小参数在环境中明确设置该值。

如果未设置,则该值为输入空间参考中输入点要素范围的宽度与高度中的较小值除以 250。

Analysis Cell Size
输入绝对障碍要素
(可选)

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

Feature Layer
带宽
(可选)

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

Double
迭代次数
(可选)

由于模型在数值上求解扩散方程,迭代计数将控制数值解的精度。迭代计数越大,预测越精确,但处理时间也将越长。障碍几何越复杂且带宽越大,精确预测所需的迭代也越多。

Long
权重字段
(可选)

用于强调某个观测。权重越大,对预测的影响就越大。对于重合的观测,为最可靠的测量值分配最大权重。

Field
输入附加障碍栅格
(可选)

基于以下公式计算出的两个相邻栅格像元之间的行程距离:

(相邻像元的平均成本值)x(像元中心间的距离)

Raster Layer
输入累积障碍栅格
(可选)

基于以下公式计算出的两个相邻栅格像元之间的行程距离:

(相邻像元的成本值之差)+(像元中心间的距离)

Raster Layer
输入流动障碍栅格
(可选)

如果想在插入数据时指明数据变化的主方向,可基于以下公式使用流动障碍:

指示符( 相邻像元的成本值 > 相邻像元的成本值)*( 相邻像元的成本值 - 相邻像元的成本值)+(像元中心间的距离),

其中,指示符(真)= 1,指示符(假)= 0。

Raster Layer

arcpy.ga.DiffusionInterpolationWithBarriers(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {in_barrier_features}, {bandwidth}, {number_iterations}, {weight_field}, {in_additive_barrier_raster}, {in_cumulative_barrier_raster}, {in_flow_barrier_raster})
名称说明数据类型
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
bandwidth
(可选)

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

Double
number_iterations
(可选)

由于模型在数值上求解扩散方程,迭代计数将控制数值解的精度。迭代计数越大,预测越精确,但处理时间也将越长。障碍几何越复杂且带宽越大,精确预测所需的迭代也越多。

Long
weight_field
(可选)

用于强调某个观测。权重越大,对预测的影响就越大。对于重合的观测,为最可靠的测量值分配最大权重。

Field
in_additive_barrier_raster
(可选)

基于以下公式计算出的两个相邻栅格像元之间的行程距离:

(相邻像元的平均成本值)x(像元中心间的距离)

Raster Layer
in_cumulative_barrier_raster
(可选)

基于以下公式计算出的两个相邻栅格像元之间的行程距离:

(相邻像元的成本值之差)+(像元中心间的距离)

Raster Layer
in_flow_barrier_raster
(可选)

如果想在插入数据时指明数据变化的主方向,可基于以下公式使用流动障碍:

指示符( 相邻像元的成本值 > 相邻像元的成本值)*( 相邻像元的成本值 - 相邻像元的成本值)+(像元中心间的距离),

其中,指示符(真)= 1,指示符(假)= 0。

Raster Layer

代码示例

DiffusionInterpolationWithBarriers 示例 2(Python 窗口)

将由障碍限制的点要素插入到矩形栅格中。

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.DiffusionInterpolationWithBarriers_ga("ca_ozone_pts", "OZONE", "outDIWB",
                                            "C:/gapyexamples/output/diwbout", "2000",
                                            "ca_outline", "", "10", "", "", "", "")
DiffusionInterpolationWithBarriers 示例 2(独立 Python 脚本)

将由障碍限制的点要素插入到矩形栅格中。

# Name: DiffusionInterpolationWithBarriers_Example_02.py
# Description: Diffusion Interpolation with Barriers uses a kernel which is 
#              based upon the heat equation and describes the variation in 
#              temperature with time in a homogeneous medium.
# 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 = "outDIWB"
outRaster = "C:/gapyexamples/output/diwbout"
cellSize = 2000.0
power = 2
inBarrier = "ca_outline.shp"
bandwidth = ""
iterations = 10
weightField = ""
addBarrier = ""
cumuBarrier = ""
flowBarrier = ""

# Execute DiffusionInterpolationWithBarriers
arcpy.DiffusionInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer,
                                            outRaster, cellSize, inBarrier,
                                            bandwidth, iterations, weightField,
                                            addBarrier, cumuBarrier, flowBarrier)

许可信息

  • Basic: 需要 Geostatistical Analyst
  • Standard: 需要 Geostatistical Analyst
  • Advanced: 需要 Geostatistical Analyst

相关主题