CN101719140B - 一种图形检索方法 - Google Patents
一种图形检索方法 Download PDFInfo
- Publication number
- CN101719140B CN101719140B CN2009102140689A CN200910214068A CN101719140B CN 101719140 B CN101719140 B CN 101719140B CN 2009102140689 A CN2009102140689 A CN 2009102140689A CN 200910214068 A CN200910214068 A CN 200910214068A CN 101719140 B CN101719140 B CN 101719140B
- Authority
- CN
- China
- Prior art keywords
- grid model
- dimensional grid
- dimensional
- reference mark
- model
- 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 47
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 239000013598 vector Substances 0.000 claims description 28
- 238000000605 extraction Methods 0.000 claims description 22
- 239000011159 matrix material Substances 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 12
- 230000008676 import Effects 0.000 claims description 10
- 239000000284 extract Substances 0.000 claims description 8
- 230000000750 progressive effect Effects 0.000 claims description 6
- 210000003954 umbilical cord Anatomy 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 description 9
- 230000008878 coupling Effects 0.000 description 8
- 238000005859 coupling reaction Methods 0.000 description 8
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007418 data mining Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明公开一种图形检索方法。该方法包括:(1)进行多边形建模,预先建立三维网格模型库;(2)三维网格模型与二维图像或图形的匹配;(3)三维网格模型骨架的提取;(4)根据骨架进行三维模型检索;(5)三维网格模型特征点提取;(6)三维网格模型控制点的计算;(7)三维网格模型控制点频谱的计算;(8)计算频谱相似度,根据相似度检索出对应的图形。本发明的技术方案可以使得检索更为方便,并且支持多模态检索。
Description
技术领域
本发明涉及图形处理技术领域,具体涉及一种图形检索方法。
背景技术
网络与计算机图形学逐渐渗入日常生活中,人们已不再满足于只能在网络上看到二维图像。另一方面,随着三维建模技术的日益成熟和计算机软硬件技术的飞速发展,三维模型的数量在最近的十年中有了飞跃性的增长。相对于二维图像,人们可以从任意角度浏览自己感兴趣部分,因此更受人喜爱,用途更广泛。网络游戏和动漫技术、网络教育技术、基于Web的信息服务关键技术及产品和数据库与数据挖掘技术等热点领域的研究均不能缺少三维图形这一媒体。
充分利用已有的三维模型数据资源,可以大大减轻设计新模型的工作量,同时也可以促进三维数据的流通和在各领域的应用。然而,如何在海量的三维模型库中快速的搜索到自己感兴趣的模型,为三维模型库建立搜索引擎是一个困难的问题。现有技术的图形检索方法是根据几何内容对三维模型进行分类检索,用户不能方便地通过检索界面表达检索要求。
发明内容
本发明要解决的技术问题是提供一种图形检索方法,能够克服现有技术的不足,实现图形的多模态检索,使得检索更为方便,完善目前的多媒体搜索技术和填补目前网络三维图形搜索引擎的空白,推动下一代智能多模态搜索引擎的实现。
本发明提供的技术方案如下:
本发明提供一种图形检索方法,包括:
1)建立三维网格模型库;
2)当用户输入的是二维图形或图像时,与三维网格模型库中的三维网格模型的轮廓进行匹配,根据匹配参数将三维网格模型投影到二维空间,得到投影的二维图像或图形,然后计算投影得到的二维图像或图形与输入的图形或图像之间的相关度,根据相关度检索得到三维网格模型;
3)当用户输入的是三维网格模型时,对输入的三维网格模型进行骨架提取,根据提取的三维网格模型骨架,在三维网格模型库中初步检索得到三维网格模型;
4)将检索得到的三维网格模型和用户输入的三维网格模型进行特征点提取,代替原始三维网格模型,再进行三角剖分,对剖分后的分割线进行分段拟合,得到原始三维网格模型的控制点,然后根据拓扑结构对控制点进行频域变换;
5)计算得到的用户输入三维网格模型的控制点频域坐标值与三维网格模型库中的三维网格模型的控制点频域坐标值之间的相似度,根据相似度检索出对应的图形;
其中,步骤4)中根据三维网格模型的拓扑结构对得到的控制点进行频域变换,包括:
(1)以控制点到网格中心的矢量的模对粗糙网格的控制点进行排序;
(2)从网格拓扑关系获得Kirchhoff矩阵
K=D-A (6)
D是对角矩阵,其对角线上的元素Dii与顶点vi的价相对应,A是网格的邻接矩阵;
对Kirchhoff矩阵进行特征值分解得到的n个特征向量wi进行升序排列,组成的n*n映射矩阵W;
(3)从先前排好序的n个控制点的空间坐标构造3个向量:
X=(x1,x2,…,xn),Y=(y1,y2,…,yn),Z=(z1,z2,…,zn) (8)
将这3个向量投影到特征向量基W上得到频域向量:
每个顶点对应的频谱的幅值Si计算公式为:
优选的,步骤1)中的三维网格模型库是对各种三维数据格式进行重新组织和多边形建模后得到的。
优选的,骨架提取过程为:首先为输入的网格模型建立渐进网格表示,然后对渐进网格不断的进行边塌缩变换,在塌缩的过程中如果一条边没有相邻三角形,则该条边标记为骨架边,并且一直保留到塌缩结束,最终获得的边构成模型的骨架。
优选的,步骤4)中对于检索得到的三维网格模型和用户输入的三维网格模型,根据其空间形状进行特征点提取,用脐带点作为特征点代替原始三维网格模型。
优选的,步骤4)中根据特征点对三维网格模型进行三角剖分,对剖分后的分割线进行分段拟合,得到分割点,将这些分割点作为原始三维网格模型的控制点。
本发明具有以下有益效果:
(1)可以根据单个二维图像检索三维图形,用户可以输入bmp、jpeg、tiff等常见格式的图片,通过本发明方法可以检索出对应的图形。
(2)支持输入的三维模型的表示方法更加广泛,对于三角形网格数据、点云数据、体数据或多边形网格模型均适用。
(3)采用了分级搜索,更快更准确。首先根据三维图形的骨架进行粗糙搜索,对搜索结果中的三维图形采用其它特征提取技术进行准确搜索,保证了实时性和准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1是用户输入二维图形或图像的处理示意图;
图2是三维与二维的匹配流程图;
图3是对于用户输入三维网格数据时的处理流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明方法主要研究图形的多模态检索,其原理是根据三维图形的特征,通过二维图像、二维图形、三维模型等信息,计算用户输入的图形或图像与三维图形之间的相似度,从而实现三维图形多模态检索。这里多模态是指支持用户以常见格式的二维图像、二维图形和三维模型进行查询。
本发明方面从三维数据中提取出很小的数据量,并将其作为对应图形或图像的主要特征,可以根据该主要特征进行检索。该特征基本不受噪声、相似变换、不同分辨率采样等因素的影响。
本发明方法主要包括以下8个环节:
(1)首先对各种数据格式的三维模型数据进行重新组织,进行多边形建模,得到三维网格模型,并建立三维网格模型库。
(2)当用户输入的是二维图形或图像时,与三维网格模型库中的三维网格模型的轮廓进行匹配,然后根据匹配参数将三维网格模型投影到二维空间,得到投影的二维图像或图形,然后计算投影得到的二维图像或图形与输入的图形或图像之间的相关度,根据相关度检索出三维网格模型。
(3)当用户输入的是三维网格模型时,对输入的三维网格模型进行骨架提取。
(4)根据提取的三维网格模型骨架,在三维网格模型库中进行快速的初步检索,得到初步检索出的三维网格模型。
(5)根据图形学理论将初步检索出的三维网格模型和用户输入的三维网格模型,根据其空间形状进行特征点提取,代替原始三维网格模型,大大减少数据量。
(6)根据提取的特征点对三维网格模型进行三角剖分,对剖分后的分割线进行分段拟合,得到分割点,将这些分割点作为原始三维网格模型的控制点。
(7)根据三维网格模型的拓扑结构对控制点进行频域变换,得到控制点的频域坐标值。
(8)计算得到的用户输入三维网格模型的控制点频域坐标值与三维网格模型库中的三维网格模型的控制点频域坐标值之间的相似度,根据相似度检索出对应的图形。
本发明的技术特点主要体现如下:
(1)***可以根据单个二维图像检索三维图形。用户可以输入bmp、jpeg、tiff等常见格式的图片,***将从图片中分割对象并提取轮廓信息,与三维网格模型库中的三维图形进行匹配,向用户返回匹配的三维网格模型。
(2)***支持输入的三维模型的表示方法更加广泛。三角形网格数据、点云数据、体数据或多边形网格模型均适用。
(3)***采用分级搜索。首先根据三维图形的骨架进行粗糙搜索,对搜索结果中的三维图形再采用特征提取技术进行准确搜索,保证了实时性和准确性。
下面对本发明做进一步详细说明。
本发明支持多模态的图形检索方法主要步骤包括:
(1)进行多边形建模,预先建立三维网格模型库;
(2)三维网格模型与二维图像或图形的匹配;
(3)三维网格模型骨架的提取;
(4)根据骨架进行三维模型检索;
(5)三维网格模型特征点提取;
(6)三维网格模型控制点的计算;
(7)三维网格模型控制点频谱的计算;
(8)计算频谱相似度,根据相似度检索出对应的图形。
以下分别对上述步骤进行详细介绍。
(1)进行多边形建模,预先建立三维网格模型库。
本发明根据三维网格模型数据,进行多边形建模,构造多边形三维网格模型。多边形建模是利用许多的多边形模拟曲面进行,多边形越多,则模型越逼近真实曲面。多边形建模是最广泛又易于实现的一种建模技术,并可以获得高精度的模型,通常采用三角形网格形式。
为了下述描述的方便和统一,利用数学符号给出三维网格模型的定义。网格模型M={V,C},由顶点集合V和连接关系集合C组成,其中集合V包含N个顶点vi且每个顶点的坐标值由(xi,yi,zi)确定,即
V={vi},i=0,1,…,N-1,vi=(xi,yi,zi) (1)
而连接关系集合C表示成
C={{ik,jk}}k=0,…m-1,0≤ik≤N-1,0≤jk≤N-1 (2)
这里{ik,jk}表示由第ik个顶点和第jk个顶点确定的第k条边。
(2)三维网格模型与二维图像或图形的匹配
图1是用户输入二维图形或图像的处理示意图。
如图1所示,对于单个图像/图形,与模型库图形进行轮廓匹配,再进行2D投影,然后进行图形匹配相关度计算。
在很多情况下,用户的检索输入为二维图像或图形,与三维网格模型的匹配具体过程如图2所示。首先,本发明采用图像的轮廓(或图形轮廓)与三维网格模型库中的三维网格模型轮廓先进行初始匹配;轮廓匹配后,将得到的三维网格模型根据匹配参数投影到2D空间,得到投影图像,然后用图像相关匹配方法进行进一步精确匹配。
三维网格模型库中的所有网格模型已完成X、Y、Z三个正方向的轮廓提取和投影计算。
对于三维网格模型,本发明方法是遍历网格中的每一个边来提取的轮廓,具体方法如下:
1.如果当前边仅与一个三角形相连接,那么它属于轮廓;
2.如果当前边与两个三角形F1和F2,则定义其法向量分别为和,当前镜头位置与当前边的一个顶点之间的向量为。如果,即和相对于镜头的轴处于不同的方向,说明F1和F2一个正对着镜头一个背对着镜头,因此当前边为轮廓边,否则,当前边不是显著的轮廓边。
对于图形,不需要轮廓提取过程。
对于图像,可以采用Sobel、Prewitt等算子提取轮廓。
本发明轮廓匹配过程可采用基于对应形状匹配方法,如Hausdoff距离;或者多边形分解匹配方法,如进行线段化处理,并且以轮廓重心为中心进行三角剖分,确定边界顶点,连接起来便成一凸多边形,用这一凸多边形近似表示原图形,将多边形重心与顶点连线,组成一系列三角形从而进行图形轮廓匹配。轮廓匹配过程简单高效,可快速实现初始检索。
本发明图像相关匹配方法是计算三维网格模型投影得到的二维图像或图形与输入的图形或图像之间的相关度,如利用归一化相关测度计算图像区域中每一对像素的相似性。对于待匹配的两幅图像I1(x,y)和I2(x,y),待检测图位置(i,j)上归一化相关测度定义为:
当相关度大于一个设定阈值,即检索到匹配的三维网格模型。
(3)三维网格模型骨架的提取
骨架是一种性质优良的图形几何特征,又称中轴(Medial Axis),是一种有效的图形描述手段。顾名思义,骨架是一种线型的几何体,居于图形的对称中心,有着与原图形相同的拓扑结构,并保留着原图形的形状信息。
图3是对于用户输入三维网格数据时的处理流程图。
图3的过程包括以下的过程,即(3)三维网格模型骨架的提取;(4)根据骨架进行三维模型检索;(5)三维网格模型特征点提取;(6)三维网格模型控制点的计算;(7)三维网格模型控制点频谱的计算;(8)计算频谱相似度,根据相似度检索出对应的图形。具体内容参加下面的描述。
本发明设计了一种快速的基于多分辨率网格的骨架提取算法:首先为三维网格模型建立渐进网格表示,然后对渐进网格不断的进行边塌缩变换,在塌缩的过程中如果一条边没有相邻三角形,那么这条边标记为骨架边,并且一直保留到塌缩结束,最终获得的边就构成了网格模型的骨架。三维网格模型库中的所有网格模型已完成骨架提取计算。用户输入网格模型时需要提取骨架。
(4)根据骨架进行三维模型检索
提取骨架是将3D图形转换为3D线段的过程,3D线段的数据量相对于原始3D图形大大减少,因此可以加快检索速度。
用户输入的三维网格模型与库中的网格模型进行骨架比较,比较骨架采用主成分PCA分析法,先将骨架定位,分段计算骨架的距离,然后直接采用欧氏距离进行排序,根据距离信息比较结果,欧氏距离差最小的即为初步检索的三维网格模型。
(5)三维网格模型特征点提取
根据图形学理论将用户输入的三维网格模型,根据其空间形状进行特征点提取,代替原始三维网格模型,大大减少数据量。三维网格模型库中的所有网格模型已完成特征点提取、控制点频谱计算。
本发明使用脐带点作为任意三角形三维网格模型的特征点。以脐带点作为三维网格模型表面的特征点,对于噪声、剪切、旋转、平移、缩放、不同分辨率采样等由于不同采集设备造成的影响具有很强的鲁棒性。由于跨过网格模型的边的曲率较大,所以曲率张量可以定义为网格模型的边的每一个点。在任意网格区域B内,定点v的曲率张量可以用下式估计:
其中,|B|是v的邻域的面积,β(e)是边e的两个邻接三角形的法向量之间的夹角,|e∩B|是区域B中的边e的长度,是沿着边e的单位法向量。v的邻接区域B是由以v为球心,以r为半径的球体与网格模型的相交圆定义的。半径r是指定曲率估计的尺度参数。
得到每一个顶点的曲率张量以后,在每个三角形网格上进行线性插值以获得连续的曲率张量场。顶点的法向量方向与曲率张量的幅值最小的特征值相对应,其它两个特征值分别对应顶点v的最小曲率和最大曲率,当这两个特征值相等时,顶点v被称为脐带点,即本发明中使用的三维网格模型的特征点。
显然曲率估计的关键是尺度参数。本发明采用不同的尺度进行曲率张量估计,以便于平滑张量场并且估计不同尺度下的脐带点对于噪声以及仿射变换等因素的鲁棒性。为了避免丢失网格模型表面的局部信息以及减小计算复杂度,尺度参数不能取得过大。本发明在寻找鲁棒性最高的尺度参数时,选用自适应遗传算法作为优化搜索技术。在鲁棒性最高的几个尺度参数中,我们选用使区域B内的平均曲率最大的点,即曲率张量的迹最大的点。
(6)三维网格模型控制点的计算
根据提取的特征点对三维网格模型进行三角剖分,对剖分后的分割线进行分段拟合,得到分割点,将这些分割点作为原始三维网格模型的控制点。
获得特征点以后,下一步任务是减小数据量,这对于特征注册以及三维网格模型的检索都非常重要。
获得特征点以后,第一步是将网格进行三角分割。在2D空间中,Delaunay三角分割能够产生形状较均匀的三角形,并且具有唯一性。在3D表面中,Delaunay三角分割不是使用欧氏距离而是使用测地距离(测地距离是指3D空间中网格表面的两个顶点沿着表面的最短距离)。本发明采用波阵面方法对整个Voronoi图及其二重Delaunay三角分割进行估计。采用该方法的好处是能够使三角分割不受采样率的大小的影响。波阵面是通过计算以种子点为球心,半径不断增长的球与三维网格表面的相交圆获得的。采用基于不同的球半径的波阵面方法比采用基于网格拓扑关系的波阵面方法得到的结果要好。
(7)三维网格模型控制点频谱的计算
根据三维网格模型的拓扑结构对控制点进行频域变换,得到控制点的频域坐标值。
由于3D网格是图形而非图像,所以每个顶点坐标没有其固有的影像函数,因此对3D网格进行变换域处理首先要构造一个影像函数,以便于将经典的变换域处理算法应用于3D网格处理。
首先,对控制点进行排序。本发明以控制点到网格中心的矢量的模对粗糙网格的控制点进行排序。
定义网格中心的坐标为
每个控制点到网格中心的矢量的模定义如下
将控制点根据模大小进行升序排列以后,开始进行频域变换。首先从网格拓扑关系中获得组合拉普拉斯算子或者叫做Kirchhoff矩阵。该矩阵定义如下:
K=D-A (6)
其中,D是对角矩阵,其对角线上的元素Dii与顶点vi的价相对应(价即是从顶点放射出的边的个数),A是网格的邻接矩阵,其元素定义如下:
对于有n个控制点的网格模型,矩阵A、D和K的尺寸均为n*n。对Kirchhoff矩阵进行特征值分解得到n个特征值λi和n个特征向量wi。将这n个特征向量进行升序排列,可以得到其对应的特征向量,这些排序后的特征向量是频率不断增大的函数基。该函数基仅取决于网格模型的拓扑结构,而与网格模型的几何特性无关。将这n个排序后的特征向量组成的n*n映射矩阵记做W。
从先前排好序的n个控制点的空间坐标构造3个向量:
X=(x1,x2,…,xn),Y=(y1,y2,…,yn),Z=(z1,z2,…,zn) (8)
将这3个向量投影到特征向量基W上,即可得到空域坐标的频域分解向量:
空间坐标也可由频率坐标还原得到:
每个顶点对应的频谱的幅值Si可以由下式计算:
(8)计算频谱相似度,根据相似度检索出对应的图形。
计算用户输入的三维网格模型的控制点频域坐标值与库中的三维网格模型的控制点频域坐标值之间的相关度,根据相关度检索出对应的图形。
比较两个网格模型的相似度时,则通过比较频域系数的波形相似度来判断。比较相似度时的准则可采用归一化相关系数:
其中为由用户输入的三维网格模型计算得到的频谱值,为库中的三维网格频谱值。NC越大,则频谱系数越相关。
给定一个阈值T,如果NC>T,则认为用户输入的三维网格与库中的三维网格相匹配,检索流程结束。
本发明具有以下有益效果:
(1)可以根据单个二维图像检索三维图形,用户可以输入bmp、jpeg、tiff等常见格式的图片,通过本发明方法可以检索出对应的图形。
(2)支持输入的三维模型的表示方法更加广泛,对于三角形网格数据、点云数据、体数据或多边形网格模型均适用。
(3)采用了分级搜索,更快更准确。首先根据三维图形的骨架进行粗糙搜索,对搜索结果中的三维图形采用其它特征提取技术进行准确搜索,保证了实时性和准确性。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器(ROM,Read Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁盘或光盘等。
以上对本发明实施例所提供的一种图形检索方法,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (5)
1.一种图形检索方法,其特征在于,包括:
1)建立三维网格模型库;
2)当用户输入的是二维图形或图像时,与三维网格模型库中的三维网格模型的轮廓进行匹配,根据匹配参数将三维网格模型投影到二维空间,得到投影的二维图像或图形,然后计算投影得到的二维图像或图形与输入的图形或图像之间的相关度,根据相关度检索得到三维网格模型;
3)当用户输入的是三维网格模型时,对输入的三维网格模型进行骨架提取,根据提取的三维网格模型骨架,在三维网格模型库中初步检索得到三维网格模型;
4)将检索得到的三维网格模型和用户输入的三维网格模型进行特征点提取,代替原始三维网格模型,再进行三角剖分,对剖分后的分割线进行分段拟合,得到原始三维网格模型的控制点,然后根据拓扑结构对控制点进行频域变换;
5)计算得到的用户输入三维网格模型的控制点频域坐标值与三维网格模型库中的三维网格模型的控制点频域坐标值之间的相似度,根据相似度检索出对应的图形;
其中,步骤4)中根据三维网格模型的拓扑结构对得到的控制点进行频域变换,包括:
(1)以控制点到网格中心的矢量的模对粗糙网格的控制点进行排序;
(2)从网格拓扑关系获得Kirchhoff矩阵
K=D-A (6)
D是对角矩阵,其对角线上的元素Dii与顶点vi的价相对应,A是网格的邻接矩阵;
对Kirchhoff矩阵进行特征值分解得到的n个特征向量wi进行升序排列,组成的n*n映射矩阵W;
(3)从先前排好序的n个控制点的空间坐标构造3个向量:
X=(x1,x2,…,xn),Y=(y1,y2,…,yn),Z=(z1,z2,…,zn) (8)
将这3个向量投影到特征向量基W上得到频域向量:
每个顶点对应的频谱的幅值Si计算公式为:
2.根据权利要求1所述的图形检索方法,其特征在于:
步骤1)中的三维网格模型库是对各种三维数据格式进行重新组织和多边形建模后得到的。
3.根据权利要求1或2所述的图形检索方法,其特征在于:
骨架提取过程为:首先为输入的网格模型建立渐进网格表示,然后对渐进网格不断的进行边塌缩变换,在塌缩的过程中如果一条边没有相邻三角形,则该条边标记为骨架边,并且一直保留到塌缩结束,最终获得的边构成模型的骨架。
4.根据权利要求1或2所述的图形检索方法,其特征在于:
步骤4)中对于检索得到的三维网格模型和用户输入的三维网格模型,根据其空间形状进行特征点提取,用脐带点作为特征点代替原始三维网格模型。
5.根据权利要求1或2所述的图形检索方法,其特征在于:
步骤4)中根据特征点对三维网格模型进行三角剖分,对剖分后的分割线进行分段拟合,得到分割点,将这些分割点作为原始三维网格模型的控制点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102140689A CN101719140B (zh) | 2009-12-23 | 2009-12-23 | 一种图形检索方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102140689A CN101719140B (zh) | 2009-12-23 | 2009-12-23 | 一种图形检索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101719140A CN101719140A (zh) | 2010-06-02 |
CN101719140B true CN101719140B (zh) | 2012-04-18 |
Family
ID=42433714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102140689A Expired - Fee Related CN101719140B (zh) | 2009-12-23 | 2009-12-23 | 一种图形检索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101719140B (zh) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102375831B (zh) * | 2010-08-13 | 2014-09-10 | 富士通株式会社 | 三维模型检索装置及其方法和模型库生成装置及其方法 |
US8437506B2 (en) * | 2010-09-07 | 2013-05-07 | Microsoft Corporation | System for fast, probabilistic skeletal tracking |
CN101957992B (zh) * | 2010-09-30 | 2012-05-02 | 清华大学 | 一种二维单射曲面数据的特征提取与匹配方法 |
CN102004922B (zh) * | 2010-12-01 | 2012-12-26 | 南京大学 | 基于骨架特征的高分辨率遥感影像飞机提取方法 |
CN102063475B (zh) * | 2010-12-22 | 2012-10-10 | 张丛喆 | 三维模型的网页用户终端呈现方式 |
CN103309938A (zh) * | 2013-04-24 | 2013-09-18 | 中国科学院遥感与数字地球研究所 | 一种多边形最大内圆的高性能实现方法 |
CN103412947A (zh) * | 2013-08-26 | 2013-11-27 | 浙江大学 | 针对空间大数据的多边形搜索方法 |
CN103544734A (zh) * | 2013-10-11 | 2014-01-29 | 深圳先进技术研究院 | 以街景为基础的三维地图建模的方法 |
TWI536186B (zh) * | 2013-12-12 | 2016-06-01 | 三緯國際立體列印科技股份有限公司 | 三維圖檔搜尋方法與三維圖檔搜尋系統 |
CN104503275B (zh) * | 2014-11-21 | 2017-03-08 | 深圳市超节点网络科技有限公司 | 基于手势的非接触式控制方法及其设备 |
CN105354866B (zh) * | 2015-10-21 | 2018-11-06 | 郑州航空工业管理学院 | 一种多边形轮廓相似度检测方法 |
CN106909552A (zh) * | 2015-12-22 | 2017-06-30 | 成都理想境界科技有限公司 | 图像检索服务器、***、相关检索及排错方法 |
CN106909539A (zh) * | 2015-12-22 | 2017-06-30 | 成都理想境界科技有限公司 | 图像检索***、服务器、数据库及相关的方法 |
CN106445981B (zh) * | 2016-02-15 | 2020-03-27 | 哈尔滨理工大学 | 一种基于小波变换的stl网格模型切片数据自适应压缩方法 |
CN107316343B (zh) * | 2016-04-26 | 2020-04-07 | 腾讯科技(深圳)有限公司 | 一种基于数据驱动的模型处理方法及设备 |
CN108596186B (zh) * | 2018-03-19 | 2021-06-22 | 西北大学 | 一种三维模型检索方法 |
CN108898128B (zh) * | 2018-07-11 | 2020-09-04 | 宁波艾腾湃智能科技有限公司 | 一种通过照片匹配数字化三维模型的防伪方法及设备 |
CN110019914B (zh) * | 2018-07-18 | 2023-06-30 | 王斌 | 一种支持三维场景交互的三维模型数据库检索方法 |
CN109034418B (zh) * | 2018-07-26 | 2021-05-28 | 国家电网公司 | 作业现场信息传输方法及*** |
CN109359534B (zh) * | 2018-09-12 | 2021-02-09 | 鲁东大学 | 一种三维物体几何特征提取方法及*** |
CN109299301B (zh) * | 2018-09-17 | 2021-09-14 | 北京工业大学 | 一种基于形状分布和曲度的三维模型检索方法 |
CN109816771B (zh) * | 2018-11-30 | 2022-11-22 | 西北大学 | 一种结合特征点拓扑和几何约束的文物碎片自动重组方法 |
CN111189440B (zh) * | 2019-12-31 | 2021-09-07 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于空间信息模型与实时图像比对的定位导航方法 |
CN111643906B (zh) * | 2020-05-29 | 2021-08-31 | 腾讯科技(深圳)有限公司 | 一种信息处理方法、装置及计算机可读存储介质 |
CN112270762A (zh) * | 2020-11-18 | 2021-01-26 | 天津大学 | 一种基于多模态融合的三维模型检索方法 |
CN112732956A (zh) * | 2020-12-24 | 2021-04-30 | 江苏智水智能科技有限责任公司 | 基于感知多模态大数据的高效查询方法 |
-
2009
- 2009-12-23 CN CN2009102140689A patent/CN101719140B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101719140A (zh) | 2010-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101719140B (zh) | 一种图形检索方法 | |
Wang et al. | Lidar point clouds to 3-D urban models $: $ A review | |
Mehra et al. | Abstraction of man-made shapes | |
De Luca et al. | Reverse engineering of architectural buildings based on a hybrid modeling approach | |
Weibel et al. | Generalising spatial data and dealing with multiple representations | |
CN108171780A (zh) | 一种基于激光雷达构建室内真实三维地图的方法 | |
CN106780458B (zh) | 一种点云骨架提取方法及装置 | |
Floriani et al. | VARIANT: a system for terrain modeling at variable resolution | |
US20240153123A1 (en) | Isogeometric Analysis Method Based on a Geometric Reconstruction Model | |
van Kreveld | Geographic information systems | |
CN104361632A (zh) | 一种基于Hermite径向基函数的三角网格补洞方法 | |
CN109344533B (zh) | 地下工井电缆网络模型的建立方法 | |
CN105631932A (zh) | 一种轮廓线指导的三维模型重构造方法 | |
CN115661374B (zh) | 一种基于空间划分和模型体素化的快速检索方法 | |
Raina et al. | Sharpness fields in point clouds using deep learning | |
Shah et al. | Simulated annealing-based fitting of CAD models to point clouds of mechanical parts’ assemblies | |
CN104715507A (zh) | 一种基于曲面片的三维地理实体自动构建方法 | |
Kamra et al. | Lightweight reconstruction of urban buildings: Data structures, algorithms, and future directions | |
Huang et al. | Meshode: A robust and scalable framework for mesh deformation | |
Kidner et al. | Multiscale terrain and topographic modelling with the implicit TIN | |
Zhou et al. | Distributed optimization based on graph filter for ensuring invariant simplification of high-volume point cloud | |
Huang et al. | Automatic CAD model reconstruction from multiple point clouds for reverse engineering | |
Pastor et al. | 3D wavelet-based multiresolution object representation | |
Gonzaga de Oliveira | A review on delaunay refinement techniques | |
Park et al. | Automatic 3-D model synthesis from measured range data |
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 |
Granted publication date: 20120418 Termination date: 20141223 |
|
EXPY | Termination of patent right or utility model |