CloudCompare——Delaunay三角剖分

在这里插入图片描述

本文由CSDN点云侠原创,CloudCompare——Delaunay三角剖分,爬虫自重。如果你不是在点云侠的博客中看到该文章,那么此处便是不要脸的爬虫与GPT生成的文章。

一、常用方法

1.分而治之算法

  分而治之算法最早是1975年由Shamos和Hoey提出的,Lewis和Rovinson在1978年利用该方法进行了三角网的剖分,随后Lee和Schachter、Dwyer等对他们的算法进行了改进和优化。
  分而治之算法的思路是将复杂问题简单化,首先将数据点分割成包含少量点的子集,如一个子集中包括三个、四个点,然后每个子集进行三角剖分,并用局部优化算法(LOP)进行优化,保证三角剖分为Delaunay三角网,最后是对每个子集的三角剖分进行合并,形成最终的整体三角网。不同的实现方法可有不同的点集划分方法、子网生成方法及合并方法。Lee和Schachter提出了分而治之的经典算法,其基本步骤如下:

  1. 把点集 V V V以横坐标为主,纵坐标为辅按升序排序,即数据点之间满足条件: ( x i , y i ) < = ( x i + 1 , y i + 1 ) (x_i,y_i)<=(x_{i+1},y_{i+1}) (xi,yi)<=(xi+1,yi+1)
  2. 把点集 V V V分为近似相等的两个子集 V L V_L VL V R V_R VR
  3. V L V_L VL V R V_R VR中生成三角网;
  4. 使用局部优化算法,使之成为Delaunay三角网;
  5. 找出连接 V L V_L VL V R V_R VR中两个凸壳的底线和顶线;
  6. 由底线至顶线合并 V L V_L VL V R V_R VR中两个三角网;
  7. 结束。

  该算法的时间复杂度为 ( N l o g N ) (NlogN) (NlogN),其中 N N N为数据点数,因此该算法执行时间效率高,但在 V L V_L VL V R V_R VR中生成三角网时,采用了递归运算,需占用较大的系统内存。

2.三角网生长算法

  三角网生长算法最早由Green和Sibson在V图中实现,稍后Brassel和Teif也发表了类似的算法。分为收缩生长算法和扩张生长算法,收缩生长算法是先生成凸壳,并以凸壳为源头,逐步缩小以形成整个三角网。收缩生长算法与数据点的分布密度有关,实际情况比较复杂。扩张生成算法与收缩算法过程刚好相反,算法是从一个三角形开始向外层扩展,最终形成覆盖整个区域的三角网,是常用的三角网生长算法。
  扩张生长算法的思想是先找出点集中最近两点连接成一条边,按照Delaunay三角网构成原则找出第三点,将这三点连接成初始三角形,再以该三角形的每一条边为基线扩展连接相邻离散点,组成新三角形,初始三角形的三边处理之后,继续以新三角形的边连接其他离散点,直到所有的离散点均包含在三角网中。该算法的基本步骤:

  1. 在区域内所有离散点任取一点作为初始三角形的第一个顶点;
  2. 找距离第一个顶点最近的点,作为初始三角形的第二个顶点,两点相连生成初始基线;
  3. 找出距离初始基线中点最近且不和已有两点在一条直线上的点作为初始三角形的第三个顶点,三点连接成三角形,即为初始三角形;
  4. 三角形的两条新边作为新的基线;
  5. 重复步骤3、4直到所有基线处理完毕。

  该算法的时间复杂度在一般情况下为 O ( N 3 / 2 ) O(N^{3/2}) O(N3/2),最坏情况下达到 O ( N 2 ) O(N^2) O(N2)。因此该算法的时间效率低,优点是占用内存空间较小。在该算法的基础上有许多的改进算法,其改进点主要集中在“第三点”的搜索和三角形重复的判断上。

3.逐点插入算法

  逐点插入算法的思想最早由Lawson(1977)提出,随后Lee和Schachter(1980)、Bowyer(1981)、Watson(1981)、Tsai(1993)等先后对其进行改进,是一种动态的构网方式。
  逐点插入算法的思路是将未处理的点加入到己经存在的Delaunay三角网中,然后判断点所在的三角形,如果点在三角形内,连接点与三角形的各顶点,利用三角剖分准则,找出影响区域,删除影响区域中的边,对影响区域利用三角剖分准则重新构网,直到区域中所有的点都加入到三角网中。该算法的基本步骤如下:

  1. 首先定义一个包含所有数据点的初始多边形,即包含所有数据点的凸闭包;
  2. 在凸闭包中构建初始Delaunay三角网;
  3. 从内部数据点中取出一点加入到三角网中;
  4. 搜寻包含该点的三角形,将该点与此三角形的三个顶点相连,形成三个三角形;
  5. 利用Delaunay三角网的剖分准则,找出影响区域,从里到外更新所有生成的三角形;
  6. 重复步骤3、4、5直到处理完点集中所有点。

  该算法一般情况下时间复杂度为 O ( N 3 / 2 ) O(N^{3/2}) O(N3/2),最坏情况下达到 O ( N 2 ) O(N^2) O(N2)。该算法的时间效率低,优点是算法思路较简单,占用内存空间较小,空间性能较好。
