克里金法 (3D Analyst)

摘要

使用克里金法将点插值成栅格表面。

了解有关克里金法工作原理的详细信息

用法

语法

Kriging_3d (in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})
参数说明数据类型
in_point_features

包含要插值到表面栅格中的 z 值的输入点要素。

Feature Layer
z_field

存放每个点的高度值或量级值的字段。

如果输入点要素包含 z 值,则该字段可以是数值型字段或者 Shape 字段。

Field
out_surface_raster

输出插值后的表面栅格。

Raster Dataset
semiVariogram_props
kriging_model

要使用的半变异函数模型。

对于克里金法、普通克里金法和泛克里金法,存在两种模型。普通克里金法模型具有五种可用的半变异函数。泛克里金法模型具有两种可用的半变异函数。每个半变异函数都具有若干可设置的可选参数。

  • 普通克里金法模型半变异函数:
    • Spherical - 球面半变异函数模型。这是默认设置。
    • Circular - 圆半变异函数模型。
    • Exponential - 指数半变异函数模型。
    • Gaussian - 高斯(或正态分布)半变异函数模型。
    • Linear - 具有基台的线性半变异函数模型。
  • 泛克里金法模型半变异函数:
    • LinearDrift - 采用一次漂移函数的泛克里金法。
    • QuadraticDrift - 采用二次漂移的泛克里金法。
  • 定义半变异函数模型之后,其余参数均是普通克里金法和泛克里金法所共有的参数。这些求解程序为:
    • Lag size - 默认值为输出栅格的像元大小。
    • MajorRange - 表示一种距离,超出此距离即认定为不相关。
    • PartialSill - 块金和基台之间的差值。
    • Nugget - 表示因过于微小而无法检测到的空间尺度下的误差和变差。块金效应被视为在原点处的不连续。

半变异函数的形式为文本字符串:

"{semivariogramType},{lagSize},{majorRange},{partialSill},{nugget}"

例如:

"Circular, 2000, 2.6, 542"
KrigingModel
cell_size
(可选)

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

如果明确设置该值,则它将是环境中的值,否则,它是输入空间参考中输入点要素范围的宽度或高度除以 250 之后得到的较小值。

Analysis Cell Size
search_radius
(可选)

定义要用来对输出栅格中各像元值进行插值的输入点。

指定搜索邻域的方式有两种:VariableFixed

Variable 使用可变的搜索半径来查找用于插值的指定数量的输入采样点。Fixed 使用指定的固定距离,将利用此距离范围内的所有输入点进行插值。Variable 默认值。

这些参数的语法为:

  • Variable, number_of_points, maximum_distance,其中:
    • number_of_points - 指定要用于执行插值的最邻近输入采样点数量的整数值。默认值为 12 个点。
    • maximum_distance - 使用地图单位指定距离,以此限制对最邻近输入采样点的搜索。默认值是范围的对角线长度。
  • Fixed, distance, minimum_number_of_points,其中
    • distance - 指定用作半径的距离,在该半径范围内的输入采样点将用于执行插值。半径值使用地图单位来表示。默认半径是输出栅格像元大小的五倍。
    • minimum_number_of_points - 定义用于插值的最小点数的整数。默认值为 0。

      如果在指定距离内没有找到所需点数,则将增加搜索距离,直至找到指定的最小点数。

      搜索半径需要增加时就会增加,直到 minimum_number_of_points 在该半径范围内,或者半径的范围越过输出栅格的下部(南)和/或上部(北)范围为止。NoData 会分配给不满足以上条件的所有位置。

Radius
out_variance_prediction_raster
(可选)

可选输出栅格中的每个像元都包含该位置的预测半方差值。

Raster Dataset

代码示例

克里金法示例 1(Python 窗口)

该示例输入一个点 shapefile,然后通过对表面插值,输出得到 Grid 栅格。

import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Kriging_3d("ca_ozone_pts.shp", "OZONE", "c:/output/krigout",
                 "Spherical", 2000, "Variable 12")
克里金法示例 2(独立脚本)

该示例输入一个点 shapefile,然后通过对表面插值,输出得到 Grid 栅格。

# Name: Kriging_3d_Ex_02.py
# Description: Interpolates a surface from points using kriging.
# Requirements: 3D Analyst Extension
# Import system modules

import arcpy
from arcpy import env

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

# Set local variables
inFeatures = "ca_ozone_pts.shp"
field = "OZONE"
outRaster = "C:/output/krigoutput02"
cellSize = 2000
outVarRaster = "C:/output/outvariance"
kModel = "CIRCULAR"
kRadius = 20000

# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")

# Execute Kriging
arcpy.Kriging_3d(inFeatures, field, outRaster, kModel, 
                 cellSize, kRadius, outVarRaster)

环境

相关主题

许可信息

ArcView: 需要 3D Analyst 或 Spatial Analyst
ArcEditor: 需要 3D Analyst 或 Spatial Analyst
ArcInfo: 需要 3D Analyst 或 Spatial Analyst

7/10/2012