gdal读取geojson

使用GDAL读取GeoJSON文件可以使用以下步骤：

1. 导入必要的模块：

import os
from osgeo import gdal
from osgeo import ogr

2. 打开GeoJSON文件：

filename = "path/to/your/geojson/file.geojson"
ds = ogr.Open(filename)

3. 获取GeoJSON中的图层：

layer = ds.GetLayer()

4. 获取图层的投影信息：

spatialRef = layer.GetSpatialRef()

5. 获取图层的要素数量：

featureCount = layer.GetFeatureCount()

6. 遍历图层中的要素并处理：

for i in range(featureCount):
    feature = layer.GetFeature(i)
    geometry = feature.GetGeometryRef()
    # 处理要素

这样就完成了使用GDAL读取GeoJSON文件的操作。