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#中 ...
随机推荐
- 布隆过滤器redis缓存
Bloom Filter布隆过滤器算法背景如果想判断一个元素是不是在一个集合里,一般想到的是将集合中所有元素保存起来,然后通过比较确定.链表.树.散列表(又叫哈希表,Hash table)等等数据结构 ...
- POJ 2376
#include<iostream>//by chengdacaizi. #include<stdio.h> #define MAXN 25005 using namespac ...
- WinForm—控制文本框只能输入整数（包括小数）
//控制控件输入的是正数(包括小数) private void txt_Weight_KeyPress(object sender, KeyPressEventArgs e) { //第一步:判断输入 ...
- 用Lingo求解线性规划问题
第一步:输入目标条件和约束条件.每行以分号隔开.然后点击工具栏上的Solve按钮,或Lingo菜单下的Solve子菜单. 第二步:检查report中的结果. 默认情况下,Lingo不进行灵敏度分析. ...
- 【转】linux下杀死进程（kill）的N种方法
转载一篇,最原始的出处已不可考,望见谅! 常规篇: 首先,用ps查看进程,方法如下: $ ps -ef ……smx 1822 1 0 11:38 ? 00:00:4 ...
- Java之BigDecimal
转载请注明源出处:http://www.cnblogs.com/lighten/p/6963836.html 1.前言 编程人员都应该知道计算机计算浮点数是不精确的,结果是近似数值,当然具体值还是和计 ...
- PL/SQL程序设计
1 PL/SQL简介 1 什么是PL/SQL? PL/SQL是 Procedure Language & Structured Query Language 的缩写.PL/SQL是对SQL语言 ...
- 【转】Hadoop vs Spark性能对比
原文地址:http://www.cnblogs.com/jerrylead/archive/2012/08/13/2636149.html 基于Spark-0.4和Hadoop-0.20.2 1. K ...
- CentOS Maven 删除 *.lastUpdated文件
find . -name "*.lastUpdated" -exec rm -rf {} \;
- vue里面computed的运用理解
computed 计算属性,是用来声明式的描述一个值依赖了其它的值,当你在模板里把数据绑定到一个计算属性上时,Vue 会在其依赖的任何值导致该计算属性改变时更新 DOM.这个功能非常强大,它可以让你的 ...