基于体半径函数的网格分割算法
- 格式:pdf
- 大小:238.54 KB
- 文档页数:3
点云常⽤分割⽅法点云分割 点云分割可谓点云处理的精髓,也是三维图像相对⼆维图像最⼤优势的体现。
点云分割的⽬的是提取点云中的不同物体,从⽽实现分⽽治之,突出重点,单独处理的⽬的。
⽽在现实点云数据中,往往对场景中的物体有⼀定先验知识。
⽐如:桌⾯墙⾯多半是⼤平⾯,桌上的罐⼦应该是圆柱体,长⽅体的盒⼦可能是⽜奶盒......对于复杂场景中的物体,其⼏何外形可以归结于简单的⼏何形状。
这为分割带来了巨⼤的便利,因为简单⼏何形状是可以⽤⽅程来描述的,或者说,可以⽤有限的参数来描述复杂的物体。
⽽⽅程则代表的物体的拓扑抽象。
于是,RanSaC算法可以很好的将此类物体分割出来。
1、RanSaC算法 RanSaC算法(随机采样⼀致)原本是⽤于数据处理的⼀种经典算法,其作⽤是在⼤量噪声情况下,提取物体中特定的成分。
下图是对RanSaC算法效果的说明。
图中有⼀些点显然是满⾜某条直线的,另外有⼀团点是纯噪声。
⽬的是在⼤量噪声的情况下找到直线⽅程,此时噪声数据量是直线的3倍。
如果⽤最⼩⼆乘法是⽆法得到这样的效果的,直线⼤约会在图中直线偏上⼀点。
关于随机采样⼀致性算法的原理参考博客: 这个算法就是从⼀堆数据⾥挑出⾃⼰最⼼仪的数据。
所谓⼼仪当然是有个标准(⽬标的形式:满⾜直线⽅程?满⾜圆⽅程?以及能容忍的误差e)。
平⾯中确定⼀条直线需要2点,确定⼀个圆则需要3点。
1. 平⾯中随机找两个点,拟合⼀条直线,并计算在容忍误差e中有多少点满⾜这条直线2. 重新随机选两点,拟合直线,看看这条直线是不是能容忍更多的点,如果是则记此直线为结果3. 重复步骤⼆(循环迭代)4. 迭代结束,记录当前结果算法的优点是噪声可以分布的任意⼴,噪声可以远⼤于模型信息。
这个算法有两个缺点,第⼀,必须先指定⼀个合适的容忍误差e。
第⼆,必须指定迭代次数作为收敛条件。
综合以上特性,本算法⾮常适合从杂乱点云中检测某些具有特殊外形的物体。
PCL中基于RanSaC的点云分割⽅法:1. //创建⼀个模型参数对象,⽤于记录结果2. pcl::ModelCoefficients::Ptr coefficients (new pcl::ModelCoefficients);3. //inliers表⽰误差能容忍的点记录的是点云的序号4. pcl::PointIndices::Ptr inliers (new pcl::PointIndices);5. // 创建⼀个分割器6. pcl::SACSegmentation<pcl::PointXYZ> seg;7. // Optional8. seg.setOptimizeCoefficients (true);9. // Mandatory-设置⽬标⼏何形状0. seg.setModelType (pcl::SACMODEL_PLANE);1. //分割⽅法:随机采样法2. seg.setMethodType (pcl::SAC_RANSAC);3. //设置误差容忍范围4. seg.setDistanceThreshold (0.01);5. //输⼊点云6. seg.setInputCloud (cloud);7. //分割点云8. seg.segment (*inliers, *coefficients);除了平⾯以外,PCL⼏乎⽀持所有的⼏何形状。
Netgen网格划分程序采用的推进前线法算法实现步骤Netgen网格划分程序采用的推进前线法advancing front method算法方法是将具体规则与规则应用程序代码分离。
算法必须检查存储在数据结构中的规则。
因此,代码复杂度与规则的数量无关。
该算法复杂,但定义良好,至少在理论上可以实现故障保护。
特别是在3D中,具体规则的选择是基于启发式的,并将其放入一个易于维护的规则描述数据库中。
我们按以下步骤进行处理。
首先,我们描述了整个算法,接着是抽象规则描述和规则应用算法。
最后,我们集中讨论了曲面和体网格生成的扩展。
基本思路首先:建立模型,构造立体几何提供常用的基本单元:立方体圆柱椭圆球圆球长方体等等复杂的模型通过对基本单元的bool操作来建立。
(对于复杂几何模型,我们需要求解三元非线性方程组,来计算用于寻找边缘的起始点)一、求解相交点图相交点示意图求解相交点的方法采用基于几何测试的二分算法,其示意图如下图所示。
图二分法求交点示意图基于几何测试的二分算法的计算流程1 建立一个立方体,将整个几何模型包含在该立方体内2 沿立方体的八条边界的中点分别建立三个相互正交的平面,三个平面将立方体划分成八个大小相等的小立方体。
3 判断交点在哪一个小立方体内4 在该小立方体内重复第二步,第三步。
5 循环的结束条件为立方体的边长小于给定的精确度值6 求出交点的坐标位置图极坐标下中的相交点图两个相贯的圆柱二、计算边缘:使用求解非线性方程组的同伦方法进行计算跟踪曲线。
1 该方法的原理是:从曲线上一个给定的点X0开始,通过预估-矫正的方法,沿着曲线小步推进,直到到达曲线的终点。
2 方法步骤:首先从选定的初始点X0出发,沿着曲线的切向量方向其次确定相邻两点之间的递归关系假设已知当前点Xk,计算曲线在改点的单位切向量Tk,预估下一点的位置,对该点的位置进行修正,采用牛顿下山法,对于平滑曲线,有这样一个关系式,它可以用于适应性步长控制。
基于谱聚类的三维⽹格分割算法(SpectralClustering) 谱聚类(Spectral Clustering)是⼀种⼴泛使⽤的数据聚类算法,[Liu et al. 2004]基于谱聚类算法⾸次提出了⼀种三维⽹格分割⽅法。
该⽅法⾸先构建⼀个相似矩阵⽤于记录⽹格上相邻⾯⽚之间的差异性,然后计算相似矩阵的前k个特征向量,这些特征向量将⽹格⾯⽚映射到k维谱空间的单位球上,最后使⽤K-means⽅法对谱空间中的数据点进⾏聚类。
具体算法过程如下:⼀.相似矩阵 ⽹格分割以⾯⽚为基本单元,为了能使算法沿着⼏何模型的凹形区域进⾏分割,⽹格相邻⾯⽚之间的距离采⽤[Katz et al. 2003]中提到的⽅法,具体形式在“”中有所解释,距离由测地距离Geod_Dist和⾓度距离Ang_Dist两部分组成,如下所⽰: 上式中有两个重要的参数δ和η,参数δ通常取值范围为[0.01, 0.05],其⽤于控制测地距离和⾓度距离之间的权重⽐例,参数η通常取值范围为[0.1, 0.2],其使得分割边界更倾向于凹形区域。
计算完相邻⾯⽚之间的距离后,相似矩阵中对应位置的值由距离的⾼斯函数得到:其中:⼆.⽹格谱聚类 谱聚类⽅法在本质上都是类似的,都是利⽤相似矩阵的特征向量将原始空间中的数据映射到谱空间,并在谱空间中进⾏聚类。
⽹格上的谱聚类⽅法如下: 1 由上述定义计算相似矩阵W; 2 计算归⼀化矩阵N:N = D-1/2WD-1/2; 3 计算矩阵N的前k个最⼤特征向量e1, e2, … , e k,以这k个特征向量为列组成矩阵U = [e1, e2, … , e k]; 4 将矩阵U的每⼀⾏单位化后得到矩阵Ȗ; 5 提取出k个初始聚类中⼼⽤于K-means聚类,该过程先计算矩阵Q = ȖȖT,然后查找矩阵Q中的最⼩元素Q rs,那么r和s点就是两个距离最远的点,然后继续查找后续点; 6 以Ȗ的⾏向量为数据样本进⾏K-means聚类。
CosmosWorks网格划分、求解器、提示与技巧一、网格划分策略网格划分,更精确地说应该称为离散化,就是将一数学模型转化为有限元模型以准备求解。
作为一种有限元方法,网格划分完成两项任务。
第一,它用一离散的模型替代连续模型。
因此,网格划分将问题简化为一系列有限多个未知域,而这些未知域符合由近似数值技术的求解结果。
第二,它用一组单元各自定义的简单多项式函数来描述我们渴望得到的解(e.g位移或温度)。
对于使用者来说,网格划分是求解问题必不可少的一步。
许多FEA初学者急切盼望格划分为全自动过程而几乎不需要自己输入什么。
随着经验的增加,就会意识到这样一个现实:网格划分常常是要求非常苛刻的任务。
商用FEA软件的发展历史见证了网格划分对F EA 用户透明的诸多尝试,然它并不是一条成功的途径。
而当网格划分过程既简单又自动执行时,它也仍旧不是一个“非手工干涉”而仅靠后台运行的任务。
作为FEA用户,我们想要有一种可以和网格划分过程交互的方法。
COSMOSW orks通过将用户从那些纯粹网格细节意义上的问题中解脱出来,找到了良好的平衡点;并使我们在需要时可以控制网格划分。
几何体准备理想情况下,我们用SolidW orks的几何体,联入COSMOSWorks环境。
在这里,我们定义分析和材料的类型,施加载荷与约束,然后为几何体划分网格并得到求解。
这种方法在简单模型下能起作用。
对于更为复杂的几何体,则要求在网格划分前作些准备。
在FEA 的几何体准备过程中,我们从特定制造,CAD 几何体出发,为分析而特地构造几何体。
我们称这个几何体为FE A 几何体。
基于两者的不同要求,我们对CAD 几何体和FE A 几何体作一区别:CAD 几何体FEA几何体必须包含机械制造所需的所有信息必须可划分网格必须允许创建能正确模拟所关心资料的网格必须允许创建能在合理时间内可求解的网格通常,CAD 几何体不能满足FEA 几何体的要求。
第 3章 ANSYS 13.0 Workbench网格划分及操作案例网格是计算机辅助工程(CAE)模拟过程中不可分割的一部分。
网格直接影响到求解精 度、求解收敛性和求解速度。
此外,建立网格模型所花费的时间往往是取得 CAE 解决方案所 耗费时间中的一个重要部分。
因此,一个越好的自动化网格工具,越能得到好的解决方案。
3.1 ANSYS 13.0 Workbench 网格划分概述ANSYS 13.0 提供了强大的自动化能力,通过实用智能的默认设置简化一个新几何体的网 格初始化,从而使得网格在第一次使用时就能生成。
此外,变化参数可以得到即时更新的网 格。
ANSYS 13.0 的网格技术提供了生成网格的灵活性,可以把正确的网格用于正确的地方, 并确保在物理模型上进行精确有效的数值模拟。
网格的节点和单元参与有限元求解,ANSYS 13.0在求解开始时会自动生成默认的网格。
可以通过预览网格,检查有限元模型是否满足要求,细化网格可以使结果更精确,但是会增 加 CPU 计算时间和需要更大的存储空间,因此需要权衡计算成本和细化网格之间的矛盾。
在 理想情况下,我们所需要的网格密度是结果随着网格细化而收敛,但要注意:细化网格不能 弥补不准确的假设和错误的输入条件。
ANSYS 13.0 的网格技术通过 ANSYS Workbench的【Mesh】组件实现。
作为下一代网格 划分平台, ANSYS 13.0 的网格技术集成 ANSYS 强大的前处理功能, 集成 ICEM CFD、 TGRID、 CFXMESH、GAMBIT网格划分功能,并计划在 ANSYS 15.0 中完全整合。
【Mesh】中可以根 据不同的物理场和求解器生成网格,物理场有流场、结构场和电磁场,流场求解可采用 【Fluent】、【CFX】、【POLYFLOW】,结构场求解可以采用显式动力算法和隐式算法。
不同的 物理场对网格的要求不一样,通常流场的网格比结构场要细密得多,因此选择不同的物理场, 也会有不同的网格划分。
基于Netgen的四面体网格剖分算法及其应用魏斌;徐华【摘要】A Netgen-based tetrahedralization algorithm was proposed,which fulfilled the tetrahedralization based on complex TIN models after close,non-manifold and surface normal vector checks of TIN model through analysis of various kinds of mo-dels.No points were added on the constraint surface of TIN models in Netgen-based tetrahedralization,leading to the high con-sistency between the original surface and the tetrahedron shape,and ensuring the accuracy of the model and boundary consistency under coplanarity.A sample application in a mineral's formation data shows that the Netgen-based tetrahedralization algorithm is stable,efficient and ensures the boundary consistency and high quality of mesh.%通过对各类模型的分析和归纳,提出一种基于Netgen的四面体剖分算法,对不规则三角网(TIN)模型进行封闭性检查、非流形检查、表面法向量检查和相交性检查,实现基于TIN复杂模型的四面体剖分.进行四面体剖分时,不在模型的约束面加点,使原始曲面与四面体形状间高度吻合,确保模型精度,保证共面情况下的边界一致性.某矿地层数据的实例应用结果表明,使用该算法能够使边界一致,得到高质量网格,兼具稳定性和高效性.【期刊名称】《计算机工程与设计》【年(卷),期】2017(038)002【总页数】6页(P368-373)【关键词】四面体剖分算法;封闭性检查;非流形检查;表面法向量检查;相交性检查;边界一致性;单元质量【作者】魏斌;徐华【作者单位】北京化工大学信息科学与技术学院,北京 100029;北京石油化工学院信息工程学院,北京 102617【正文语种】中文【中图分类】TP391.41如何解决复杂模型剖分时的边界一致性和生成高质量单元是四面体剖分的难点之一。
一种基于热核的三维人体网格模型分割方法【摘要】提出了一种基于热核的三维人体模型分割方法。
首先计算三维人体模型顶点的热核签名值;然后采用阀值分割法利用设定的分割阀值对模型进行初始分割;最后对阀值分割结果中大于阀值的部分进行凝聚聚类,得到最终的分割结果。
实验结果表明,分割结果不仅符合视觉最小值原则,还具有很好的语义信息;同时分割结果具有与姿势无关和鲁棒性好的特点。
【关键词】热核签名;人体模型;模型分割;谱聚类1.引言随着三维人体扫描技术的发展,三维人体模型已经无处不在,广泛应用于人机工程、游戏、动画、服装等领域。
由于人体模型越来越精细,复杂度越来越高,为了简化人体结构的复杂性,同时更好地适用于其它应用领域的需求,通常要求将三维人体模型分割成头部、躯干、上肢和下肢等部分,为人体模型的高级表示和进一步的语义特征分析提供基础。
模型分割是形状特征分析的基础[1]。
由于人体外形受年龄、性别、姿态等的影响,人体模型的自动分割成为一个具有挑战性的问题,国内外学者对此问题进行了比较深入的研究。
Nurre等[2]应用一种六杆的棒形人体模板,通过整体几何的匹配来达到细节特征识别,进而把人体分成与模板对应的头、左臂、躯干、右臂、左腿、右腿六个功能部分;Wang等[3] 通过对水平面与人体模型相交的截面形状,应用模糊规则确定拐点位置来识别腋窝和胯特征点,然后把人体模型分割为与文献[2]相同的六部分;Werghi等[4]提出了一种基于Reeb图的三维人体分割算法。
通过选取测地距离函数作为莫尔斯函数建立Reeb图,该方法对体形、姿态等变化具有鲁棒性。
其缺点是计算量大,算法效率不高;Katz等[5]提出了一种通用的三维网格模型模糊聚类分割算法,该算法可以处理任意姿态、不同部件比例的模型,并且应用了层次分割的思想,在避免了过分割和边界锯齿的同时,能够得到不同粒度层次分割结果。
此算法应用与人体模型可以得到较好的分割结果。
本文以三维人体测量为背景,采用坐标无关、几何不变的模型特征描述子,实现一种基于热核的人体模型分割方法。
ANSYSWorkbenchMesh⽹格划分(⾃⼰总结)Workbench Mesh⽹格划分分析步骤⽹格划分⼯具平台就是为ANSYS软件的不同物理场和求解器提供相应的⽹格⽂件,Workbench中集成了很多⽹格划分软件/应⽤程序,有ICEM CFD,TGrid,CFX,GAMBIT,ANSYS Prep/Post等。
⽹格⽂件有两类:①有限元分析(FEM)的结构⽹格:结构动⼒学分析,电磁场仿真,显⽰动⼒学分析(AUTODYN,ANSYS LS DYNA);②计算流体⼒学(CFD 分析)分析的⽹格:⽤于ANSYS CFX,ANSYS FLUENT,Polyflow;这两类⽹格的具体要求如下:(1)结构⽹格:①细化⽹格来捕捉关⼼部位的梯度,例如温度、应变能、应⼒能、位移等;②⼤部分可划分为四⾯体⽹格,但六⾯体单元仍然是⾸选;③有些显⽰有限元求解器需要六⾯体⽹格;④结构⽹格的四⾯体单元通常是⼆阶的(单元边上包含中节点);(2)CFD⽹格:①细化⽹格来捕捉关⼼的梯度,例如速度、压⼒、温度等;②由于是流体分析,⽹格的质量和平滑度对结果的精确度⾄关重要,这导致较⼤的⽹格数量,经常数百万的单元;③⼤部分可划分为四⾯体⽹格,但六⾯体单元仍然是⾸选,流体分析中,同样的求解精度,六⾯体节点数少于四⾯体⽹格的⼀半。
④CFD⽹格的四⾯体单元通常是⼀阶的(单元边上不包含中节点)⼀般⽽⾔,针对不同分析类型有不同的⽹格划分要求:①结构分析:使⽤⾼阶单元划分较为粗糙的⽹格;②CFD:好的,平滑过渡的⽹格,边界层转化(不同CFD 求解器也有不同的要求);③显⽰动⼒学分析:需要均匀尺⼨的⽹格;注:上⾯的⼏项分别对应Advanced中的Element Midside Nodes,以及Sizeing中的Relevance Center,Smoothing,Transition。
⽹格划分的⽬的是对CFD (流体) 和FEM (结构) 模型实现离散化,把求解域分解成可得到精确解的适当数量的单元。
基于凹凸性方法的杂乱场景点云分割算法黄镇;韩慧妍;韩燮【摘要】为了处理杂乱场景中多个物体的分割问题,提出了一种基于RGB-D点云数据的分割方法;该方法先将场景点云超体聚类分解为基于体素网格的邻接图,然后对邻接图的边缘进行分类创建凸度图,再通过区域生长合并具有凸关系的分块从而得到未知物体.此外,提出用欧几里得算法对区域生长进行改进,发现对于碗和杯子这类具有内部凹面的物体有较好地分割效果.在对象分割数据库和手动提取场景中的实验结果,表明该方法可以在杂乱的桌面场景中分割各种形状的对象.%In order to deal with the problem of segmentation of multiple objects in clutter scene , a segmentation method based on RGB-D point cloud data was proposed .The method firstly decomposed the scene cloud superclass cluster into the adjacency graph based on the voxel mesh , then classified the edge of the adjacent graph to create the convexity graph , and then merged the convexity with the convexity through the regional growth to get the un-known object .In addition , it was proposed to improve the regional growth with the Euclidean algorithm and found that the objects with internal concave surfaces such as bowl and cup had a better segmentation effect .The results in the object segmentation database and the manual extraction of the scene indicate that the method can segment ob -jects of various shapes in cluttered desktop scenes .【期刊名称】《科学技术与工程》【年(卷),期】2018(018)014【总页数】5页(P43-47)【关键词】计算机视觉;超体聚类;凹凸性;欧几里得;点云分割【作者】黄镇;韩慧妍;韩燮【作者单位】中北大学大数据学院,太原030051;中北大学大数据学院,太原030051;中北大学大数据学院,太原030051【正文语种】中文【中图分类】TP391.41尽管进行了很长时间的研究,将场景分割成物体仍是计算机视觉中比较有挑战性的课题之一[1]。