移动窗口克里金法 (地统计分析)
用法
-
地统计模型源是一个地统计图层或地统计模型 (XML)。
-
输入数据集的数据点数必须大于 10。但是,建议将其与具有非静态趋势的大型数据集结合使用。
-
对于支持空值的数据格式(如文件和个人地理数据库要素类),使用空值表示无法对该位置做出预测或显示值作为输入时被忽略。对于不支持空值的数据格式(如 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 [in_dataset,...] |
输入数据集显示两项内容:
| 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)
相关主题
7/10/2012