使用GDAL/GEOS求面特征的并集

存在这样一个示例的矢量文件，包含了两个重叠的面特征：

一个很常见的需求是求取这个矢量中所有面元素的并集，通过GDAL/GEOS很容易实现这个功能，具体代码如下：

#include <iostream>#include <gdal/ogrsf_frmts.h>using namespace std;bool WritePolygon(string filePath, OGRPolygon *pOgrMerged)
{//创建GDALDriver* driver = GetGDALDriverManager()->GetDriverByName("ESRI Shapefile");if (!driver){printf("Get Driver ESRI Shapefile Error！\n");return false;}GDALDataset* dataset = driver->Create(filePath.c_str(), 0, 0, 0, GDT_Unknown, NULL);OGRLayer* poLayer = dataset->CreateLayer("houseType", NULL, wkbPolygon, NULL);//创建特征{OGRFeature *poFeature = new OGRFeature(poLayer->GetLayerDefn());poFeature->SetGeometry(pOgrMerged);if (poLayer->CreateFeature(poFeature) != OGRERR_NONE){printf("Failed to create feature in shapefile.\n");return false;}}//释放GDALClose(dataset);dataset = nullptr;//GDALDestroyDriverManager();return true;
}int main()
{GDALAllRegister();CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");  //支持中文路径CPLSetConfigOption("SHAPE_ENCODING", "");  //解决中文乱码问题string filePath = "D:/Work/OSGWork/shpTest/test/src.shp";GDALDataset *poDS = (GDALDataset*)GDALOpenEx(filePath.c_str(), GDAL_OF_VECTOR, NULL, NULL, NULL);if (!poDS){printf("无法读取该文件，试检查格式是否正确！");return 1;}if (poDS->GetLayerCount()<1){printf("该文件的层数小于1，试检查格式是否正确！");return 1;}OGRLayer  *poLayer = poDS->GetLayer(0); //读取层poLayer->ResetReading();std::unique_ptr<OGRPolygon> pOgrMerged(new OGRPolygon());OGRFeature *poFeature;while ((poFeature = poLayer->GetNextFeature()) != NULL){//OGRGeometry *pGeo = poFeature->GetGeometryRef();OGRwkbGeometryType pGeoType = pGeo->getGeometryType();if (pGeoType == wkbPolygon){OGRPolygon  *pPolygon = (OGRPolygon*)pGeo;if (!pPolygon){continue;}OGRPolygon* pTemp = static_cast<OGRPolygon*>(pOgrMerged->Union(pPolygon));if (pTemp){pOgrMerged.reset(pTemp);}        }OGRFeature::DestroyFeature(poFeature);}GDALClose(poDS);poDS = nullptr;if (pOgrMerged && pOgrMerged->IsValid() && pOgrMerged->getExteriorRing()){string path = "D:/Work/OSGWork/shpTest/test/dst.shp";WritePolygon(path, pOgrMerged.get());}return 0;
}

在这段代码中，遍历了示例矢量文件中的每个面元素，求取了所有面元素的并集，得到最终一个面元素，并将这个面元素保存成新的矢量文件。这个矢量文件用ArcMap打开显示如下：

使用GDAL/GEOS求面特征的并集相关推荐

[NA]Lab5:反幂法求矩阵特征根
任务概述给定一个n*n矩阵A以及A的某一个特征根的近似值p,从一个给定向量 x⃗ (0) {\vec{x}}^{(0)}出发,求该特征根的高精度近似值,并且求出其对应的无穷范数下的单位特征向量. 函 ...
JAVA Hashset求集合的交集并集差集
JAVA Hashset求集合的交集并集差集 Hashset是Set接口较为常见的一个子类,该子类的最大特点是不允许保存重复的元素,并且所有的内容都采用散列(无序)的方式进行存储. package c ...
开源GIS库GDAL/GEOS在Windows下的编译与 Codeblocks配置
最近在研究开源GIS库GDAL,我是采用自己下载源码编译的方式,主要进行以下几个步骤: 编译GEOS,参考自GEOS的wiki 下载GEOS源代码,目前的版本为3.5.0,解压到C盘根目录.打开VS2 ...
Java中求集合交集、并集、差集
前言如果是自定义类,需要重写equals方法,如果使用Set集合,还需要重写hashCode方法:如下以User为例子: public class User{private Integer id;p ...
系统稳定性（基于matlab求传递函数特征根）
本文包含以下内容: 一.系统稳定性概念二.系统稳定性判断一.系统稳定性概念系统加入扰动后偏离了原来的状态,当把干扰去掉后,系统如果能恢复到原来的状态,则说明系统稳定. 二.系统稳定性判断脉冲信 ...
matlab 用古典雅可比方法求矩阵特征根（仅使用基础函数）
我们先看看<数值计算方法>这本书上关于古典雅可比方法的例题: %author FoddcusL FAFU %The Jacobi For root%输入:目标矩阵(input):目标接近的 ...
编译好的GDAL库,支持ECW格式,支持proj,支持geos
这几天因为工作需要,把GDAL重新编译了一下,现支持ECW格式,即可以用GDALOPEN读取JPEG2000格式的影像数据;编译时也添加了对geos和proj的支持. 内含3个文件夹,分别是gdal, ...
Java实验——定义一个类，该类中包含以下几个方法（静态）：实现两个字符串数组的逆序排序，输出结果为字符串数组；求两个整形数组的交集；求两个浮点型数组的并集；
目录一.实现两个字符串数组的逆序排序,输出结果为字符串数组解题思路: 具体代码: 运行截图: 二.求两个整形数组的交集具体代码: 运行截图: 三.求两个浮点型数组的并集解题思路: 具体代码: ...
C/C++语言实现一百个数组之间求任意数组的交集和并集
1:背景要求 1.C++ 2.给定一百个数组,求任意两个集合之间的交集和并集(最好能按照学生的方式做,不被发现的那种)撰写设计书,主要内容包括(500-1000字): ①背景介绍或实验目的: ②功能设 ...

使用GDAL/GEOS求面特征的并集

使用GDAL/GEOS求面特征的并集相关推荐

最新文章

热门文章