基于网格模型三维区域面积测量.docVIP

基于网格模型三维区域面积测量 　　摘要摘要：提出一种基于网格模型的物体区域面积计算方法，首先通过人机交互选取区域特征点，根据特征点生成整个区域的边界，采用平面与网格求交线的方法生成两点间最直的测地线，用测地线圈划整个区域的边界。重建被切割后的区域网格，通过区域生长法得到区域内所有三角面片的拓扑信息和坐标信息，将所有小面片的面积相加得到整个所选区域的面积。创新之处在于运用改进的网格切割算法，在切割过程中通过不断修正切割面使生成的切割线达到最优。实验证明，所提方法可行，可有效计算三维物体表面的区域面积，相对误差较小。 　　关键词关键词：三维模型； 网格切割；测地线；网格重建；区域生长；面积计算 　　DOIDOI：10.11907/rjdk.1511425 　　中图分类号：TP319 　　文献标识码：A文章编号文章编号2016）002009804 　　0引言 　　三维物体的区域面积测量在人体烧伤[1]、医学整形、纹理映射[2]、地理测量[3]、文物修复等领域有着重要应用。传统的测量方法有光学测量法[4]、超声波法[5]等，但这些传统方法对设备要求高，而且价格昂贵、不方便操作。近年来，由于数字获取设备和计算机技术的发展，通过三维数字彩色扫描[6，7]获取三维物体的几何和纹理信息数据，并构建三维网格模型越发便捷。基于此，本文提出基于网格模型的三维物体兴趣区域划分、提取及面积测量方法。对于兴趣区域的提取一般可分为基于颜色的自动提取和基于边界划分区域生成。文献[8]基于颜色分割的区域面积提取方法可以通过颜色和空间坐标信息将不同颜色区域分割成不同部分，基于颜色的兴趣区域自动划分虽然方便，但存在许多不足，如当颜色边界信息比较模糊时，提取的区域和需要选取的区域相差较大，在要求比较精确的领域（如人体烧伤面积计算），该方法存在明显不足。文献[9]提出了基于交互方式的区域生成方法，能准确生成兴趣区域，但该方法对噪声敏感，且模型点数量过大会导致运行速度过慢。为此，本文采用基于平面切割的兴趣区域边界划分方法，通过目测颜色信息指定边缘的若干关键点，并由平面切割算法自动生成兴趣区域的边界，然后采用区域生长方法[10]快速、精确提取指定区域，进而计算其面积。本文的创新点在于运用改进的网格切割算法生成区域边界，在切割过程中通过不断修正切割面使切割边界达到最优，从而使选取的区域边界更准确，实验证明，通过本文方法计算的面积更准确、有效。 　　1基于测地线的区域边界划分 　　测地线是光滑曲面上的最短路线，如地球表面上的最短路程。在网格模型上表现为最短路径，又称为近似测地曲线。Chen[11]基于Dijkstra算法提出了网格上最短路径生成的方法，但所给出的时间复杂度为O（n2），当网格点达到105数量级时速度很慢，而且该方法对网格的要求也比较严格，网格必须是连通的。Lee等[12]以基于边界条件的平面网格切割法生成最直路线，该方法时间复杂度为O（n），速度快，但该方法只处理带边界的网格。本文提出用不带边界的平面切割算法来生成三维网格测地线。 　　1.1基于平面切割算法的曲面测地线 　　在经典空间几何中，空间曲面的测地线有如下性质：①测地线上每一点的密切平面总是垂直于此点的切平面；②空间曲面上两点之间的测地线是两点之间最短的距离。本文提出的基于平面切割的测地线生成方法正是基于此思想。在求测地线过程中，本文使用一个尽可能垂直于网格曲面的平面去切割网格曲面，所生成的交线就是两点间的测地线。 　　理想情况下，切割面应垂直于网格曲面所对应的平面，也就是说切割面的法向量与测地线上每一个点主法向量内积之和尽可能小。定义曲面上测地线的平均法向量N=1l∫l0Nt（s）ds。其中，Nt（s）为以s为参数的法向量函数，l为测地线的长度。由于网格上的最短路径由一段段直线段组成，因此平均法向量可以离散化表示为N=1l∑n1liNi，li和Ni分别表示切割路径所经过的第i个三角面片的长度和法向量。设起点与终点的方向向量Vbe=Ve-Vb（Ve，Vb分别为终点和起点的方向向量），则切割平面的法向量N=1l∑n1liNi×Vbe，由于确定切割面之前，无法得知切割路径，为得到最优的测地线，本文采用基于修正的切割平面生成算法。其基本思想是：如果当前切割点为起始点时切平面的法向量为N=Vbe×Nbe，Vbe=Vb-Ve，Nbe=Nb+Ne，其中Vb，Ve为起点和终点的方向向量，Nb，Ne为起点和终点的法向量。如果当前切割点不是起始点，将该点作为新的起点V′b，该点的法向量为N′b，更新切割平面的方向向量V′be=V′b-Ve，更新基平面的法向量N′be=N′b+Ne，更新切平面的法向量N=V′be×N′be。关于切割点N′b法向量的计算可分两种情况：当点在三角