CN112017302A - 基于cad模型的投影标记和机器视觉的实时注册方法 - Google Patents

基于cad模型的投影标记和机器视觉的实时注册方法 Download PDF

Info

Publication number
CN112017302A
CN112017302A CN202010890907.5A CN202010890907A CN112017302A CN 112017302 A CN112017302 A CN 112017302A CN 202010890907 A CN202010890907 A CN 202010890907A CN 112017302 A CN112017302 A CN 112017302A
Authority
CN
China
Prior art keywords
real
virtual
image
camera
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.)
Pending
Application number
CN202010890907.5A
Other languages
English (en)
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.)
Nanjing Aoxiang Intelligent Manufacturing Technology Co ltd
Original Assignee
Nanjing Aoxiang Intelligent Manufacturing Technology Co ltd
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 Nanjing Aoxiang Intelligent Manufacturing Technology Co ltd filed Critical Nanjing Aoxiang Intelligent Manufacturing Technology Co ltd
Priority to CN202010890907.5A priority Critical patent/CN112017302A/zh
Publication of CN112017302A publication Critical patent/CN112017302A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/006Mixed reality
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Architecture (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及AR领域,公开了一种基于CAD模型的投影标记和机器视觉的实时注册方法,步骤为:根据虚拟模型图像和真实物体图像的相似度计算得出对应关系;对真实相机进行标定,确定虚拟模型和真实物体的坐标转换关系;在虚拟模型四周选定投影标记点,通过系列计算获得投影仪坐标系下的投影坐标,再将标记投影于真实四周;计算出自由相机和投影标记点位置的相对位置关系,计算虚拟模型与投影标记点之间的相对位姿,利用虚拟相机拍照获得虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚实注册,能够在不对真实物体张贴标记的情况下,就能够实现投影标记和视觉的虚实注册,并且精度高,实施起来更加方便实用。

Description

基于CAD模型的投影标记和机器视觉的实时注册方法
技术领域
本发明涉及AR领域,更具体地说,它涉及一种基于CAD模型的投影标记和机器视觉的实时注册方法。
背景技术
随着AR技术的发展,虚实注册技术成为了AR领域内的重要支撑技术与研究焦点,通过虚实注册,使虚拟模型与真实物体实时的虚实融合具有重要的应用意义,这就需要实时的识别物体并计算物体的位姿,使虚拟模型可以根据真实物体位姿,实时调整位姿。现有的基于标记的AR注册方法,必须限定物体与标记的相对位置,且在工业背景下,物体不允许张贴标记;基于三维模型各姿态的模板识别的方法,注册精度与模板数量正相关,想达到较好的精度,计算量大,无法保证实时性。基于深度相机采集点云匹配的方法,需要额外的深度相机设备,且计算量极大,同样无法保证实时性。
发明内容
本发明的上述技术目的是通过以下技术方案得以实现的:基于CAD模型的投影标记和机器视觉的实时注册方法,包括以下步骤:
采用虚拟相机对真实尺寸的CAD虚拟模型进行拍摄,得到虚拟模型图像;采用真实相机对真实物体进行拍摄,得到真实物体图像,并根据虚拟模型图像和真实物体图像的相似度计算得出两者之间的对应关系;
在得出虚拟模型图像和真实物体图像之间的对应关系后,对真实相机进行标定,通过虚拟模型图像所在像素坐标系和真实物体图像所在世界坐标系的转换关系,确定虚拟模型和真实物体坐标的转换关系;
在虚拟模型四周选定投影标记点,通过虚拟模型与真实模型之间的坐标转换关系,计算得到投影标记点在真实模型四周的坐标,通过像素坐标系与世界坐标系的转换关系,将投影标记点的世界坐标转化为真实相机的像素坐标;通过相机投影仪标定方法,得到相机像素坐标到投影仪坐标的转换关系,并由此获得投影仪坐标系下的投影坐标,再按投影坐标将标记投影于真实物体四周;
计算出自由相机和投影标记点位置的相对位置关系,从投影标记点位姿得到自由相机位姿;根据自由相机位姿,在虚拟模型所在的虚拟空间坐标系中设置相同位姿的虚拟相机;设置想要叠加于真实物体的虚拟模型相对于真实模型的坐标,计算虚拟模型与投影标记点之间的相对位姿,将虚拟模型按照此相对位姿放置在虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
作为本发明的一种优选技术方案,虚拟模型图像和真实物体图像的获得方法是:将真实尺寸的CAD模型放置在虚拟空间中,并建立虚拟空间坐标系,通过虚拟相机对CAD模型进行拍摄得到虚拟模型图像;将真实相机按照虚拟相机和虚拟模型之间的相对位置对真实物体进行拍摄得到真实物体图像。
作为本发明的一种优选技术方案,虚拟空间坐标系的建立方法为:将虚拟模型底部中心作为虚拟空间坐标原点,以虚拟模型底面为X-Y平面,垂直于模型底面为Z轴。
作为本发明的一种优选技术方案,虚拟模型和真实物体之间对应关系的计算方法为:
A1、将视角向下的虚拟相机设置在虚拟模型的中心上方,使得虚拟相机中心与虚拟模型中心一致;
A2、在X-Y平面上,在[0,2π)范围内,均匀划分为m个角度,每个角度为α,在每一角度上,通过虚拟相机拍摄虚拟模型,获得每一角度姿态下的模型顶层图像模板;设图像模板集合为UA,Um为各角度的图像模板,有:
UA={Um};
所有图像模板对应旋转角为0的基准模板U0的旋转矩阵为:
Figure BDA0002656948900000021
A3、将每张图像模板通过拉普拉斯算子进行处理,取拉普拉斯过零点作为边缘,使得图像模板转换为边缘二值图像;
A4、将真实相机设于真实物体顶部,平行与底面向下实时视频采集真实物体图像,设置真实相机与虚拟相机高度、分辨率、视野大小完全一致,使得真实相机拍摄真实物体的尺寸与虚拟相机拍摄虚拟模型的尺寸完全一致,并且使世界坐标原点为真实相机光轴与桌面交点,此时真实世界坐标与虚拟世界坐标完全对应;
通过A3中的相同方法,对采集的每一帧真实物体图像进行边缘二值处理;
A5、将A4中实时采集的每一帧真实物体图像的边缘二值图像与虚拟模型的图像模板的边缘二值图像进行倒角距离匹配;设Um={ui}与V={vj}设定为模板与真实尺寸的边缘图上的各自边缘点的集合;在Um与V之间的倒角距离由每个ui∈Um的点与V中边缘的最近距离的平均值来决定,其中n为U集合中边缘点的数量:
Figure BDA0002656948900000031
A6、在A5的基础上,对模板U进行一定范围内的旋转平移,迭代缩小倒角距离范围;将U0模板放置在采集图像的中心位置与采集的真实物体图像进行匹配;设W函数为与参数s相关的2D的欧几里得变换函数,其中s∈SE(2),s=(θ,tx,ty),这里tx,ty分别表示为沿着x,y轴方向的平移变换t,θ为平面内的旋转变换;图像上的点x通过W函数执行变换操作:
Figure BDA0002656948900000032
m模板的最佳匹配变换系数sm,通过下式计算获得,得到当前模板与待识别图像之间的最佳旋转平移关系:
Figure BDA0002656948900000033
这里W(Um;s)={W(ui,s)};
A7、根据A6倒角匹配计算公式,遍历各角度的模型图像模板Um∈UA,获取最佳匹配模板
Figure BDA0002656948900000034
与对应的变换系数
Figure BDA0002656948900000035
满足公式如下:
Figure BDA0002656948900000036
Figure BDA0002656948900000037
A8、确定模板坐标与最佳匹配之间的变换矩阵,根据相对于基准模板U0旋转变换矩阵
Figure BDA0002656948900000038
并与
Figure BDA0002656948900000039
与输入图像V最佳匹配的旋转平移变换矩阵
Figure BDA00026569489000000310
中的旋转变换
Figure BDA00026569489000000311
Figure BDA00026569489000000312
求解最佳匹配状态与基准模板U0的旋转变换矩阵
Figure BDA00026569489000000313
Figure BDA0002656948900000041
Figure BDA0002656948900000042
通过变换矩阵
Figure BDA0002656948900000043
实现将虚拟模型模板与真实模型实时图像变换一致,采集图像上的真实模型坐标相对于放置在图像中心的U0模板的像素坐标V为:
Figure BDA0002656948900000044
A9、对真实相机进行标定,选定张正友标定法,计算得到世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,设U0模板像素坐标对应的世界坐标为UW,V真实物体像素坐标对应的世界坐标VW,有公式:
U0=MWCUW
V=MWCVW
且根据步骤A8中像素坐标U0与V之间的转换关系,所以VW与UW之间的转换关系为:
Figure BDA0002656948900000045
Figure BDA0002656948900000046
此时虚拟模型通过
Figure BDA0002656948900000047
转换矩阵,与真实物体的位姿完全一致。
作为本发明的一种优选技术方案,确定投影仪投影标记点的方法为:
B1、在相对于虚拟模型的+x,-x,+y,-y方向上各投影一个不同的人工标记,设偏移量为O,则各标记中心与物体中心的偏移矩阵为:
Figure BDA0002656948900000048
Figure BDA0002656948900000049
设i为投影标记序号,则标记偏移矩阵集合
Figure BDA00026569489000000410
可以表示为:
Figure BDA00026569489000000411
标记坐标Pm相对于虚拟模型UW的坐标为:
Figure BDA00026569489000000412
则标记坐标Pm的像素坐标Pmc,根据步骤A9计算的世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,所以标记对于虚拟模型的像素坐标Pm为:
Pmc=MWCPm;并且,根据步骤A8中,模板像素坐标与采集图像坐标的转换关系为
Figure BDA00026569489000000413
所以标记在真实物体周围的像素坐标Pc可以根据下式计算得到:
Figure BDA0002656948900000051
所以真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW为:
Figure BDA0002656948900000052
B2、通过相机投影仪标定算法,计算得到相机坐标系与投影仪坐标系转换矩阵Mcp,以及投影仪坐标与世界坐标的转换矩阵Mpw,通过以下式子,即可得出投影仪坐标下的标记位置Pp与标记真实位置坐标Pw
Pp=McpPc
Pw=McpMpwPc
投影仪根据Pp坐标将标记投影到真实物体四周;
B3,确定投影标记真实投影位置Pw与真实物体位置VW间的转换关系,根据步骤B1中说明的真实标记像素坐标与虚拟标记像素坐标之间的转换关系,以及B2中,有如下关系:
Figure BDA0002656948900000053
并且根据A9中真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW的转换关系,带入上式子可知真实模型坐标VW与真实物体四周投影标记Pw之间的转换关系为:
Figure BDA0002656948900000054
作为本发明的一种优选技术方案,得到自由相机拍摄下的真实物体投影图像的方法是:
C1、使用ORB特征构建每个标记的特征以及各个特征相对于自身标记的位置,建立标记特征模板特征点集合世界坐标集合
Figure BDA0002656948900000055
通过自由相机拍摄标记,将采集到的图片使用ORB进行处理,得到相机拍摄的标记特征点像素坐标集合
Figure BDA0002656948900000056
将每个标记模板特征与采集图像中的采集到的特征,通过RANSAC随机采样一致性算法,将对应位置的特征点匹配起来;
C2、根据投影标记确定相机位姿,首先通过张正友标定法,计算出自由相机内参MC,设投影的真实标记Pw与自由相机坐标系C之间的变换矩阵为矩阵Mcw,求解公式如下:
Pfc=MCMcwPw
根据Pfc与Pw的特征点之间对应关系,可列n个等式方程,n=特征点数量且满足≥4,即可求解Mcw
C3、计算自由相机C相对于投影标记Pw的位姿,根据C2,可知:
C=McwPw
将虚拟相机按照C位姿放置在虚拟空间中;
C4、设定要叠加在真实模型上的虚拟模型位姿与真实模型位姿VW之间的转换关系,设转换矩阵为Mvr,转换关系如下:
F=MvrVW
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置V间的转换关系,与C4虚拟模型坐标F与真实模型V之间的位置关系,确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure BDA0002656948900000061
Figure BDA0002656948900000062
将虚拟物体按照F位姿放置虚拟空间中;
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置VW之间的转换关系:
Figure BDA0002656948900000063
和C4中虚拟模型坐标F与真实模型VW之间的位置关系:
F=MvrVW
确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure BDA0002656948900000064
Figure BDA0002656948900000065
将虚拟物体按照F位姿放置虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
本发明的有益效果为:解决了当前虚实注册方法中,注册精度与实时性的矛盾,解决基于标记的注册方法的相对位置约束与工业实际约束、以及基于多姿态模板匹配注册方法的计算量大,实时性差的约束,能够在不对真实物体张贴标记的情况下,就能够实现投影标记和视觉的虚实注册,并且精度高,实施起来更加方便实用。
附图说明
图1是本发明的方法流程图。
具体实施方式
以下结合附图对本发明作进一步详细说明。
本发明提供一种基于CAD模型的投影标记和机器视觉的实时注册方法,包括以下步骤:
采用虚拟相机对真实尺寸的CAD虚拟模型进行拍摄,得到虚拟模型图像;采用真实相机对真实物体进行拍摄,得到真实物体图像,并根据虚拟模型图像和真实物体图像的相似度计算得出两者之间的对应关系;
在得出虚拟模型图像和真实物体图像之间的对应关系后,对真实相机进行标定,通过虚拟模型图像所在像素坐标系和真实物体图像所在世界坐标系的转换关系,确定虚拟模型和真实物体坐标的转换关系;
在虚拟模型四周选定投影标记点,通过虚拟模型与真实模型之间的坐标转换关系,计算得到投影标记点在真实模型四周的坐标,通过像素坐标系与世界坐标系的转换关系,将投影标记点的世界坐标转化为真实相机的像素坐标;通过相机投影仪标定方法,得到相机像素坐标到投影仪坐标的转换关系,并由此获得投影仪坐标系下的投影坐标,再按投影坐标将标记投影于真实物体四周;
计算出自由相机和投影标记点位置的相对位置关系,从投影标记点位姿得到自由相机位姿;根据自由相机位姿,在虚拟模型所在的虚拟空间坐标系中设置相同位姿的虚拟相机;设置想要叠加于真实物体的虚拟模型相对于真实模型的坐标,计算虚拟模型与投影标记点之间的相对位姿,将虚拟模型按照此相对位姿放置在虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
具体的,虚拟模型图像和真实物体图像的获得方法是:将真实尺寸的CAD模型放置在虚拟空间中,并建立虚拟空间坐标系,通过虚拟相机对CAD模型进行拍摄得到虚拟模型图像;将真实相机按照虚拟相机和虚拟模型之间的相对位置对真实物体进行拍摄得到真实物体图像。
具体的,虚拟空间坐标系的建立方法为:将虚拟模型底部中心作为虚拟空间坐标原点,以虚拟模型底面为X-Y平面,垂直于模型底面为Z轴。
具体的,虚拟模型和真实物体之间对应关系的计算方法为:
A1、将视角向下的虚拟相机设置在虚拟模型的中心上方,使得虚拟相机中心与虚拟模型中心一致;
A2、在X-Y平面上,在[0,2π)范围内,均匀划分为m个角度,每个角度为α,在每一角度上,通过虚拟相机拍摄虚拟模型,获得每一角度姿态下的模型顶层图像模板;设图像模板集合为UA,Um为各角度的图像模板,有:
UA={Um};
所有图像模板对应旋转角为0的基准模板U0的旋转矩阵为:
Figure BDA0002656948900000081
A3、将每张图像模板通过拉普拉斯算子进行处理,取拉普拉斯过零点作为边缘,使得图像模板转换为边缘二值图像;
A4、将真实相机设于真实物体顶部,平行与底面向下实时视频采集真实物体图像,设置真实相机与虚拟相机高度、分辨率、视野大小完全一致,使得真实相机拍摄真实物体的尺寸与虚拟相机拍摄虚拟模型的尺寸完全一致,并且使世界坐标原点为真实相机光轴与桌面交点,此时真实世界坐标与虚拟世界坐标完全对应;
通过A3中的相同方法,对采集的每一帧真实物体图像进行边缘二值处理;
A5、将A4中实时采集的每一帧真实物体图像的边缘二值图像与虚拟模型的图像模板的边缘二值图像进行倒角距离匹配;设Um={ui}与V={vj}设定为模板与真实尺寸的边缘图上的各自边缘点的集合;在Um与V之间的倒角距离由每个ui∈Um的点与V中边缘的最近距离的平均值来决定,其中n为U集合中边缘点的数量:
Figure BDA0002656948900000091
A6、在A5的基础上,对模板U进行一定范围内的旋转平移,迭代缩小倒角距离范围;将U0模板放置在采集图像的中心位置与采集的真实物体图像进行匹配;设W函数为与参数s相关的2D的欧几里得变换函数,其中s∈SE(2),s=(θ,tx,ty),这里tx,ty分别表示为沿着x,y轴方向的平移变换t,θ为平面内的旋转变换;图像上的点x通过W函数执行变换操作:
Figure BDA0002656948900000092
m模板的最佳匹配变换系数sm,通过下式计算获得,得到当前模板与待识别图像之间的最佳旋转平移关系:
Figure BDA0002656948900000093
这里W(Um;s)={W(ui,s)};
A7、根据A6倒角匹配计算公式,遍历各角度的模型图像模板Um∈UA,获取最佳匹配模板
Figure BDA0002656948900000094
与对应的变换系数
Figure BDA0002656948900000095
满足公式如下:
Figure BDA0002656948900000096
Figure BDA0002656948900000097
A8、确定模板坐标与最佳匹配之间的变换矩阵,根据相对于基准模板U0旋转变换矩阵
Figure BDA0002656948900000098
并与
Figure BDA0002656948900000099
与输入图像V最佳匹配的旋转平移变换矩阵
Figure BDA00026569489000000910
中的旋转变换
Figure BDA00026569489000000911
Figure BDA00026569489000000912
求解最佳匹配状态与基准模板U0的旋转变换矩阵
Figure BDA00026569489000000913
Figure BDA00026569489000000914
Figure BDA00026569489000000915
通过变换矩阵
Figure BDA00026569489000000916
实现将虚拟模型模板与真实模型实时图像变换一致,采集图像上的真实模型坐标相对于放置在图像中心的U0模板的像素坐标V为:
Figure BDA00026569489000000917
A9、对真实相机进行标定,选定张正友标定法,计算得到世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,设U0模板像素坐标对应的世界坐标为UW,V真实物体像素坐标对应的世界坐标VW,有公式:
U0=MWCUW
V=MWCVW
且根据步骤A8中像素坐标U0与V之间的转换关系,所以VW与UW之间的转换关系为:
Figure BDA0002656948900000101
Figure BDA0002656948900000102
此时虚拟模型通过
Figure BDA0002656948900000103
转换矩阵,与真实物体的位姿完全一致。
具体的,确定投影仪投影标记点的方法为:
B1、在相对于虚拟模型的+x,-x,+y,-y方向上各投影一个不同的人工标记,设偏移量为O,则各标记中心与物体中心的偏移矩阵为:
Figure BDA0002656948900000104
Figure BDA0002656948900000105
设i为投影标记序号,则标记偏移矩阵集合
Figure BDA0002656948900000106
可以表示为:
Figure BDA0002656948900000107
标记坐标Pm相对于虚拟模型UW的坐标为:
Figure BDA0002656948900000108
则标记坐标Pm的像素坐标Pmc,根据步骤A9计算的世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,所以标记对于虚拟模型的像素坐标Pm为:
Pmc=MWCPm;并且,根据步骤A8中,模板像素坐标与采集图像坐标的转换关系为
Figure BDA0002656948900000109
所以标记在真实物体周围的像素坐标Pc可以根据下式计算得到:
Figure BDA00026569489000001010
所以真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW为:
Figure BDA00026569489000001011
B2、通过相机投影仪标定算法,计算得到相机坐标系与投影仪坐标系转换矩阵Mcp,以及投影仪坐标与世界坐标的转换矩阵Mpw,通过以下式子,即可得出投影仪坐标下的标记位置Pp与标记真实位置坐标Pw
Pp=McpPc
Pw=McpMpwPc
投影仪根据Pp坐标将标记投影到真实物体四周;
B3,确定投影标记真实投影位置Pw与真实物体位置VW间的转换关系,根据步骤B1中说明的真实标记像素坐标与虚拟标记像素坐标之间的转换关系,以及B2中,有如下关系:
Figure BDA0002656948900000111
并且根据A9中真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW的转换关系,带入上式子可知真实模型坐标VW与真实物体四周投影标记Pw之间的转换关系为:
Figure BDA0002656948900000112
具体的,得到自由相机拍摄下的真实物体投影图像的方法是:
C1、使用ORB特征构建每个标记的特征以及各个特征相对于自身标记的位置,建立标记特征模板特征点集合世界坐标集合
Figure BDA0002656948900000113
通过自由相机拍摄标记,将采集到的图片使用ORB进行处理,得到相机拍摄的标记特征点像素坐标集合
Figure BDA0002656948900000114
将每个标记模板特征与采集图像中的采集到的特征,通过RANSAC随机采样一致性算法,将对应位置的特征点匹配起来;
C2、根据投影标记确定相机位姿,首先通过张正友标定法,计算出自由相机内参MC,设投影的真实标记Pw与自由相机坐标系C之间的变换矩阵为矩阵Mcw,求解公式如下:
Pfc=MCMcwPw
根据Pfc与Pw的特征点之间对应关系,可列n个等式方程,n=特征点数量且满足≥4,即可求解Mcw
C3、计算自由相机C相对于投影标记Pw的位姿,根据C2,可知:
C=McwPw
将虚拟相机按照C位姿放置在虚拟空间中;
C4、设定要叠加在真实模型上的虚拟模型位姿与真实模型位姿VW之间的转换关系,设转换矩阵为Mvr,转换关系如下:
F=MvrVW
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置V间的转换关系,与C4虚拟模型坐标F与真实模型V之间的位置关系,确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure BDA0002656948900000121
Figure BDA0002656948900000122
将虚拟物体按照F位姿放置虚拟空间中;
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置VW之间的转换关系:
Figure BDA0002656948900000123
和C4中虚拟模型坐标F与真实模型VW之间的位置关系:
F=MvrVW
确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure BDA0002656948900000124
Figure BDA0002656948900000125
将虚拟物体按照F位姿放置虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:包括以下步骤:
采用虚拟相机对真实尺寸的CAD虚拟模型进行拍摄,得到虚拟模型图像;采用真实相机对真实物体进行拍摄,得到真实物体图像,并根据虚拟模型图像和真实物体图像的相似度计算得出两者之间的对应关系;
在得出虚拟模型图像和真实物体图像之间的对应关系后,对真实相机进行标定,通过虚拟模型图像所在像素坐标系和真实物体图像所在世界坐标系的转换关系,确定虚拟模型和真实物体坐标的转换关系;
在虚拟模型四周选定投影标记点,通过虚拟模型与真实模型之间的坐标转换关系,计算得到投影标记点在真实模型四周的坐标,通过像素坐标系与世界坐标系的转换关系,将投影标记点的世界坐标转化为真实相机的像素坐标;通过相机投影仪标定方法,得到相机像素坐标到投影仪坐标的转换关系,并由此获得投影仪坐标系下的投影坐标,再按投影坐标将标记投影于真实物体四周;
计算出自由相机和投影标记点位置的相对位置关系,从投影标记点位姿得到自由相机位姿;根据自由相机位姿,在虚拟模型所在的虚拟空间坐标系中设置相同位姿的虚拟相机;设置想要叠加于真实物体的虚拟模型相对于真实模型的坐标,计算虚拟模型与投影标记点之间的相对位姿,将虚拟模型按照此相对位姿放置在虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
2.根据权利要求1所述的基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:虚拟模型图像和真实物体图像的获得方法是:将真实尺寸的CAD模型放置在虚拟空间中,并建立虚拟空间坐标系,通过虚拟相机对CAD模型进行拍摄得到虚拟模型图像;将真实相机按照虚拟相机和虚拟模型之间的相对位置对真实物体进行拍摄得到真实物体图像。
3.根据权利要求2所述的基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:虚拟空间坐标系的建立方法为:将虚拟模型底部中心作为虚拟空间坐标原点,以虚拟模型底面为X-Y平面,垂直于模型底面为Z轴。
4.根据权利要求3所述的基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:虚拟模型和真实物体之间对应关系的计算方法为:
A1、将视角向下的虚拟相机设置在虚拟模型的中心上方,使得虚拟相机中心与虚拟模型中心一致;
A2、在X-Y平面上,在[0,2π)范围内,均匀划分为m个角度,每个角度为α,在每一角度上,通过虚拟相机拍摄虚拟模型,获得每一角度姿态下的模型顶层图像模板;设图像模板集合为UA,Um为各角度的图像模板,有:
UA={Um};
所有图像模板对应旋转角为0的基准模板U0的旋转矩阵为:
Figure FDA0002656948890000021
A3、将每张图像模板通过拉普拉斯算子进行处理,取拉普拉斯过零点作为边缘,使得图像模板转换为边缘二值图像;
A4、将真实相机设于真实物体顶部,平行与底面向下实时视频采集真实物体图像,设置真实相机与虚拟相机高度、分辨率、视野大小完全一致,使得真实相机拍摄真实物体的尺寸与虚拟相机拍摄虚拟模型的尺寸完全一致,并且使世界坐标原点为真实相机光轴与桌面交点,此时真实世界坐标与虚拟世界坐标完全对应;
通过A3中的相同方法,对采集的每一帧真实物体图像进行边缘二值处理;
A5、将A4中实时采集的每一帧真实物体图像的边缘二值图像与虚拟模型的图像模板的边缘二值图像进行倒角距离匹配;设Um={ui}与V={vj}设定为模板与真实尺寸的边缘图上的各自边缘点的集合;在Um与V之间的倒角距离由每个ui∈Um的点与V中边缘的最近距离的平均值来决定,其中n为U集合中边缘点的数量:
Figure FDA0002656948890000022
A6、在A5的基础上,对模板U进行一定范围内的旋转平移,迭代缩小倒角距离范围;将U0模板放置在采集图像的中心位置与采集的真实物体图像进行匹配;设W函数为与参数s相关的2D的欧几里得变换函数,其中s∈SE(2),s=(θ,tx,ty),这里tx,ty分别表示为沿着x,y轴方向的平移变换t,θ为平面内的旋转变换;图像上的点x通过W函数执行变换操作:
Figure FDA0002656948890000031
m模板的最佳匹配变换系数sm,通过下式计算获得,得到当前模板与待识别图像之间的最佳旋转平移关系:
Figure FDA0002656948890000032
这里W(Um;s)={W(ui,s)};
A7、根据A6倒角匹配计算公式,遍历各角度的模型图像模板Um∈UA,获取最佳匹配模板
Figure FDA0002656948890000033
与对应的变换系数
Figure FDA0002656948890000034
满足公式如下:
Figure FDA0002656948890000035
Figure FDA0002656948890000036
A8、确定模板坐标与最佳匹配之间的变换矩阵,根据相对于基准模板U0旋转变换矩阵
Figure FDA0002656948890000037
并与
Figure FDA0002656948890000038
与输入图像V最佳匹配的旋转平移变换矩阵
Figure FDA0002656948890000039
中的旋转变换
Figure FDA00026569488900000310
Figure FDA00026569488900000311
求解最佳匹配状态与基准模板U0的旋转变换矩阵
Figure FDA00026569488900000312
Figure FDA00026569488900000313
Figure FDA00026569488900000314
通过变换矩阵
Figure FDA00026569488900000315
实现将虚拟模型模板与真实模型实时图像变换一致,采集图像上的真实模型坐标相对于放置在图像中心的U0模板的像素坐标V为:
Figure FDA00026569488900000316
A9、对真实相机进行标定,选定张正友标定法,计算得到世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,设U0模板像素坐标对应的世界坐标为UW,V真实物体像素坐标对应的世界坐标VW,有公式:
U0=MWCUW
V=MWCVW
且根据步骤A8中像素坐标U0与V之间的转换关系,所以VW与UW之间的转换关系为:
Figure FDA00026569488900000317
Figure FDA0002656948890000041
此时虚拟模型通过
Figure FDA0002656948890000042
转换矩阵,与真实物体的位姿完全一致。
5.根据权利要求4所述的基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:确定投影仪投影标记点的方法为:
B1、在相对于虚拟模型的+x,-x,+y,-y方向上各投影一个不同的人工标记,设偏移量为O,则各标记中心与物体中心的偏移矩阵为:
Figure FDA0002656948890000043
Figure FDA0002656948890000044
设i为投影标记序号,则标记偏移矩阵集合
Figure FDA0002656948890000045
可以表示为:
Figure FDA0002656948890000046
标记坐标Pm相对于虚拟模型UW的坐标为:
Figure FDA0002656948890000047
则标记坐标Pm的像素坐标Pmc,根据步骤A9计算的世界坐标系到真实相机像素坐标系之间的转换矩阵MWC,所以标记对于虚拟模型的像素坐标Pm为:Pmc=MWCPm;并且,根据步骤A8中,模板像素坐标与采集图像坐标的转换关系为
Figure FDA0002656948890000048
所以标记在真实物体周围的像素坐标Pc可以根据下式计算得到:
Figure FDA0002656948890000049
所以真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW为:
Figure FDA00026569488900000410
B2、通过相机投影仪标定算法,计算得到相机坐标系与投影仪坐标系转换矩阵Mcp,以及投影仪坐标与世界坐标的转换矩阵Mpw,通过以下式子,即可得出投影仪坐标下的标记位置Pp与标记真实位置坐标Pw
Pp=McpPc
Pw=McpMpwPc
投影仪根据Pp坐标将标记投影到真实物体四周;
B3,确定投影标记真实投影位置Pw与真实物体位置VW间的转换关系,根据步骤B1中说明的真实标记像素坐标与虚拟标记像素坐标之间的转换关系,以及B2中,有如下关系:
Figure FDA0002656948890000051
并且根据A9中真实物体四周的标记像素坐标Pc相对于虚拟模型的坐标UW的转换关系,带入上式子可知真实模型坐标VW与真实物体四周投影标记Pw之间的转换关系为:
Figure FDA0002656948890000052
6.根据权利要求5所述的基于CAD模型的投影标记和机器视觉的实时注册方法,其特征是:得到自由相机拍摄下的真实物体投影图像的方法是:
C1、使用ORB特征构建每个标记的特征以及各个特征相对于自身标记的位置,建立标记特征模板特征点集合世界坐标集合
Figure FDA0002656948890000053
通过自由相机拍摄标记,将采集到的图片使用ORB进行处理,得到相机拍摄的标记特征点像素坐标集合
Figure FDA0002656948890000054
将每个标记模板特征与采集图像中的采集到的特征,通过RANSAC随机采样一致性算法,将对应位置的特征点匹配起来;
C2、根据投影标记确定相机位姿,首先通过张正友标定法,计算出自由相机内参MC,设投影的真实标记Pw与自由相机坐标系C之间的变换矩阵为矩阵Mcw,求解公式如下:
Pfc=MCMcwPw
根据Pfc与Pw的特征点之间对应关系,可列n个等式方程,n=特征点数量且满足≥4,即可求解Mcw
C3、计算自由相机C相对于投影标记Pw的位姿,根据C2,可知:
C=McwPw
将虚拟相机按照C位姿放置在虚拟空间中;
C4、设定要叠加在真实模型上的虚拟模型位姿与真实模型位姿VW之间的转换关系,设转换矩阵为Mvr,转换关系如下:
F=MvrVW
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置V间的转换关系,与C4虚拟模型坐标F与真实模型V之间的位置关系,确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure FDA0002656948890000061
Figure FDA0002656948890000062
将虚拟物体按照F位姿放置虚拟空间中;
C5、以投影标记位姿Pw,作为虚拟坐标系基准,根据B3中得到的投影标记真实投影位置Pw与真实物体位置VW之间的转换关系:
Figure FDA0002656948890000063
和C4中虚拟模型坐标F与真实模型VW之间的位置关系:
F=MvrVW
确定虚拟模型位姿F与标记真实投影位置Pw的相对位置,转换关系如下:
Figure FDA0002656948890000064
Figure FDA0002656948890000065
将虚拟物体按照F位姿放置虚拟空间中;利用虚拟相机拍照虚拟模型,获得此位姿下虚拟模型图像;将虚拟模型图像叠加于自由相机拍摄的图像上,完成虚拟模型与真实物体的虚实注册。
CN202010890907.5A 2020-08-29 2020-08-29 基于cad模型的投影标记和机器视觉的实时注册方法 Pending CN112017302A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010890907.5A CN112017302A (zh) 2020-08-29 2020-08-29 基于cad模型的投影标记和机器视觉的实时注册方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010890907.5A CN112017302A (zh) 2020-08-29 2020-08-29 基于cad模型的投影标记和机器视觉的实时注册方法

Publications (1)

Publication Number Publication Date
CN112017302A true CN112017302A (zh) 2020-12-01

Family

ID=73503242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010890907.5A Pending CN112017302A (zh) 2020-08-29 2020-08-29 基于cad模型的投影标记和机器视觉的实时注册方法

Country Status (1)

Country Link
CN (1) CN112017302A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113538578A (zh) * 2021-06-22 2021-10-22 恒睿(重庆)人工智能技术研究院有限公司 目标定位方法、装置、计算机设备和存储介质
CN113610986A (zh) * 2021-07-06 2021-11-05 同济大学 一种基于生物不变特征的数字物理耦合方法和系统
CN114048541A (zh) * 2022-01-12 2022-02-15 深圳前海贾维斯数据咨询有限公司 基于数字孪生的资产空间标记方法及系统
CN114820808A (zh) * 2022-03-31 2022-07-29 上海理工大学 基于光学追踪器和ar眼镜的虚实融合标定方法、系统
CN116541923A (zh) * 2023-04-07 2023-08-04 中国民用航空飞行学院 一种基于vr的带支架设备室内安装基础定位方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113538578A (zh) * 2021-06-22 2021-10-22 恒睿(重庆)人工智能技术研究院有限公司 目标定位方法、装置、计算机设备和存储介质
CN113538578B (zh) * 2021-06-22 2023-07-25 恒睿(重庆)人工智能技术研究院有限公司 目标定位方法、装置、计算机设备和存储介质
CN113610986A (zh) * 2021-07-06 2021-11-05 同济大学 一种基于生物不变特征的数字物理耦合方法和系统
CN113610986B (zh) * 2021-07-06 2024-04-23 同济大学 一种基于生物不变特征的数字物理耦合方法和系统
CN114048541A (zh) * 2022-01-12 2022-02-15 深圳前海贾维斯数据咨询有限公司 基于数字孪生的资产空间标记方法及系统
CN114048541B (zh) * 2022-01-12 2022-04-22 深圳前海贾维斯数据咨询有限公司 基于数字孪生的资产空间标记方法及系统
CN114820808A (zh) * 2022-03-31 2022-07-29 上海理工大学 基于光学追踪器和ar眼镜的虚实融合标定方法、系统
CN116541923A (zh) * 2023-04-07 2023-08-04 中国民用航空飞行学院 一种基于vr的带支架设备室内安装基础定位方法
CN116541923B (zh) * 2023-04-07 2023-12-19 中国民用航空飞行学院 一种基于vr的带支架设备室内安装基础定位方法

Similar Documents

Publication Publication Date Title
CN112017302A (zh) 基于cad模型的投影标记和机器视觉的实时注册方法
CN106553195B (zh) 工业机器人抓取过程中的物体6自由度定位方法及系统
CN108648240B (zh) 基于点云特征地图配准的无重叠视场相机姿态标定方法
WO2014024579A1 (ja) 光学データ処理装置、光学データ処理システム、光学データ処理方法、および光学データ処理用プログラム
CN111121655B (zh) 一种共面等大多孔型工件位姿与孔径视觉检测方法
CN107767456A (zh) 一种基于rgb‑d相机的物体三维重建方法
CN108288294A (zh) 一种3d相机群的外参标定方法
CN110717936B (zh) 一种基于相机姿态估计的图像拼接方法
CN113034593B (zh) 6d位姿标注方法、系统及存储介质
Pathak et al. Dense 3D reconstruction from two spherical images via optical flow-based equirectangular epipolar rectification
CN109613974B (zh) 一种大场景下的ar家居体验方法
CN113112545B (zh) 基于计算机视觉的手持移动打印装置定位方法
KR20240089161A (ko) 촬영 측정 방법, 장치, 기기 및 저장 매체
CN115564842A (zh) 双目鱼眼相机的参数标定方法、装置、设备及存储介质
CN109785373A (zh) 一种基于散斑的六自由度位姿估计系统及方法
CN114022542A (zh) 一种基于三维重建的3d数据库制作方法
WO2022156044A1 (zh) 一种即时获得线上物体实际尺寸测量方法
CN116468609A (zh) 基于SuperGlue的两段式变焦相机多图像拼接方法和系统
CN110363801B (zh) 工件实物与工件三维cad模型的对应点匹配方法
CN115423863A (zh) 相机位姿估计方法、装置及计算机可读存储介质
WO2020133080A1 (zh) 物体定位方法、装置、计算机设备及存储介质
Pathak et al. Virtual reality with motion parallax by dense optical flow-based depth generation from two spherical images
KR101673144B1 (ko) 부분 선형화 기반의 3차원 영상 정합 방법
Ventura et al. Structure and motion in urban environments using upright panoramas
EP2879090B1 (en) Aligning ground based images and aerial imagery

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20201201

WD01 Invention patent application deemed withdrawn after publication