需要 Geostatistical Analyst 许可。
使用方法
此工具的主要用途是将地统计图层导出为多种输出类型。例如,只需运行一次此工具,即可将通过经验贝叶斯克里金法生成的地统计图层导出为预测栅格、标准误差栅格、分位数栅格及概率栅格。由此可快速导出各种输出类型。
只有使用此工具才能通过 EBK 回归预测工具创建标准误差图、分位数图或概率图。要获得上述输出,请使用 EBK 回归预测创建地统计图层,并将该地统计图层指定为此工具的输入。
如果输入地统计图层表示 Indicator、Probability 或 Disjunctive 克里金方法,则只能使用用于创建地统计图层的阈值。对于上述各种方法,半变异函数仅适用于计算某一特定的阈值,而不能用于重新计算任何其他阈值。
如果输入地统计图层表示 3D 插值模型,则必须使用输出高程参数来指定要导出的等值线的高程。像元中的点数(水平)和像元中的点数(垂直)参数不适用于 3D 模型。无论这些参数提供的值如何,该工具在水平和垂直方向均会使用一个点。
语法
arcpy.ga.GALayerToRasters(in_geostat_layer, out_raster, {output_type}, {quantile_probability_value}, {cell_size}, {points_per_block_horz}, {points_per_block_vert}, {additional_rasters}, {out_elevation})
参数 | 说明 | 数据类型 |
in_geostat_layer | 要分析的地统计图层。 | Geostatistical Layer |
out_raster | 要创建的主要输出栅格。其他栅格可以通过附加输出栅格参数进行创建。 | Raster Dataset |
output_type (可选) | 输出栅格的表面类型。 有关详细信息,请参阅插值模型可以生成何种类型的输出表面?
| String |
quantile_probability_value (可选) | 如果输出表面类型设置为分位数,则使用此参数输入请求的分位数。如果输出表面类型设置为概率,则使用此参数输入请求的阈值,然后即可计算超出此阈值的概率。 注:如果输入地统计图层是通过未超出地图选项创建的指示器概率图或标准误差图,则将计算不超出阈值的概率。该功能适用于此工具的所有概率栅格输出。 | Double |
cell_size (可选) | 输出栅格的像元大小。该值将由输出栅格和附加输出栅格参数共享。 | Analysis Cell Size |
points_per_block_horz (可选) | 针对水平方向上的各像元用于分块内插的预测数。 | Long |
points_per_block_vert (可选) | 针对垂直方向上的各像元用于分块内插的预测数。 | Long |
additional_rasters [[out_raster, output_type, quantile_probability_value],...] (可选) | 用于为每一个附加输出栅格提供名称、输出类型及分位数或概率值。有关详细信息,请参阅上述参数的描述。这些附加栅格将与输出栅格保存在同一位置。 | Value Table |
out_elevation (可选) | 对于 3D 插值模型,可以导出任何高程处的栅格。可以使用此参数来指定要导出的高程。如果留空,则将从输入图层继承高程。单位将默认为输入图层的相同单位。 | Linear Unit |
派生输出
名称 | 说明 | 数据类型 |
out_additional_rasters | 适用于附加输出栅格的容器。 | 多值(栅格) |
代码示例
将地统计图层转换为栅格。
import arcpy
arcpy.GALayerToRasters_ga("ebk", "C:/gapyexamples/output/ras_predict",
"PREDICTION", None, "5000", 1, 1, "")
将一个地统计图层转换为多个栅格。
# Name: GALayerToRaster_Example_02.py
# Description: Exports a geostatistical layer to a
# - prediction raster
# - standard error of prediction raster
# - quantile raster
# Requirements: Geostatistical Analyst Extension
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/output"
# Set local variables
inLayer = "C:/gapyexamples/data/ebk.lyr"
outPredict = "ras_predict"
cellSize = 5000
cellptsHor = 1
cellptsVer = 1
# Execute GALayerToRasters
arcpy.GALayerToRasters_ga(inLayer, outPredict, "PREDICTION", None, cellSize,
cellptsHor, cellptsVer,
"ras_se PREDICTION_STANDARD_ERROR #;ras_quantile QUANTILE 0.67")
许可信息
- Basic: 需要 地统计分析
- Standard: 需要 地统计分析
- Advanced: 需要 地统计分析