使用OGR与GDAL读写矢量数据格式详解

"GDAL之OGR入门，用于读写矢量数据格式" GDAL（Geospatial Data Abstraction Library）是一个开源的库，它提供了对多种栅格和矢量地理数据格式的支持。OGR是GDAL的一部分，专注于处理矢量数据，提供了读取、写入和操作各种矢量数据格式的能力。这篇文档主要介绍了如何使用OGR的C++ API来读取和写入矢量数据。 1. OGR注册所有格式 在使用OGR读取数据前，首先要注册所有支持的数据格式。这可以通过调用`OGRRegisterAll()`函数完成，该函数会注册GDAL/OGR支持的所有矢量数据格式，如ESRI Shapefile、GeoJSON、GML、KML等。 2. 打开数据源 要读取数据，需要创建一个`OGRDataSource`对象来代表数据源。数据源可以是本地文件、数据库或网络服务。例如，打开一个Shapefile文件，可以使用`OGRSFDriverRegistrar::Open()`函数，传入文件名和一个布尔值表示是否需要写入权限。 ```cpp OGRDataSource* poDS = OGRSFDriverRegistrar::Open("point.shp", FALSE); ``` 3. 获取图层 `OGRDataSource`对象可以包含多个图层（Layer）。可以使用`GetLayerCount()`获取图层数量，`GetLayer()`或`GetLayerByName()`根据索引或名称获取特定图层。例如，获取名为"point"的图层： ```cpp OGRLayer* poLayer = poDS->GetLayerByName("point"); ``` 4. 读取要素（Features） 图层中的数据是以要素的形式存储，每个要素通常包含几何和属性信息。遍历图层中的所有要素，可以使用`ResetReading()`重置读取指针，然后连续调用`GetNextFeature()`来获取下一个要素。 ```cpp poLayer->ResetReading(); OGRFeature* poFeature; while ((poFeature = poLayer->GetNextFeature()) != NULL) { // 处理poFeature... } ``` 5. 设置过滤器 在遍历要素前，可以设置属性过滤器（Attribute Filter）或空间过滤器（Spatial Filter）来限制读取的要素。属性过滤器基于要素的属性值，而空间过滤器则基于要素的空间位置。 6. 访问几何和属性 每个`OGRFeature`对象包含一个`OGRGeometry`实例，表示其几何形状，以及一组属性字段。可以通过`GetGeometryRef()`获取几何，`GetFieldCount()`和`GetFieldDefnRef()`获取字段数量和定义，`GetFieldAsXXX()`系列方法获取字段值。 7. 写入数据 写入矢量数据时，首先创建一个`OGRDataSource`，然后创建新的`OGRLayer`，定义其属性字段和几何类型。接着创建`OGRFeature`，设置其几何和属性值，最后添加到图层中。 8. 保存和关闭数据源 完成读写操作后，记得使用`FlushCache()`将内存中的更改写入磁盘，并使用`DeleteDataSource()`或`Close()`关闭数据源。 GDAL/OGR提供了一套强大的接口，使得开发者能够轻松地处理各种矢量数据格式。无论是读取现有数据，还是创建新的数据文件，OGR都能提供灵活且高效的解决方案。通过学习和掌握这些基本操作，可以构建复杂的应用程序，处理大规模的地理空间数据。