基于SuperMap iObjects C++之地形凹凸点提取

最新推荐文章于 2025-05-09 15:16:33 发布

supermapsupport

最新推荐文章于 2025-05-09 15:16:33 发布

阅读量578

点赞数

CC 4.0 BY-SA版权

分类专栏：组件GIS 文章标签： c++ 开发语言

本文为博主原创文章，欢迎转载，SuperMap技术问答社区 http://ask.supermap.com

本文链接：https://blog.youkuaiyun.com/supermapsupport/article/details/127749825

组件GIS 专栏收录该内容

177 篇文章

订阅专栏

在实际项目中，有时候会遇到这样的需求，根据DEM地形数据，我们需要根据某两点的连线，找到经过的所有山顶点和凹陷点。

山顶点指哪些在特定邻域分析范围内，该点都比周围点高的区域。凹陷点指哪些在特定邻域分析范围内，该点都比周围点低的区域区域。DEM数据提供了丰富的地形信息，基于DEM数据如何正确有效的提取这些地形特征信息，在数字地形分析中具有重要意义。

对于客户所遇到的这种需求，下面跟着我一起来看看在SuperMap iObjects c++中，是如何实现吧。

首先，我们还是需要明确的一点，前期的数据处理，在这里包括dem数据的入库，在程序开发中不是每一次都要进行的操作，也不是具体的业务需求，所以所以我们通常都是使用SuperMap提供的专门用于前期处理数据的iDesktop软件工作进行数据的入库，处理。先在桌面软件中，将数据制作完成，地图配图，保存为工作空间文件（*.smwu）,组件开发直接打开数据即可，程序开发只实现业务需求。

数据准备

这里，我们首先需要将拿到的DEM地形数据先导入到SuperMap的数据源（*.udb或者*.udbx）中，不管是tif还是其他格式的数据都一样.

注意：导入的时候，数据类型需要修改为“栅格”，不能按照默认的影像方式导入。

关于如何SuperMap iObjects c++环境的部署，以及如何打开数据不清楚的小伙伴，这里可以参考之前的博客：

环境配置（windows环境VS + Qt）：SuperMap iObjects for C++ 入门详解（VS + Qt）_supermapsupport的博客-优快云博客

环境配置（windows+MFC）：SuperMap iObjects C++之MFC快速入门_supermapsupport的博客-优快云博客_sumpermap c++配置

环境配置（Qt Creator）：https://blog.youkuaiyun.com/supermapsupport/article/details/52609945

数据源管理：SuperMap iObjects C++之数据源管理_supermapsupport的博客-优快云博客_supermap怎么打开udb数据源

好了，程序也建好了，数据也准备好了，下面就一起来看看功能是如何实现的吧：

、计算山顶点（凸点）：

第一步：

拿到对应的栅格数据，在SuperMap iObjects c++中，栅格数据的对象类型为UGDatasetRaster ；

第二步：确定需要分析的目标线，我们需要分析的就是从起点到终点的所有途径点，所以我们需要把这个起始点确定好。

第三步：确定一个分析系数：想要结果的精细度是更精细还是略微粗糙，越精细，分析的点越多，得到的结果越密集，反之，这里，我们把这个分析系数命名为interval；

完整实现步骤：

传入参数：

demDataset：栅格数据集；

Geoline：待分析的路线；

Interval：分析系数；

实现代码参考：

UGPoint2Ds* MMapControl::FindHeightPoint(UGDatasetRaster *demDataset, UGGeoLine *geoline, UGdouble interval)

{

    try {

        UGPoint2Ds *heightPoints=new UGPoint2Ds();

        if(interval==0.0)

        {

            interval=demDataset->GetBounds().Width()/demDataset->GetWidth();

        }

        UGdouble num=geoline->GetLength()/interval;

        for (int i=1;i<num-1;i++) {



            UGdouble distanceBack = num * (i-1);

            UGdouble distanceNow = num * i;

            UGdouble distanceNext = num * (i+1);



            UGPoint2D point2dBack = geoline->FindPointOnLineByDistance(distanceBack);

            UGPoint2D point2dNow = geoline->FindPointOnLineByDistance(distanceNow);

            UGPoint2D point2dNext = geoline->FindPointOnLineByDistance(distanceNext);



            UGPoint pointBack;

            demDataset->XYToGrid(point2dBack,pointBack);

            UGPoint pointNow;

            demDataset->XYToGrid(point2dNow,pointNow);

            UGPoint pointNext;

            demDataset->XYToGrid(point2dNext,pointNext);



            UGdouble heightBack = demDataset->GetValue(pointBack.x, pointBack.y);

            UGdouble heightNow = demDataset->GetValue(pointNow.x, pointNow.y);

            UGdouble heightNext = demDataset->GetValue(pointNext.x, pointNext.y);



            double BN = heightBack - heightNow;

            double NN = heightNow - heightNext;



            if (BN * NN < 0)

            {

                if (BN < 0)

                {

                    heightPoints->Add(point2dNow);

                }

            }

        }

        return heightPoints;

    } catch (double d) {



    }

}

关键代码分析：

1、UGdouble num=geoline->GetLength()/interval;

根据分析系数，求出我们将在线上等距离提取出多少个点进行分析，系数interval越小，分析精细越高，反之。

2、根据在线上按照等距离查找，找到相连的三个点坐标；

UGPoint2D point2dBack = geoline->FindPointOnLineByDistance(distanceBack);

UGPoint2D point2dNow = geoline->FindPointOnLineByDistance(distanceNow);

UGPoint2D point2dNext = geoline->FindPointOnLineByDistance(distanceNext);

3、根据demDataset->XYToGrid拿到对应点位在栅格数据集中所处的行列号，结果UGPoint 的x就是所处的行号，y就是所处的列号；

UGPoint pointBack;

demDataset->XYToGrid(point2dBack,pointBack);

UGPoint pointNow;

demDataset->XYToGrid(point2dNow,pointNow);