发明名称 密集点云模型快速精确配准方法
摘要 本发明公开了一种密集点云模型快速精确配准方法,通过对齐的方法实现涡轮叶片密集点云模型和CAD模型的初始配准;对得到的初始配准点云数据进行适当的简化,计算SVD_ICP算法中每步迭代的旋转矩阵和平移矩阵,最终采用SVD_ICP算法实现控制点集{P}和CAD模型的精确配准。本发明在实现点云模型与CAD模型预配准的前提下,通过对密集点云数据进行简化,以简化后的数据作为配准控制点集,在SVD_ICP算法的基础上计算求得总的旋转矩阵和平移矩阵,然后把求得的总变换矩阵应用于简化前的原始密集点云,实现其与CAD模型的快速配准,兼顾了配准精度和配准速度。
申请公布号 CN102169579A 申请公布日期 2011.08.31
申请号 CN201110080192.8 申请日期 2011.03.31
申请人 西北工业大学 发明人 程云勇;卜昆;张定华;黄胜利;董一巍;黄魁东;张现东;窦杨柳;刘金钢;杨小宁;袁帅;傅蒋威;周丽敏;乔燕;张亮
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.一种密集点云模型快速精确配准方法,其特征在于包括下述步骤:步骤1分别选取测量点云模型与CAD模型上任意三组对应点对(P<sub>1</sub>,Q<sub>1</sub>),(P<sub>2</sub>,Q<sub>2</sub>),(P<sub>3</sub>,Q<sub>3</sub>),其中P<sub>1</sub>,P<sub>2</sub>,P<sub>3</sub>不在同一平面;令<img file="FSA00000463496700011.GIF" wi="267" he="167" /><img file="FSA00000463496700012.GIF" wi="308" he="167" />作矢量v<sub>3</sub>和w<sub>3</sub>,<img file="FSA00000463496700013.GIF" wi="434" he="324" />作矢量矢量v<sub>2</sub>和w<sub>2</sub>,<img file="FSA00000463496700014.GIF" wi="303" he="137" />v<sub>1</sub>、v<sub>2</sub>和v<sub>3</sub>构成正交坐标系,表示为单位矢量<img file="FSA00000463496700015.GIF" wi="218" he="216" />同理得<img file="FSA00000463496700016.GIF" wi="235" he="215" />测量点云模型中的任一点P<sub>i</sub>变换到CAD模型上对应点Q<sub>i</sub>的关系式为:Q<sub>i</sub>=RP<sub>i</sub>+T,其中R=[w][v]<sup>-1</sup>,T=Q<sub>1</sub>-[w][v]<sup>-1</sup>P<sub>1</sub>;步骤2利用均匀采样法将步骤1得到的密集点云模型简化至10000-50000个点,以简化后的点云作为控制点集{P},寻找CAD模型上与之对应的点集{P’};步骤3以步骤2得到的控制点集{P}作为配准控制点集,计算SVD_ICP算法中每步迭代的旋转矩阵和平移矩阵:1)利用公式<img file="FSA00000463496700017.GIF" wi="193" he="100" />和<img file="FSA00000463496700018.GIF" wi="221" he="108" />求取两点集{P}和{P’}的质心p和p’;2)根据公式q<sub>i</sub>=p<sub>i</sub>-p和q′<sub>i</sub>=p′<sub>i</sub>-p′计算两点集{P}和{P’}中每个点相对于其质心的位移;3)根据公式<img file="FSA00000463496700019.GIF" wi="225" he="110" />求取矩阵H,并由公式H=UΛV<sup>T</sup>对矩阵H进行奇异值分解,其中U为H的左奇异矩阵,V为H的右奇异矩阵,二者均为3×3的正交矩阵;4)根据公式R=VU<sup>T</sup>和T=p′-Rp计算第i步迭代的旋转矩阵R<sub>i</sub>和平移矩阵T<sub>i</sub>;步骤4采用SVD_ICP算法实现控制点集{P}和CAD模型的精确配准,方法如下:根据步骤3求得每迭代一次的旋转矩阵R<sub>i</sub>和平移矩阵T<sub>i</sub>,从而利用公式R<sub>z</sub>=R<sub>i</sub>*R<sub>i-1</sub>*.....*R<sub>1</sub>和T<sub>z</sub>=R<sub>i</sub>*T<sub>i-1</sub>+T’<sub>i</sub>求得控制点集{P}相对于CAD模型总的旋转矩阵R<sub>z</sub>和平移矩阵T<sub>z</sub>,其中R<sub>z</sub>是3×3的矩阵,T<sub>z</sub>是3×1的矩阵;由公式Q<sub>i</sub>=R<sub>z</sub>P<sub>i</sub>+T<sub>z</sub>,其中P<sub>i</sub>为点云模型中的任一点,Q<sub>i</sub>为变换后与P<sub>i</sub>对应的点;把R<sub>z</sub>和T<sub>z</sub>应用于简化前的原始点云模型,实现其与CAD模型的快速精确配准。
地址 710072 陕西省西安市友谊西路127号