[发明专利]一种三维模型形状特征二叉树的构建方法

摘要

本发明一种三维模型形状特征二叉树的构建方法，现有的三维模型形状特征计算方法主要是将三维模型投影到各个方向上的一个平面上，获得三维模型的投影图像，通过计算投影图像的形状特征来获得三维模型的形状特征。本发明方法的具体步骤是：首先计算三维模型M的每个三角片的面积及总面积并计算三维模型所有顶点与中心的距离，对三维模型进行比例变化。其次划分三角片，映射所有元素到球面，把PRO(r_min，r_max)作为二叉树T根结点的球面投影图像。再次分别计算K结点的左右孩子结点的球面投影图像PRO(r₁，r)和PRO(r，r₂)。最后对二叉树T的所有结点包含的球面投影图像计算球面调和特征向量，得到三维模型M的特征二叉树。本发明可以对三维模型进行整体到局部的相似性比较，具有较好的准确性。

主权项

1.一种三维模型形状特征二叉树的构建方法，其特征在于该方法的具体步骤是：步骤(1).计算三维模型M的每个三角片的面积S_i，i＝1，2，3…N，N为三维模型的三角片数量，并计算M的总面积S=Σi=1NSi;]]>计算M所有三角片的中心坐标(x_ic，y_ic，z_ic)，按公式1计算M的中心O_C，其坐标值为(x_oc，y_oc，z_oc)；xoc=Σi=1N(xic×Si)/N,]]>yoc=Σi=1N(yic×Si)/N,]]>zoc=Σi=1N(zic×Si)/N]]>步骤(2).以M的中心O_c为三维坐标系的坐标原点，计算三维模型所有顶点与O_c的距离，记与三维模型中心最大的距离为D_max；对三维模型的所有顶点的坐标(x_i，y_i，z_i)变换为(x_i/D_max，y_i/D_max，z_i/D_max)；三维模型的每个三角片面积变换为S′_i＝S_i/D_max；步骤(3).如果M存在面积S′_i＞1/2000的三角片，那么这个三角片划分为6个子三角片，划分的具体方法是：三角片的三个顶点分别与三角形的中心相连，其连线的延长线与对边相交，三角片内的所有连线与延长线将三角片划分为6个子三角片；如果划分后存在面积S″＞1/2000的子三角片，那么采用相同的方法对这个子三角片进行三角片划分，直到进行三角片划分后的三维模型的所有三角片的面积都小于1/2000为止；设定进行三角片划分后的三维模型为M′，计算M′所有的三角片的中心点坐标，将三维模型M′的所有三角片的中心点和三角片的顶点一起作为三维模型M′的采样点集合P；P＝{p|p是M的三角片的顶点|P是M的三角片的中心点}P中的采样点p到O_c的距离称为顶点半径；步骤(4).记L(p，O_c)是采样点集合P中的元素p到O_c的距离；P(r₁，r₂)是如下公式所示的集合：P(r₁，r₂)＝{p|r₁＜L(p，O_c)≤r₂，r₁＜r₂}，r₁、r₂为大于零的实数记半径为r₂球面为S_r2，把S_r2球面按经纬线均匀划分，经度线的划分单位是5°，纬度线的划分单位是5°；划分后，S_r2球面被划分为一系列的经纬度网格，设定所有经纬网格的权值为0，然后把P(r₁，r₂)的所有元素按极坐标方向映射到球面S_r2上面，并且投影后的每个经纬网格的权值等于P(r₁，r₂)中映射到相应经纬网格内的所有顶点半径的最大值；把带有权值S_r2经纬度网格称为球面投影图像PRO(r₁，r₂)；步骤(5).设r_min和r_max分别是三维模型M′的采样点集合的最小顶点半径和最大顶点半径，把PRO(r_min，r_max)作为二叉树T根结点的球面投影图像；步骤(6).设二叉树T的某个结点K包含集合P(r₁，r₂)，它的球面投影图像为PRO(r₁，r₂)，设r＝(r₁+r₂)/2，把P(r₁，r₂)划分为两个子集P(r₁，r)和P(r，r₂)，把P(r₁，r)作为K结点的左孩子结点的集合，把P(r，r₂)作为K结点的右孩子结点的集合，分别计算K结点的左右孩子结点的球面投影图像PRO(r₁，r)和PRO(r，r₂)；从根节点开始，按步骤(5)进行递归集合划分；步骤(7).对二叉树T的所有结点包含的球面投影图像计算球面调和特征向量SHD；球面调和特征向量SHD的计算如下列公式所示：SHn(θ,φ)=Σm=-nnan,mYnm(θ,φ)]]>其中an,m=∫02π∫0πsin(θ)f(θ,φ)Ynm(θ,φ)dθdφ]]>Ynm(θ,φ)=(2n+1)(n-m)!4π(n+m)!Pnm(cos(θ))eimφ]]>SHD={||SH1(θ,φ)||||SH0(θ,φ)||,||SH2(θ,φ)||||SH0(θ,φ)||,||SH3(θ,φ)||||SH0(θ,φ)||,···,||SHn(θ,φ)||||SH0(θ,φ)||}]]>其中，θ是球面投影图像的经度数，φ是球面投影图像的纬度数，n＝1，2，3…，m＝±1，±2，±3…±n，P_n^m(x)是拉格朗日多项式；二叉树T的每个结点的球面投影图像采用步骤(7)的方法计算SHD，得到三维模型M的特征二叉树。