高斯地统计模拟 (地统计分析)
用法
-
输入地统计图层必须是对数据集执行了简单克里金法后所生成的结果。其他类型的克里金法所产生的地统计图层不适用于此工具。
此外：
- 建议对数据进行“正态得分”变换以确保输入数据符合标准正态分布。
- 应对已聚类的数据进行去聚操作（使用裁剪轮廓法处理像元或面）以使输入直方图能准确地表示采样总体。在实现中将重新生成（均分）此直方图。
-
要生成条件实现，条件数据应与用于构建模拟所基于的简单克里金模型的数据相同。但是，其他数据集可用作实现的条件。
-
此工具所生成的输出可采用以下识别方式：
- 当已选择保存模拟栅格选项时，将使用前缀加 s0 至 sN（其中 N 是实现的数量）来命名模拟栅格。
- 当已选择后处理选项时，将相应地使用前缀加 MIN、MAX、MEAN、STDDEV、QUARTILE1、MEDIAN、QUARTILE3、QUANTILE 或 P_THRSHLD 来命名输出栅格。
- 当已选择对感兴趣区（统计面）进行后处理时，将使用前缀加面要素类名称来命名输出面要素类。
-
使用不同的前缀识别来自不同模拟运行的输出。如果使用相同的前缀，则在创建新结果之前会清除所有以该前缀开头的现有结果。或者，也可以将来自不同模拟运行的输出存储在不同的文件夹中。
-
如果选中，则输出面要素类将含有每个面内模拟值的汇总统计数据。要了解有关这些汇总统计数据的详细信息，请参阅“高斯地统计模拟”的工作原理。
-
表示感兴趣区的面必须完全包含在模拟栅格范围之内。如果面有任何部分被模拟栅格中的 NoData 值覆盖，则面属性表将包含无效结果。在这种情况下，CELL_COUNT 字段将显示面内含有模拟值的像元数量，并且该数量将以负值表示。
-
“种子”值（在“环境”变量中设置 - 请参阅随机数生成器）指定在模拟中使用的随机数的序列。默认情况下，“种子”值被设为 0，从而使每个模拟都将使用一个新的随机数序列。如果将“种子”设为大于 0 的值，将不断重复模拟结果直到“种子”值发生变化。当需要复制模拟研究的结果时，将种子设置为非 0 值将非常有用。
-
如果已选择保存模拟栅格，则只会向 ArcMap 中的内容列表添加前两个。但是，您可以浏览至输出工作空间并添加其余的模拟栅格。
-
处理重合点的“环境”设置（在 Geostatistical Analyst 设置下）不会影响条件或非条件模拟的结果。构建好简单克里金图层（将用作模拟的输入）后，可在“地统计向导”中管理带有重合点的输入数据集。
-
对于条件模拟，将对落入相同像元内部的条件数据集的点进行平均化，并且实现的条件将遵循此平均值。如果输入像元很大，则会有许多点落入各个像元内并被平均化，并且实现的条件将遵循这些数量（相对）较少的平均值。
-
当前软件限制如下：
- 最大栅格大小被限制为 2049^2 像元（即：对于方形栅格为 2,049 行及 2,049 列）。
- 单次运行中可以请求的最大实现数为 4,500。请注意，工作空间中可包含的最大栅格数为 4,999。
- 基于周期半变异函数模型（J-Bessel 和孔洞效应）的模拟可能不准确。
-
没有足够内存执行请求的操作错误可能说明以所请求像元大小生成的输出栅格过大。
-
对于支持空值的数据格式（如文件和个人地理数据库要素类），使用空值表示无法对该位置做出预测或显示值作为输入时被忽略。对于不支持空值的数据格式（如 shapefile），使用值 -1.7976931348623158e+308（C++ 定义的常量 DBL_MAX 的负数）表示无法对该位置做出预测。
语法
| 参数 | 说明 | 数据类型 |
in_geostat_layer |
输入由简单克里金模型生成的地统计图层。 | Geostatistical Layer |
number_of_realizations |
要执行的模拟数量。 | Long |
output_workspace |
存储所有模拟结果。 | Workspace |
output_simulation_prefix |
自动添加到输出数据库名称中的字母数字前缀（包含 1 至 3 个字符）。 | String |
in_conditioning_features （可选） |
用作实现条件的要素。如果留空，将生成非条件实现。 | Feature Layer |
conditioning_field （可选） |
用作实现条件的字段。如果留空，将生成非条件实现。 | Field |
conditioning_measurement_error_field （可选） |
在输入半变异函数模型中为所有输入数据指定恒定的测量误差。如果每个采样位置的测量误差不同，请使用此字段。将应用输入的测量单位。 如果没有测量误差值，请将此字段留空。 | Field |
cell_size （可选） |
要创建的输出栅格的像元大小。 可在“环境设置”的“栅格分析”下显式设置此值。如果未设置，则该值为输入空间参考中输入点要素范围的宽度与高度中的较小值除以 250。 | Analysis Cell Size |
in_bounding_dataset （可选） |
将分析限制在这些要素的边界面的范围之内。如果输入为点要素，则会自动创建凸包。然后将在该面的内部执行实现。 | Feature Layer |
save_simulated_rasters （可选） |
确定是否将模拟栅格保存到磁盘。 | Boolean |
quantile （可选） |
要生成的输出栅格的分位数。 | Double |
threshold （可选） |
要基于每个像元生成的输出栅格的阈值（以超出所设阈值次数的百分比形式）。 | Double |
in_stats_polygons （可选） |
这些面表示要计算汇总统计数据的感兴趣区域。 | Feature Layer |
raster_stat_type [raster_stat_type,...] （可选） |
将逐个像元地对模拟栅格进行后处理，计算每个所选统计类型并在输出栅格中给出结果。
| String |
代码示例
执行非条件模拟。
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GaussianGeostatisticalSimulations_ga("C:/gapyexamples/data/kriging.lyr", "10",
"C:/gapyexamples/output", "ggs", "", "",
"2000", "", "", "", "", "", "MEAN")
执行非条件模拟。
# Name: GaussianGeostatisticalSimulations_Example_02.py
# Description: This tool performs conditional or unconditional geostatistical
# simulation based on a Simple Kriging model.
# 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"
numRealizations = 10
outWorkspace = "C:/gapyexamples/output"
cellSize = 2000
prefix = "ggs"
rasstatType = "MEAN"
conFeatures = ""
conField = ""
boundingData = ""
savesimRasters = ""
quantile = ""
threshold = ""
statsPolygons = ""
errorField = ""
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute GaussianGeostatisticalSimulations
arcpy.GaussianGeostatisticalSimulations_ga(inLayer, numRealizations, outWorkspace, prefix,
conFeatures, conField, cellSize, boundingData,
savesimRasters, quantile, threshold, statsPolygons,
rasstatType, errorField)