Python利用GDAL库：统计“志远县”图斑个数与“罗丰县”面积总和

在Python中，我们可以使用GDAL库来读取和操作Shapefile数据。以下是一个示例代码，假设你已经安装了GDAL库： ```python from osgeo import ogr # 1. 读取name字段中“志远县”图斑的个数 count = 0 driver = ogr.GetDriverByName("ESRI Shapefile") data_folder = "data1" shp_file = os.path.join(data_folder, "township.shp") # 打开shapefile dataset = driver.Open(shp_file, 0) if dataset is None: print("无法打开shapefile") else: layer = dataset.GetLayer() # 遍历图层中的每个feature for feature in layer: # 获取name字段的值 name = feature.GetField("name") # 检查是否为"志远县" if name == "志远县": count += 1 print(f"名称为'志远县'的图斑个数: {count}") # 2. 读取area字段（面积），并计算“罗丰县”的面积之和 total_area = 0 for feature in layer: area = feature.GetField("area") # 检查是否为"罗丰县" if feature.GetField("name") == "罗丰县": total_area += area print(f"罗丰县的面积之和: {total_area}") ``` 请确保将`data1`和`township.shp`替换为你的实际文件路径。这段代码首先打开shapefile，然后遍历图层中的每个图斑，检查name字段的值是否为"志远县"或"罗丰县"，并进行相应的计数或面积累加。 ######[AI写代码神器 | 389点数解答 | 2024-06-30 05:13:27]