标注 | 说明 | 数据类型 |
输入表 | 包含用于计算统计数据的字段的输入表。 | Table View; Raster Layer |
输出表 | 输出表用于存储计算的统计数据。 | Table |
统计字段 (可选) | 指定包含用于计算指定统计数据的属性值的一个或多个数值字段。 可以指定多项统计和字段组合。 空值将被排除在所有统计计算之外。 可使用第一种和最后一种统计来对文本属性字段进行汇总。 可使用任何一种统计来对数值属性字段进行汇总。 可用统计类型如下:
| Value Table |
案例分组字段 (可选) | “输入”中用于为每个唯一属性值(如果指定多个字段,则为属性值组合)单独计算统计数据的一个或多个字段。 | Field |
摘要
为表中字段计算汇总统计数据。
使用情况
输出表将由包含统计运算结果的字段组成。
使用此工具可执行以下统计运算:总和、平均值、最小值、最大值、范围、标准差、计数、第一个、最后一个、中值、方差和唯一值。
将使用以下命名约定为每种统计类型创建字段:SUM_<field>、MEAN_<field>、MIN_<field>、MAX_<field>、RANGE_<field>、STD_<field>、COUNT_<field>、FIRST_<field>、LAST_<field>、MEDIAN_<field>、VARIANCE_<field> 和 UNIQUE_<field>(其中 <field> 是计算统计数据的输入字段的名称)。 当输出表是 dBASE 表时,字段名称会被截断为 10 个字符。
如果已指定案例分组字段,则单独为每个唯一属性值计算统计数据。 如果未指定案例分组字段,则输出表中将仅包含一条记录。 如果已指定一个案例分组字段,则每个案例分组字段值均有一条对应的记录。
空值将被排除在所有统计计算之外。 例如,10、5 和空值的平均值为 7.5 ((10+5)/2)。
使用图层时,仅使用当前所选要素计算统计数据。
参数
arcpy.analysis.Statistics(in_table, out_table, {statistics_fields}, {case_field})
名称 | 说明 | 数据类型 |
in_table | 包含用于计算统计数据的字段的输入表。 | Table View; Raster Layer |
out_table | 输出表用于存储计算的统计数据。 | Table |
statistics_fields [[field, {statistic_type}],...] (可选) | 指定包含用于计算指定统计数据的属性值的一个或多个数值字段。 可以指定多项统计和字段组合。 空值将被排除在所有统计计算之外。 可使用第一种和最后一种统计来对文本属性字段进行汇总。 可使用任何一种统计来对数值属性字段进行汇总。 可用统计类型如下:
| Value Table |
case_field [case_field,...] (可选) | “输入”中用于为每个唯一属性值(如果指定多个字段,则为属性值组合)单独计算统计数据的一个或多个字段。 | Field |
代码示例
以下 Python 窗口脚本演示了如何在即时模式下使用 Statistics 工具。
import arcpy
arcpy.env.workspace = "C:/data/Habitat_Analysis.gdb"
arcpy.Statistics_analysis("futrds", "C:/output/output.gdb/stats", [["Shape_Length", "SUM"]], "NM")
以下独立脚本汇总了主要道路 150 英尺范围内的植被区域。
# Description: Summarize the vegetation by area within 150 feet of major roads
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/data"
# Set local variables
inRoads = "majorrds.shp"
outBuffer = "C:/output/output.gdb/buffer_out"
bufferDistance = "250 feet"
inVegetation = "Habitat_Analysis.gdb/vegtype"
outClip = "C:/output/output.gdb/clip_out"
joinField = "HOLLAND95"
joinTable = "c:/data/vegtable.dbf"
joinedField = "HABITAT"
outStatsTable = "C:/output/output.gdb/stats_out"
statsFields = [["Shape_Area", "SUM"]]
# Execute Buffer to get a buffer of major roads
arcpy.Buffer_analysis(inRoads, outBuffer, bufferDistance, dissolve_option="ALL")
# Execute Clip using the buffer output to get a clipped feature class of
# vegetation
arcpy.Clip_analysis(inVegetation, outBuffer, outClip)
# Execute JoinField to add the vegetation type
arcpy.JoinField_management(outClip, joinField, joinTable, joinField, joinedField)
# Execute Statistics to get the area of each vegetation type within the
# clipped buffer.
arcpy.Statistics_analysis(outClip, outStatsTable, statsFields, joinedField)
以下独立脚本循环遍历数据集的属性字段并构造 statistics_fields 参数,从而计算各个数值字段的 SUM 统计量。
# Description: Script that runs the Summary Statistic tool to calculate the
# Sum statistic for every numeric field based on a unique case
# field.
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/data/f.gdb"
# Set local variables
intable = "intable"
outtable = "sumstats"
casefield = "Name"
stats = []
# Loop through all fields in the Input Table
for field in arcpy.ListFields(intable):
# Just find the fields that have a numeric type
if field.type in ("Double", "Integer", "Single", "SmallInteger"):
# Add the field name and Sum statistic type to the list of fields to
# summarize
stats.append([field.name, "Sum"])
# Correct formatting of stats [["Field1", "Sum"], ["Field2", "Sum"], ...]
# Run the Summary Statistics tool with the stats list
arcpy.Statistics_analysis(intable, outtable, stats, casefield)
以下脚本使用 pandas DataFrame 来访问和显示 Statistics 工具的表格结果。
import arcpy
import pandas
import os
arcpy.env.overwriteOutput = True
in_table = r"d:\data\states.shp"
out_table = r"in_memory\stats_table"
stat_fields = [['POP1990', 'SUM'], ['POP1997', 'SUM']]
stats = arcpy.Statistics_analysis(in_table, out_table, stat_fields,
case_field='SUB_REGION')
# Get a list of field names to display
field_names = [i.name for i in arcpy.ListFields(out_table) if i.type != 'OID']
# Open a cursor to extract results from stats table
cursor = arcpy.da.SearchCursor(out_table, field_names)
# Create a pandas dataframe to display results
df = pandas.DataFrame(data=[row for row in cursor],
columns=field_names)
print(df)
许可信息
- Basic: 是
- Standard: 是
- Advanced: 是