CN103729510B - 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 - Google Patents
基于内蕴变换的三维复杂模型精确镜像对称性计算方法 Download PDFInfo
- Publication number
- CN103729510B CN103729510B CN201310733139.2A CN201310733139A CN103729510B CN 103729510 B CN103729510 B CN 103729510B CN 201310733139 A CN201310733139 A CN 201310733139A CN 103729510 B CN103729510 B CN 103729510B
- Authority
- CN
- China
- Prior art keywords
- symmetry
- model
- plane
- conversion
- cos
- 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.)
- Active
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于内蕴变换的三维复杂模型精确镜像对称性计算方法。首先交互设定或计算模型的粗略对称面,然后根据粗略对称面,变换给定模型,使粗略对称面与一坐标平面重合,接着根据镜像对称模型的内蕴变换,计算精确对称面;最后变换模型及精确对称面,得到模型原始位置的精确对称面。本发明基于镜像对称模型的严格内蕴变换关系,精度可控,具有较强泛化性能,广泛适用于工程应用领域。
Description
技术领域
本发明涉及三维几何建模与计算领域,具体为一种基于内蕴变换的三维复杂模型精确镜像对称性计算方法。
背景技术
对称是自然界普遍存在的现象,无论人造的还是天然的物体,其外形大都存在不同程度的对称性,其中以镜像对称为主。本发明主要针对工程中经过人为加工的具有镜像对称性的模型,特别是根据这些实物模型上采集的数据进行三维重建或误差分析中的对称性分析计算问题。
同天然物体对比,人造物体的对称性精度要高很多,对这类模型进行对称性分析计算的精度要求一般较高。已有的对称性分析计算算法主要针对计算机动画、计算机视觉等领域应用,精度普遍无法满足工程应用需求。
随着科学技术的发展,具有对称性的复杂外形零部件在工程中应用越来越广泛,特别是在飞机、船舶、汽车、高速列车等制造业中,产品的外形都具有镜像对称性。这种镜像对称性影响了产品的外观质量,对飞机、轮船、跑车、高速列车等类型产品的操控性能还具有决定性的影响。在对这类产品或实物模型的外形进行三维重建时,保证重建模型的对称性是非常重要的;另外,目前还缺乏有效手段对这类产品外形进行对称误差的分析计算。因此,研究开发一种精确镜像对称性计算方法是非常必要的。
发明内容
本发明的目的是提供一种基于内蕴变换的三维复杂模型精确镜像对称性计算方法,用于基于实物模型的三维重建或产品外形的对称误差检测等精度要求较高的工程应用领域。
为了达到上述目的,本发明所采用的技术方案为:
基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:包括以下步骤:
(1)、交互设定或计算给定模型的粗略对称面;
(2)、根据步骤(1)得到的粗略对称面,变换给定模型,使粗略对称面与一个坐标平面重合,得到变换后模型;
(3)、根据镜像对称模型的内蕴变换,计算步骤(2)中得到的变换后模型的精确对称面;
(4)、再次变换步骤(2)得到的变换后模型及精确对称面,得到模型原始位置的精确对称面。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:所述步骤(1)中,交互设定或计算粗略对称面的方法为:将视图旋转至近似垂直于理论对称面位置,交互设定一条直线,从而粗略对称面的位置,粗略对称面由一个点P=(x,y,z)T和一个法矢V=(a,b,c)T定义,其中x、y、z分别为P的三个坐标值,a、b、c分别为矢量V的三个分量值,a2+b2+c2=1。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:所述步骤(2)中,变换给定模型用到齐次变换矩阵,齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:所述步骤(3)中,设给定模型为M,给定模型M与其关于x-z坐标平面的镜像对称模型Mx-z的内蕴变换为:M=RzRyRxRxR-yRzMx-z+(RzRyRxRxR-yRz+I)T,其中
令Rx-z=RzRyRxRxR-yRz,其元素满足以下关系:
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:步骤(3)按以下过程实现:
(a)、采用模型注册算法,得到Rx-Z和(Rx-z+I)T的值;
(b)、根据S301中Rx-z的计算结果,得到一个非线性方程组,采用非线性方程组求解算法求解,得到sin(θx),cos(θx),sin(θy),cos(θy),sin(θz),cos(θz)的值;
(c)、根据步骤(a)和步骤(b)的计算结果,得到一个线性方程组,求解该线性方程组得到Tx,Ty,Tz的值;
(d)、令y=(0,1,0)T,模型M的精确对称面的法矢为RzRyRxy,其上一点为(Tx,Ty,Tz)T。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:所述步骤(4)中,变换步骤(2)得到的变换后模型及精确对称面用到齐次变换矩阵,齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值,将模型变换回原始位置并将对称面变换回模型原始位置对应的对称面。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:通过设定粗略对称面,提高了步骤(3)中模型注册算法计算的可靠性和稳定性。
所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:通过设定模型注册算法中的收敛条件,可以控制镜像对称性计算的精度,满足高精度工程应用的要求。
本发明用于基于实物模型的三维重建或产品外形的对称误差检测等应用领域,具有计算稳定、可靠,结果精度可控的优点,能够满足高精度工程应用的要求。
附图说明
图1为本发明的流程图。
图2为针对一实测采集的汽车前盖点云数据,采用本发明所述的三维复杂模型的精确镜像对称性计算方法,得到的镜像对称面计算结果示意图。
图3为采用发明所述的三维复杂模型的精确镜像对称性计算方法,计算轴对称模型的对称轴的示意图。
具体实施方式
基于内蕴变换的三维复杂模型精确镜像对称性计算方法,包括以下步骤:
(1)、交互设定或计算给定模型的粗略对称面;
(2)、根据步骤(1)得到的粗略对称面,变换给定模型,使粗略对称面与一个坐标平面重合,得到变换后模型;
(3)、根据镜像对称模型的内蕴变换,计算步骤(2)中得到的变换后模型的精确对称面;
(4)、再次变换步骤(2)得到的变换后模型及精确对称面,得到模型原始位置的精确对称面。
步骤(1)中,交互设定或计算粗略对称面的方法为:将视图旋转至近似垂直于理论对称面位置,交互设定一条直线,从而粗略对称面的位置,粗略对称面由一个点P=(x,y,z)T和一个法矢V=(a,b,c)T定义,其中x、y、z分别为P的三个坐标值,a、b、c分别为矢量V的三个分量值,a2+b2+c2=1。
步骤(2)中,变换给定模型用到齐次变换矩阵,齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值。
步骤(3)中,设给定模型为M,给定模型M与其关于x-z坐标平面的镜像对称模型Mx-z的内蕴变换为:
M=RzRyRxRxR-yRzMx-z+(RzRyRxRxR-yRz+I)T,其中
令Rx-z=RzRyRxRxR-yRz,其元素满足以下关系:
步骤(3)按以下过程实现:
(a)、采用模型注册算法,得到Rx-z和(Rx-z+I)T的值;
(b)、根据S301中Rx-z的计算结果,得到一个非线性方程组,采用非线性方程组求解算法求解,得到sin(θx),cos(θx),sin(θy),cos(θy),sin(θz),cos(θz)的值;
(c)、根据步骤(a)和步骤(b)的计算结果,得到一个线性方程组,求解该线性方程组得到Tx,Ty,Tz的值;
(d)、令y=(0,1,0)T,模型M的精确对称面的法矢为RzRyRxy,其上一点为(Tx,Ty,Tz)T。
步骤(4)中,变换步骤(2)得到的变换后模型及精确对称面用到齐次变换矩阵,齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值,将模型变换回原始位置并将对称面变换回模型原始位置对应的对称面。
通过设定粗略对称面,提高了步骤(3)中模型注册算法计算的可靠性和稳定性。
通过设定模型注册算法中的收敛条件,可以控制镜像对称性计算的精度,满足高精度工程应用的要求。
下面结合附图1至3和实施例,对发明的具体实施方式做进一步描述。以下实施例仅用于说明本发明,但不用来限制本发明的范围。
流程图如图1所示的一种三维复杂模型的精确镜像对称性计算方法,主要包括步骤:
S1.交互设定或计算模型的粗略对称面。交互设定粗略对称面的方法为:将视图旋转至近似垂直于理论对称面位置,交互设定一条直线,从而粗略对称面的位置。粗略对称面由一个点P=(x,y,z)T和一个法矢V=(a,b,c)T定义,其中a2+b2+c2=1。模型的粗略对称面的也可采用已有的针对其他应用的对称面计算方法,但是从工程应用角度看,这些算法一般比较复杂、计算费时。
S2.根据粗略对称面,变换给定模型,使粗略对称面与一坐标平面重合。给定模型的齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别P的三个坐标值,a、b、c分别为法矢V的三个分量值。假设模型上点的坐标为(xi,yi,zi)T,i=O,…,n,变换后的齐次坐标为(x'i,y'i,z'i,1)T=Tm(xi,yi,zi,1)T,i=0,…,n。坐标变换的目的是为了保证S301中模型注册时具有良好的初值,提高算法的收敛性。
S3.根据镜像对称模型的内蕴变换,计算精确对称面。模型M由(x′i,y′i,z′i)T,i=0,…,n定义,M与其关于x-z坐标平面的镜像对称模型Mx-z的内蕴变换为:
M=RzRyRxRxR-yRzMx-z+(RZRyRxRxR-yRz+I)T
其中 T=(Tx,Ty,Tz)T。Rx、Ry、Rz分别为绕x、y、z轴的旋转变换矩阵,θx、θy、θz分别为绕x、y、z轴的旋转角度,T为平移变换矩阵,均为未知量。
本实施例中该步骤主要包括:
S301:令Rx-z=RzRyRxRxR-yRz,针对M与Mx-z,采用模型注册算法,得到Rx-z和(Rx-z+I)T的值。模型注册常用的算法就是ICP(Iterative Closest Points)算法及其改进算法,在计算中采用点到平面距离模式,在这种模式下ICP算法的收敛性较好。根据具体应用需要,也可以采用其他模型注册或模型匹配算法。
S302:Rx-z为一3x3的矩阵,其元素满足以下关系:
因此得到一个变量为sin(θx),cos(θx),sin(θy),cos(θy),sin(θz),cos(θz)的6个方程的非线性方程组。为了保证变量之间的关系,需再添加三个方程,即sin(θx)2+cos(θx)2=1、sin(θy)2+cos(θy)2=1、sin(θz)2+cos(θz)2=1。采用Levenberg-Marquardt算法或其他非线性方程组求解算法,得到sin(θx),cos(θx),sin(θy),cos(θy),sin(θz),cos(θz)的值;
S303:根据S301和S302的计算结果,得到变量为Tx,Ty,Tz的一个三元齐次线性方程组,求解该线性方程组得到Tx,Ty,Tz的值;
S304:令y=(0,1,0),模型M的精确对称面的法矢为
(vx,vy,vz)=RzRyRxy,其上一点为(Tx,Ty,Tz)T。
S4.再次变换步骤(2)变换后模型及精确对称面,得到模型原始位置的精确对称面。模型及精确对称面的齐次变换矩阵为:
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别P的三个坐标值,a、b、c分别为法矢V的三个分量值。将模型变换回原始位置,即 同时将对称面变换回模型原始位置对应的对称面位置。即,给定模型的对称面由一点P=(x,y,z)T和一个法矢V=(a,b,c)T定义,其中 图2给出了一个采用本发明所述的对称性计算方法,对实测采集的汽车前盖点云数据计算的镜像对称面计算结果。
以上实施方式仅用于说明本发明,而非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型(例如可以采用本发明计算轴对称模型的两个对称面,其交即为模型的对称轴,如图3所示),因此所有等同的技术方案也属于本发明的保护范畴。
Claims (5)
1.基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:包括以下步骤:
(1)、交互设定或计算给定模型的粗略对称面;
(2)、根据步骤(1)得到的粗略对称面,变换给定模型,使粗略对称面与一个坐标平面重合,得到变换后模型;
(3)、根据镜像对称模型的内蕴变换,计算步骤(2)中得到的变换后模型的精确对称面;
(4)、再次变换步骤(2)得到的变换后模型及精确对称面,得到模型原始位置的精确对称面;
所述步骤(2)中,变换给定模型用到齐次变换矩阵,齐次变换矩阵为:
,
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值;
所述步骤(3)中,设给定模型为M,给定模型M与其关于x-z坐标平面的镜像对称模型Mx-z的内蕴变换为:M=RzRyRxRxR-yRzMx-z+(RzRyRxR-yRz+1)T,
其中 T=(Tx,Ty,Tz)T,Rx、Ry、Rz分别为绕x、y、z轴的旋转变换矩阵,θx、θy、θz分别为绕x、y、z轴的旋转角度,T为平移变换矩阵,均为未知量;
令Rx-z=RzRyRxR-yRz,其元素满足以下关系:
所述步骤(4)中,变换步骤(2)得到的变换后模型及精确对称面用到齐次变换矩阵,齐次变换矩阵为:
,
其中α为粗略对称面的法矢V与Y轴之间的夹角,x、y、z分别为点P的三个坐标值,a、b、c分别为矢量V的三个分量值,将模型变换回原始位置并将对称面变换回模型原始位置对应的对称面。
2.根据权利要求1所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:所述步骤(1)中,交互设定或计算粗略对称面的方法为:将视图旋转至近似垂直于理论对称面位置,交互设定一条直线,从而计算粗略对称面的位置,粗略对称面由一个点P=(x,y,z)T和一个法矢V=(a,b,c)T定义,其中x、y、z分别为P的三个坐标值,a、b、c分别为矢量V的三个分量值,a2+b2+c2=1。
3.根据权利要求1所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:步骤(3)按以下过程实现:
(a)、采用模型注册算法,得到Rx-z和(Rx-z+1)T的值;
(b)、根据Rx-z的计算结果,得到一个非线性方程组,采用非线性方程组求解算法求解,得到sin(θx),cos(θx),sin(θy),cos(θy),sin(θz),cos(θz)的值;
(c)、根据步骤(a)和步骤(b)的计算结果,得到一个线性方程组,求解该线性方程组得到Tx,Ty,Tz的值;
(d)、令y=(0,1,0)T,模型M的精确对称面的法矢为RzRyRxy,其上一点为(Tx,Ty,Tz)T。
4.根据权利要求3所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:通过设定粗略对称面,提高了步骤(3)中模型注册算法计算的可靠性和稳定性。
5.根据权利要求3所述的基于内蕴变换的三维复杂模型精确镜像对称性计算方法,其特征在于:通过设定模型注册算法中的收敛条件,可以控制镜像对称性计算的精度,满足高精度工程应用的要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310733139.2A CN103729510B (zh) | 2013-12-25 | 2013-12-25 | 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310733139.2A CN103729510B (zh) | 2013-12-25 | 2013-12-25 | 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103729510A CN103729510A (zh) | 2014-04-16 |
CN103729510B true CN103729510B (zh) | 2016-09-14 |
Family
ID=50453584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310733139.2A Active CN103729510B (zh) | 2013-12-25 | 2013-12-25 | 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103729510B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204740B (zh) * | 2016-07-27 | 2020-02-07 | 扬州大学 | 基于切片数据对称性分析的三维缺损颜面模型重建方法 |
CN106370100A (zh) * | 2016-08-17 | 2017-02-01 | 北京汽车股份有限公司 | 车身对称偏差检测方法及*** |
CN106484965B (zh) * | 2016-09-21 | 2018-05-29 | 南京信息工程大学 | 基于点云数据的工业构件仿真安装方法 |
CN108305289B (zh) * | 2018-01-25 | 2020-06-30 | 山东师范大学 | 基于最小二乘法的三维模型对称性特征检测方法及*** |
CN110360944B (zh) * | 2019-06-25 | 2021-03-26 | 江苏省特种设备安全监督检验研究院 | 一种基于三维点云的吊钩形变监测与显示方法 |
CN113012294B (zh) * | 2021-03-30 | 2022-12-13 | 上海汇像信息技术有限公司 | 三维模型的对称重建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101131730A (zh) * | 2007-09-25 | 2008-02-27 | 浙江大学 | 一种弱化表情形变影响的三维人脸识别方法 |
CN101315665A (zh) * | 2008-06-27 | 2008-12-03 | 浙江大学 | 一种基于三维鼻形的身份识别方法 |
CN101937579A (zh) * | 2010-09-20 | 2011-01-05 | 南京大学 | 一种利用透视草图创建三维曲面模型的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070080967A1 (en) * | 2005-10-11 | 2007-04-12 | Animetrics Inc. | Generation of normalized 2D imagery and ID systems via 2D to 3D lifting of multifeatured objects |
-
2013
- 2013-12-25 CN CN201310733139.2A patent/CN103729510B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101131730A (zh) * | 2007-09-25 | 2008-02-27 | 浙江大学 | 一种弱化表情形变影响的三维人脸识别方法 |
CN101315665A (zh) * | 2008-06-27 | 2008-12-03 | 浙江大学 | 一种基于三维鼻形的身份识别方法 |
CN101937579A (zh) * | 2010-09-20 | 2011-01-05 | 南京大学 | 一种利用透视草图创建三维曲面模型的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103729510A (zh) | 2014-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103729510B (zh) | 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 | |
Wu et al. | Hand-eye calibration: 4-D procrustes analysis approach | |
CN106248035A (zh) | 一种基于点云模型的面轮廓度精确检测的方法及*** | |
CN101329764B (zh) | 采用两个任意共面圆进行摄像机标定的方法 | |
Ren et al. | Coordinate transformation uncertainty analysis in large-scale metrology | |
CN106568381B (zh) | 一种基于标准平面的线激光测量***标定方法 | |
US20170371022A1 (en) | System and method for determining geo location of a target using a cone coordinate system | |
CN106248014A (zh) | 一种基于单相片的三维坐标测量方法及装置 | |
CN105205824A (zh) | 基于高精度辅助相机和球靶标的多摄像机全局校准方法 | |
CN103591949A (zh) | 三轴姿态测量***非正交性误差的正交补偿方法 | |
CN108759665A (zh) | 一种基于坐标转换的空间目标三维重建精度分析方法 | |
CN109753549A (zh) | 一种大圆航线距离计算方法 | |
CN104517291A (zh) | 基于目标同轴圆特征的位姿测量方法 | |
CN102721375B (zh) | 一种强反光金属结构件的在位测量中多次反光抑制方法 | |
CN103175470A (zh) | 基于线结构光视觉传感器的基准球定位测量方法 | |
CN103810313B (zh) | 一种stl模型到空间分割模型的转换方法 | |
Ravishankar et al. | Automated inspection of aircraft parts using a modified ICP algorithm | |
CN102087117B (zh) | 飞船交会对接用测距敏感器精度的地面测量方法 | |
CN105300382B (zh) | 一种基于惯性测量单元的大角度船体变形测量方法 | |
CN107576275A (zh) | 一种利用摄影测量技术对充气结构进行应变场测量的方法 | |
CN104634248A (zh) | 一种双目视觉下的转轴标定方法 | |
CN107091627A (zh) | 一种多表面***的综合测量与评估方法 | |
Liu et al. | Precise pose and radius estimation of circular target based on binocular vision | |
CN102445171A (zh) | 一种确定自由曲面上任意点处法向矢量的方法 | |
Lu et al. | Sphere-based calibration method for trinocular vision sensor |
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 |