GDAL矢量(.shp)文件读写与创建

本文主要是介绍GDAL矢量(.shp)文件读写与创建，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

gdal18版本读写矢量与创建，具体李明录老师的书本更详细
以下代码是求两个矢量的交集：

    GDALAllRegister();OGRRegisterAll();CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");OGRDataSource *podata = OGRSFDriverRegistrar::Open(str_checkshp.c_str(),1);///str_checkshp为矢量路径，字符串类if (podata == NULL){cout <<"file:" <<str_checkshp << " open checkshp fail！！" << endl;}OGRLayer *poSourceLayer = podata->GetLayer(0);if (poSourceLayer == NULL){cout << "file:" << str_checkshp << "open checkshp_layer fail！！" << endl;}poSourceLayer->ResetReading();OGRFeature* pFea = poSourceLayer->GetNextFeature();if (pFea == NULL){cout << "file:" << str_checkshp << " feature is not exist！！" << endl;}OGRGeometry *poGeometry= poFeature->GetGeometryRef();
///获取矢量2的geo
OGRDataSource *podata2 = OGRSFDriverRegistrar::Open(str_checkshp2.c_str(),1);
///str_checkshp2为矢量路径，字符串类型if (podata2 == NULL){cout <<"file:" <<str_checkshp2 << " open checkshp fail！！" << endl;}OGRLayer *poSourceLayer2 = podata2->GetLayer(0);if (poSourceLayer2 == NULL){cout << "file:" << str_checkshp2 << "open checkshp_layer fail！！" << endl;}poSourceLayer2->ResetReading();OGRFeature* pFea2 = poSourceLayer2->GetNextFeature();if (pFea2 == NULL){cout << "file:" << str_checkshp2 << " feature is not exist！！" << endl;}OGRGeometry *poGeometry2= poFeature2->GetGeometryRef();///创建矢量OGRSFDriver *poDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName(pszDriverName);///获取驱动if (poDriver == NULL){cout << "驱动不可用" << endl;return FALSE;}OGRDataSource *poNewDS = poDriver->CreateDataSource(strShpNewfliedir.c_str(), NULL);OGRSpatialReference *sp = poSourceLayer2->GetSpatialRef();OGRSpatialReference* oSRS = sp;///获取空间参考OGRLayer *poNewLayer = poNewDS->CreateLayer("newlayer", oSRS, wkbPolygon, NULL);
OGRFeatureDefn *poDefn = poLayer->GetLayerDefn();OGRFeature *poFeatureIntersection = OGRFeature::CreateFeature(poDefn);OGRGeometry *poGeometry3 = poGeometry1->Difference(poGeometry2);poFeatureIntersection->SetGeometry(poGeometry3);poLayer->CreateFeature(poFeatureIntersection);OGRFeature::DestroyFeature(poFeatureIntersection);OGRFeature::DestroyFeature(pFea );OGRFeature::DestroyFeature(pFea2);GDALClose(poNewDS);GDALClose(podata2);GDALClose(podata);return 0;

GDAL201矢量读写与创建
暂时没有中文书籍，就记录了下来，方便拷贝，提高码代码效率

if (OGRGeometryFactory::haveGEOS() == false){cout << "GDAL库未包含GEOS库" << endl;return 0;}// 打开栅格文件GDALDataset* poSrcDS1 = (GDALDataset*)GDALOpenEx(pszSrc1File, GDAL_OF_VECTOR, NULL, NULL, NULL);if (poSrcDS1 == NULL){return 0;}OGRLayer* poLayer1 = poSrcDS1->GetLayer(0);poLayer1->ResetReading();OGRFeature *poFeature1 = poLayer1->GetNextFeature();OGRGeometry *poGeometry1 = poFeature1->GetGeometryRef();GDALDataset* poSrcDS2 = (GDALDataset*)GDALOpenEx(pszSrc21File, GDAL_OF_VECTOR, NULL, NULL, NULL);if (poSrcDS2 == NULL){return 0;}OGRLayer* poLayer2 = poSrcDS2->GetLayer(0);poLayer2->ResetReading();OGRFeature *poFeature2 = poLayer2->GetNextFeature();OGRGeometry *poGeometry2= poFeature2->GetGeometryRef();OGRGeometry *poGeometry3 = poGeometry1->Difference(poGeometry2);// 创建输出矢量文件GDALDriver *poDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName("ESRI Shapefile");if (poDriver == NULL){printf("%s驱动不可用！\n", "ESRI Shapefile");GDALClose((GDALDatasetH)poSrcDS1);GDALClose((GDALDatasetH)poSrcDS2);return 0;}//根据文件名创建输出矢量文件GDALDataset* poDstDS = poDriver->Create(pszDstFile, 0, 0, 0, GDT_Unknown, NULL);if (poDstDS == NULL){GDALClose((GDALDatasetH)poSrcDS1);GDALClose((GDALDatasetH)poSrcDS2);return 0;}// 定义空间参考，与输入图像相同OGRSpatialReference *poSpatialRef = new OGRSpatialReference(poSrcDS1->GetProjectionRef());if (poSpatialRef==NULL){return 0;}OGRLayer* poLayer = poDstDS->CreateLayer("Result", poSpatialRef, wkbPolygon, NULL);if (poDstDS == NULL){GDALClose((GDALDatasetH)poSrcDS1);GDALClose((GDALDatasetH)poSrcDS2);delete poSpatialRef;poSpatialRef = NULL;return 0;}OGRFeatureDefn *poDefn = poLayer->GetLayerDefn();OGRFeature *poFeatureIntersection = OGRFeature::CreateFeature(poDefn);poFeatureIntersection->SetGeometry(poGeometry3);poLayer->CreateFeature(poFeatureIntersection);OGRFeature::DestroyFeature(poFeatureIntersection);OGRFeature::DestroyFeature(poFeature1);OGRFeature::DestroyFeature(poFeature2);GDALClose(poSrcDS1);GDALClose(poSrcDS2);GDALClose(poDstDS);return 0;