CN101201845B - 基于中轴点集分层螺旋信息的三维模型检索方法 - Google Patents
基于中轴点集分层螺旋信息的三维模型检索方法 Download PDFInfo
- Publication number
- CN101201845B CN101201845B CN2007101950618A CN200710195061A CN101201845B CN 101201845 B CN101201845 B CN 101201845B CN 2007101950618 A CN2007101950618 A CN 2007101950618A CN 200710195061 A CN200710195061 A CN 200710195061A CN 101201845 B CN101201845 B CN 101201845B
- Authority
- CN
- China
- Prior art keywords
- dimensional model
- information
- model
- point
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 124
- 239000013598 vector Substances 0.000 claims abstract description 18
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims description 16
- 238000012360 testing method Methods 0.000 claims description 14
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- 230000007935 neutral effect Effects 0.000 claims 3
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 3
- 239000000284 extract Substances 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 8
- 230000009466 transformation Effects 0.000 description 8
- 238000005070 sampling Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000000844 transformation Methods 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 3
- 241000238421 Arthropoda Species 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011524 similarity measure Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012733 comparative method Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明提出一种基于中轴点集分层螺旋信息的三维模型检索方法,该方法提取三维模型的中轴点集分层螺旋信息作为特征描述符,用一种基于Hausdorff距离的多级加权相似性匹配方法度量三维模型之间的相似性。首先对三维模型数据库中的三维模型进行坐标标准化及体素化预处理;然后提取所述三维模型的特征并转化为特征描述符,生成特征数据库;检索时将用户提交的三维模型按照上述步骤方法依次进行标准化、体素化及特征提取以得到其三维模型的特征描述符;然后再将该特征描述符与特征数据库中的特征进行匹配,生成检索结果。该方法的优点是:形状特征描述更加完整;形状特征描述符的比较被简化为不等长向量之间的距离计算问题;方法性能好,效率高。
Description
技术领域
本发明涉及一种三维模型检索方法,特别是涉及一种基于内容的三维模型检索方法。
背景技术
随着互联网的发展,三维模型数量正在成倍的增长。过去十年中,在医药、化学、建筑和娱乐等许多领域,三维模型的获取和可视化工具已经成为数据处理的一个完整部分,这些工具的发展同样使得可用的三维模型呈爆炸式的增长。面对庞大的三维模型海洋,大家面临的问题已经从“如何产生三维模型”渐变为“如何查找三维模型”。从大规模数据库中检索模型的需求正在快速增长,而形状分析所关心的核心问题是设计有效和健壮的匹配方法。
三维模型的相似性度量是基于形状的识别、检索、聚类和分类的基本问题,主要应用于计算机视觉、机械工程和分子生物学等领域。但是近期以来,一是计算机软件和硬件快速发展,桌面计算机的三维图形硬件和CPU已经能够快速处理和显示三维数据,这促使了三维模型向多领域发展的需求;二是建模和数字技术的发展使构建三维模型变得更加简便;三是互联网对全球用户建模的三维模型都是开放的,它使三维模型得到了广泛的分布和共享。因此,无论是在因特网或是其它地方,三维模型的数量在飞速增长的同时,三维形状分析和匹配的应用也已经扩展到广泛的领域之中。以上这些更加强调了从海量三维模型数据库中精确的检索出所需要的三维目标的技术发展的需求。
处理二维图像时可以用形状的属性识别侧面影像或轮廓线,如使用曲率或图像的一些性质:颜色、纹理或小波系数。这些技术过去常用在搜索二维图像数据库中。
过去已有一些处理三维形状匹配和特征提取的方法。例如:(1)基于全局特征的方法,其中一种是将球面函数集合作为一种反射对称描述符;(2)基于统计特征的方法,其中一种是将模型的特征表示为一种形状分布,源于一个测量三维目标的全局几何属性的形状函数,该方法准确度的有效性取决于采样点的数量,更多的采样点得到更高的准确度,相反,效率和采样点的数量呈反比;(3)基于频域空间映射的方法,其中一种是一种新的基于球面调和的数学工具,以此来获得模型的旋转不变表示并可以用来辨别模型的差异,该种方法的缺点是不支持局部匹配,因为它们没有反映特征和物体各部分的相互关系;(4)基于视图的方法,其中一种是实现一种二维草图查询接口,并且开发了匹配二维草图和三维目标的方法,另外还提出了用一种二维形状采样统计函数的傅里叶描述符作为提取三维模型的特征的方法,这些方法的优点是接口友好,任何用户均能够使用二维草图工具检索相似的三维模型。上述四类方法的缺点是它们不能得到形状的具体细节信息,并且健壮性不佳。(5)基于局部特征的方法,其中一种是一种上下文的方法用来度量三维形状的相似性,该类方法可用于局部匹配,但必须首先对准模型,且效率较低。
还有一类是基于图的度量三维形状相似性的方法。其中基于模型图的方法主要用于CAD/CAM领域,很难应用于类似人和动物这样的具有自然形状的模型。在基于Reeb图的方法中,Hilaga等(Hilaga,M.,Shinagawa,Y.,Kohmura,T.,Kunii,T.L.:Topology matching for fullyautomatic similarity estimation of 3D shapes.Proceedings of the 28thannual conference on Computer graphics and interactive techniques(SIGGRAPH 2001),Los Angeles,CA,USA(2001)203-212)通过比较三维模型的多分辨率Reeb图(MRGs,Multiresolutional Reeb Graphs)来计算模型之间的相似性。由测地距离定义Reeb图适合匹配关节式物体,但是对拓扑变化比较敏感。基于Reeb图方法的不足主要有以下几点:(1)受表面连通性的影响;(2)对几何属性比拓扑属性更敏感;(3)在子图匹配时不总是能够经受住检验;(4)不总是能够表示骨架;(5)产生了不同密度的顶点。总之,由测地距离定义的Reeb图适合匹配关节目标,但是对拓扑改变敏感。另外,他们不能被应用于任意的网格结构,因为类似于丢失面片这样的拓扑问题会对测地距离的计算造成干扰。基于骨架的方法计算模型的“骨架”并将它转换为骨架图,以此作为该模型的形状描述符。骨架的概念由Blum(Blum,H.:Atransformation for extracting new descriptors of shape.In:Whaten-Dunn,W.(ed.):Proceedings of the Symposium on Models for Perception ofSpeech and Visual Form.MIT Press,Cambridge,MA(1967)362-380)提出。二维中的骨架是中轴,三维中的骨架是中间曲面等。基于骨架图的方法有四类:细化和边界传播、基于距离场、几何和广义场函数。细化方法能产生连通的曲线,但是该曲线需要平滑处理且对噪声敏感;距离场和几何方法不能直接生成曲线骨架,且最后产生的骨架点需要被修剪,然后再连通为曲线骨架;势场方法能够产生最细和最清晰的骨架,但是耗费时间长。总之,拓扑一般被表示为如图和树这样的相关数据结构,实际上绝大多数模型只能表示为图。因此这类方法的相似性评估比其它类型方法的更为复杂。另外,基于图的相似性度量具有一定的局限性:图比较的计算量和图的大小成正比。
发明内容
鉴于以上缺陷,本发明的主要目的在于提供一种基于中轴点集分层螺旋信息的三维模型检索方法,该方法选择的检索项和匹配方法包括精确的形状表示,能够很好的描述形状的特征,能够自动、快速、健壮和有效地实现。
为了达到以上目的,本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法,其提取三维模型的中轴点集分层螺旋信息作为特征描述符,用一种基于Hausdorff距离的多级加权相似性匹配方法度量三维模型之间的相似性。该方法适用于多边形模型也适用于体元模型。中轴点是一种向量结构,它存储骨架化过程中产生的各种实体。基于中轴点的方法的优点是:中轴点能够用做子图同构,并且计算成本较低;另外,存储的子拓扑局部属性能够被用做更加精确的比较。
本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法,其主要通过以下步骤实现:
(1)对三维模型数据库中的每个模型进行坐标标准化预处理;
(2)对坐标标准化预处理后的除体元模型之外的三维模型进行体素化预处理;
(3)对每个体素化的三维模型进行特征提取,得到每个三维模型的特征描述符,生成特征数据库;
(4)将用户提交的三维模型首先按照步骤(1)中的方法进行坐标标准化预处理、然后按照步骤(2)中的方法体素化其三维模型,再按照步骤(3)中的方法进行特征提取得到该三维模型的特征描述符;
(5)将用户提交的三维模型特征描述符与特征数据库中的特征进行匹配,计算相似距离并排序,最终生成检索结果。
其中所述步骤(1)对三维模型数据库中的每个模型进行坐标标准化预处理,获取标准化后的三维模型,具体步骤如下:确定三维模型包围盒的三个主轴;确定三维模型包围盒的中心和三个半长;确定三维模型包围盒主轴的正方向;根据包围盒的主轴、中心和正方向确定三维模型的本征坐标系。
所述步骤(2)对坐标标准化预处理后的除体元模型之外的三维模型进行体素化预处理,获取体素化的三维模型,具体步骤如下:面片顶点体素化;面片边的体素化;面片的体素化;三维模型内部体素化。
所述步骤(3)对体素化的三维模型进行特征提取,得到每个三维模型的特征描述符,生成特征数据库,具体步骤如下:对每个体素化的三维模型通过排斥力场函数提取模型的中轴点;对每个体素化的三维模型建立中轴点集分层螺旋信息特征描述符,生成特征数据库。
所述步骤(4)将用户提交的查询三维模型首先进行坐标标准化预处理、然后体素化为体元模型,再进行特征提取得到该三维模型的特征描述符,具体步骤如下:对用户提交的查询三维模型按照步骤(1)中的方法进行坐标标准化预处理,获取标准化后的三维模型;对坐标标准化预处理后的三维模型按照步骤(2)中的方法进行体素化预处理,获取体素化的三维模型;对体素化的三维模型按照步骤(3)的方法进行特征提取,得到用户提交的查询三维模型的特征描述符。
所述步骤(5)将用户提交的三维模型特征描述符与特征数据库中的特征进行匹配,计算相似距离并排序,最终生成检索结果。具体步骤如下:将用户提交的查询模型特征和特征数据库中的特征用一种Hausdorff距离的变种进行特征匹配;将特征匹配得到的所有数据依据距离大小进行相似度排序;根据相似度排序得到最终检索结果。
本发明的优点是:
(1)本发明提出的方法能够有效地区别和检索不同分组的三维模型,该方法在模型分类和相似性匹配时达到了比较高的准确度;
(2)本发明提出的方法能够保证三维模型的变换不变性和健壮性,其算法对刚体具有平移、旋转和缩放不变性,并且基本不受输入模型的噪声和网格化简的干扰;
(3)本发明提出的方法效率较高,该方法检索具有成千上万个多边形面片的三维模型平均耗时十多秒,更重要的是,如果是在上一次检索结果中二次检索三维模型或在模型库内检索三维模型时,只需要计算描述符的匹配时间,大约二毫秒即可返回查询结果。
综上所述,本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法对三维模型的分类和检索性能好,特征提取时间快,另外,本发明还可用于CAD、医学、分子生物学等领域。
附图说明
图1为本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法的流程图;
图2为本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法与基于体素化和球面调和的方法的总体查准-查全曲线比较图;
图3为本发明提供的基于中轴点集分层螺旋信息的三维模型检索方法与基于体素化和球面调和的方法在部分类上的查准-查全曲线比较图。
具体实施方式
三维模型匹配的关键是提取一个能够有效匹配三维形状的描述符。该描述符应具有以下特性:(1)相似变换不变性;(2)快速计算;(3)具有简明的存储结构。
本发明使用中轴点集代替了中轴结构,且将中轴点集和它们的统计信息组合起来作为特征描述符。这样做的好处是:降低了计算复杂度,这是因为没有将中轴点集连接为骨架;保留了关键的骨架结构,因为通过基于势场方法提取的中轴点集是最重要的骨架信息;组合了模型的曲面信息;其中相似性度量方法减少了噪声的影响。
本发明提供一种基于中轴点集分层螺旋信息的三维模型检索方法,其流程如图1所示,包括以下步骤:
步骤一、对三维模型数据库中的每个三维模型进行坐标标准化预处理;
本实施例进行位姿标准化的方式是找到三维模型的唯一的包围盒,其采用一种法向量分布法来计算三维模型唯一的包围盒。该方法选取最大法向量分布作为一个主轴,首先计算出每个面片的面积,然后累加具有同样或相反法线方向的面片的面积。其中具有相同方向的法线被认为是具有相同的分布,具有最大面积的法线被定义为第一个主轴,剩余的法线分布中满足下列条件的法线方向被定义为第二个主轴:(1)具有最大的面积;(2)与第一个主轴正交。通过前面两个主轴的叉积运算得到最后一个主轴,本实施例采用通过面片的法线分布产生模型包围盒主轴方向的方法实现三维模型的位姿对准和标准化。
步骤二、对坐标标准化预处理后的除体元模型之外的三维模型进行体素化预处理;
位姿标准化后再将除体元模型之外的三维模型体素化为L×L×L的体元模型,其具体步骤为:首先将面片顶点体素化;再体素化面片的边,使用Bresenham算法解决平行于坐标平面的特殊边的体素化问题;然后体素化面片,体素化时分别将三角面片投影至XOZ、ZOY、YOX平面以得到面片的完整体素化,最后将三维模型内部体素化。
步骤三、对每个体素化的三维模型进行特征提取,得到每个三维模型的特征描述符,生成特征数据库;
对于体素化的三维模型,首先通过排斥力场函数提取三维模型的中轴点集,具体描述如下:
排斥力场方法的关键思想是给三维物体的边界充电以在物体内部产生排斥场。本步骤中,首先假设边界体元为点电荷,边界点被定义为一个物体体元距离坐标原点最近的立方体的角,内部点是在每个面上都邻接三维物体体元的模型体格中距离坐标原点最近的角,坐标原点被定义为三维模型包围盒的一个角。一个内部点的排斥力可以通过式(1)计算:
其中,CiP/Ri m是点电荷Ci产生的力,FP是在内部点P的排斥合力,CiP是从Ci到P的标准化向量,它指明了力的方向,R是P和电荷Ci之间的距离,幂m称为该力函数的阶。
本步骤中,临界点是矢量力消失的点,只有一部分中轴点被当作临界点。矢量场中的临界点很难定位,尤其是因为它们不一定位于给定的采样位置,但是它们一般都位于采样点之间。如果力矢量(x,y,z)的三个分量都消失,矢量场就为零。也就是说,当某个区域中每一个矢量分量的符号都改变时,该区域就是包含临界点的候选区域。具体的说,首先假设体单元是最小的区域,并计算格子单元八个角的力场值,对于该单元所有的矢量分量(x,y,z),如果即有正值又有负值,那么该单元就是潜在的候选单元,然后再将候选单元递归地分为八个子单元,并且对每个子单元重复做上述的候选性测试。其中,本实施例定义了下述规则:如果有两个或两个以上的候选子单元,则比较它们的力值并且选择具有最小力值的子单元替换原先的候选单元作为新的候选单元。处理过程终止的条件是:该单元未能通过候选性测试,或者该单元太小(用户指定大小的阈值),此时该单元仍然作为候选单元保留。如果是后一种情况,临界点的坐标就被定义为该单元的坐标。但是如果所有的八个子单元都未能通过候选性测试,本实施例定义该单元的坐标为候选临界点。本实施例把从上述算法所提取的临界点和候选临界点称为中轴点。
接着,提取了三维模型的所有中轴点后,本实施例建立三维模型的中轴点集分层螺旋信息特征描述符,具体描述如下:
从三维模型提取的所有中轴点的集合被称为中轴点集。中轴点集描述的特征仅仅表示了物体中轴的大致轮廓,它不能评估两个拓扑相同但形状却不同的三维物体。总体来说,本实施例是将三维物体沿着X轴均匀的分层,并判断中轴点所在的层,然后从中轴点沿着平行于YOZ平面的方向均匀的发出射线,使从每个中轴点所发出的射线数量和两两相邻射线之间的夹角相等。这样,360度被均匀的分配到每一层,并保证任何不同层的中轴点具有不同的射线方向,任何同层的中轴点具有相同的射线方向。
首先本实施例做出以下几个定义:
定义1.分层螺旋射线集合
首先定义全局射线总数,全局射线总数是指从某中轴点的发射的射线数量与存在中轴点的层的数量的乘积。设全局射线总数为R,分层总数为N,且N可被R整除。则,每一层分布R/N条射线。以沿着三维模型的X轴分层为例,在YOZ平面上等角度(Δθ)发射R条射线,其中Δθ=360/R。则第k层的射线与Y轴的夹角集合为θk={θk,j|0≤j<R/N},通过式(2)计算:
θk,j=(k+j·N)·Δθ(2)
在这里,θk,j表示第k层中任一个中轴点的第j+1条射线与Y轴的夹角。则
θk={kΔθ,(k+N)Δθ,(k+2N)Δθ,...,(k+(R/N-1)gN)Δθ},N层射线集合从全局上形成螺旋射线集合θ={θk|0≤k<N}。
定义2.中轴点信息
第k层的第i个中轴点的信息定义为向量Mi=(Ci,Di),其中,坐标信息Ci=(Xi,Yi,Zi)代表中轴点的坐标位置;距离信息为Di={Dij| 0≤j<R/N},在这里,Dij表示中轴点i的第j+1条射线。从中轴点i放射状地发射θk,j角度的射线Rij,定义Rij与模型表面首次相交的点为Interij,则Dij为中轴点i与Interij之间的距离。
定义3中轴点集分层螺旋信息
组合所有中轴点信息Mi的集合就是中轴点集分层螺旋信息M=[M0,M1,...,Mn-1]-1,即模型的特征描述符。
提取中轴点集分层螺旋信息的步骤如下。首先,均匀地沿着X轴将一个三维物体分为N层,并构造分层螺旋射线集合。第二,判断中轴点i属于的层,并使点i沿着该层的射线方向放射状地等角度发出R/N条射线,在这里,R是分辨率参数和射线的分布总数(即全局射线总数)。第三,我们沿着第j+1条射线用射线跟踪算法计算Interij的坐标值(与射线相交的最近的体格单元的坐标)。接着计算每个中轴点i和Interij之间的距离Dij,并结合坐标分量距离构造中轴点i的信息。最后,得到模型的特征描述符M=[M0,M1,...,Mm-1]-1,即该模型的中轴点集分层螺旋信息。
步骤四,将用户提交的三维模型首先进行预处理、然后体素化为体元模型,再进行特征提取得到该三维模型的特征描述符;
将用户提交的查询三维模型首先进行坐标标准化预处理、然后体素化为体元模型,再进行特征提取得到该三维模型的特征描述符,具体步骤如下:对用户提交的查询三维模型按照步骤一中的方法进行坐标标准化预处理,获取标准化后的三维模型;对坐标标准化后的除体元模型之外的三维模型按照步骤二中的方法进行体素化预处理,获取体素化的三维模型;对体素化的三维模型按照步骤三中的方法进行特征提取,得到用户提交的查询三维模型的特征描述符。
步骤五、将用户提交的三维模型特征描述符与特征数据库中的所有特征进行匹配,计算相似距离并排序,最终生成检索结果。
构建完用户提交的三维模型的中轴点集分层螺旋信息特征描述符之后,本实施例用一种匹配方法将该特征描述符与特征数据库中的特征进行匹配。有许多相似性匹配方法,如Euclidean距离、Manhattan距离、Hausdorff距离等。本实施例提出了一种基于Hausdorff距离的多级加权相似性匹配方法来进行三维模型之间的匹配。首先定义如下几个概念。
定义4中轴点加权信息
在中轴点信息的基础上,本实施例定义中轴点加权信息为:
其中,Ci +是坐标加权信息,Di +是距离加权信息,ω0是中轴点类型的权值。
定义5中轴点加权信息的距离
任意两个加权信息向量Mi +和Mi +的距离被定义为:
其中, ω是层的权值,
ω1是点的坐标信息或距离信息的权值。
本实施例采用一种基于Hausdorff距离的匹配器匹配用户提交的三维模型A的特征描述符与特征数据库中的三维模型B的特征的步骤为:
第一步,根据中轴点的类型对其赋予不同的权值,如式(3)。若为临界点,权值为ω0;若为候选临界点,权值为1-ω0,亦为第1级加权。这样,A和B各自描述为加权的中轴点集分层螺旋信息A=[a0,...,am-1]-1和B=[b0,...,bn-1]-1,其中,au和bv是中轴点加权信息,u∈[1,m],v∈[1,n]。
第二步,分别计算A和B的所有中轴点加权信息中Ci +和Di +的加权Euclidean距离,权值ω1调节坐标加权信息的Euclidean距离和距离加权信息的Euclidean距离所占的比重,亦为第2级加权。另外,根据式(4)计算两个中轴点的距离,判断两个中轴点是否同层,用对应的权值ω或1-ω加权这两个中轴点加权信息的距离,即第3级加权。A和B的任意两个中轴点加权信息au和bv之间的距离为d(au,bv)。
第三步,m可能不等于n,本实施例用Hausdorff距离度量模型A和模型B的差异值H(A,B):
H(A,B)=max{h(A,B),h(B,A)}
其中,
综上所述,将用户提交的查询模型和模型库中所有模型的特征按照上述匹配步骤计算,得出所有的特征匹配距离,然后将这些距离排序后按照相似度大小生成检索结果。
以下对本实施例提供的中轴点集分层螺旋信息相似性度量方法对三维模型的变换和干扰的健壮性进行相应测试。首先选择5个具有代表性的三维模型,并对所述模型进行五种变换。在此基础上,再建立了一个有30个模型的数据库,它包含每个模型的6个不同版本,其中包括1个原型及5个变换。所做变换如下:
·缩放:X、Y、Z每一维放大10倍;
·各相异性缩放:Y维放大5%,Z维放大10%;
·旋转:分别绕X维、Y维和Z维旋转45度;
·剪切:分别将Y维和Z维放大X维的5%和10%;
·插入:随机插入5%的多边形复件。
在所述数据库上,对本实施例方法区分不同类的能力进行测试,具体是计算30个模型描述符的差异性度量,结果证明了本实施例的方法具有很好的鉴别力,对刚体具有平移、旋转和缩放不变性。
此后再测试具有不同数量的多边形网格的相似三维物体的差异性。本实施例用简化封套(Simplification Envelopes)软件生成了Stanford Bunny的8个版本,其面片数量被剖分为600到70000不等。提取每个版本的中轴点集并计算它们之间的差异性度量值,实验结果可以明显地观察到每个化简模型版本得到的中轴点集大都非常相似,进一步从量化的数值来看,每个版本的平均差异性度量值都非常接近,说明它们之间的距离几乎相同。该部分实验证实了本实施例方法对改变网格、模型化简或模型退化具有良好的鲁棒性。
上述两个测试表明本实施例方法能够保证三维模型的变换不变性和健壮性。
接下来再测试本实施例的匹配方法能否在一个更大、更多样的数据库中表现出较好的效果。先从互联网上下载430个不同种类的模型,它们被分为35个类。在表2中列出了部分类的名称和所包含的模型数量。
该数据库具有一些显著的特点。首先,每个类包含任意数量的三维模型,即为随机的来自互联网;第二,类之间的相似性有很大差异;第三,每个类中模型的相似性也有差异。
为了对比测试本实施例的三维模型检索方法(EFS)与基于体素化和球面调和的三维模型检索方法(VOXSH)(Funkhouser,T.,Min,P.,Kazhdan,M.,Chen,J.,Halderman,A.,Dobkin,D.,Jacobs,D.:A searchengine for 3D models.ACM Transactions on Graphics 22(2003)83-105)对模型分类和检索的效果和效率,本实施例进行了分类和检索实验:比较数据库中每个模型的描述符,并根据归类与查询模型最相近的所有模型的原则进行分类。表1列出了反映总体分类性能的各项实验结果。另外,本实施例还在每个类上做了性能测试,随机抽取16个基本类,将这些基本类的分类实验结果在表2中列出。同时,本实施例通过生成部分典型模型的检索结果和全部数据库的查准-查全图来做测试。如图2和图3所示,方框曲线表示本实施例的EFS方法的查准-查全曲线,三角形曲线表示VOXSH方法的查准-查全曲线,图2是EFS、VOXSH的总体查准-查全曲线比较图,图3是EFS、VOXSH在部分类上的查准-查全曲线比较图。
表1 EFS和VOXSH的总体性能比较
方法(Method) | 最近邻(NearestNeighbor) | 第一级(Firsttier) | 第二级(Secondtier) | E度量(E-Measure) | 折衷累积增益(DCG) | 匹配时间(MatchTime(ms)) |
EFS | 0.205 | 0.113 | 0.175 | 0.100 | 0.363 | 2.378 |
VOXSH | 0.095 | 0.056 | 0.094 | 0.065 | 0.306 | 4.061 |
表2 EFS和VOXSH在部分类上的性能比较
分类和它的数量(Class & Itsnumber) | 方法(Method) | 最近邻(NearestNeighbor) | 第一级(Firsttier) | 第二级(Secondtier) | E度量(E-Measure) | 折衷累积增益(DCG) | 匹配时间(MatchTime(ms)) |
冲翼艇(winged_vehicle)89 | EFS | 0.685 | 0.489 | 0.767 | 0.297 | 0.826 | 1.241625 |
VOXSH | 0.528 | 0.308 | 0.549 | 0.184 | 0.736 | 4.044447 | |
直升机(Helicopter)17 | EFS | 0.176 | 0.092 | 0.151 | 0.100 | 0.430 | 1.04487 |
VOXSH | 0.118 | 0.044 | 0.092 | 0.061 | 0.378 | 4.056771 | |
节肢动物(arthropod)11 | EFS | 0.273 | 0.155 | 0.236 | 0.182 | 0.462 | 1.321564 |
VOXSH | 0.000 | 0.055 | 0.164 | 0.113 | 0.344 | 4.053065 | |
四足动物(Quadruped)10 | EFS | 0.300 | 0.122 | 0.167 | 0.107 | 0.403 | 5.122093 |
VOXSH | 0.100 | 0.044 | 0.067 | 0.059 | 0.291 | 4.05186 | |
水下动物(underwater_creature) | EFS | 0.100 | 0.111 | 0.167 | 0.117 | 0.406 | 2.605814 |
VOXSH | 0.000 | 0.067 | 0.122 | 0.078 | 0.328 | 4.037209 |
10 | |||||||
头(head)9 | EFS | 0.444 | 0.153 | 0.208 | 0.094 | 0.396 | 3.418863 |
VOXSH | 0.667 | 0.222 | 0.250 | 0.117 | 0.485 | 4.049612 | |
桥梁(Bridge)10 | EFS | 0.200 | 0.044 | 0.056 | 0.039 | 0.297 | 1.646977 |
VOXSH | 0.000 | 0.011 | 0.033 | 0.044 | 0.264 | 4.059069 | |
建筑物(Building)26 | EFS | 0.154 | 0.074 | 0.172 | 0.084 | 0.439 | 1.84186 |
VOXSH | 0.000 | 0.080 | 0.148 | 0.090 | 0.421 | 4.04991 | |
座位(Seat)12 | EFS | 0.333 | 0.152 | 0.189 | 0.120 | 0.423 | 0.503682 |
VOXSH | 0.000 | 0.038 | 0.068 | 0.062 | 0.328 | 4.042054 | |
桌子(Table)12 | EFS | 0.250 | 0.045 | 0.068 | 0.050 | 0.336 | 1.37093 |
VOXSH | 0.000 | 0.030 | 0.053 | 0.035 | 0.297 | 4.051744 | |
枪械(Gun)8 | EFS | 0.000 | 0.179 | 0.304 | 0.141 | 0.376 | 5.048547 |
VOXSH | 0.000 | 0.036 | 0.054 | 0.083 | 0.292 | 4.052325 | |
容器(Container)26 | EFS | 0.577 | 0.266 | 0.446 | 0.286 | 0.641 | 0.23873 |
VOXSH | 0.308 | 0.123 | 0.180 | 0.123 | 0.502 | 4.074061 | |
植物(Plant)36 | EFS | 0.306 | 0.161 | 0.306 | 0.158 | 0.566 | 1.937984 |
VOXSH | 0.306 | 0.121 | 0.221 | 0.119 | 0.529 | 4.05969 | |
器皿(Vessel)13 | EFS | 0.000 | 0.109 | 0.141 | 0.098 | 0.391 | 2.20805 |
VOXSH | 0.000 | 0.064 | 0.141 | 0.094 | 0.356 | 4.069588 | |
手掌(Handheld)27 | EFS | 0.148 | 0.094 | 0.174 | 0.101 | 0.456 | 2.260723 |
VOXSH | 0.000 | 0.073 | 0.128 | 0.074 | 0.417 | 4.177519 | |
车(Car)25 | EFS | 0.640 | 0.275 | 0.362 | 0.269 | 0.621 | 8.963349 |
VOXSH | 0.160 | 0.098 | 0.185 | 0.106 | 0.468 | 3.995721 |
首先在数据库中检索到最相似的9个结果,所有的结果在大约17秒的时间返回。在该实验中,构造查询模型的特征描述符占用了绝大部分的检索时间,另外描述符的平均比较时间(2.38毫秒)比VOXSH方法的平均比较时间(4.06毫秒)快近1倍,如表1中匹配时间列所示,本实施例提出的EFS方法在各个类上的平均相似性度量时间绝大多数都比VOXSH方法更快。如此一来,当用户检索那些描述符已存在于特征库的模型时,速度将得到显著的加快(如从检索结果中进行二次检索)。
然后再对比VOXSH方法和本实施例提出的EFS方法在测试数据库中的分类性能指标、总体查准-查全曲线图和在部分类上的查准-查全曲线图。如表1所示,用于对比方法性能的折衷累积增益(DCG)指标表明,本实施例提出的EFS方法比VOXSH方法高出5.7个百分点,而且其余的各项指标都要优于VOXSH方法,其中最近邻(NearestNeighbor)高出了11个百分点。如图2所示,从总体查准-查全曲线比较图可以看出,代表本实施例EFS方法的方框曲线位于代表VOXSH方法的三角形曲线的上方。上述比较说明,本实施例方法的分类和检索效果优于VOXSH方法。其次,对各个类的分类和检索性能较好。为了分析方法对不同种类模型的适应性,表2列出了本实施例EFS方法和VOXSH方法在部分类上的性能统计,数据显示EFS方法各个类上的性能指标都明显优于VOXSH方法。另外,不同的方法对于不同类别的检索效果也有所不同,从图3可以看出,本实施例的EFS方法对于绝大多数类别的查准-查全率要高于VOXSH方法。从上述实验结果可以看出,本实施例提出的方法对于分类和检索该测试库中的三维模型,性能优于VOXSH方法。
以下再测试多级加权系数的影响,对两个模型采用基于Hausdorff距离的多级加权相似性匹配方法计算分层螺旋信息的相似度时,本实施例定义了三种权值,ω0、ω1和ω分别代表加权中轴点类型、中轴点的位置信息或距离信息以及是否同层。在此测试中,通过设定不同的权值进行了多次测试,测试实验结果表明,本发明方法的检索性能均优于VOXSH方法,并可以通过调整权值,满足对查准率或查全率的不同需要。
综上所述,本发明通过采用中轴点集分层螺旋信息作为三维形状的特征,并且提出一种基于Hausdorff距离的多级加权相似性匹配方法度量三维模型的相似性。本发明的关键是提供了一种框架,在该框架中,三维模型被转换为自然参数的函数,从而通过本发明提供的度量方法得到健壮高效和较好的三维模型相似性度量和检索结果。
Claims (2)
1.基于中轴点集分层螺旋信息的三维模型检索方法,其特征在于,包括以下步骤:
(1)对三维模型数据库中的每个三维模型进行坐标标准化预处理;
所述步骤(1)对三维模型数据库中的每个模型进行坐标标准化预处理,获取标准化后的三维模型,具体步骤如下:
a.确定三维模型包围盒的三个主轴;
b.确定三维模型包围盒的中心和三个半长;
c.确定三维模型包围盒主轴的正方向;
d.根据包围盒的主轴、中心和正方向确定三维模型的本征坐标系;
(2)对坐标标准化预处理后的除体元模型之外的三维模型进行体素化预处理;
所述步骤(2)对坐标标准化预处理后的除体元模型之外的三维模型进行体素化预处理,获取体素化的三维模型,具体步骤如下:
A.面片顶点体素化;
B.面片边的体素化;
C.面片的体素化;
D.三维模型内部体素化;
(3)对每个体素化的三维模型进行特征提取,得到每个三维模型的特征描述符,生成特征数据库;
所述步骤(3)对体素化的三维模型进行特征提取,得到每个三维 模型的特征描述符,生成特征数据库,具体步骤如下:
a.对每个体素化的三维模型通过排斥力场函数提取模型的中轴点;
b.对每个体素化的三维模型建立中轴点集分层螺旋信息特征描述符,生成特征数据库;
对于体素化的三维模型,首先通过排斥力场函数提取三维模型的中轴点集,具体步骤如下:
首先假设边界体元为点电荷,边界点被定义为一个物体体元距离坐标原点最近的立方体的角,内部点是在每个面上都邻接三维物体体元的模型体格中距离坐标原点最近的角,坐标原点被定义为三维模型包围盒的一个角,一个内部点的排斥力可以通过式(1)计算:
如果力矢量(x,y,z)的三个分量都消失,矢量场就为零,也就是说,当某个区域中每一个矢量分量的符号都改变时,该区域就是包含临界点的候选区域,具体的说,首先假设体单元是最小的区域,并计算格子单元八个角的力场值,对于该单元所有的矢量分量(x,y,z),如果既有正值又有负值,那么该单元就是潜在的候选单元,然后再将候选单元递归地分为八个子单元,并且对每个子单元重复做上述的候选性测试,其中,定义了下述规则:如果有两个或两个以上的候 选子单元,则比较它们的力值并且选择具有最小力值的子单元替换原先的候选单元作为新的候选单元,处理过程终止的条件是:该单元未能通过候选性测试,或者该单元太小,此时该单元仍然作为候选单元保留,如果是后一种情况,临界点的坐标就被定义为该单元的坐标。但是如果所有的八个子单元都未能通过候选性测试,定义该单元的坐标为候选临界点,把从上述算法所提取的临界点和候选临界点称为中轴点;
(4)将用户提交的三维模型首先按照步骤(1)中的方法进行坐标标准化预处理、然后按照步骤(2)中的方法体素化其三维模型,再按照步骤(3)中的方法进行特征提取得到该三维模型的特征描述符;
(5)将用户提交的三维模型特征描述符与特征数据库中的特征进行匹配,计算相似距离并排序,最终生成检索结果。
2.根据权利要求1所述的基于中轴点集分层螺旋信息的三维模型检索方法,其特征在于:所述步骤(5)将用户提交的三维模型特征描述符与特征数据库中的特征进行匹配,计算相似距离并排序,最终生成检索结果,具体步骤如下:
(1)采用一种基于Hausdorff距离的多级加权相似性匹配方法对用户提交的查询模型特征和特征数据库中的每一个三维模型特征进行匹配,将特征匹配得到的所有数据依据距离大小进行相似度排序;
首先定义如下几个概念;
定义中轴点加权信息
在中轴点信息的基础上,定义中轴点加权信息为:
定义中轴点加权信息的距离
一种基于Hausdorff距离的匹配器匹配用户提交的三维模型A的特征描述符与特征数据库中的三维模型B的特征的步骤为:
第一步,根据中轴点的类型对其赋予不同的权值,如式(2),若为临界点,权值为ω0;若为候选临界点,权值为1-ω0,亦为第1级加权,这样,A和B各自描述为加权的中轴点集分层螺旋信息A=[a0,...,am-1]-1和B=[b0,...,bn-1]-1,其中,au和bv是中轴点加权信息,u∈[1,m],v∈[1,n],
第二步,分别计算A和B的所有中轴点加权信息中 和 的加权Euclidean距离,权值ω1调节坐标加权信息的Euclidean距离和距离加权信息的Euclidean距离所占的比重,亦为第2级加权,另外,根据式(3)计算两个中轴点的距离,判断两个中轴点是否同层,用对应的权值ω或1-ω加权这两个中轴点加权信息的距离,即第3级加 权,A和B的任意两个中轴点加权信息au和bv之间的距离为d(au,bv),
第三步,用Hausdorff距离度量模型A和模型B的差异值H(A,B):
H(A,B)=max{h(A,B),h(B,A)}
(2)将用户提交的查询模型和模型库中所有模型的特征按照上述匹配步骤计算,得出所有的特征匹配距离,然后将这些距离排序后按照相似度大小生成检索结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101950618A CN101201845B (zh) | 2007-09-24 | 2007-12-11 | 基于中轴点集分层螺旋信息的三维模型检索方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710122370 | 2007-09-24 | ||
CN200710122370.2 | 2007-09-24 | ||
CN2007101950618A CN101201845B (zh) | 2007-09-24 | 2007-12-11 | 基于中轴点集分层螺旋信息的三维模型检索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101201845A CN101201845A (zh) | 2008-06-18 |
CN101201845B true CN101201845B (zh) | 2012-01-11 |
Family
ID=39517012
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101950618A Expired - Fee Related CN101201845B (zh) | 2007-09-24 | 2007-12-11 | 基于中轴点集分层螺旋信息的三维模型检索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101201845B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103262120B (zh) * | 2010-12-09 | 2017-03-22 | 皇家飞利浦电子股份有限公司 | 图像数据的体积描绘 |
CN102495926B (zh) * | 2011-12-01 | 2014-06-25 | 北京数码大方科技股份有限公司 | 三维原始模型的检验方法及装置 |
CN102693285A (zh) * | 2012-05-11 | 2012-09-26 | 北京师范大学 | 一种基于形状认知的三维建筑模型匹配与检索的方法 |
CN104462163B (zh) * | 2014-03-06 | 2018-01-16 | 北京工商大学 | 一种三维模型表征方法、检索方法及检索系统 |
CN105354593B (zh) * | 2015-10-22 | 2018-10-30 | 南京大学 | 一种基于nmf的三维模型分类方法 |
CN106650747B (zh) * | 2016-11-16 | 2019-08-02 | 佛山科学技术学院 | 一种基于压缩感知的三维模型特征提取方法 |
CN107688599B (zh) * | 2017-07-05 | 2018-11-16 | 国网山东省电力公司 | 一种快速检索三维模型的方法 |
CN107578470B (zh) * | 2017-10-16 | 2020-11-06 | 北京锐安科技有限公司 | 一种三维模型处理方法、装置、服务器及存储介质 |
CN108921045B (zh) * | 2018-06-11 | 2021-08-03 | 佛山科学技术学院 | 一种三维模型的空间特征提取和匹配方法及其装置 |
CN112001044B (zh) * | 2020-07-31 | 2024-04-02 | 宁波智讯联科科技有限公司 | 零件cad模型相似度估算方法与系统 |
CN113155453B (zh) * | 2021-01-08 | 2023-11-07 | 沈阳大学 | 一种多特征融合评估预测轴承性能寿命方法 |
CN114055781B (zh) * | 2021-10-24 | 2023-12-29 | 扬州大学 | 基于点体素相关场的燃油箱焊接机械臂自适应校正方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1527251A (zh) * | 2002-12-05 | 2004-09-08 | ���ǵ�����ʽ���� | 用感性三维形状描述检索三维图形模型数据库的方法和设备 |
EP1530157A1 (en) * | 2002-07-10 | 2005-05-11 | NEC Corporation | Image matching system using 3-dimensional object model, image matching method, and image matching program |
CN101004748A (zh) * | 2006-10-27 | 2007-07-25 | 北京航空航天大学 | 基于二维草图的三维模型检索方法 |
-
2007
- 2007-12-11 CN CN2007101950618A patent/CN101201845B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1530157A1 (en) * | 2002-07-10 | 2005-05-11 | NEC Corporation | Image matching system using 3-dimensional object model, image matching method, and image matching program |
CN1527251A (zh) * | 2002-12-05 | 2004-09-08 | ���ǵ�����ʽ���� | 用感性三维形状描述检索三维图形模型数据库的方法和设备 |
CN101004748A (zh) * | 2006-10-27 | 2007-07-25 | 北京航空航天大学 | 基于二维草图的三维模型检索方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101201845A (zh) | 2008-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101201845B (zh) | 基于中轴点集分层螺旋信息的三维模型检索方法 | |
Liu et al. | Construction of iso-contours, bisectors, and Voronoi diagrams on triangulated surfaces | |
Caelli et al. | An eigenspace projection clustering method for inexact graph matching | |
US8015125B2 (en) | Multi-scale segmentation and partial matching 3D models | |
Zhong | Intrinsic shape signatures: A shape descriptor for 3D object recognition | |
Iyer et al. | Three-dimensional shape searching: state-of-the-art review and future trends | |
Yang et al. | Content-based 3-D model retrieval: A survey | |
Kazmi et al. | A survey of 2D and 3D shape descriptors | |
Papadakis et al. | PANORAMA: A 3D shape descriptor based on panoramic views for unsupervised 3D object retrieval | |
JP4516957B2 (ja) | 3次元オブジェクトについて検索を行なうための方法、システムおよびデータ構造 | |
CN103914571B (zh) | 一种基于网格分割的三维模型检索方法 | |
Bespalov et al. | Scale-space representation of 3d models and topological matching | |
Qie et al. | Enhanced invariance class partitioning using discrete curvatures and conformal geometry | |
CN111597367B (zh) | 基于视图和哈希算法的三维模型检索方法 | |
Wang et al. | Geometric decomposition of 3D surface meshes using Morse theory and region growing | |
Boyer et al. | Saliency sequential surface organization for free-form object recognition | |
Adán et al. | Modeling wave set: Definition and application of a new topological organization for 3d object modeling | |
Greenspan | Geometric probing of dense range data | |
Memon et al. | 3D shape retrieval using bag of word approaches | |
Jiang et al. | Partial shape matching of 3D models based on the Laplace-Beltrami operator eigenfunction | |
Mahmoudi et al. | 3D objects retrieval using curvature scale space and Zernike moments | |
Martin et al. | Turning shape decision problems into measures | |
Mademlis et al. | Ellipsoidal harmonics for 3-D shape description and retrieval | |
Wang et al. | Shape complexity from image similarity | |
Otu | Methods for 3D Shape Description, Indexing, Matching and Retrieval |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120111 Termination date: 20211211 |