需要 Spatial Analyst 许可。
描述
汇总另一个数据集区域内的栅格数据值并以表的形式显示结果。
插图
使用
区域定义为输入中具有相同值的所有区域。各区域无需相连。栅格和要素都可用于区域输入。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)的像元大小与输入值栅格(Python 中的 in_value_raster)不同,则输出像元大小将为输入的最大值,且输入值栅格将在内部用作捕捉栅格。如果像元大小相同,但是像元未对齐,则输入值栅格将在内部用作捕捉栅格。在执行区域操作之前,上述任一情况都会触发内部重采样。
如果区域输入和值输入均为具有相同像元大小的栅格,并且像元已对齐,则可将其直接应用于工具中,且工具执行期间不会对其进行内部重采样。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)为栅格,则其必须为整型栅格。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)为要素,则将使用输入值栅格(Python 中的 in_value_raster)的像元大小和像元对齐方式在内部将其转换为栅格。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)为要素,对于不与值栅格的任何像元中心重叠的任何区域要素,这些区域将不会转换为内部区域栅格。因此,这些区域将不会在输出中表示。您可以通过为保留要素区域所需细节层次的像元大小环境指定适当的值来管理这些区域并在分析环境中指定这些区域。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)为点要素,则其可能具有值输入栅格的任何特定像元内所含的多个点。对于这样的像元,区域值将通过具有最低 ObjectID 字段的点(例如 OID 或 FID)来确定。
如果输入栅格或要素区域数据(Python 中的 in_zone_data)具有重叠面,则无法针对每个单独的面执行区域分析。因为要素输入已转换为栅格,所以每个位置只能拥有一个值。
另一种方法就是要为每个面区域反复进行区域操作并核对结果。
在指定输入栅格或要素区域数据(Python 中的 in_zone_data)时,默认区域字段将为第一个可用的整型或文本字段。如果不存在其他有效字段,则 ObjectID 字段(如 OID 或 FID)将成为默认字段。
输入值栅格(Python 中的 in_value_raster)可以为整型,也可以为浮点型。但是,当它是浮点型时,将无法计算众数、中值、少数和变异性。
对于众数和少数计算,当存在平局时,区域的输出将基于平局值的最低值。有关详细信息,请参阅分区统计的工作原理。
将根据统计类型的设置在输出表中创建一个字段或一系列字段。值输入为整型时,可对所有统计数据(最小值、最大值、范围、平均值、标准差、总和、变异度、众数、少数和中值)进行计算。如果值输入为浮点型,将不会计算众数、少数、中值和变异度统计数据。
输出表中各项下的各个值的数据类型取决于所执行的分区计算。有关任何统计的特定行为的详细信息,请参阅“分区统计”的工作原理。
输出表中的行数与区域数相同。
默认情况下,此工具会利用多核处理器。可供使用的最大核数为四。
如果希望工具使用较少的核,请使用 parallelProcessingFactor 环境设置。
有关适用于此工具的地理处理环境的详细信息,请参阅分析环境和 Spatial Analyst。
语法
ZonalStatisticsAsTable(in_zone_data, zone_field, in_value_raster, out_table, {ignore_nodata}, {statistics_type})
参数 | 说明 | 数据类型 |
in_zone_data | 定义区域的数据集。 可通过整型栅格或要素图层来定义区域。 | Raster Layer; Feature Layer |
zone_field | 保存定义每个区域的值的字段。 该字段可以是区域数据集的整型字段或字符串型字段。 | Field |
in_value_raster | 含有要计算统计数据的值的栅格。 | Raster Layer |
out_table | 将包含每个区域中值的汇总的输出表。 表的格式由输出位置和路径确定。默认情况下,输出为一张地理数据库表。如果路径不在地理数据库中,则格式将由表达式确定。如果扩展名为 .dbf,则将采用 dBASE 格式。如果未指定扩展名,则输出将为 INFO 表。 | Table |
ignore_nodata (可选) | 指示值输入中的 NoData 值是否会影响其所落入区域的结果。
| Boolean |
statistics_type (可选) | 要计算的统计类型。
| String |
代码示例
ZonalStatisticsAsTable 示例 1(Python 窗口)
本示例对面要素定义的区域内的栅格值进行了汇总,并将结果记录在表中。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outZSaT = ZonalStatisticsAsTable("zones.shp", "Classes", "valueforzone",
"zonalstattblout", "NODATA", "SUM")
ZonalStatisticsAsTable 示例 2(独立脚本)
本示例对面要素定义的区域内的栅格值进行了汇总,并将结果记录在 .dbf 文件中。
# Name: ZonalStatisticsAsTable_Ex_02.py
# Description: Summarizes values of a raster within the zones of
# another dataset and reports the results to a table.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inZoneData = "zones.shp"
zoneField = "Classes"
inValueRaster = "valueforzone"
outTable = "zonalstattblout02.dbf"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ZonalStatisticsAsTable
outZSaT = ZonalStatisticsAsTable(inZoneData, zoneField, inValueRaster,
outTable, "NODATA", "MEAN")
环境
许可信息
- Basic: 需要 Spatial Analyst
- Standard: 需要 Spatial Analyst
- Advanced: 需要 Spatial Analyst