CN103065353A - 三维模型特征提取方法和系统、三维模型检索方法和系统 - Google Patents

三维模型特征提取方法和系统、三维模型检索方法和系统 Download PDF

Info

Publication number
CN103065353A
CN103065353A CN2012105630431A CN201210563043A CN103065353A CN 103065353 A CN103065353 A CN 103065353A CN 2012105630431 A CN2012105630431 A CN 2012105630431A CN 201210563043 A CN201210563043 A CN 201210563043A CN 103065353 A CN103065353 A CN 103065353A
Authority
CN
China
Prior art keywords
voxel
model
distribution characteristics
coordinate
space
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.)
Granted
Application number
CN2012105630431A
Other languages
English (en)
Other versions
CN103065353B (zh
Inventor
张宝运
李晓明
胡金星
修文群
彭大为
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhuhai Zhongke advanced technology industry Co.,Ltd.
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201210563043.1A priority Critical patent/CN103065353B/zh
Publication of CN103065353A publication Critical patent/CN103065353A/zh
Application granted granted Critical
Publication of CN103065353B publication Critical patent/CN103065353B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Generation (AREA)

Abstract

一种三维模型特征提取方法,包括以下步骤:对三维模型进行体素化处理,生成体素模型;获取体素模型的重心;计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离;生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,所述同心球面的球心为所述重心;计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值;根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。此外,还提供一种三维模型特征提取系统。上述方法和系统提取的特征的精度不受三维模型的缩放和旋转的影响。

Description

三维模型特征提取方法和系统、三维模型检索方法和系统
【技术领域】
本发明涉及三维模型处理领域,特别涉及一种三维模型特征提取方法和系统,以及一种三维模型检索方法和系统。
【背景技术】
如今,三维模型的应用的非常广泛,工业产品设计、虚拟现实、虚拟人、三维游戏、教育、影视动画等都广泛使用三维模型。
为了便于进行三维模型的数字化存储、三维模型的匹配检索等处理,需要提取三维模型的特征。传统的方法提取三维模型的轮廓特征,而轮廓特征相对比较简单,无法得到模型的总体形状特征,而且同一三维模型在不同的缩放比例或旋转角度状态下,所提取到的轮廓特征有可能是不一样的,即提取的轮廓特征会受到三维模型的缩放和旋转的影响。
【发明内容】
基于此,有必要提供一种提取的特征不受三维模型的缩放和旋转影响的三维模型特征提取方法和系统。
一种三维模型特征提取方法,包括以下步骤:
对三维模型进行体素化处理,生成体素模型;
获取体素模型的重心;
计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离;
生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,所述同心球面的球心为所述重心;
计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值;
根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
在其中一个实施例中,所述对三维模型进行体素化处理,生成体素模型的步骤包括:
获取三维模型的网格模型;
对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部;
获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值;
生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值;
将生成的包围盒均匀切割成预设大小的体素;
根据网格模型的网格向量获取网格模型内部的体素,并由所述网格模型内部的体素组成体素模型。
在其中一个实施例中,根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征的步骤包括:
根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(r1,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量;
根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,frN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
一种三维模型特征提取系统,包括:
体素化模块,用于对三维模型进行体素化处理,生成体素模型;
重心获取模块,用于获取体素模型的重心;
体素距离计算模块,用于计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离;
区域划分模块,用于生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,所述同心球面的球心为所述重心;
体素数量获取模块,用于计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值;
分布特征计算模块,用于根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
在其中一个实施例中,所述体素化模块包括:
网格模型获取单元,用于获取三维模型的网格模型;
平移变换单元,用于对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部;
坐标极值获取单元,用于获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值;
包围盒生成单元,用于生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值;
体素切割单元,用于将生成的包围盒均匀切割成预设大小的体素;
体素模型生成单元,用于根据网格模型的网格向量获取网格模型内部的体素,并由所述网格模型内部的体素组成体素模型。
在其中一个实施例中,所述分布特征计算模块包括:
相对密度计算单元,用于根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(r1,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量;
数值差分计算单元,用于根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,drN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
分布特征获取单元,用于获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
上述三维模型特征提取方法和系统,统计三维模型的体素模型中体素在不同空间区域的分布情况,得到体素模型的体素密度分布特征,由于三维模型的缩放和旋转不影响三维模型的体素在空间区域中的分布情况,因此,上述三维模型特征提取方法和系统提取的特征的精度不受三维模型的缩放和旋转的影响。另外,两个形状相似的三维模型的体素模型中体素在空间区域中的分布也是相似的,因此,上述方法提取的特征能准确反映三维模型的总体形状。而且,上述方法和系统提取的特征的数据量小、便于进行存储和计算处理。
一种三维模型检索方法,包括以下步骤:
提取待检索模型的体素密度分布特征;
计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离;
获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型;
其中,所述提取待检索模型的体素密度分布特征的步骤按照上述任一实施例中的三维模型特征提取方法提取待检索模型的体素密度分布特征。
在其中一个实施例中,在提取待检索模型的体素密度分布特征的步骤之前,所述方法还包括:
按照所述提取待检索模型的体素密度分布特征的方法提取三维模型库中三维模型的体素密度分布特征,将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
在其中一个实施例中,计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离的步骤包括:
分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离;
将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
一种三维模型检索系统,包括:
体素密度分布特征提取模块,用于提取待检索模型的体素密度分布特征,所述体素密度分布特征提取模块包括权利要求7至9任一所述的三维模型特征提取系统中的模块;
分布特征距离计算模块,用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离;
匹配模型获取模块,用于获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型。
在其中一个实施例中,所述体素密度分布特征提取模块还用于提取三维模型库中三维模型的体素密度分布特征;
所述系统还包括体素密度分布特征库生成模块,用于将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
在其中一个实施例中,所述分布特征距离计算模块包括:
欧氏距离计算单元,用于分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离;
归一化单元,用于将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
分布特征距离计算单元,用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
上述三维模型检索方法和系统,根据提取三维模型的体素密度分布特征来进行匹配检索,由于体素密度分布特征能准确反映三维模型的总体形状且数据量小,因此,上述方法和系统能快速准确地检索到与待检索三维模型最匹配的三维模型,提取检索的速度的准确度。
【附图说明】
图1为一个实施例中的三维模型特征提取方法的流程示意图;
图2为一个实施例中图1的步骤S10的流程示意图;
图3为将体素模型划分空间区域的示意图;
图4为一个实施例中的三维模型检索方法的流程示意图;
图5为一个实施例中图4的步骤S80的流程示意示意图;
图6为一个实施例中的三维模型特征提取系统的结构示意图;
图7为一个实施例中体素化模块的结构示意图;
图8为一个实施例中分布特征计算模块的结构示意图;
图9为一个实施例中的三维模型检索系统的结构示意图;
图10为一个实施例中分布特征距离计算模块的结构示意图。
【具体实施方式】
如图1所示,在一个实施例中,一种三维模型特征提取方法,包括以下步骤:
步骤S10,对三维模型进行体素化处理,生成体素模型。
如图2所示,在一个实施例中,步骤S10包括:
步骤S101,获取三维模型的网格模型。
具体的,可按照预设的网格模型生成方法生成三维模型的网络模型,在此不作赘述。
步骤S102,对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部。
优选的,可将网络模型的中心平移到坐标系原点。
步骤S103,获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值。
步骤S104,生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值。
步骤S105,将生成的包围盒均匀切割成预设大小的体素。
步骤S106,根据网格模型的网格向量获取网格模型内部的体素,并由网格模型内部的体素组成体素模型。
步骤S20,获取体素模型的重心。
具体的,可计算体素模型中所有体素的x坐标值的平均值、y坐标值的平均值、z坐标值的平均值,坐标(x坐标值的平均值,y坐标值的平均值,z坐标值的平均值)对应的点即为体素模型的重心。
步骤S30,计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离。
步骤S40,生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,同心球面的球心为体素模型的重心。
图3为将体素模型划分空间区域的示意图。如图3所示,将球体的一条半径等分为N段,分别以球体的球心为球心、以球心到等分点的距离为半径作球面,可作出N-1个球面,这N-1个球面将球体划分为N个空间区域。其中,N为预设值。
步骤S50,计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值。
步骤S60,根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
在一个实施例中,步骤S60包括以下步骤:
根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(r1,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量;
进一步的,根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,drN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
进一步的,获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
上述三维模型特征提取方法,统计三维模型的体素模型中体素在不同空间区域的分布情况,得到体素模型的体素密度分布特征,由于三维模型的缩放和旋转不影响三维模型的体素在空间区域中的分布情况,因此,上述三维模型特征提取方法提取的特征的精度不受三维模型的缩放和旋转的影响。另外,两个形状相似的三维模型的体素模型中体素在空间区域中的分布也是相似的,因此,上述方法提取的特征能准确反映三维模型的总体形状。而且,上述方法提取的特征的数据量小、便于进行存储和计算处理。
如图4所示,在一个实施例中,一种三维模型检索方法,包括以下步骤:
步骤S70,提取待检索模型的体素密度分布特征。
具体的,可按照上述任一实施例中的三维模型特征提取方法提取待检索模型的体素密度分布特征。
步骤S80,计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离。
如图5所示,在一个实施例中,步骤S80包括以下步骤:
步骤S801,分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离。
步骤S802,将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
步骤S803,计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
例如,两个体素密度分布特征的归一化后的相对密度向量的欧氏距离为dr,归一化后的一阶数值差分向量的欧氏距离为ddr,则可计算两个体素密度分布特征的距离为a*dr+b*ddr,其中,a和b为预设值,a与b的和为1。
步骤S90,获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型。
步骤S90获取的三维模型即为与待检索模型最匹配的三维模型。
在一个实施例中,在步骤S70之前,上述三维模型检索方法还包括步骤:按照上述提取待检索模型的体素密度分布特征的方法提取三维模型库中三维模型的体素密度分布特征,将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
上述三维模型检索方法,根据提取三维模型的体素密度分布特征来进行匹配检索,由于体素密度分布特征能准确反映三维模型的总体形状且数据量小,因此,上述方法能快速准确地检索到与待检索三维模型最匹配的三维模型,提取检索的速度的准确度。
如图6所示,在一个实施例中,一种三维模型特征提取系统,包括体素化模块10、重心获取模块20、体素距离计算模块30、区域划分模块40、体素数量获取模块50和分布特征计算模块60,其中:
体素化模块10用于对三维模型进行体素化处理,生成体素模型。
如图7所示,在一个实施例中,体素化模块10包括网格模型获取单元101、平移变换单元102、坐标极值获取单元103、包围盒生成单元104、体素切割单元105和体素模型生成单元106,其中:
网格模型获取单元101用于获取三维模型的网格模型。
具体的,网格模型获取单元101可按照预设的网格模型生成方法生成三维模型的网络模型,在此不作赘述。
平移变换单元102用于对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部。
优选的,网格模型获取单元101可将网络模型的中心平移到坐标系原点。
坐标极值获取单元103用于获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值。
包围盒生成单元104用于生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值。
体素切割单元105用于将生成的包围盒均匀切割成预设大小的体素。
体素模型生成单元106用于根据网格模型的网格向量获取网格模型内部的体素,并由所述网格模型内部的体素组成体素模型。
重心获取模块20用于获取体素模型的重心。
具体的,重心获取模块20可计算体素模型中所有体素的x坐标值的平均值、y坐标值的平均值、z坐标值的平均值,坐标(x坐标值的平均值,y坐标值的平均值,z坐标值的平均值)对应的点即为体素模型的重心。
体素距离计算模块30用于计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离。
区域划分模块40用于生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,同心球面的球心为体素模型的重心。
区域划分模块40可将球体的一条半径等分为N段,分别以球体的球心为球心、以球心到等分点的距离为半径作球面,可作出N-1个球面,这N-1个球面将球体划分为N个空间区域。其中,N为预设值。
体素数量获取模块50用于计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值。
分布特征计算模块60用于根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
如图8所示,在一个实施例中,分布特征计算模块60包括相对密度计算单元601、数值差分计算单元602和分布特征获取单元603,其中:
相对密度计算单元601用于根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(r1,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量。
数值差分计算单元602用于根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,drN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
分布特征获取单元603用于获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
上述三维模型特征提取系统,统计三维模型的体素模型中体素在不同空间区域的分布情况,得到体素模型的体素密度分布特征,由于三维模型的缩放和旋转不影响三维模型的体素在空间区域中的分布情况,因此,上述三维模型特征提取系统提取的特征的精度不受三维模型的缩放和旋转的影响。另外,两个形状相似的三维模型的体素模型中体素在空间区域中的分布也是相似的,因此,上述系统提取的特征能准确反映三维模型的总体形状。而且,上述系统提取的特征的数据量小、便于进行存储和计算处理。
如图9所示,在一个实施例中,一种三维模型检索系统,包括体素密度分布特征提取模块70、分布特征距离计算模块80和匹配模型获取模块90,其中:
体素密度分布特征提取模块70用于提取待检索模型的体素密度分布特征。具体的,体素密度分布特征提取模块70包括上述任一实施例中的三维模型特征提取系统中的模块。
分布特征距离计算模块80用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离。
如图10所示,在一个实施例中,分布特征距离计算模块80包括欧氏距离计算单元801、归一化单元802和分布特征距离计算单元803,其中:
欧氏距离计算单元801用于分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离;
归一化单元802用于将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
分布特征距离计算单元803用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
例如,两个体素密度分布特征的归一化后的相对密度向量的欧氏距离为dr,归一化后的一阶数值差分向量的欧氏距离为ddr,则分布特征距离计算单元803可计算两个体素密度分布特征的距离为a*dr+b*ddr,其中,a和b为预设值,a与b的和为1。
匹配模型获取模块90用于获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型。
匹配模型获取模块90获取的三维模型即为与待检索模型最匹配的三维模型。
在一个实施例中,体素密度分布特征提取模块70还用于提取三维模型库中三维模型的体素密度分布特征;上述三维模型检索系统还包括体素密度分布特征库生成模块(图中未示出),用于将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
上述三维模型检索系统,根据提取三维模型的体素密度分布特征来进行匹配检索,由于体素密度分布特征能准确反映三维模型的总体形状且数据量小,因此,上述系统能快速准确地检索到与待检索三维模型最匹配的三维模型,提取检索的速度的准确度。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (12)

1.一种三维模型特征提取方法,包括以下步骤:
对三维模型进行体素化处理,生成体素模型;
获取体素模型的重心;
计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离;
生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,所述同心球面的球心为所述重心;
计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值;
根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
2.根据权利要求1所述的三维模型特征提取方法,其特征在于,所述对三维模型进行体素化处理,生成体素模型的步骤包括:
获取三维模型的网格模型;
对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部;
获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值;
生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值;
将生成的包围盒均匀切割成预设大小的体素;
根据网格模型的网格向量获取网格模型内部的体素,并由所述网格模型内部的体素组成体素模型。
3.根据权利要求1所述的三维模型特征提取方法,其特征在于,根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征的步骤包括:
根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(ri,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量;
根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,drN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
4.一种三维模型检索方法,包括以下步骤:
提取待检索模型的体素密度分布特征;
计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离;
获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型;
其中,所述提取待检索模型的体素密度分布特征的步骤按照权利要求1至3中任一所述的三维模型特征提取方法提取待检索模型的体素密度分布特征。
5.根据权利要求4所述的三维模型检索方法,其特征在于,在提取待检索模型的体素密度分布特征的步骤之前,所述方法还包括:
按照所述提取待检索模型的体素密度分布特征的方法提取三维模型库中三维模型的体素密度分布特征,将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
6.根据权利要求4所述的三维模型检索方法,其特征在于,计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离的步骤包括:
分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离;
将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
7.一种三维模型特征提取系统,其特征在于,包括:
体素化模块,用于对三维模型进行体素化处理,生成体素模型;
重心获取模块,用于获取体素模型的重心;
体素距离计算模块,用于计算体素模型中每个体素到重心的距离,获取计算得到的距离中的最大距离;
区域划分模块,用于生成以重心为球心、以最大距离为半径的球体,将球体用等距离的预设数量个同心球面划分为多个空间区域,所述同心球面的球心为所述重心;
体素数量获取模块,用于计算各空间区域中体素的数量,获取各空间区域中体素数量的最大值;
分布特征计算模块,用于根据各空间区域中体素的数量以及各空间区域中体素数量的最大值计算体素模型的体素密度分布特征。
8.根据权利要求7所述的三维模型特征提取系统,其特征在于,所述体素化模块包括:
网格模型获取单元,用于获取三维模型的网格模型;
平移变换单元,用于对网格模型进行平移变换,使得平移变换后坐标系原点位于网格模型内部;
坐标极值获取单元,用于获取网格模型的所有顶点的x坐标值中的最大x坐标值和最小x坐标值、y坐标值中的最大y坐标值和最小y坐标值、z坐标值中的最大z坐标值和最小z坐标值;
包围盒生成单元,用于生成以第一点和第二点为对角线的包围盒,所述第一点的x坐标为最大x坐标值、y坐标为最大y坐标值、z坐标为最大z坐标值,所述第二点的组成的x坐标为最小x坐标值、y坐标为最小y坐标值、z坐标为最小z坐标值;
体素切割单元,用于将生成的包围盒均匀切割成预设大小的体素;
体素模型生成单元,用于根据网格模型的网格向量获取网格模型内部的体素,并由所述网格模型内部的体素组成体素模型。
9.根据权利要求7所述的三维模型特征提取系统,其特征在于,所述分布特征计算模块包括:
相对密度计算单元,用于根据以下公式计算各空间区域中体素的相对密度,得到体素密度分布特征的相对密度向量(r1,...,rN):
r i = S i S max ( i = 1 , . . . , N ) ,
其中,ri为划分的空间区域中第i个空间区域的体素的相对密度,Si为第i个空间区域中体素的数量,Smax为划分的空间区域的体素数量的最大值,N为划分的空间区域的数量;
数值差分计算单元,用于根据以下公式计算各空间区域的体素在极坐标方向上的一阶数值差分,得到体素密度分布特征的一阶数值差分向量(dri,...,drN):
dr i = | r 1 - r 2 | , i = 1 | 2 r 1 - r i - 1 - r i + 1 | , 1 < i < N | r N - r N - 1 | , i = N
其中,dri为划分的空间区域中第i个空间区域的体素在极坐标方向上的一阶数值差分,ri为上述计算得到的第i个空间区域的体素的相对密度,N为划分的空间区域的数量;
分布特征获取单元,用于获取体素密度分布特征为:(r1,...,rN;dri,...,drN),其中,r1,...,rN分别为上述计算得到的第1,…,N个空间区域的体素的相对密度,dri,...,drN分别为上述计算得到的第1,…,N个空间区域的体素在极坐标方向上的一阶数值差分。
10.一种三维模型检索系统,其特征在于,包括:
体素密度分布特征提取模块,用于提取待检索模型的体素密度分布特征,所述体素密度分布特征提取模块包括权利要求7至9任一所述的三维模型特征提取系统中的模块;
分布特征距离计算模块,用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离;
匹配模型获取模块,用于获取三维模型体素密度分布特征库中与待检索模型的体素密度分布特征距离最小的体素密度分布特征对应的三维模型。
11.根据权利要求10所述的三维模型检索系统,其特征在于,所述体素密度分布特征提取模块还用于提取三维模型库中三维模型的体素密度分布特征;
所述系统还包括体素密度分布特征库生成模块,用于将提取的三维模型的体素密度分布特征组成三维模型体素密度分布特征库。
12.根据权利要求10所述的三维模型检索系统,其特征在于,所述分布特征距离计算模块包括:
欧氏距离计算单元,用于分别计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离;
归一化单元,用于将计算得到的相对密度向量的欧氏距离和一阶数值差分向量的欧氏距离采用高斯模型进行归一化;
分布特征距离计算单元,用于计算待检索模型的体素密度分布特征与三维模型体素密度分布特征库中的体素密度分布特征的距离为归一化后的相对密度向量的欧氏距离与一阶数值差分向量的欧氏距离的加权值,相对密度向量的欧氏距离的权值系数与一阶数值差分向量的欧氏距离的权值系统之和等于一。
CN201210563043.1A 2012-12-22 2012-12-22 三维模型特征提取方法和系统、三维模型检索方法和系统 Active CN103065353B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210563043.1A CN103065353B (zh) 2012-12-22 2012-12-22 三维模型特征提取方法和系统、三维模型检索方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210563043.1A CN103065353B (zh) 2012-12-22 2012-12-22 三维模型特征提取方法和系统、三维模型检索方法和系统

Publications (2)

Publication Number Publication Date
CN103065353A true CN103065353A (zh) 2013-04-24
CN103065353B CN103065353B (zh) 2015-09-09

Family

ID=48107969

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210563043.1A Active CN103065353B (zh) 2012-12-22 2012-12-22 三维模型特征提取方法和系统、三维模型检索方法和系统

Country Status (1)

Country Link
CN (1) CN103065353B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617603A (zh) * 2013-12-06 2014-03-05 南京大学 一种三维数字几何网格模型结构的自动修复方法
CN103914571A (zh) * 2014-04-25 2014-07-09 南京大学 一种基于网格分割的三维模型检索方法
CN105787445A (zh) * 2016-02-24 2016-07-20 武汉迈步科技有限公司 一种车载激光扫描数据中杆状物自动提取方法及系统
CN106126085A (zh) * 2016-06-13 2016-11-16 华南理工大学 一种面向移动设备的多点触控手势交互三维建模方法
CN106599053A (zh) * 2016-11-16 2017-04-26 佛山科学技术学院 一种三维模型检索方法
CN106650747A (zh) * 2016-11-16 2017-05-10 佛山科学技术学院 一种基于压缩感知的三维模型特征提取方法
CN107154076A (zh) * 2017-06-02 2017-09-12 王征 一种可动态自定义3d模型尺寸的系统和方法
CN108921045A (zh) * 2018-06-11 2018-11-30 佛山科学技术学院 一种三维模型的空间特征提取和匹配方法及其装置
CN109814125A (zh) * 2017-11-21 2019-05-28 通用汽车环球科技运作有限责任公司 用于确定激光雷达点的速度的系统和方法
CN109918524A (zh) * 2019-01-16 2019-06-21 中德(珠海)人工智能研究院有限公司 一种三维模型检索方法和系统
CN110889903A (zh) * 2019-11-26 2020-03-17 南京理工大学 一种基于切层策略的三维模型快速体素化方法
CN111080518A (zh) * 2019-11-27 2020-04-28 贝壳技术有限公司 点云图像的拼接方法、装置及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477529A (zh) * 2008-12-01 2009-07-08 清华大学 一种三维对象的检索方法和装置
CN102110159A (zh) * 2011-02-21 2011-06-29 清华大学 一种cad三维模型检索方法及系统
CN102254353A (zh) * 2011-06-10 2011-11-23 南京大学 一种三维模型的特征描述方法
CN102375831A (zh) * 2010-08-13 2012-03-14 富士通株式会社 三维模型检索装置及其方法和模型库生成装置及其方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477529A (zh) * 2008-12-01 2009-07-08 清华大学 一种三维对象的检索方法和装置
CN102375831A (zh) * 2010-08-13 2012-03-14 富士通株式会社 三维模型检索装置及其方法和模型库生成装置及其方法
CN102110159A (zh) * 2011-02-21 2011-06-29 清华大学 一种cad三维模型检索方法及系统
CN102254353A (zh) * 2011-06-10 2011-11-23 南京大学 一种三维模型的特征描述方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐鹏捷: "基于内容的三维模型检索技术", 《中国优秀硕士学位论文全文数据库》 *
王枫: "球面调和在模型检索中的应用", 《中国优秀硕士学位论文全文数据库》 *

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617603A (zh) * 2013-12-06 2014-03-05 南京大学 一种三维数字几何网格模型结构的自动修复方法
CN103617603B (zh) * 2013-12-06 2016-02-24 南京大学 一种三维数字几何网格模型结构的自动修复方法
CN103914571A (zh) * 2014-04-25 2014-07-09 南京大学 一种基于网格分割的三维模型检索方法
CN103914571B (zh) * 2014-04-25 2017-02-15 南京大学 一种基于网格分割的三维模型检索方法
CN105787445A (zh) * 2016-02-24 2016-07-20 武汉迈步科技有限公司 一种车载激光扫描数据中杆状物自动提取方法及系统
CN106126085A (zh) * 2016-06-13 2016-11-16 华南理工大学 一种面向移动设备的多点触控手势交互三维建模方法
CN106650747B (zh) * 2016-11-16 2019-08-02 佛山科学技术学院 一种基于压缩感知的三维模型特征提取方法
CN106650747A (zh) * 2016-11-16 2017-05-10 佛山科学技术学院 一种基于压缩感知的三维模型特征提取方法
CN106599053A (zh) * 2016-11-16 2017-04-26 佛山科学技术学院 一种三维模型检索方法
CN106599053B (zh) * 2016-11-16 2020-02-04 佛山科学技术学院 一种三维模型检索方法
CN107154076A (zh) * 2017-06-02 2017-09-12 王征 一种可动态自定义3d模型尺寸的系统和方法
CN109814125A (zh) * 2017-11-21 2019-05-28 通用汽车环球科技运作有限责任公司 用于确定激光雷达点的速度的系统和方法
CN108921045A (zh) * 2018-06-11 2018-11-30 佛山科学技术学院 一种三维模型的空间特征提取和匹配方法及其装置
CN109918524A (zh) * 2019-01-16 2019-06-21 中德(珠海)人工智能研究院有限公司 一种三维模型检索方法和系统
CN110889903A (zh) * 2019-11-26 2020-03-17 南京理工大学 一种基于切层策略的三维模型快速体素化方法
CN110889903B (zh) * 2019-11-26 2020-10-30 南京理工大学 一种基于切层策略的三维模型快速体素化方法
CN111080518A (zh) * 2019-11-27 2020-04-28 贝壳技术有限公司 点云图像的拼接方法、装置及存储介质
CN111080518B (zh) * 2019-11-27 2023-10-20 贝壳技术有限公司 点云图像的拼接方法、装置及存储介质

Also Published As

Publication number Publication date
CN103065353B (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN103065353B (zh) 三维模型特征提取方法和系统、三维模型检索方法和系统
CN102855663B (zh) 一种根据激光雷达栅格点云构建csg模型的方法
CN102254343B (zh) 一种基于凸包与obb的三维网格模型骨架提取方法
CN101872488B (zh) 曲面渲染系统及方法
CN103247041A (zh) 一种基于局部采样的多几何特征点云数据的分割方法
CN101526611B (zh) 针对舰船的雷达散射截面检测方法
CN102609992A (zh) 基于三角网格变形体的自碰撞检测方法
CN101084498A (zh) 三维模型的检索方法和装置
CN103839292A (zh) 对等值面进行采样并生成高质量三角网格的方法
CN107680168A (zh) 三维重建中基于平面拟合的网格简化方法
Weller et al. A brief overview of collision detection
CN105160344A (zh) 一种三维点云的局部特征提取方法及装置
CN102567465B (zh) 复眼视觉基础上的三维模型几何相似性搜索方法
CN103473811A (zh) 基于二维手绘线画图的三维实体模型便捷生成方法
CN103279621A (zh) 一种装配体的自动拆卸方法
CN104299255A (zh) 一种三维地形模型的渲染方法
CN105785919B (zh) 点云五轴无全局干涉刀轴矢量可行域计算方法
Jiang et al. Efficient Booleans algorithms for triangulated meshes of geometric modeling
Wang et al. Simplified representation for 3D point cloud data
Wang et al. Quick segmentation for complex curved surface with local geometric feature based on IGES and Open CASCADE
Wu et al. Variational mannequin approximation using spheres and capsules
CN107610231A (zh) 一种动态碰撞检测方法
Wang et al. A collision detection algorithm using AABB and octree space division
CN103236078A (zh) 一种基于面元的复杂场景三维重建方法
Su et al. View frustum culling algorithm based on optimized scene management structure

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200630

Address after: Office building of Shenzhen Institute of advanced technology A-207 518000 in Guangdong city of Shenzhen province Nanshan District City Road No. 1068 Chinese Academy of Shenzhen University Academy of Sciences

Patentee after: Shenzhen advanced science and technology Cci Capital Ltd.

Address before: 1068 No. 518055 Guangdong city in Shenzhen Province, Nanshan District City Xili University School Avenue

Patentee before: SHENZHEN INSTITUTES OF ADVANCED TECHNOLOGY CHINESE ACADEMY OF SCIENCES

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201029

Address after: 519000 the 12 floor of A District, No. 1, Port Road, Tangjia Bay, Zhuhai High-tech Zone, Guangdong.

Patentee after: ZHUHAI INSTITUTE OF ADVANCED TECHNOLOGY CHINESE ACADEMY OF SCIENCES Co.,Ltd.

Address before: Office building of Shenzhen Institute of advanced technology A-207 518000 in Guangdong city of Shenzhen province Nanshan District City Road No. 1068 Chinese Academy of Shenzhen University Academy of Sciences

Patentee before: Shenzhen advanced science and technology Cci Capital Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210420

Address after: Room 1102, building 5, Longyuan intelligent industrial park, No.2, hagongda Road, Tangjiawan Town, high tech Zone, Zhuhai City, Guangdong Province, 519000

Patentee after: Zhuhai Zhongke advanced technology industry Co.,Ltd.

Address before: 519000 the 12 floor of A District, No. 1, Port Road, Tangjia Bay, Zhuhai High-tech Zone, Guangdong.

Patentee before: ZHUHAI INSTITUTE OF ADVANCED TECHNOLOGY CHINESE ACADEMY OF SCIENCES Co.,Ltd.