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

当识别到环状面时，您可以使用`ogr`库中的方法来获取环状面的坐标，并将其打印出来。以下是修改后的代码：

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:
            print("Found a closed ring polygon:", feature.GetFID())
            # 获取环状面的坐标
            coordinates = []
            for i in range(ring.GetPointCount()):
                point = ring.GetPoint(i)
                coordinates.append((point[0], point[1]))
            print("Coordinates:", coordinates)

# 关闭数据源
data_source = None

在上述代码中，当找到一个闭合环的多边形时，我们首先打印出找到的多边形的要素ID，然后使用循环遍历多边形的每个点，并将其坐标存储在`coordinates`列表中。最后，将`coordinates`列表打印出来。

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

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