Python用ogr实现遍历LINESTRING中的点

使用Python的GDAL/OGR库可以轻松实现遍历LINESTRING中的点。以下是一个简单的示例代码：

from osgeo import ogr

# 打开shapefile文件
driver = ogr.GetDriverByName('ESRI Shapefile')
dataSource = driver.Open('path/to/shapefile.shp', 0)
layer = dataSource.GetLayer()

# 循环遍历feature
for feature in layer:
    # 获取geometry
    geometry = feature.GetGeometryRef()
    # 如果是LINESTRING类型
    if geometry.GetGeometryName() == 'LINESTRING':
        # 遍历线上的点
        for i in range(geometry.GetPointCount()):
            # 获取点坐标
            point = geometry.GetPoint(i)
            print(point)

在这个示例中，我们打开一个shapefile文件并获取其第一个图层。然后，我们循环遍历每个feature，并检查其geometry类型是否为LINESTRING。如果是，我们遍历线上的每个点并输出其坐标。