通常,分而治之算法的效率是最高的,但是当数据量比较大,分块比较多时,块之间的合并算法复杂,且花时也多,并且递归执行,内存消耗大;三角网生长算法和逐点插入法的执行效率差不多,但是在构网的过程中逐点插入法能保持较好的空间特性。

二、软件操作

1.Delaunay 2.5D(XY plane)

  通过菜单栏的'Edit > Mesh > Delaunay 2.5D (XY plane)'找到该功能。
在这里插入图片描述
其实现流程为

  1. 点云简单地在XY平面上进行2D投影。
  2. 对相应的2D点进行三角剖分。
  3. 将网格结构应用于3D点。

  默认情况下,2D Delaunay三角剖分是在点云凸包上完成。因此,CloudCompare将要求用户指定三角形边的最大长度(如下图所示)。这允许移除不一定有意义的最大三角形(通常在边界上)。如果该值保持为零,则会保留Delaunay三角剖分输出的所有三角形。
在这里插入图片描述
  这适用于表面较为平坦和具有适当定向(即Z维度为垂直的)的点云。如果点云没有正确定向但在某个方向上仍然相当平坦,则应考虑使用此方法的另一个版本'Mesh> Delaunay 2.5D(best fitting plane)
以下是某道路机载点云的Delaunay 2.5D模型生成结果
在这里插入图片描述
  CloudCompare默认以面片的形式显示mesh网格模型,因此放大之后也看不到所谓的“网格”
在这里插入图片描述
  在属性选项中,选中以线的形式显示格网,即可清楚看到Delaunay三角网的网形状,操作如下:
在这里插入图片描述
Delaunay三角网的构网形状如下:
在这里插入图片描述
  此外,还有以点的形式显示模型。
在这里插入图片描述
  以点的形式显示结果如下,可以看出,路面起伏特征及模型的纹理细节更加清晰。所以Delaunay 算法也是可以用于地面滤波与变化监测的???

在这里插入图片描述

2.Delaunay 2.5D(XY best fitting plane)

  该方法与 Delaunay 2.5D (XY plane)功能基本相同。
其实现流程为

  1. 点云投影在最小二乘拟合的最佳平面。
  2. 对相应的2D点进行三角剖分。
  3. 将网格结构应用于3D点。

  其他操作与Delaunay 2.5D相同。
在这里插入图片描述

三、相关代码

  • 10
    点赞
  • 47
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
点云Delaunay网格化是一种常见的点云处理方法,通过将点云数据转化为Delaunay三角网格来进行分析和处理。Delaunay三角化的概念最早由Delaunay于1934年提出,其基本原理是通过点云数据中的一组点,构造出最优三角剖分的方法。 在点云Delaunay网格化过程中,首先需要针对点云数据进行三维凸包的计算,将点云包围在一个凸多面体内。凸包计算可以使用相关算法,例如QuickHull、Jarvis等算法。凸包计算完成后,将凸多面体的外表面取出,并去除冗余点(例如凸多面体的顶点),得到一组简化的点数据。 接下来,通过在简化的点数据上进行Delaunay三角网格的计算。Delaunay三角剖分是一种无重叠的、无长细比的三角剖分方法,使得任意两个三角形的外接圆不包含其他点。常见的Delaunay三角剖分算法有Bowyer-Watson算法、法向量算法等。Delaunay网格化可以使用这些算法构建三角网格。 最后,通过对Delaunay三角网格进行加密和平滑,可以获得更加精细和平滑的点云网格。加密可以通过插入新的点进行,使得Delaunay三角剖分的三角形边界更加接近点云表面。平滑可以通过迭代优化的方法,例如Laplacian平滑、Moving Least Squares等方法,使得点云网格的形状更加平滑。 点云Delaunay网格化在计算机图形学、计算机视觉、三维重建等领域有着广泛的应用。通过将点云数据转化为Delaunay三角网格,可以方便地进行点云的分析、可视化、特征提取等操作,为后续处理提供了基础。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

点云侠

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值