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

2024-09-01 09:18
文章标签 创建 读写 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;

最后附上gdal201版本的编译完成的库,包含了geos。欢迎下载
http://download.csdn.net/download/u013230291/10000023

GDAL18下载地址:http://download.csdn.net/download/u013230291/9920097

这篇关于GDAL矢量(.shp)文件读写与创建的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1126531

相关文章

Python中使用uv创建环境及原理举例详解

《Python中使用uv创建环境及原理举例详解》uv是Astral团队开发的高性能Python工具,整合包管理、虚拟环境、Python版本控制等功能,:本文主要介绍Python中使用uv创建环境及... 目录一、uv工具简介核心特点:二、安装uv1. 通过pip安装2. 通过脚本安装验证安装:配置镜像源(可

Java中实现线程的创建和启动的方法

《Java中实现线程的创建和启动的方法》在Java中,实现线程的创建和启动是两个不同但紧密相关的概念,理解为什么要启动线程(调用start()方法)而非直接调用run()方法,是掌握多线程编程的关键,... 目录1. 线程的生命周期2. start() vs run() 的本质区别3. 为什么必须通过 st

Macos创建python虚拟环境的详细步骤教学

《Macos创建python虚拟环境的详细步骤教学》在macOS上创建Python虚拟环境主要通过Python内置的venv模块实现,也可使用第三方工具如virtualenv,下面小编来和大家简单聊聊... 目录一、使用 python 内置 venv 模块(推荐)二、使用 virtualenv(兼容旧版 P

使用Python和Matplotlib实现可视化字体轮廓(从路径数据到矢量图形)

《使用Python和Matplotlib实现可视化字体轮廓(从路径数据到矢量图形)》字体设计和矢量图形处理是编程中一个有趣且实用的领域,通过Python的matplotlib库,我们可以轻松将字体轮廓... 目录背景知识字体轮廓的表示实现步骤1. 安装依赖库2. 准备数据3. 解析路径指令4. 绘制图形关键

Linux lvm实例之如何创建一个专用于MySQL数据存储的LVM卷组

《Linuxlvm实例之如何创建一个专用于MySQL数据存储的LVM卷组》:本文主要介绍使用Linux创建一个专用于MySQL数据存储的LVM卷组的实例,具有很好的参考价值,希望对大家有所帮助,... 目录在Centos 7上创建卷China编程组并配置mysql数据目录1. 检查现有磁盘2. 创建物理卷3. 创

ShardingSphere之读写分离方式

《ShardingSphere之读写分离方式》:本文主要介绍ShardingSphere之读写分离方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录ShardingSphere-读写分离读写分离mysql主从集群创建 user 表主节点执行见表语句项目代码读写分

Java 如何创建和使用ExecutorService

《Java如何创建和使用ExecutorService》ExecutorService是Java中用来管理和执行多线程任务的一种高级工具,可以有效地管理线程的生命周期和任务的执行过程,特别是在需要处... 目录一、什么是ExecutorService?二、ExecutorService的核心功能三、如何创建

使用Python创建一个功能完整的Windows风格计算器程序

《使用Python创建一个功能完整的Windows风格计算器程序》:本文主要介绍如何使用Python和Tkinter创建一个功能完整的Windows风格计算器程序,包括基本运算、高级科学计算(如三... 目录python实现Windows系统计算器程序(含高级功能)1. 使用Tkinter实现基础计算器2.

CentOS和Ubuntu系统使用shell脚本创建用户和设置密码

《CentOS和Ubuntu系统使用shell脚本创建用户和设置密码》在Linux系统中,你可以使用useradd命令来创建新用户,使用echo和chpasswd命令来设置密码,本文写了一个shell... 在linux系统中,你可以使用useradd命令来创建新用户,使用echo和chpasswd命令来设

使用Python和Pyecharts创建交互式地图

《使用Python和Pyecharts创建交互式地图》在数据可视化领域,创建交互式地图是一种强大的方式,可以使受众能够以引人入胜且信息丰富的方式探索地理数据,下面我们看看如何使用Python和Pyec... 目录简介Pyecharts 简介创建上海地图代码说明运行结果总结简介在数据可视化领域,创建交互式地