(三)建筑物多边形化简系列——去除冗余点
发布日期:2022-02-10 11:37:01 浏览次数:50 分类:技术文章

本文共 3088 字,大约阅读时间需要 10 分钟。

制图综合和建筑物数据处理等都涉及到建筑物多边形的化简。制图综合中,由于比例尺的变小,建筑物在小比例尺地图上所占面积变小,这意味着建筑物图形的形状精度也有一定的损失,为了更好地表示原有建筑物的特征(面积、图形和方向),需要对建筑物多边形进行化简。另外,从遥感影像提取的建筑物矢量数据需要经过图形化简等一系列操作才能作为可使用的数据。

首先明确,矢量数据结构中,建筑物多边形是由一些列环(circle)组成,环则由一串首尾相同的点构造。这意味着对建筑物多边形化简其实是对这些环(线状)处理,更具体地说,是对环的各个点进行处理(删除、移位)。

对建筑物多边形的化简到最后得到能较好拟合建筑物的成果多边形,这中间需要一系列操作。首先要进行去除冗余点的操作。

1)冗余点

对于建筑物多边形边界上的邻近三点,若它们所组成两条直线的夹角为平角(或零角),则该三点共线,且位于中间位置那个点的取舍对建筑物多边形的 几 何形状没有影响,因此该三点中位于中间位置那个点是冗余的。

需要注意的是由于数据采集误差、制图不规范等原因,使得地图上建筑物多边形边界上出现冗余点的情况很少,出于方法的合理性考虑,选取5°作为缓冲角度,对冗余点的定义做出以下修正定义:对于建筑物多边形边界上的邻近三点,若它们所组成两条直线的夹角与平角(或零角)相差的角度小于5°时,称该三点中位于中间位置的那个点为冗余点。

2)基本思想

取环上邻近三点作为一个单元,分析夹角:若夹角小于阈值,则删除中间点,加入下一个点构成新的分析单元;若夹角大于阈值,则保留中间点,继续遍历。

3)代码实现

// 去除冗余点,zf,0717void CGeoPolygon::RemoveRedundant(void){    for(int i = 0;i
temp; //用于存放冗余点的下标 bool odinary = true; //设定状态,以确定三点单元的首点 int count = 0; //记录删除的点数,一旦有正常情况恢复为0 int ptsNum = circles[i]->pts.size(); for(int n = 0;n
pts[n]; middle = circles[i]->pts[n+1]; last = circles[i]->pts[n+2]; } else //不正常,存在冗余点的情况 { first = circles[i]->pts[n-count]; //存在冗余点,删除后判断冗余点之后的点是否为冗余点还需要之前三点单元的首点 middle = circles[i]->pts[n+1]; last = circles[i]->pts[n+2]; } double angle = calAngle(first,middle,last); //计算角度 if(angle >= 10&&angle<=170) //角度阈值设为10°(可根据需要调整),当角度大于阈值,情况正常,不需要去除点 { odinary = true; count = 0; //记录删除的点数,一旦有正常情况恢复为0 } else // 存在冗余点 { temp.push_back(n+1); //记录冗余点在circles[i]->pts中下标 odinary = false; //情况不正常 count++; //连续不正常时都要记录 } }
vector
tempPoints; //定义临时点集 if(tempPoints.size()!=0) vector
().swap(tempPoints); //防止点集不为空 for(int j = 0;j
pts[j]); } tempPoints.swap(circles[i]->pts); //交换,得到去除冗余点之后的环circles[i]->pts }}

  其中calAngle函数是计算三点的夹角,具体代码如下:

// 计算三点之间的角度,zf,0717,要#include
double CGeoPolygon::calAngle(CMyPoint* p1, CMyPoint* p2, CMyPoint* p3){ double x1 = p1->Getx(); double y1 = p1->Gety(); double x2 = p2->Getx(); double y2 = p2->Gety(); double x3 = p3->Getx(); double y3 = p3->Gety(); double angle,p1p2,p2p3,p1p3; p1p2 = sqrt(((x1 - x2) * (x1 - x2) + (y1 -y2) * (y1 -y2))); p2p3 = sqrt(((x3 - x2) * (x3 - x2) + (y3 -y2) * (y3 -y2))); p1p3 = sqrt(((x1 - x3) * (x1 - x3) + (y1 -y3) * (y1 -y3))); angle = acos((p1p2 * p1p2 + p2p3 * p2p3 - p1p3 * p1p3) / (2 * p1p2 * p2p3)) * 180.0 / 3.1415926; if(angle<0) angle = angle*(-1); return angle;}

  isContained函数判断某个数是不是在数组中,见下

// 判断某个数是不是在数组中,zf,0717bool CGeoPolygon::isContained(vector
temp, int i){ bool answer = false; for(int j = 0;j

4)效果展示

去除冗余点之前

 

去除冗余点之后

5)小结

去除冗余点是建筑物化简的第一步,虽然看起来前后差别不大,但是细细观察还是有一定区别,更近一步的化简将在下一篇博客中讲述。

转载地址:https://blog.csdn.net/Fan_z_0802/article/details/81090580 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!

上一篇:计算机知识记录本
下一篇:(二)建筑物多边形化简系列——多边形点数化简

发表评论

最新留言

留言是一种美德,欢迎回访!
[***.207.175.100]2024年03月19日 23时23分35秒

关于作者

    喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!

推荐文章

qpython3安装lxml_在python的lxml中使用xml目录? 2019-04-21
java 幂取模_快速幂取模算法 2019-04-21
java build path jre_java-如何在安装了jre 7后为Jre 6设置路径? 2019-04-21
java上传下载源码_javaweb简单实现文件上传与下载源代码 2019-04-21
java socket udp 广播_1.Java 的屏幕广播(基于UDP),2.多线程下载器 2019-04-21
java控制热敏打印机的例子.rar_stm32控制热敏打印机 2019-04-21
java clone equals_(原)java中对象复制、==、equals 2019-04-21
java滚动字幕实训报告_Java实习报告 (7000字).doc 2019-04-21
php7 memcached.exe,PHP7 下安装 memcache 和 memcached 扩展 2019-04-21
计算机二级java技巧,计算机二级报java难考吗 2019-04-21
php foreach 数据库,php – 使用foreach将数据库检索的数据排列在HTML表中 2019-04-21
拉格朗日matlab编程例题,Matlab习题讲解.doc 2019-04-21
case是不是php语言关键字,PHP语言 switch 的一个注意点 2019-04-21
linux php mkdir失败,linux – mkdir错误:参数无效 2019-04-21
config.php渗透,phpMyAdmin 渗透利用总结 2019-04-21
java list 合并 重复的数据_Java ArrayList合并并删除重复数据3种方法 2019-04-21
android volley 上传图片 和参数,android - 使用android中的volley将图像上传到multipart中的服务器 - 堆栈内存溢出... 2019-04-21
android开发的取消清空按钮,Android开发实现带清空按钮的EditText示例 2019-04-21
android gp服务,ArcGIS Runtime SDK for Android开发之调用GP服务(异步调用) 2019-04-21
mysql整体会滚_滚mysql 2019-04-21