Arcgis小技巧【13】——数据统计(Statistics)相关工具汇总

欢迎关注同名微信公众号，更多文章推送：

在Arcgis中可以通过属性表中字段的【统计】功能或使用统计相关的工具对属性表进行数据统计。

在Arcgis工具箱中有一组【统计分析】工具集，不仅包含对属性数据执行标准统计分析（例如平均值、最小值、最大值和标准差）的工具，也包含对重叠和相邻要素计算面积、长度和计数统计的工具。

这里以统计分析工具集为基础，聊聊Arcgis中的数据统计工具。

---

一、属性表中字段的【统计】功能

属性表中字段的【统计】功能是最简单的数据统计方式，虽然功能简单，但胜在方便，所以使用率算是挺高的。

打开属性表，在需要统计的字段上右击，即可选择【统计】功能：

---

二、【汇总统计数据】工具（Statistics）

【汇总统计数据】工具是最常用的一个统计工具，通常用于按字段统计面积等。

下图就是一个典型的按【JQDLBM】字段来统计用地面积的例子：

生成的结果简单明了：

生成表中的【FREQUENCY】是按分类统计的图斑个数，第一行的【36】即【JQDLBM】值为【0101】的图斑有36个的意思。

实际上【统计分析】工具集里的另一个【频数】工具，就是实现这么个功能……，而且在统计类型里也有这么一个选项【COUNT】可以得到同样的结果。其它选项如下：

另一个参数【案例分组字段】就按某组字段来分组的意思，上面的例子以单个字段【JQDLBM】来分。

如果你的图斑是包括多个村庄，你想统计不同村庄的不同用地类型，就可以用2个字段作为分组字段：

工具运行结果如下：

---

三、【交集制表】工具（TabulateIntersection）

【交集制表】工具和上面的【汇总统计数据】工具有点类似，但应用场景有点不一样。

例如一个村庄可以分为几个片区，但这个分区并不在村庄用地的相关字段里，而是在另一个要素类中定义。

你想统计上图中这两个片区的用地面积，就可以使用【交集制表】工具。需要注意的是，这个分区要素并不一定要和统计要素完全重叠，工具运行的时候只会统计两个要素之间完全重叠的区域。

工具运行结果如下：

这里也会生成一个额外的【PERCENTAGE】字段，意思是该分区统计的面积占分区面积之比。例如南村的用地并没有布满区域，所以这个比值就只有38%。

---

四、【范围内汇总】工具（SummarizeWithin）

需要说明的是，【范围内汇总】工具是ArcGIS Pro中独有的，从功能上看，相当于【交集制表】工具的加强版，不仅能生成汇总表，还能复制一个分区范围的要素，并且将汇总数据添加到分区范围要素上。

工具位置也在【统计分析】工具集下。

从工具参数可以看出，有2个输出，一个是汇总表，还有一个是输出的分区要素，上面有添加的汇总信息。还有一些汇总参数，像少数、众数什么的，不过感觉用处不大。

【范围内汇总】工具生成的汇总表和【交集制表】工具差不多：

输出的分区要素：

和【交集制表】工具相比，感觉也就是锦上添花吧，说不上有什么特别优秀的地方，可以看情况选用。

---

五、【数据透视表】工具（PivotTable）

此工具通过在“输入表”中减少记录中的冗余并简化一对多关系来创建表。可以对生成的汇总表进一步处理，特别是用多个分组字段生成的汇总表。

工具位置在【数据管理工具】的【表】工具集下。

以下图的汇总表为例：

这张表的表达不够直观，特别是2个村庄之间的数据对比比较差，我们需要将【所属边界】放到横坐标（字段名）上（相当于Excel里的列）。

工具参数如下：

这几个参数经常搞混，多用几次，熟悉就好了。生成结果如下：

转换后的数据透视表就直观多了，导出Excel后也便于后续处理。

---

六、【字段统计数据转表】工具（FieldStatisticsToTable）

【字段统计数据转表】工具也是ArcGIS Pro独有的，工具位置在【数据管理工具】的【字段】工具集下。

这个工具怎么说呢，就是花里胡哨……。基本功能也就是【汇总统计数据】工具做的事，只是它能汇总的内容更多更花哨（无聊？）。

统计结果如下：

---

以上就是ArcGIS中和数据统计相关的内容，最常用的实际上也就是属性表字段的【统计】功能和【汇总统计数据】工具，其它工具也就是在特定场景下才用得上。不过可以看出ArcGIS给我们准备的功能或工具是很齐全的，也可见是个很有深度的软件。

---

七、自制脚本工具汇总

如果上面的工具都不能满足你的要求，还可以尝试用脚本来实现。

