深圳珠宝品牌网站设计,做网站的边框,sem竞价托管,php网站开发入门统计分析
现在数据已经清洗完成#xff0c;接下来我们对特定区域的土地利用类型进行统计分析。假设我们有一个城市中心区域的边界 shapefile 文件#xff0c;我们想要统计该区域内每种土地利用类型的面积。
# 设置城市中心区域边界
city_center_boundary rC:\GIS_Da…统计分析
现在数据已经清洗完成接下来我们对特定区域的土地利用类型进行统计分析。假设我们有一个城市中心区域的边界 shapefile 文件我们想要统计该区域内每种土地利用类型的面积。
# 设置城市中心区域边界
city_center_boundary rC:\GIS_Data\City_Center_Boundary.shp# 创建一个临时的要素类用于存储城市中心区域内的土地利用数据
temp_feature_class rC:\GIS_Data\Land_Use\Temp_Land_Use.shp
arcpy.Clip_analysis(feature_classes[0], city_center_boundary, temp_feature_class)# 统计每种土地利用类型的面积
statistics_table rC:\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 对裁剪后的数据进行统计分析生成一个统计表格其中包含每种土地利用类型的总面积。