移动窗口克里金法 (地统计分析)

需要 Geostatistical Analyst 许可。

描述

基于较小的邻域,经过所有位置点重新计算变程、块金和偏基台半变异函数参数。

了解有关“移动窗口克里金法”工作原理的详细信息

使用方法

  • 地统计模型源是应用克里金模型,而非经验贝叶斯克里金模型地的统计图层或地统计模型 (XML)。

  • 输入数据集必须至少包含 10 个点该工具才能执行。但是,该工具在处理含有非静态趋势的大型数据集时效果最佳。

  • 在 Python 脚本中,GeostatisticalDatasets ArcPy 类将用于填充输入数据集参数。

  • 对于支持空值的数据格式(如文件地理数据库要素类),使用空值表示无法对该位置做出预测或该值用作输入时应被忽略。对于不支持空值的数据格式(如 shapefile),使用值 -1.7976931348623158e+308(C++ 定义的常量 DBL_MAX 的负数)表示无法对该位置做出预测。

语法

GAMovingWindowKriging(in_ga_model_source, in_datasets, in_locations, neighbors_max, out_featureclass, {cell_size}, {out_surface_grid})
参数说明数据类型
in_ga_model_source

要分析的地统计模型源。

File; Geostatistical Layer
in_datasets

GeostatisticalDatasets 对象。

也可以是以分号分隔的元素字符串。各元素包含以下组件:

  • 数据集的目录路径与名称或当前内容表中的图层名称,后跟一个空格。
  • 一系列字段名称,各字段名称以空格分隔。如果是栅格,则使用像元值。
Geostatistical Value Table
in_locations

将执行预测的点位置。

Feature Layer
neighbors_max

要在移动窗口中使用的相邻要素的数目。

Long
out_featureclass

存储结果的要素类。

Feature Class
cell_size
(可选)

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

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

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

Analysis Cell Size
out_surface_grid
(可选)

使用局部多项式插值法将输出要素类中的预测值插到栅格。

Raster Dataset

代码示例

MovingWindowKriging 示例 1(Python 窗口)

选定点位置处的预测值。

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GAMovingWindowKriging_ga(
    "C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
    "C:/gapyexamples/data/obs_pts.shp", "10", "C:/gapyexamples/output/outMWK", "", "")
MovingWindowKriging 示例 2(独立脚本)

选定点位置处的预测值。

# Name: MovingWindowKriging_Example_02.py
# Description: The kriging model is automatically estimated for each neighborhood
#              as the kriging interpolation moves through all the location points.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

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

# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
inPoints = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
obsPoints = "C:/gapyexamples/data/obs_pts.shp"
maxNeighbors = 10
outPoints = "C:/gapyexamples/output/outMWK"

# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
                               outPoints)

许可信息

  • Basic: 需要 地统计分析
  • Standard: 需要 地统计分析
  • Advanced: 需要 地统计分析

相关主题