CN108305289B - 基于最小二乘法的三维模型对称性特征检测方法及*** - Google Patents

基于最小二乘法的三维模型对称性特征检测方法及*** Download PDF

Info

Publication number
CN108305289B
CN108305289B CN201810073992.9A CN201810073992A CN108305289B CN 108305289 B CN108305289 B CN 108305289B CN 201810073992 A CN201810073992 A CN 201810073992A CN 108305289 B CN108305289 B CN 108305289B
Authority
CN
China
Prior art keywords
plane
dimensional model
symmetry
symmetrical
points
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
Application number
CN201810073992.9A
Other languages
English (en)
Other versions
CN108305289A (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.)
Shandong Normal University
Original Assignee
Shandong Normal University
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 Shandong Normal University filed Critical Shandong Normal University
Priority to CN201810073992.9A priority Critical patent/CN108305289B/zh
Publication of CN108305289A publication Critical patent/CN108305289A/zh
Application granted granted Critical
Publication of CN108305289B publication Critical patent/CN108305289B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/68Analysis of geometric attributes of symmetry
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/757Matching configurations of points or features

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Geometry (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了基于最小二乘法的三维模型对称性特征检测方法及***,步骤:建立三维模型的最小包围盒,建立新的坐标系;根据空间分辨率划分最小包围盒的平面,将划分的六个面中相对的面设置为一组,在最小包围盒的每个平面均取若干个等间距标记点,将每组标记点投影到三维模型表面得到投影点,计算标记点与投影点的匹配度,根据匹配度筛选对称数据组;将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面;检验校正初始对称面,不断细化初始对称面直到得到准确的对称面。本发明由计算最小包围盒、筛选对称数据组、平面反射对称性检测、检验校正对称面四部分组成,检测三维模型的对称面。

Description

基于最小二乘法的三维模型对称性特征检测方法及***
技术领域
本发明涉及基于最小二乘法的三维模型对称性特征检测方法及***。
背景技术
计算机群体仿真对于现代社会研究具有重要的意义,随着关注度的提高,群体仿真应用的领域也越来越广,相关技术也逐渐发展成熟。因此,人们对群体仿真中应用到的模型要求也越来越高,既想保持模型的精度,又不能占用太多的内存空间,三维网格简化算法也由最初的全局简化逐渐过渡为局部特征保留的简化方法。考虑到这些三维模型本身大部分都具有对称性的特征,比如动物、家具、人物、汽车等,如何利用物体的对称性减少简化步骤中的计算量,对提高三维模型简化算法的性能具有十分重要的意义。但在实际应用步骤中,并不是所有的三维模型都是完美对称的。对于非完美对称的三维模型,如何高效、准确的检测到模型的对称面,具有十分重要的研究意义。
发明内容
为了解决现有技术的不足,本发明提供了基于最小二乘法的三维模型对称性特征检测方法及***,本发明由计算最小包围盒、筛选对称数据组、平面反射对称性检测、检验校正对称面四部分组成,检测三维模型的对称面。
为了实现上述目的,本发明采用如下方案:
基于最小二乘法的三维模型对称性特征检测方法,包括以下步骤:
对三维模型的最小包围盒的六个平面的每个平面,依据空间分辨率进行单元格划分,在划分的最小包围盒的每个平面上选取若干个等间距标记点,将划分的六个平面的每两个相对的面设置为一组,将每组标记点投影到三维模型表面得到投影点,计算标记点与投影点的匹配度,根据匹配度筛选对称数据组;
将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面。
优选的,在对三维模型的最小包围盒的六个平面的每个平面,依据空间分辨率进行划分步骤之前,还包括:建立三维模型的最小包围盒,建立新的坐标系。
优选的,在找到初始对称面步骤之后,还包括:
检验校正初始对称面,不断细化初始对称面直到得到准确的对称面。
优选的,所述建立新的坐标系,是以最小包围盒的中心点为坐标轴原点,建立新的坐标系。
优选的,所述根据空间分辨率划分最小包围盒的平面,包括步骤如下:
计算三维模型顶点个数与最小包围盒的体积的比值,将比值与设定阈值进行比较,如果比值大于等于设定阈值,则将最小包围盒的六个平面按照比值等间距划分为若干个单元格;如果比值小于设定阈值,则将最小包围盒的六个平面按照像素划分为若干个单元格。
优选的,计算标记点与投影点的匹配度,包括步骤如下:
两个对称的标记点都有两个投影点,如果所述两个投影点的横坐标、纵坐标与竖坐标三项中有两项的数值一致,只有一项不一致,那么判断不一致的坐标数值偏差是否在设定的偏差测度范围内,如果是,就判断两个标记点匹配,否则两个标记点不匹配。
如果两个对称的标记点,只有一个投影点或者没有投影点,那么直接判断为两个标记点不匹配。
如果两个标记点匹配,那么计算匹配的标记点对数与该组标记点对数的比值,得到的比值就是标记点与投影点的匹配度。
优选的,根据匹配度筛选对称数据组的步骤如下:
匹配度值最高的标记点组对应的投影点集合,即为对称数据组。
优选的,将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面的步骤如下:
将每份对称数据组内的离散数据点经过拟合后得出对应的直线,将所有的直线连接成面得到一条初始对称面。
优选的,所述的检验校正初始对称面的步骤如下:
将初始对称面与预设对称平面进行比对,初始对称面与预设对称平面有部分平行或重合,则该平行或重合部分定义为完全对称部分,其余部分定义为非完全对称部分;
对于非完全对称部分,则将等间距缩小一半,重复寻找初始对称面的步骤。
基于最小二乘法的三维模型对称性特征检测***,包括:存储器、处理器以及存储在存储器上并在处理器上运行的计算机程序,所述计算机程序被处理器运行时,完成上述任一方法所述的步骤。
一种计算机可读存储介质,其上运行有计算机指令,所述计算机指令被处理器运行时,完成上述任一方法所述的步骤。
本发明的有益效果:
(1)本发明的基于最小二乘法的三维模型对称性特征检测方法,在预处理阶段,根据模型分布重新定义坐标轴方向建立新的坐标系,使得对应点的坐标只有某项坐标不同,为后续的计算提供了便利、节省了时间;
(2)本发明的基于最小二乘法的三维模型对称性特征检测方法,根据模型的空间分辨率等间距缩减计算范围,从原始的遍历整个模型的点进行计算缩减为计算小范围数据,不仅计算量上大大减少,也没丢失模型的对称特性;
(3)本发明的基于最小二乘法的三维模型对称性特征检测方法,利用微积分的思想,先将整体细分为局部,再由局部细节汇总为整体,既降低了算法的复杂度避免重复迭代,又将对称的细节特征得到较好的保留,最大程度上依托模型的对称性找到了模型的对称曲面,在实际应用步骤中,为后续的三维模型的简化工作提供了便利。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1是本发明的模块示意图;
图2(a)和图2(b)是本发明中包围盒概念示意图;
图3是本发明中网格划分示意图;
图4(a)和图4(b)是本发明中某模型包围盒及坐标系变化对比示意图;
图5是本发明中某模型的网格划分示意图;
图6是本发明中某个面向模型投影的效果示意图;
图7是本发明中某模型的子集平面示意图;
图8是本发明中某模型的某子集平面拟合直线的示意图;
图9(a)、图9(b)和图9(c)是本发明中某模型的对称检测示意图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面结合附图与实施例对本发明作进一步说明。
下面结合附图与实施例对本发明做进一步说明:
本发明可以应用到3D打印领域,还可以应用到三维动画设计领域。
如图1所示,基于最小二乘法的三维模型对称性特征检测方法,包括以下步骤:
步骤(1):建立三维模型的最小包围盒,建立新的坐标系;
步骤(2):根据空间分辨率划分最小包围盒的平面,在最小包围盒的每个平面均取若干个等间距标记点,将标记点投影到三维模型表面得到对称特征点,计算标记点与对称特征点的匹配度,根据匹配度筛选对称数据组;
步骤(3):将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面;
步骤(4):检验校正初始对称面,不断细化初始对称面直到得到准确的对称面。
优选的,所述建立新的坐标系,是以最小包围盒的中心点为坐标轴原点,建立新的坐标系。
优选的,所述步骤(1)中,所述的建立三维模型的最小包围盒的步骤如下:
步骤(1-1):假设三维模型的顶点集合P={p1,p2,...,pn},其中,pi=(xi,yi,zi)∈R,i=1,2,...,n组成,令
Figure BDA0001558971680000041
为三维模型的质心;
步骤(1-2):根据最小包围盒内的点的性质,将最小包围盒的点的最大值设为正无穷大,将最小包围盒的点的最小值设为负无穷大;
步骤(1-3):遍历全部点,并扩展最小包围盒的边界框直到它包含三维模型所有点为止。
如图2(a)和图2(b)所示,优选的,所述步骤(1-2)中,最小包围盒内的点的性质为:
(1-2-1):最小包围盒内的点满足下列不等式:
Xmin≤X≤Xmax;(1-1)
Ymin≤Y≤Ymax;(1-2)
Zmin≤Z≤Zmax;(1-3)
其中,
X表示最小包围盒内的三维模型的点的X坐标,
Y表示最小包围盒内的三维模型的点的Y坐标,
Z表示最小包围盒内的三维模型的点的X坐标,
Xmin表示最小包围盒内的三维模型的点的X坐标的最小值,
Ymin表示最小包围盒内的三维模型的点的Y坐标的最小值,
Zmin表示最小包围盒内的三维模型的点的Z坐标的最小值,
Xmax表示最小包围盒内的三维模型的点的X坐标的最大值,
Ymax表示最小包围盒内的三维模型的点的Y坐标的最大值,
Zmax表示最小包围盒内的三维模型的点的Z坐标的最大值。
(1-2-2):最小包围盒内的两个点:
Pmin=[Xmin Ymin Zmin];(2-1)
Pmax=[Xmax Ymax Zmax];(2-2)
其中,
Pmin表示最小包围盒内的三维模型的点的各坐标值取最小值的点,也是最小包围盒的八个顶点中各坐标值最小的一个顶点,
Pmax表示最小包围盒内的三维模型的点的各坐标值取最大值的点,也是最小包围盒的八个顶点中各坐标值最大的一个顶点,由Pmin和Pmax两个顶点定义最小包围盒。
(1-2-3):最小包围盒的中心点C为:
C=(Pmin+Pmax)/2;(3)
(1-2-4):尺寸向量S是从Pmin指向Pmax的向量,尺寸向量S包含最小包围盒矩形边界的长、宽和高;
S=Pmax-Pmin;(4)
(1-2-5):计算质心点pm与中心点C的距离作为偏差测度E。
优选的,所述步骤(1)中,建立新的坐标系的方法如下:
步骤(1-4):将坐标轴的原点位置移到计算出的中心点C;
步骤(1-5):采用霍特林变换调整三维模型的坐标,使得三维模型的几何形状不变,形成新的坐标系,同时更新Pmin和Pmax的坐标值:
Pmin=(xmin,ymin,zmin);
Pmax=(xmax,ymax,zmax);
其中,xmin表示新坐标系下最小包围盒内的三维模型的点的X坐标的最小值,ymin表示新坐标系下最小包围盒内的三维模型的点的Y坐标的最小值;zmin表示新坐标系下最小包围盒内的三维模型的点的Z坐标的最小值,xmax表示新坐标系下最小包围盒内的三维模型的点的X坐标的最大值,ymax表示新坐标系下最小包围盒内的三维模型的点的Y坐标的最大值,zmax表示新坐标系下最小包围盒内的三维模型的点的Z坐标的最大值;
最小包围盒各个平面的表达式为:
F1:x=a;F2:y=b;F3:z=c;
F4:x=-a;F5:y=-b;F6:z=-c,
其中,F1和F4、F2和F5、F3和F6分别是最小包围盒中相对的两个平面,a、b和c均为常数;本示例中某模型包围盒及坐标系变化对比如图4(a)和图4(b)所示。
优选的,所述步骤(2)中,根据空间分辨率划分最小包围盒的步骤为:
计算三维模型顶点个数n与最小包围盒的体积M之比值,依据比值对最小包围盒进行划分:
Figure BDA0001558971680000061
时,将最小包围盒的六个面依据比值等间距划分为若干单元格;
Figure BDA0001558971680000062
时,将最小包围盒的六个面根据像素划分为若干单元格,如图3所示。
将相对的面分为一组,得到三组分别有2*X'*Y'、2*Y'*Z'、2*X'*Z'个标记点的标记点集合A、B、C,其中A=A1∪A2,B=B1∪B2,C=C1∪C2
其中,
集合A1是划分平面F1后取得标记点的集合,
集合A2是划分平面F4后取得标记点的集合,
集合B1是划分平面F2后取得标记点的集合,
集合B2是划分平面F5后取得标记点的集合,
集合C1是划分平面F3后取得标记点的集合,
集合C2是划分平面F6后取得标记点的集合,
A1={aij(x,y,z)|x=a,-b≤y≤b,-c≤z≤c},
A2={aij'(x,y,z)|x=-a,-b≤y≤b,-c≤z≤c},
其中,aij是集合A1中的点,aij'是集合A2中的点,0≤i≤X'-1,0≤j≤Y'-1;
B1={bmn(x,y,z)|y=b,-c≤z≤c,-a≤x≤a},
B2{bmn'(x,y,z)|y=-b,-c≤z≤c,-a≤x≤a},
其中,bmn是集合B1中的点,bmn'是集合B2中的点,0≤m≤Y'-1,0≤n≤Z'-1;
C1={cpq(x,y,z)|z=c,-a≤x≤a,-b≤y≤b},
C2={cpq'(x,y,z)|z=-c,-a≤x≤a,-b≤y≤b},
其中,cpq是集合C1中的点,cpq'是集合C2中的点,0≤p≤X'-1,0≤q≤Z'-1;
每组的两个集合中,下标相同的点为一对对称标记点,即
集合A1中的aij与集合A2中的aij',当ij相同时,aij与aij'是一对对称标记点;
集合B1中的bmn与集合B2中的bmn',当mn相同时,bmn与bmn'是一对对称标记点;
集合C1中的cpq与集合C2中的cpq',当pq相同时,cpq与cpq'是一对对称标记点。
实施例:在所述步骤(2)中,根据空间分辨率划分包围盒的步骤为:
计算三维模型顶点个数N=5844126与包围盒的体积M=520之比N/M=12365>100,将包围盒的六个面划分为单元格,本示例划分效果如图5所示,将划分的六个面中相对的面设置为一组,即可得到三组分别有2*30*20、2*20*11、2*30*11个标记点的数据集A=A1∪A2,B=B1∪B2,C=C1∪C2,
本示例其中一面的投影示意图如图6所示;
优选的,所述步骤(2)中,计算标记点与对称特征点的匹配度的步骤为:
步骤(2-1):所有标记点做投影映射到三维模型表面,三组数据集合A、B、C根据投影点落在三维模型的位置,分别取投影点构成三组投影点集合Ffa、Ffb、Ffc
A→Ffa(A1→Ffa1,A2→Ffa2);
B→Ffb(B1→Ffb1,B2→Ffb2);
C→Ffc(C1→Ffc1,C2→Ffc2);
Ffa=Ffa1∪Ffa2
Ffb=Ffb1∪Ffb2
Ffc=Ffc1∪Ffc2
其中,
Ffa是集合A映射到三维模型表面的投影点构成的集合,
Ffb是集合B映射到三维模型表面的投影点构成的集合,
Ffc是集合C映射到三维模型表面的投影点构成的集合,
Ffa1是集合A1映射到三维模型表面的投影点构成的集合,
Ffa2是集合A2映射到三维模型表面的投影点构成的集合,
Ffb1是集合B1映射到三维模型表面的投影点构成的集合,
Ffb2是集合B2映射到三维模型表面的投影点构成的集合,
Ffc1是集合C1映射到三维模型表面的投影点构成的集合,
Ffc2是集合C2映射到三维模型表面的投影点构成的集合。
Ffa1={faij(x,y,z)|x≤a,-b≤y≤b,-c≤z≤c},
Ffa2={faij'(x,y,z)|-a≤x,-b≤y≤b,-c≤z≤c},
其中,faij是集合Ffa1中的点,aij'是集合Ffa2中的点,0≤i≤X'-1,0≤j≤Y'-1;
Ffb1={fbmn(x,y,z)|y≤b,-a≤x≤a,-c≤z≤c},
Ffb2={fbmn'(x,y,z)|-b≤y,-a≤x≤a,-c≤z≤c},
其中fbmn是集合Ffb1中的点,fbmn'是集合Ffb2中的点,0≤m≤Y'-1,0≤n≤Z'-1;
Ffc1={fcpq(x,y,z)|z≤c,-a≤x≤a,-b≤y≤b},
Ffc2={fcpq'(x,y,z)|-c≤z,-a≤x≤a,-b≤y≤b},
其中fcpq是集合Ffc1中的点,fcpq'是集合Ffc2中的点,0≤p≤X'-1,0≤q≤Z'-1;
每组的两个集合中,下标相同的点为一对对称投影点,即:
集合Ffa1中的faij与集合Ffa2中的faij',当ij相同时,faij与faij'是一对对称投影点;
集合Ffb1中的fbmn与集合Ffb2中的fbmn',当mn相同时,fbmn与fbmn'是一对对称投影点;
集合Ffc1中的fcpq与集合Ffc2中的fcpq',当pq相同时,fcpq与fcpq'是一对对称投影点。
步骤(2-2):判断三组标记点集合中每组的每一对标记点是否都在三维模型表面有投影:
若每一对标记点都有一对投影点时,这一对投影点的坐标(x1,y1,z1)与(x2,y2,z2)只有一对数值不同,假设是z1和z2不同(A组标记点集合是x1和x2不同,B组标记点集合是y1和y2不同,C组标记点集合是z1和z2不同),且数值偏差e=|z1-z2|小于等于偏差测度E,则记为匹配;
反之,若数值偏差e=|z1-z2|大于偏差测度E,则记为不匹配;
若这一对标记点至少有一个在三维模型上没有投影点时,该投影点坐标记为(0,0,0)并
,认为其不匹配。
步骤(2-3):将匹配的标记点对数R与该组标记点对数G的比值作为该组匹配度H,即
Figure BDA0001558971680000091
其中,
Figure BDA0001558971680000092
Ha表示标记点集合A的匹配度,
Figure BDA0001558971680000093
Hb表示标记点集合B的匹配度,
Figure BDA0001558971680000094
Hc表示标记点集合C的匹配度;
匹配度H越高说明该组数据的对称性越明显,取匹配度H最高的组max(Ha,Hb,Hc)对应的投影点集合为对称数据组T,
当max(Ha,Hb,Hc)=Ha时,T=Ffa
当max(Ha,Hb,Hc)=Hb时,T=Ffb
当max(Ha,Hb,Hc)=Hc时,T=Ffc
步骤(2-4):过质心做与对称数据组所在平面的平行面为预设对称平面。
在所述步骤(2)中,计算匹配度的步骤为:
第一步:若两个相对的标记点都在三维模型上有投影点,即点坐标不等于(0,0,0)时,以Fa为例:Ffa1中的点faij的下标ij与Ffa2中的点faij'的下标ij相等时,这两个点即为一组相对的点,两者的坐标(x1,y1,z1)与(x2,y2,z2)只有z数值不同,而这个数值偏差e=|z1-z2|在误差测度E以内,则记为匹配的对称点,反之不是匹配对;若两个相对的标记点至少有一个在三维模型上没有投影点,即至少其中一个坐标为(0,0,0),则不认为其是匹配对。
第二步:匹配的对称点对数与该组标记点对数的比值就作为该组匹配度P,匹配度越高说明该组数据的对称性越明显,于是取匹配度P最高的对称点集作为对称数据组;在本示例中,T1的匹配度为421,T2的匹配度为95,T3的匹配度为113,所以取集合T1作为对称数据组。
优选的,所述步骤(3)步骤如下:
步骤(3-1):对称数据组T的两个数据点集Ffr1=fruv(x,y,z)和Ffr2=fruv'(x,y,z)
当T=Ffa时,r=a,u=i,v=j,Ffa1=faij(x,y,z),Ffa2=faij'(x,y,z);
当T=Ffb时,r=b,u=m,v=n,Ffb1=fbmn(x,y,z),Ffb2=fbmn'(x,y,z);
当T=Ffc时,r=c,u=p,v=q,Ffc1=fcpq(x,y,z),Ffc2=fcpq'(x,y,z)。
根据点的下标uv,取u和v中较大的数为子集个数,将T分成若干份子集,T={T0',T1',...,TW'}例如:
当u>v时,W=0,...u-1,
T0'=fr0v(x,y,z)∪fr0v'(x,y,z);
T1'=fr1v(x,y,z)∪fr1v'(x,y,z);
……
Tu-1'=fr(u-1)v(x,y,z)∪fr(u-1)v'(x,y,z);
当u<v时,W=0,...v-1,
T0'=fru0(x,y,z)∪fru0'(x,y,z);
T1'=fru1(x,y,z)∪fru1'(x,y,z);
……
Tv-1'=fru(v-1)(x,y,z)∪fru(v-1)'(x,y,z);
即Tw'=Tw'1∪Tw'2,Tw'1表示每一份子集中来自集合Ffr1=fruv(x,y,z)中的点构成的子集,Tw'2表示每一份子集中来自集合Ffr2=fruv(x,y,z)中的点构成的子集;每份子集Tw'中的点都分布在同一平面上,该平面分别与得到对称数据组T的两个相对面垂直,本示例中某模型划分的若干子集平面示意图如图7,即在某一平面内有若干离散数据点(xv,yv,zv),v=0,1,...,n-1,该平面的方程根据所在区间表示为
A1x+B1y+C1z+D1=0,(5)
x,y,z为变量,a',b',c',d'为常数;
步骤(3-2):三维平面中一条直线由两个相交的平面表示,设另一个平面方程为:
A'x+B'y+C'z+D'=0,(6)
其中,A'2+B'2+C'2≠0,x,y,z为变量,A',B',C',D'为常数;
化简可得:
Figure BDA0001558971680000111
记:
Figure BDA0001558971680000112
代入(7)
得:z=a0x+a1y+a2, (8)
步骤(3-3):对于每一份子集Tw'中的离散的n个数据点,n≥3,(xv,yv,zv)v=0,1,...,n-1;要用点(xv,yv,zv)v=0,1,...,n-1拟合计算公式(7)的平面方程,从几何意义上来说,就是寻求与给出的数据点的距离的平方和最小的拟合直线,则使S最小:
Figure BDA0001558971680000113
其中,S表示直线到给出的数据点的距离的平方和;要使得S最小,应满足
Figure BDA0001558971680000114
k=0,1,2,
Figure BDA0001558971680000115
表示对ak求S的偏导数;
即:
Figure BDA0001558971680000121
有:
Figure BDA0001558971680000122
或,
Figure BDA0001558971680000123
解公式(10)-(12)线性方程组,得:a0,a1,a2,代入(8)得
z=a0x+a1y+a2,(13)
所得的拟合直线可由方程(5)和方程(13)表示为:
Figure BDA0001558971680000124
A1,B1,C1,D1,a0,a1,a2为常数,x,y,z为变量。
步骤(3-4):若干份对称数据组内的离散数据点经过拟合后得出若干条直线,连线成面得到一条初始对称面。
本示例中,某模型的某子集平面拟合直线的示意图如图8所示;
30份数据组内的离散数据点经过拟合后得出30条直线,连线成面可得到一条初始对称面,本示例中,某模型的初始对称面如图9(a)所示;
优选的,在所述步骤(4)中,所述的检验校正初始对称面的步骤如下:
步骤(4-1):将初始对称面与预设对称平面进行比对,初始对称面与预设对称平面有部分平行或重合,则该平行或重合部分定义为完全对称部分,其余部分定义为非完全对称部分;
步骤(4-2):对于非完全对称部分,则将步骤(2)中的等间距缩小一半,重复执行步骤(2)-(4)。本示例中某模型的对称检测效果对比图9(b)和图9(c)所示。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (10)

1.基于最小二乘法的三维模型对称性特征检测方法,其特征是,包括以下步骤:
对三维模型的最小包围盒的六个平面的每个平面,依据空间分辨率进行单元格划分,在划分的最小包围盒的每个平面上选取若干个等间距标记点,将划分的六个平面的每两个相对的面设置为一组,将每组标记点投影到三维模型表面得到投影点,根据投影点是否构成匹配的对称点计算出两个标记点的匹配度,根据匹配度筛选对称数据组;
将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面。
2.如权利要求1所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,在对三维模型的最小包围盒的六个平面的每个平面,依据空间分辨率进行划分步骤之前,还包括:建立三维模型的最小包围盒,建立新的坐标系。
3.如权利要求1所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,在找到初始对称面步骤之后,还包括:
检验校正初始对称面,不断细化初始对称面直到得到准确的对称面。
4.如权利要求2所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,所述建立新的坐标系,是以最小包围盒的中心点为坐标轴原点,建立新的坐标系。
5.如权利要求1所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,根据空间分辨率进行单元格划分,包括步骤如下:
计算三维模型顶点个数与最小包围盒的体积的比值,将比值与设定阈值进行比较,如果比值大于等于设定阈值,则将最小包围盒的六个平面按照比值等间距划分为若干个单元格;如果比值小于设定阈值,则将最小包围盒的六个平面按照像素划分为若干个单元格。
6.如权利要求1所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,计算标记点与投影点的匹配度,包括步骤如下:
两个对称的标记点都有两个投影点,如果所述两个投影点的横坐标、纵坐标与竖坐标三项中有两项的数值一致,只有一项不一致,那么判断不一致的坐标数值偏差是否在设定的偏差测度范围内,如果是,就判断两个标记点匹配,否则两个标记点不匹配;
如果两个对称的标记点,只有一个投影点或者没有投影点,那么直接判断为两个标记点不匹配;
如果两个标记点匹配,那么计算匹配的标记点对数与该组标记点对数的比值,得到的比值就是标记点与投影点的匹配度。
7.如权利要求1所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,根据匹配度筛选对称数据组的步骤如下:
匹配度值最高的标记点组对应的投影点集合,即为对称数据组;
或者,
将对称数据组根据最小二乘法进行平面反射对称性的检测,找到初始对称面的步骤如下:
将每份对称数据组内的离散数据点经过拟合后得出对应的直线,将所有的直线连接成面得到一条初始对称面。
8.如权利要求3所述的基于最小二乘法的三维模型对称性特征检测方法,其特征是,所述的检验校正初始对称面的步骤如下:
将初始对称面与预设对称平面进行比对,初始对称面与预设对称平面有部分平行或重合,则该平行或重合部分定义为完全对称部分,其余部分定义为非完全对称部分;
对于非完全对称部分,则将等间距缩小一半,重复寻找初始对称面的步骤。
9.基于最小二乘法的三维模型对称性特征检测***,其特征是,包括:存储器、处理器以及存储在存储器上并在处理器上运行的计算机程序,所述计算机程序被处理器运行时,完成上述权利要求1-8任一所述方法的步骤。
10.一种计算机可读存储介质,其特征是,其上运行有计算机指令,所述计算机指令被处理器运行时,完成上述权利要求1-8任一所述方法的步骤。
CN201810073992.9A 2018-01-25 2018-01-25 基于最小二乘法的三维模型对称性特征检测方法及*** Expired - Fee Related CN108305289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810073992.9A CN108305289B (zh) 2018-01-25 2018-01-25 基于最小二乘法的三维模型对称性特征检测方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810073992.9A CN108305289B (zh) 2018-01-25 2018-01-25 基于最小二乘法的三维模型对称性特征检测方法及***

Publications (2)

Publication Number Publication Date
CN108305289A CN108305289A (zh) 2018-07-20
CN108305289B true CN108305289B (zh) 2020-06-30

Family

ID=62866464

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810073992.9A Expired - Fee Related CN108305289B (zh) 2018-01-25 2018-01-25 基于最小二乘法的三维模型对称性特征检测方法及***

Country Status (1)

Country Link
CN (1) CN108305289B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4086849A1 (en) * 2021-12-24 2022-11-09 Univerza v Mariboru Procedure for detecting the reflectional symmetry of geometrical objects

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110377997B (zh) * 2019-07-10 2022-11-18 武汉理工大学 一种建筑铝模板配模设计漏板检查方法及漏板检查***
CN110706276B (zh) * 2019-09-27 2022-05-10 中国科学院软件研究所 一种基于骨架引导的三维模型内蕴对称性检测方法
CN112802201B (zh) * 2019-10-25 2023-12-29 北京博超时代软件有限公司 实体模型间并行最近距离获取方法和装置
CN110986788B (zh) * 2019-11-15 2020-10-23 华南农业大学 一种基于三维点云牲畜表型体尺数据的自动测量方法
CN112652075B (zh) * 2021-01-30 2022-08-09 上海汇像信息技术有限公司 对称物体3d模型的对称拟合方法
CN116597117B (zh) * 2023-07-18 2023-10-13 中国石油大学(华东) 一种基于物体对称性的六面体网格生成方法
CN117876554B (zh) * 2024-03-12 2024-05-28 中南建筑设计院股份有限公司 一种基于凸包的板件最小包围盒计算方法和***

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102945569A (zh) * 2012-10-23 2013-02-27 西北工业大学 一种基于热核信号的三维模型对称性分析方法
WO2013075339A1 (en) * 2011-11-25 2013-05-30 Thomson Licensing Methods and apparatus for reflective symmetry based 3d model compression
CN103729510A (zh) * 2013-12-25 2014-04-16 合肥工业大学 基于内蕴变换的三维复杂模型精确镜像对称性计算方法
CN105631932A (zh) * 2015-12-22 2016-06-01 北京航空航天大学 一种轮廓线指导的三维模型重构造方法
CN107146251A (zh) * 2017-03-28 2017-09-08 浙江大学 一种三维人脸模型的对称性分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013075339A1 (en) * 2011-11-25 2013-05-30 Thomson Licensing Methods and apparatus for reflective symmetry based 3d model compression
CN102945569A (zh) * 2012-10-23 2013-02-27 西北工业大学 一种基于热核信号的三维模型对称性分析方法
CN103729510A (zh) * 2013-12-25 2014-04-16 合肥工业大学 基于内蕴变换的三维复杂模型精确镜像对称性计算方法
CN105631932A (zh) * 2015-12-22 2016-06-01 北京航空航天大学 一种轮廓线指导的三维模型重构造方法
CN107146251A (zh) * 2017-03-28 2017-09-08 浙江大学 一种三维人脸模型的对称性分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Efficient 3D reflection symmetry detection: A view-based approach;Bo Li 等;《Graphical Models》;20150926;2-14 *
一种结构对称的破碎青铜文物的数字修复方法;税午阳 等;《文物保护与考古科学》;20151130;第27卷(第4期);113-117 *
医疗三维几何模型的对称面定位研究;程小富;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180115(第01期);I138-1246 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4086849A1 (en) * 2021-12-24 2022-11-09 Univerza v Mariboru Procedure for detecting the reflectional symmetry of geometrical objects

Also Published As

Publication number Publication date
CN108305289A (zh) 2018-07-20

Similar Documents

Publication Publication Date Title
CN108305289B (zh) 基于最小二乘法的三维模型对称性特征检测方法及***
CN106909877B (zh) 一种基于点线综合特征的视觉同时建图与定位方法
CN107610061B (zh) 一种基于二维投影的保边点云孔洞修补方法
CN102930246B (zh) 一种基于点云片段分割的室内场景识别方法
CN108171780A (zh) 一种基于激光雷达构建室内真实三维地图的方法
Mokhtarian et al. Multi-scale free-form 3D object recognition using 3D models
CN112347550B (zh) 耦合式室内三维语义建图及建模方法
JP2007502473A (ja) 有界三次元対象物の分類と空間ローカライゼーションのためのコンピュータ視覚システム
CN107622530B (zh) 一种高效鲁棒的三角网切割方法
US20240153123A1 (en) Isogeometric Analysis Method Based on a Geometric Reconstruction Model
US20160138914A1 (en) System and method for analyzing data
CN105678252A (zh) 依据人脸三角网格自适应细分和高斯小波的迭代插值方法
CN116662600B (zh) 一种基于轻量结构化线地图的视觉定位方法
Lee et al. Region matching and depth finding for 3D objects in stereo aerial photographs
CN111681300B (zh) 轮廓素描线组成的目标区域获取方法
Lee et al. A new hybrid method for 3D object recovery from 2D drawings and its validation against the cubic corner method and the optimisation-based method
CN115861247A (zh) 一种高分辨率遥感影像轮廓多级正则化方法、***及应用
CN115661374A (zh) 一种基于空间划分和模型体素化的快速检索方法
CN111914832A (zh) 一种rgb-d相机在动态场景下的slam方法
Ruan et al. Efficient exact collision detection between ellipsoids and superquadrics via closed-form Minkowski sums
CN113744324A (zh) 一种多相似性测度结合的立体匹配方法
CN106909874A (zh) 一种鼻尖定位方法和装置
CN114926536B (zh) 一种基于语义的定位与建图方法、***及智能机器人
CN106023314A (zh) 一种基于回转轴方向映射的b样条母曲线拟合方法
CN113313629B (zh) 交叉路口自动识别方法、***及其模型保存方法、***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
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: 20200630

Termination date: 20220125