[发明专利]一种基于凸包与OBB的三维网格模型骨架提取方法

摘要

本发明公开了一种利用凸包与OBB的三维网格模型骨架提取方法。本发明通过对三维网格模型进行分割，生成多个子网格模型，再通过对各子网格的点集求最小凸包，使用凸包顶点构造子网格的凸包近似，以及计算每个凸包点集的形心，构造原始骨架点，然后用OBB包围盒进行重叠计算，求出相交点集，生成关节骨架点，最终对原始骨架点与关节骨架点进行连接，最终形成完整骨架。本发明的优点是数据准确、操作方便，计算效率高等优点。本发明在计算机动画、物体识别和匹配、模型重建、虚拟导航等领域也有着重要的应用价值。

主权项

1.一种基于凸包与OBB的三维网格模型骨架提取方法，其特征在于包括下述步骤：(1)网格分割对一整体网格模型，基于特征的区域生长算法对网格进行分割处理；首先输入一个三维网格模型，在三维网络模型上标注出前景曲线与背景曲线，然后计算从视点发出的射线通过前后背景曲线与网格的交点，定义为C＝{C_F，C_B}，其中C_F、C_B分别为前景曲线与背景曲线的投影点集；再根据Isophotic度量距离进行区域增长；Isophotic度量考虑的因素有：曲面上边的长度、曲面法向量、曲率的度量；Isophotic度量的公式化离散形式表述为：dΓ(p,q)=|p-q|+ω|np-nq|+ω*f(kpq‾),]]>其中d_Γ(p，q)为p，q两点的Isophotic度量距离；n_p，n_q分别为两点的法向量；为曲面沿线方向的曲率；ω、ω^*为权重，设为5；f运算定义如下：f(kD)=kD,kD≥0e|kD|-1,kD<0]]>以C_F、C_B为“种子点”，根据Isophotic度量应用区域生长策略，将全部网格顶点划分为两类，即前景点和背景点，也就将网格分为两部分；在边界处，应用基于蛇形曲线的最小化能量方法进行边界优化，最终生成所要的分割结果；(a)将网格分为头部与身体两个部分；(b)再上一步基础上进一步将四肢划分出来；(c)为最后分割效果，相邻子网格之间用不同颜色表示，加以区分；(2)凸包近似对各子网格进行凸包近似处理：凸包是包含点集合S中所有对象的最小凸集，凸包的构造涉及到两个因素，即从点集S中筛选凸壳顶点以及建立各顶点之间的拓扑关系，采用改进的增量算法来构建三维凸包，具体如下：遍历点集S，取出x坐标最大与最小的两点，与y坐标最大和z坐标最大的点构建一初始四面体，然后将多面体的每个面的三个顶点均逆时针排列，然后每次从点集S中取出一个点，对该点进行可见性判断；若该点对于所有面都不可见，则该点位于凸包内，不做处理；若有可见面，则遍历棱边，若该棱边的两个邻接面均可见，则删除该棱边，若该棱边只有一个可见面，则用新加入的点与该棱边创建新面，并创建两条新边，最后删除所有可见面，这样依次加入新点，直至点集S为空；其中可见性判断用面法向量与待检测点到面上一端点向量的点积的正负号判断；(3)确定初始骨架接下来构造每个子网格的包围盒：OBB是三维空间中的一个轴向任意的长方体包围盒；将三角形顶点集合看作是三变量的概率分布函数，在对OBB包围盒的进行计算时，利用三角形顶点的均值和协方差矩阵来计算包围盒的中心位置和轴向；设三角形集的第i个三角形的顶点矢量为p_i，q_i和r_i，包围盒包围的三角面片数，即凸壳面片数为n，则凸包顶点的均值为：Cm=13nΣi=1n(pi+qi+ri).]]>其协方差矩阵为：D=Σi=1n(pi‾x2+qi‾x2+ri‾x2)Σi=1n(pi‾xpi‾y+qi‾xqi‾y+ri‾xri‾y)Σi=1n(pi‾x,pi‾z+qi‾xqi‾z+ri‾xri‾z)Σi=1n(pi‾xpi‾y+qi‾xqi‾y+ri‾xri‾y)Σi=1n(pi‾y2+qi‾y2+ri‾y2)Σi=1n(pi‾ypi‾z+qi‾yqi‾z+ri‾yri‾z)Σi=1n(pi‾xpi‾z+qi‾xqi‾z+ri‾xri‾z)Σi=1n(pi‾ypi‾z+qi‾yqi‾z+ri‾yri‾z)Σi=1n(pi‾z2+qi‾z2+ri‾z2),]]>其中向量分别为的x，y，z分量，点q与r相关变量定义类似；协方差矩阵D_j，k是一个实对称矩阵，所以它的三个特征向量是正交的，正规化后可作为一个基底，这就确定了OBB的三个轴向，分别计算各个顶点在该基底的三个轴向上的最大值和最小值，以确定该OBB的大小；对各子网格求得到凸包顶点的OBB包围盒后，将各个均值C_m作为初始骨架点，连接各相邻部分的初始骨架点，就得到了初始骨架；(4)生成完整骨架上述步骤中的初始骨架只是将各个包围盒的中心点简单连接起来，若两部分的轴向成一定的角度，则连线就会远离中心位置，不满足中心性要求，需对初始骨架进行优化处理，增加一些重要关节骨架点；该处理过程分两步，首先确定两相邻子网格点集的OBB交集内所有点，然后再求取这些点的中心点作为关节骨架点，这样连接后就生成了完整骨架；(4.1)求取相邻OBB交集内的点集已知点集A的中心点C_A、三个轴向单位向量A_xA_yA_z和在三个轴向上的长度中心点C_A到相邻点集B中的任意点的向量为p→A→Bi=pBi-CA.]]>点在局部坐标系A(轴向为A_xA_yA_z，原点为C_A)中的坐标为：pBAxi=p‾A→BigAxpBAyi=p→A→BigAypBAzi=p‾A→BigAz,]]>其中为点在局部坐标系A中的A_x轴向上的坐标，其余依次类推；若满足以下不等式方程组-d(-Ax)≤pBAxi≤d(Ax)-d(-Ay)≤pBAyi≤d(Ay)-d(-Az)≤pBAzi≤d(Az),]]>则该点位于点集A的OBB包围盒内；依次对B内所有点进行判断，确定点集B内所有位于点集A的OBB包围盒内的点；然后反向操作，求出点集A内所有位于点集B的OBB包围盒内的点，两者的并集即为我们要求的OBB交集内所有点的集合；(4.2)确定关节点，生成骨架按上述步骤求出所有相邻OBB包围盒相交区域内的点集后，利用下式Cj=1nΣi=1n(pi),]]>其中n为OBB交集内的点集总数，p_i为该点集中的第i个点；求取OBB交集内所有点的均值坐标，将之作为关节骨架点，并与两个相关联的骨架点相连接，从而形成完整的骨架。