移动窗口克里金法 (地统计分析)
用法
-
地统计模型源是一个地统计图层或地统计模型 (XML)。
-
输入数据集必须至少包含 10 个点该工具才能执行。但是,该工具在处理含有非静态趋势的大型数据集时效果最佳。
在 Python 脚本中,GeostatisticalDatasets ArcPy 类将用于填充输入数据集参数。
-
对于支持空值的数据格式(如文件和个人地理数据库要素类),使用空值表示无法对该位置做出预测或该值用作输入时应被忽略。对于不支持空值的数据格式(如 shapefile),使用值 -1.7976931348623158e+308(C++ 定义的常量 DBL_MAX 的负数)表示无法对该位置做出预测。
语法
GAMovingWindowKriging_ga (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 | 也可以是以分号分隔的元素字符串。各元素包含以下组件:
| 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"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
outPoints)
相关主题
许可信息
ArcGIS for Desktop Basic: 需要 Geostatistical Analyst
ArcGIS for Desktop Standard: 需要 Geostatistical Analyst
ArcGIS for Desktop Advanced: 需要 Geostatistical Analyst
5/10/2014