在GIS 工作流中实现数据处理(3)

统计分析

现在数据已经清洗完成，接下来我们对特定区域的土地利用类型进行统计分析。假设我们有一个城市中心区域的边界 shapefile 文件，我们想要统计该区域内每种土地利用类型的面积。

# 设置城市中心区域边界
city_center_boundary = r"C:\GIS_Data\City_Center_Boundary.shp"# 创建一个临时的要素类，用于存储城市中心区域内的土地利用数据
temp_feature_class = r"C:\GIS_Data\Land_Use\Temp_Land_Use.shp"
arcpy.Clip_analysis(feature_classes[0], city_center_boundary, temp_feature_class)# 统计每种土地利用类型的面积
statistics_table = r"C:\GIS_Data\Land_Use\Land_Use_Statistics.dbf"
arcpy.Statistics_analysis(temp_feature_class, statistics_table, [["Shape_Area", "SUM"]], "Land_Use_Type")

通过 arcpy.Clip_analysis 将城市中心区域边界内的土地利用数据裁剪出来，然后使用 arcpy.Statistics_analysis 对裁剪后的数据进行统计分析，生成一个统计表格，其中包含每种土地利用类型的总面积。