OGR中空间叠加函数Union

最新推荐文章于 2023-07-25 19:00:21 发布

weixin_30716725

最新推荐文章于 2023-07-25 19:00:21 发布

阅读量2.7k

点赞数 1

CC 4.0 BY-SA版权

文章标签： python

原文链接：http://www.cnblogs.com/abella/p/9473121.html

本文介绍如何利用OGR库中的Union函数实现地理信息系统中两个图层的合并操作，并详细解析了函数参数的意义及使用方法。同时提供了具体的代码实例，展示了如何设置参数以达到期望的合并效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在OGRLayer的对象中调用Union函数了。先看一下Union函数的原型：

Union(OGRLayer pLayerMethod, OGRLayer pLayerResult, char **papszOptions=NULL, GDALProgressFunc pfnProgress=NULL, void *pProgressArg=NULL)，根据OGR提供的解释来看，前三个参数是比较重要的：

pLayerMethod：要跟当前图层做Union的那个图层，如果我要对layer1、layer2做Union操作，通过layer1对象调用Union函数，那么layer2就是pLayerMethod，不能为空。

pLayerResult：很明显，存储Union结果的图层，也不能为空。这里可以是用ArcGIS创建好的没有要素的图层，也可以是在程序中用OGR创建的新图层，个人倾向与在程序中用代码创建一个，简单的代码就能实现。

papszOptions：这个二维指针中有其实包括了四个参数——SKIP_FAILURES=YES/NO、PROMOTE_TO_MULTI=YES/NO、INPUT_PREFIX=string、METHOD_PREFIX=string。当SKIP_FAILURES=YES的时候，如果某个要素的合并出现错误，则会跳过该要素继续后续要素的合并。PROMOTE_TO_MULTI=YES时，可以将单空间对象转换为多空间对象，如将Polygons转换为MultiPolygons， LineStrings转换为MultiLineStrings。INPUT_PREFIX=string，输入图层的属性字段在结果图层中的前置标记。METHOD_PREFIX=string，操作图层的属性字段在结果图层中的前置标记。还是拿layer1和layer2来举例，layer1、layer2合并后的结果图层中包含有layer1和layer2的所有属性字段，那么如果我们设置INPUT_PREFIX=1，METHOD_PREFIX=2，那么在结果图层的字段中，来自layer1的所有字段将在字段名称前加上前缀"1"，来自layer2的将加前缀"2"。

后面的两个参数没做尝试。

下面看一个具体的代码示例

	char *filePath = "D:\\CeShi_Data\\CESHI_NEW";
	char *layerName1 = "layer1";
	char *layerName2 = "layer2";
	OGRLayer *pLayer1 = NULL;
	OGRLayer *pLayer2 = NULL;
	OGRDataSource *pODS = NULL;

	OGRRegisterAll();
	pODS = OGRSFDriverRegistrar::Open(filePath,TRUE);

	//读取取要进行Union的两个图层
	pLayer1 = pODS->GetLayerByName(layerName1);
	pLayer2 = pODS->GetLayerByName(layerName2);
	//创建结果图层
	OGRLayer *pResultLayer = NULL;
	pResultLayer = pODS->CreateLayer("result",pLayer2->GetSpatialRef(),wkbMultiPolygon,NULL);
	//配置Union函数中的第三个参数
	char **p = new char *[4];
	p[0] = "SKIP_FAILURES=YES";
	p[1] = "PROMOTE_TO_MULTI=YES";
	p[2] = "INPUT_PREFIX=1";
	p[3] = "METHOD_PREFIX=2";
	pLayer2->Union(pLayer1,pResultLayer,p,NULL,NULL);  

	//将对pResultLayer的编辑写入文件，如果不加这句,result文件中将没有记录
	pResultLayer->SyncToDisk();
	OGRDataSource::DestroyDataSource(pODS);

如下图：

layer1

layer2:

结果图层result:

如果layer1、layer2中的属性字段相同，如示例中，在结果图层中想保存一份属性字段，则在创建结果图层的时给结果图层创建同layer1一样的属性字段集即可，在代码中做如下修改：

	//创建结果图层
	OGRLayer *pResultLayer = NULL;
	pResultLayer = pODS->CreateLayer("result",pLayer2->GetSpatialRef(),wkbMultiPolygon,NULL);
	//为结果图层创建跟输入图层layer1相同的属性字段
	OGRFeatureDefn *pOGRFeatureDefn = pLayer1->GetLayerDefn();
	for (int i = 0 ; i < pOGRFeatureDefn->GetFieldCount(); i++)
		pResultLayer->CreateField( pOGRFeatureDefn->GetFieldDefn(i));
	pLayer2->Union(pLayer1,pResultLayer,p,NULL,NULL);

生成的结果图层的属性表如下：