需要 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 | 也可以是以分号分隔的元素字符串。各元素包含以下组件:
| 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 |
代码示例
选定点位置处的预测值。
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", "", "")
选定点位置处的预测值。
# 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: 需要 地统计分析