GDAL线面互转换(2)
在上一个文章中介绍了线转化为面和面转化为线,其主要的实现思路就是把面中的点取出来构成线,把线中的点取出来构成面,实际上就是一个硬拷贝,无奈客户的实际需求并非如此,客户想要线转面的时候几条相交线构成面,面转线的时候相同的线去除,所以又重新对功能进行了调整。
关于线构面,这个过程中也是挺曲折的,在GDAL的群里问了好久,大家给的答案一致是需要自己写,线构面的算法需要使用左转或者右转算法,并且在网上查了相关的资料,发现蒋波涛编著的《插件式GIS应用框架的设计与实现:基于C#和AE 9.2》一书中有ArcGIS实现GDAL的完整代码,没办法,埋头去干吧,先做了一部分前期的验证,验证的时候,公司一大牛问我OGRGeometry中的Polygonize ()方法是干啥的?仔细查阅相关资料,就是要做线转面的啊!赶紧验证,发现失败,后来又仔细分析,似乎该方法需要OGRGeometry的空间拓扑关系正确,如何保证这点儿呢?考虑用Union来把一根根线合并了拓扑关系应该没有问题吧,实际测试了一下,的确正确,心中窃喜。后来拿实际的省界线来合成省面,发现效率超级慢,跟踪代码发现,OGRGeometry合并的时候随着合并线段的增加时间也在快速的增加……问题找到,上网查找解决方案,看到两篇文章(http://www.docin.com/p-1155026547.html)(http://www.docin.com/p-1238395230.html)介绍说可以考虑使用UnionCascaded(级联求并)可以大大的加快效率,实际验证确不知道怎么使用,考虑到项目的进度要求,这个问题暂且搁置了,后续再提高吧。
对于面构线,就相对简单了不少,先把一个个的面转化成线(硬拷贝,参考上一博文中的方法),再把转化成的线Union了,生成一个拓扑关系正确的OGRMultiPolygon,再调用Simplify方法,得到的结果即为想要的线,不过当线特别多还很复杂的时候效率也高,问题和线构面的时候一样。
在此补充上源代码,在此留个记录,后续考虑怎么解决这个问题吧……
/*
* @brief ConvertPolygonToPolyline 面图层转换为线图层
* @param[in] tString polylinePath 转换后线图层文件路径
* @param[in] OGRLayer* pLayer 要转换的面图层文件
* @param[in] Envelope envelope 要转换数据的范围
* @param[in] vector<long> vecFIDs 选中的要素ID列表
* @param[in] pOGRSpatialReference 要转换的数据的空间参考(如果为空表示坐标系信息不变)
* @return bool 是否成功
* @author
* @date
* @note 2015年11月04日 小八创建;
*/
bool FeatureLayerOperator::ConvertPolygonToPolylineEx(tString polylinePath,OGRLayer* pLayer,Envelope envelope,vector<long> vecFIDs,OGRSpatialReference* pOGRSpatialReference)
{
// 判断
if(pLayer==NULL) return false;
if(true==polylinePath.empty()) return false; // 坐标系读取
OGRSpatialReference* pOGRSpatialReference_Source=pLayer->GetSpatialRef();
bool isSameCoordSystem = false;
if(pOGRSpatialReference == NULL)
{
pOGRSpatialReference=pOGRSpatialReference_Source;
isSameCoordSystem=true;
}
else if(pOGRSpatialReference!=NULL && pOGRSpatialReference_Source!=NULL)
{
isSameCoordSystem=pOGRSpatialReference_Source->IsSame(pOGRSpatialReference);
} // 创建Shape文件
OGRDataSource* pOGRDataSource=CreateShapeFile(polylinePath,pOGRSpatialReference,wkbLineString);
if(pOGRDataSource==NULL) return false; OGRLayer* pOGRLayer=pOGRDataSource->GetLayer();
if(pOGRLayer==NULL) return false; // 面转线再合并
OGRFeature* pOGRFeature_Old;
OGRGeometry* pTempGeometry=NULL;
OGRGeometry* pTempGeometryUnion=NULL; // 当前选择导出
if(false==vecFIDs.empty()&&vecFIDs.size()>)
{
for(int i=;i<vecFIDs.size();i++)
{
pOGRFeature_Old=pLayer->GetFeature(vecFIDs[i]);
pTempGeometry=ConvertPolygonToPolylineGeo(pOGRFeature_Old->GetGeometryRef());
pTempGeometry->assignSpatialReference(pOGRSpatialReference_Source);
if(false == isSameCoordSystem)pTempGeometry->transformTo(pOGRSpatialReference); if(pTempGeometryUnion==NULL) pTempGeometryUnion=pTempGeometry;
else pTempGeometryUnion=pTempGeometryUnion->Union(pTempGeometry);
}
}
else
{
if(false!=envelope.isNull()&&envelope.getMaxX()!=envelope.getMinX())
{
pLayer->SetSpatialFilterRect(envelope.getMinX(),envelope.getMinY(),envelope.getMaxX(),envelope.getMaxY());
} pOGRFeature_Old=pLayer->GetNextFeature();
while(NULL!= pOGRFeature_Old)
{
pTempGeometry=ConvertPolygonToPolylineGeo(pOGRFeature_Old->GetGeometryRef());
pTempGeometry->assignSpatialReference(pOGRSpatialReference_Source);
if(false == isSameCoordSystem)pTempGeometry->transformTo(pOGRSpatialReference); if(pTempGeometryUnion==NULL) pTempGeometryUnion=pTempGeometry;
else pTempGeometryUnion=pTempGeometryUnion->Union(pTempGeometry); pOGRFeature_Old=pLayer->GetNextFeature();
}
} // 获得Simply的距离
double distanceValue=0.0;
if(NULL==pOGRSpatialReference) // 如果为空
{
OGREnvelope pTempOGREnvelope ;
pTempGeometryUnion->getEnvelope(&pTempOGREnvelope);
if(pTempOGREnvelope.MaxX<) distanceValue=0.00000001;
else distanceValue=0.01;
}
else if(true==pOGRSpatialReference->IsProjected()) // 如果是Project的
{
distanceValue=0.01;
}
else // 如果是Geo的
{
distanceValue=0.00000001;
}
pTempGeometryUnion=pTempGeometryUnion->Simplify(distanceValue); // 在ShapeFile文件中添加数据行
OGRFeature* pOGRFeature_New;
OGRGeometry* pOGRGeometry;
OGRFeatureDefn* pOGRFeatureDefn=NULL;
pOGRFeatureDefn=pOGRLayer->GetLayerDefn(); OGRwkbGeometryType ogrGeometryType=pTempGeometryUnion->getGeometryType();
ogrGeometryType=wkbFlatten(ogrGeometryType); if(ogrGeometryType==OGRwkbGeometryType::wkbMultiLineString)
{
OGRGeometryCollection* pOGRGeometryCollectionTarget=(OGRGeometryCollection*) pTempGeometryUnion;
int geometryCount=pOGRGeometryCollectionTarget->getNumGeometries();
for(int i=;i<geometryCount;i++)
{
pOGRFeature_New=OGRFeature::CreateFeature(pOGRFeatureDefn);
pOGRGeometry=pOGRGeometryCollectionTarget->getGeometryRef(i);
pOGRFeature_New->SetGeometry(pOGRGeometry);
pOGRLayer->CreateFeature(pOGRFeature_New); OGRFeature::DestroyFeature(pOGRFeature_New);
pOGRFeature_New=NULL;
}
}
else if(ogrGeometryType==OGRwkbGeometryType::wkbLineString)
{
pOGRFeature_New=OGRFeature::CreateFeature(pOGRFeatureDefn);
pOGRGeometry=pTempGeometryUnion;
pOGRFeature_New->SetGeometry(pOGRGeometry);
pOGRLayer->CreateFeature(pOGRFeature_New); OGRFeature::DestroyFeature(pOGRFeature_New);
pOGRFeature_New=NULL; }
OGRDataSource::DestroyDataSource(pOGRDataSource); // 销毁pTargetGeometrys
OGRGeometryFactory::destroyGeometry(pTempGeometryUnion);
pTempGeometryUnion=NULL; return true;
} /*
* @brief ConvertPolygonToPolyline 线图层转换为面图层
* @param[in] tString polylinePath 转换后面图层文件路径
* @param[in] OGRLayer* pLayer 要转换的线图层文件
* @param[in] Envelope envelope 要转换数据的范围
* @param[in] vector<long> vecFIDs 选中的要素ID列表
* @param[in] pOGRSpatialReference 要转换的数据的空间参考(如果为空表示坐标系信息不变)
* @return bool 是否成功
* @author
* @date
* @note 2015年11月04日 小八创建;
*/
bool FeatureLayerOperator::ConvertPolylineToPolygonEx(tString polylinePath,OGRLayer* pLayer,Envelope envelope,vector<long> vecFIDs,OGRSpatialReference* pOGRSpatialReference)
{
// 判断
if(pLayer==NULL) return false;
if(true==polylinePath.empty()) return false; // 坐标系读取
OGRSpatialReference* pOGRSpatialReference_Source=pLayer->GetSpatialRef();
bool isSameCoordSystem = false;
if(pOGRSpatialReference == NULL)
{
pOGRSpatialReference=pOGRSpatialReference_Source;
isSameCoordSystem=true;
}
else if(pOGRSpatialReference!=NULL && pOGRSpatialReference_Source!=NULL)
{
isSameCoordSystem=pOGRSpatialReference_Source->IsSame(pOGRSpatialReference);
} // 创建Shape文件
OGRDataSource* pOGRDataSource=CreateShapeFile(polylinePath,pOGRSpatialReference,wkbPolygon);
if(pOGRDataSource==NULL) return false; OGRLayer* pOGRLayer=pOGRDataSource->GetLayer();
if(pOGRLayer==NULL) return false; // 面合并
OGRFeature* pOGRFeature_Old;
OGRGeometry* pTempGeometry=NULL;
OGRGeometry* pTempGeometryUnion=NULL; // 当前选择导出
if(false==vecFIDs.empty()&&vecFIDs.size()>)
{
for(int i=;i<vecFIDs.size();i++)
{
pOGRFeature_Old=pLayer->GetFeature(vecFIDs[i]);
pTempGeometry=pOGRFeature_Old->GetGeometryRef();
if(false == isSameCoordSystem)pTempGeometry->transformTo(pOGRSpatialReference); if(pTempGeometryUnion==NULL) pTempGeometryUnion=pTempGeometry;
else pTempGeometryUnion=pTempGeometryUnion->Union(pTempGeometry);
}
}
else
{
if(false!=envelope.isNull()&&envelope.getMaxX()!=envelope.getMinX())
{
pLayer->SetSpatialFilterRect(envelope.getMinX(),envelope.getMinY(),envelope.getMaxX(),envelope.getMaxY());
} pOGRFeature_Old=pLayer->GetNextFeature();
while(NULL!= pOGRFeature_Old)
{
pTempGeometry=pOGRFeature_Old->GetGeometryRef();
if(false == isSameCoordSystem)pTempGeometry->transformTo(pOGRSpatialReference); if(pTempGeometryUnion==NULL) pTempGeometryUnion=pTempGeometry;
else pTempGeometryUnion=pTempGeometryUnion->Union(pTempGeometry); pOGRFeature_Old=pLayer->GetNextFeature();
}
} OGRGeometry* pOGRGeometryUnion = pTempGeometryUnion->Polygonize(); // 在ShapeFile文件中添加数据行
OGRFeature* pOGRFeature_New;
OGRGeometry* pOGRGeometry;
OGRFeatureDefn* pOGRFeatureDefn=NULL;
pOGRFeatureDefn=pOGRLayer->GetLayerDefn(); OGRwkbGeometryType ogrGeometryType=pOGRGeometryUnion->getGeometryType();
ogrGeometryType=wkbFlatten(ogrGeometryType); if(ogrGeometryType==OGRwkbGeometryType::wkbGeometryCollection||ogrGeometryType==OGRwkbGeometryType::wkbMultiPolygon)
{
OGRGeometryCollection* pOGRGeometryCollectionTarget=(OGRGeometryCollection*) pOGRGeometryUnion;
int geometryCount=pOGRGeometryCollectionTarget->getNumGeometries();
for(int i=;i<geometryCount;i++)
{
pOGRFeature_New=OGRFeature::CreateFeature(pOGRFeatureDefn);
pOGRGeometry=pOGRGeometryCollectionTarget->getGeometryRef(i);
pOGRFeature_New->SetGeometry(pOGRGeometry);
pOGRLayer->CreateFeature(pOGRFeature_New); OGRFeature::DestroyFeature(pOGRFeature_New);
pOGRFeature_New=NULL;
}
}
else if(ogrGeometryType==OGRwkbGeometryType::wkbPolygon)
{
pOGRFeature_New=OGRFeature::CreateFeature(pOGRFeatureDefn);
pOGRGeometry=pOGRGeometryUnion;
pOGRFeature_New->SetGeometry(pOGRGeometry);
pOGRLayer->CreateFeature(pOGRFeature_New); OGRFeature::DestroyFeature(pOGRFeature_New);
pOGRFeature_New=NULL; }
OGRDataSource::DestroyDataSource(pOGRDataSource); // 销毁pTargetGeometrys
OGRGeometryFactory::destroyGeometry(pTempGeometryUnion);
pTempGeometryUnion=NULL; return true;
}
GDAL线面互转换(2)的更多相关文章
- GDAL线面互转换
最近因为项目需要,需做GDAL线面互转的功能,查阅部分资料完成,下面把核心部分的代码贴出来,留个记录,也欢迎大家提问题指正完善. /* * @brief ConvertPolygonToPolylin ...
- opengl中场景变换|2D与3D互转换(转)
opengl中场景变换|2D与3D互转换 我们生活在一个三维的世界——如果要观察一个物体,我们可以: 1.从不同的位置去观察它.(视图变换) 2.移动或者旋转它,当然了,如果它只是计算机里面的物体,我 ...
- PHP时间戳和日期互转换
在php中我们要把时间戳转换日期可以直接使用date函数来实现,如果要把日期转换成时间戳可以使用strtotime()函数实现,下面我来给大家举例说明. 1.php中时间转换函数 strtotime ...
- Javascript Array和String的互转换。
Array类可以如下定义: var aValues = new Array(); 如果预先知道数组的长度,可以用参数传递长度 var aValues = new Array(20); -------- ...
- Javascript Array和String的互转换
Array类可以如下定义: var aValues = new Array(); 如果预先知道数组的长度,可以用参数传递长度 var aValues = new Array(20); -------- ...
- C#二进制与字符串互转换,十六进制转换为字符串、float、int
/// <summary> /// 将 字符串 转成 二进制 “10011100000000011100011111111101” /// </summary> /// ...
- Image与Base64String的互转换
public Image Base64ToImage(string base64String) { // Convert Base64 String to byte[] byte[] imageByt ...
- Gdal 1.11.0 添加 Postgresql 9.1 sqlite3 支持
OS环境Ubuntu12.04 32bit 因为公司一个功能要用到gdal 的ogr2ogr命令转换shp数据,需要能往postgis和sqlite 中插入数据. 用gdal1.11.0的源码默认安装 ...
- c#中将IP地址转换成无符号整形数的方法与逆变换方法
我们知道 IP地址就是给每个连接在Internet上的主机分配的一个32bit地址. 按照TCP/IP协议规定,IP地址用二进制来表示,每个IP地址长32bit,比特换算成字节,就是4个字节.而c#中 ...
随机推荐
- Java爬虫——Gecco简单入门程序(根据下一页一直爬数据)
为了完成作业,所以学习了一下爬虫Gecco,这个爬虫集合了以往所有的爬虫的特点,但是官方教程中关于Gecco的教程介绍的过于简单,本篇博客是根据原博客的地址修改的,原博客中只有程序的截图,而没有给出一 ...
- 使用python+selenium对web进行自动化测试
想用python代码,对web网页进行自动化测试 web自动化测试和手动测试的区别: 手动测试:通过手动去对网页的功能进行点点点 web自动化:可以通过代码,自动对网页点点点 首先,将python+s ...
- POJ 2390
import java.util.*; public class Main { public static void main(String args[]){ double interest; Sca ...
- Spark安装过程
Precondition:jdk.Scala安装,/etc/profile文件部分内容如下: JAVA_HOME=/home/Spark/husor/jdk CLASSPATH=.:$JAVA_HOM ...
- Centos6.7安装Pycharm及升级JDK
首先到pycharm官网下载pycharm压缩包 wget https://www.jetbrains.com/pycharm/download/download-thanks.html?platfo ...
- spring boot整合RabbitMQ(Direct模式)
springboot集成RabbitMQ非常简单,如果只是简单的使用配置非常少,springboot提供了spring-boot-starter-amqp项目对消息各种支持. Direct Excha ...
- Python -- Gui编程 -- Tkinter的使用 -- 基本控件
1.按钮 tkBtton.py import tkinter root = tkinter.Tk() btn1 = tkinter.Button(root, anchor=tkinter.E,\ te ...
- wap尝试调取app(网易新闻为例)
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- redis-springboot-redistemplate更改序列化方式
redisTemplate 默认的序列化方式为 jdkSerializeable, StringRedisTemplate的默认序列化方式为StringRedisSerializer 可以通过手动配置 ...
- 从入门到不放弃系列之Koa2
一.Koa2入门 本来是想Express入门的,但是既然都是要学,干嘛不学最新的呢? 其实我想说,我本来只是想学个小程序开发,现在已经陆陆续续开了好多坑了.. 本文参考廖雪峰教程 二.Async 最新 ...