如上图所示，要分组统计的字段有2个，但和汇总工具不同的是，汇总结果并非同时满足2个字段值，而是分开统计，即按字段1、字段2分别统计，再将表合并起来。在实际工作中，你可以得到一级地类和二级地类的统计表格。

用脚本工具是很自由的，你可以再增加一点其它功能，比如说，生成一个【合计值】，即所有地块的总面积，再增加一个字段，统计各类用地占总用地面积的【比例】。

甚至还可以对【面积单位进行转换】，如转成公顷、平方公里等，以及保留【固定的小数位数】。

生成的示例结果如下：

脚本工具参数设置如下：

限于篇幅，具体代码不详细说了，直接贴上：

# -*- coding: utf-8 -*-
import arcpy
import os

input_table = arcpy.GetParameterAsText(0)       # 输入表
input_field = arcpy.GetParameterAsText(1)       # 输入要统计的字段
input_sta = arcpy.GetParameterAsText(2)         # 输入合计字段值(可选)
is_ha = arcpy.GetParameter(3)                   # 是否转成公顷
is_round = arcpy.GetParameter(4)                # 是否保留2位小数
output_table = arcpy.GetParameterAsText(5)      # 输出表

# 获取文档位置
default_path = os.path.dirname(os.path.dirname(os.path.dirname(__file__)))  # D:\【软件资料】\GIS相关\GisPro工具箱
tem_gdb = default_path + r'\3-默认数据库\GisPro工具箱.gdb'  # 临时数据库

# 汇总多个字段
def MultyStatic():
    # 解析字段列表
    list_field = input_field.split(';')
    # 汇总多个字段
    list_table = []
    for i in range(len(list_field)):
        area_static = arcpy.Statistics_analysis(input_table, tem_gdb + r'\area_static_' + str(i), [['Shape_Area', "SUM"]], [list_field[i]])
        arcpy.AddField_management(area_static, '统计字段', "TEXT")   # 添加公共字段
        arcpy.CalculateField_management(area_static, '统计字段',expression='!' + list_field[i] + '!')    # 计算公共字段
        list_table.append(area_static)
    # 汇总总量
    area_total = arcpy.Statistics_analysis(input_table, tem_gdb + r'\area_static_total', [['Shape_Area', "SUM"]])
    arcpy.AddField_management(area_total, '统计字段', "TEXT")  # 添加公共字段
    arcpy.CalculateField_management(area_total, '统计字段', expression='\'' + input_sta + '\'')  # 计算公共字段
    list_table.append(area_total)
    # 合并表格
    area_all = arcpy.Merge_management(list_table, tem_gdb + r'\yd_all')
    # 去除空值，输出table
    arcpy.TableSelect_analysis(area_all, output_table, where_clause="统计字段 IS NOT NULL")
    # 添加统计值和比例字段
    arcpy.AddField_management(output_table, '统计值', 'FLOAT')
    arcpy.AddField_management(output_table, '比例', 'FLOAT')
    # 计算统计值和比例字段
    arcpy.CalculateField_management(output_table, '统计值', expression='!SUM_SHAPE_Area!')
    total = float(str(GetFieldValue(output_table, '统计值', "统计字段 = '" + input_sta + "'")))
    fields = ['统计值', '比例']
    with arcpy.da.UpdateCursor(output_table, fields) as cursor:
        for row in cursor:
            row[1] = row[0]/total*100
            if is_ha is True:               # 转换成公顷
                row[0] = row[0]/10000
            cursor.updateRow(row)
    # 判断是否保留2位小数
    with arcpy.da.UpdateCursor(output_table, fields) as cursor:
        for row in cursor:
            if is_round is True:
                row[0] = round(row[0], 2)
                row[1] = round(row[1], 2)
            cursor.updateRow(row)
    # 删除中间字段
    arcpy.DeleteField_management(output_table, list_field)
    arcpy.DeleteField_management(output_table, ['FREQUENCY', 'SUM_SHAPE_Area'])
    # 删除中间要素
    arcpy.Delete_management(area_static)
    arcpy.Delete_management(area_total)
    arcpy.Delete_management(area_all)

# 获取要素或表的某一个字段的第一个值，默认字段为Shape_Area。筛选表达式为：keyField=key
def GetFieldValue(input_table, value_field, sql):
    # 按sql语句筛选出所要的表格
    table = arcpy.TableSelect_analysis(input_table, r'memory\table_sel', sql)
    # 定义一个空值，用于储存返回的结果
    filedValue = '0'
    # 读取值【如果筛选值有多个，只取第一个值】
    with arcpy.da.SearchCursor(table, [value_field]) as cursor:
        for row in cursor:
            filedValue = row[0]
            break
    arcpy.SetParameterAsText(3, str(filedValue))


if __name__ == '__main__':
    MultyStatic()