在上一个文章中介绍了线转化为面和面转化为线,其主要的实现思路就是把面中的点取出来构成线,把线中的点取出来构成面,实际上就是一个硬拷贝,无奈客户的实际需求并非如此,客户想要线转面的时候几条相交线构成面,面转线的时候相同的线去除,所以又重新对功能进行了调整。

关于线构面,这个过程中也是挺曲折的,在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)的更多相关文章

  1. GDAL线面互转换

    最近因为项目需要,需做GDAL线面互转的功能,查阅部分资料完成,下面把核心部分的代码贴出来,留个记录,也欢迎大家提问题指正完善. /* * @brief ConvertPolygonToPolylin ...

  2. opengl中场景变换|2D与3D互转换(转)

    opengl中场景变换|2D与3D互转换 我们生活在一个三维的世界——如果要观察一个物体,我们可以: 1.从不同的位置去观察它.(视图变换) 2.移动或者旋转它,当然了,如果它只是计算机里面的物体,我 ...

  3. PHP时间戳和日期互转换

    在php中我们要把时间戳转换日期可以直接使用date函数来实现,如果要把日期转换成时间戳可以使用strtotime()函数实现,下面我来给大家举例说明. 1.php中时间转换函数 strtotime ...

  4. Javascript Array和String的互转换。

    Array类可以如下定义: var aValues = new Array(); 如果预先知道数组的长度,可以用参数传递长度 var aValues = new Array(20); -------- ...

  5. Javascript Array和String的互转换

    Array类可以如下定义: var aValues = new Array(); 如果预先知道数组的长度,可以用参数传递长度 var aValues = new Array(20); -------- ...

  6. C#二进制与字符串互转换,十六进制转换为字符串、float、int

    /// <summary>     /// 将 字符串 转成 二进制 “10011100000000011100011111111101” /// </summary> /// ...

  7. Image与Base64String的互转换

    public Image Base64ToImage(string base64String) { // Convert Base64 String to byte[] byte[] imageByt ...

  8. Gdal 1.11.0 添加 Postgresql 9.1 sqlite3 支持

    OS环境Ubuntu12.04 32bit 因为公司一个功能要用到gdal 的ogr2ogr命令转换shp数据,需要能往postgis和sqlite 中插入数据. 用gdal1.11.0的源码默认安装 ...

  9. c#中将IP地址转换成无符号整形数的方法与逆变换方法

    我们知道 IP地址就是给每个连接在Internet上的主机分配的一个32bit地址. 按照TCP/IP协议规定,IP地址用二进制来表示,每个IP地址长32bit,比特换算成字节,就是4个字节.而c#中 ...

随机推荐

  1. Java爬虫——Gecco简单入门程序(根据下一页一直爬数据)

    为了完成作业,所以学习了一下爬虫Gecco,这个爬虫集合了以往所有的爬虫的特点,但是官方教程中关于Gecco的教程介绍的过于简单,本篇博客是根据原博客的地址修改的,原博客中只有程序的截图,而没有给出一 ...

  2. 使用python+selenium对web进行自动化测试

    想用python代码,对web网页进行自动化测试 web自动化测试和手动测试的区别: 手动测试:通过手动去对网页的功能进行点点点 web自动化:可以通过代码,自动对网页点点点 首先,将python+s ...

  3. POJ 2390

    import java.util.*; public class Main { public static void main(String args[]){ double interest; Sca ...

  4. Spark安装过程

    Precondition:jdk.Scala安装,/etc/profile文件部分内容如下: JAVA_HOME=/home/Spark/husor/jdk CLASSPATH=.:$JAVA_HOM ...

  5. Centos6.7安装Pycharm及升级JDK

    首先到pycharm官网下载pycharm压缩包 wget https://www.jetbrains.com/pycharm/download/download-thanks.html?platfo ...

  6. spring boot整合RabbitMQ(Direct模式)

    springboot集成RabbitMQ非常简单,如果只是简单的使用配置非常少,springboot提供了spring-boot-starter-amqp项目对消息各种支持. Direct Excha ...

  7. Python -- Gui编程 -- Tkinter的使用 -- 基本控件

    1.按钮 tkBtton.py import tkinter root = tkinter.Tk() btn1 = tkinter.Button(root, anchor=tkinter.E,\ te ...

  8. wap尝试调取app(网易新闻为例)

    <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...

  9. redis-springboot-redistemplate更改序列化方式

    redisTemplate 默认的序列化方式为 jdkSerializeable, StringRedisTemplate的默认序列化方式为StringRedisSerializer 可以通过手动配置 ...

  10. 从入门到不放弃系列之Koa2

    一.Koa2入门 本来是想Express入门的,但是既然都是要学,干嘛不学最新的呢? 其实我想说,我本来只是想学个小程序开发,现在已经陆陆续续开了好多坑了.. 本文参考廖雪峰教程 二.Async 最新 ...