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(0);
    if(pOGRLayer==NULL) return false;

    // 面转线再合并
    OGRFeature* pOGRFeature_Old;
    OGRGeometry* pTempGeometry=NULL;
    OGRGeometry* pTempGeometryUnion=NULL;

    // 当前选择导出
    if(false==vecFIDs.empty()&&vecFIDs.size()>0)
    {
        for(int i=0;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<180) 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=0;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(0);
    if(pOGRLayer==NULL) return false;

    // 面合并
    OGRFeature* pOGRFeature_Old;
    OGRGeometry* pTempGeometry=NULL;
    OGRGeometry* pTempGeometryUnion=NULL;

    // 当前选择导出
    if(false==vecFIDs.empty()&&vecFIDs.size()>0)
    {
        for(int i=0;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=0;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;
}