如何在C++中使用GDAL库将SHP文件转换为GeoJSON格式？

在C++中，你可以通过GDAL（Geospatial Data Abstraction Library）库来处理地理空间数据，包括将SHP（Shapefile）格式转换为GeoJSON。以下是一个简化的步骤说明：

1. **包含头文件**：
首先，你需要在你的代码中包含必要的GDAL库头文件：

   #include <gdal.h>
   #include <ogr ogr_api.h>

2. **初始化GDAL**：
在使用GDAL之前，需要初始化GDAL环境：

   GDALAllRegister();
   OGRSFDriverManager::GetDriverByName("ESRI Shapefile");

3. **读取SHP文件**：
使用`OGRDataSource`来打开SHP文件，并获取数据源对象：

   OGRDataSource* dataSource = OGROpen("input.shp", FALSE);

4. **获取layer（层）**：
从数据源中获取对应的数据层：

   OGRLayer* layer = dataSource->GetLayer(0);

5. **创建GeoJSON writer**：
使用GDAL的内置支持，如`OGRGeoJSONWriter`，创建一个GeoJSON文件的对象并设置输出路径：

   OGRGeoJSONWriter oGeoJSONWriter;
   oGeoJSONWriter.SetOutputFile("output.geojson");

6. **遍历features（特征）**：
对于数据层中的每个特征，将其转换为GeoJSON并写入：

   for (OGRFeature* feature = layer->GetNextFeature(); feature; feature = layer->GetNextFeature()) {
       OGRGeoJSONFeatureBuilder builder;
       builder.AddGeometry(feature->GetGeometryRef());
       OGRGeoJSONFeature* geoJsonFeature = builder.BuildFeature(nullptr);
       oGeoJSONWriter.WriteFeature(geoJsonFeature);
   }

7. **清理资源**：
最后别忘了关闭和释放所有已分配的资源：

   delete geoJsonFeature;
   layer->ReleaseResultSet(NULL);
   dataSource->ReleaseResultSet(layer);
   dataSource->Close();
   CPLFree(dataSource);