Python ogr 读取shp，识别镂空面，并把镂空面坐标打印出来

要识别镂空面并打印其坐标，您可以使用Python中的ogr库来读取shp文件并进行处理。以下是一个示例代码：

from osgeo import ogr

# 打开shp文件
shp_file = "path/to/your/shp/file.shp"
driver = ogr.GetDriverByName("ESRI Shapefile")
data_source = driver.Open(shp_file, 0)

# 获取第一个图层
layer = data_source.GetLayer(0)

# 遍历要素
for feature in layer:
    geometry = feature.GetGeometryRef()

    # 判断是否为镂空面
    if geometry.GetGeometryName() == "POLYGON":
        ring = geometry.GetGeometryRef(0)
        if ring.IsClosed() and ring.GetPointCount() > 3:
            # 判断是否有内环
            if geometry.GetGeometryCount() > 1:
                print("Found a polygon with holes:", feature.GetFID())

                # 获取外环的坐标
                exterior_ring = geometry.GetGeometryRef(0)
                exterior_coordinates = []
                for i in range(exterior_ring.GetPointCount()):
                    point = exterior_ring.GetPoint(i)
                    exterior_coordinates.append((point[0], point[1]))
                print("Exterior Coordinates:", exterior_coordinates)

                # 获取内环的坐标
                for i in range(1, geometry.GetGeometryCount()):
                    interior_ring = geometry.GetGeometryRef(i)
                    interior_coordinates = []
                    for j in range(interior_ring.GetPointCount()):
                        point = interior_ring.GetPoint(j)
                        interior_coordinates.append((point[0], point[1]))
                    print("Interior Coordinates:", interior_coordinates)

# 关闭数据源
data_source = None

在上述代码中，我们首先打开shp文件并获取第一个图层。然后，遍历图层中的要素，判断每个要素的几何类型是否为多边形（即镂空面）。如果满足条件，则打印出找到的镂空面的要素ID。接下来，我们获取外环（多边形的边界）的坐标，并将其打印出来。然后，通过循环获取内环（镂空部分）的坐标，并将其分别打印出来。

请确保将代码中的`"path/to/your/shp/file.shp"`替换为实际的shp文件路径。

希望这可以满足您的需求！如果还有其他问题，请随时提问。