CN102721375A - 一种强反光金属结构件的在位测量中多次反光抑制方法 - Google Patents
一种强反光金属结构件的在位测量中多次反光抑制方法 Download PDFInfo
- Publication number
- CN102721375A CN102721375A CN2012102107428A CN201210210742A CN102721375A CN 102721375 A CN102721375 A CN 102721375A CN 2012102107428 A CN2012102107428 A CN 2012102107428A CN 201210210742 A CN201210210742 A CN 201210210742A CN 102721375 A CN102721375 A CN 102721375A
- Authority
- CN
- China
- Prior art keywords
- phase
- coordinate
- projection
- camera
- projector
- 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.)
- Granted
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种强反光金属结构件的在位测量中多次反光抑制方法,该方法能够实现对强反光金属结构件的在位测量,它有八大步骤。该方法首先对左右相机和投射器进行标定;再利用在位测量机床提供的参数求取投射器与工件坐标系的相对位姿;然后通过零件CAD的模型的各个面法线方向进行测量规划;进而工件上需要测量的点投影到数字投射器上,确定出需要投射的区域后选择性投射;最后通过对相机拍的图片的相位解算完成物体的三维复现。本发明除了具有速度快、精度高和三维点云稠密等优点外,还实现了无人工干预下的全自动在位测量,可操作性强,能够满足现场强反光金属结构件在位快速测量的需求。
Description
技术领域
本发明涉及一种强反光金属结构件的在位测量中多次反光抑制方法,本方法可用于对强反光金属结构件的在位测量。本发明属于机器视觉技术领域。
背景技术
大型金属结构件在航空、航天、船舶、汽车等高端装备制造领域应用十分广泛,其加工精度对装备的性能具有重要作用。传统金属整体结构件的材料主要为铝合金或钛合金,经过精密铣削加工之后结构件表面极为光滑,会出现强反光现象,使现有商用光学三维视觉测量仪器失效,无法正确测量金属等强反光表面的三维形貌。尤其在基于自适应快速条纹投射的在线三维视觉检测方法中,金属结构件表面的强反光导致相机采集到的条纹图像饱和或过暗,无法正确解相,出现大面积三维点云缺失,而且即便是测量出的三维点云,其测量精度也较低。解决问题的关键在于提高测量系统的动态范围,使其与金属表面反射率范围一致。而且,现场测量环境下,环境光干扰不易控制,测量系统必须能够适应现场环境光的要求,大幅提高了强反光金属结构件在位测量的难度。
因此,在强反光金属结构件在位测量中,金属表面强放光问题和多次反光问题是亟待解决的关键瓶颈,测量基础理论和关键技术亟待突破。研究强反光表面三维形貌测量理论和方法,建立多次反光条件下新的测量理论,突破强反光金属结构件在位测量中的技术瓶颈,建立完整的强反光金属结构件在位测量中多次反光抑制理论和方法体系具有重要的应用和推广价值。
发明内容
本发明的技术解决问题是:提出一种强反光金属结构件的在位测量中多次反光抑制方法,将选择性投射多次反光抑制技术、强反光表面的高动态范围条纹图像采集技术与高速、高动态范围和选择性投射的数字条纹投射器相结合,通过选择性投射多次反光抑制技术实现对被测表面条纹的选择性投射,通过强反光表面的高动态范围条纹图像采集技术实现提高合成高动态范围条纹图像,通过高速、高动态范围和选择性投射的数字条纹投射器实现投射的条纹具有较大的亮度变化范围、较高的响应速度和刷新频率,从而实现强反光金属结构件的在位测量。
本发明的技术解决方案为:一种强反光金属结构件的在位测量中多次反光抑制方法,其特征在于,它包括以下步骤:
步骤一:利用基于平面靶标的相机标定算法,标定左相机、右相机和投射器内参数,标定出左相机与右相机之间的外参数和左相机与投射器之间的外参数。
首先将平面靶标至于投射器和相机前,通过相机拍摄靶标计算圆心点坐标;然后投射横纵条纹,计算圆心点坐标对应的相位,并由相位计算对应的数字投射器坐标,从而建立其靶标圆心三维坐标与投射器像平面上二维图像坐标的关系;最后利用基于平面靶标的相机标定算法,标定投影仪内参数。
利用靶标作为中介标定出投影仪与左相机之间的外参数和左相机与右相机之间的外参数。
步骤二:通过零件的CAD模型进行测量规划。
首先根据高斯球原理,确定测量传感器的测量姿态(方位角);然后根据测量传感器单视场的测量范围,确定每个姿态下测量传感器的空间坐标。再根据坐标转换原理,将传感器的空间坐标转换为机床下的坐标值。最后,就划分出来了需要选择性投射的各个面。
步骤三:计算被测工件和数字投射器之间的坐标关系。
即需要确定旋转矩阵R和平移向量T。数字投射器与左相机坐标系之间的相对位姿可由数字投射器标定给出,左相机坐标系与拼接测量坐标系之间的相对位姿可由拼接测量系统给出,拼接测量坐标系与工件坐标系之间的相对位姿可由坐标对齐系统给出,由坐标系转换关系可得到数字投射器与工件坐标系之间的相对位姿。
步骤四:将工件上需要测量的点投影到数字投射器上,确定出需要投射的区域。
将工件上需要测量的点的坐标转换成数字投射器坐标,由此确定出需要投射的区域,为选择性投射做准备。
步骤五:投射的区域进行腐蚀及膨胀运算后,选择性投射。
根据步骤二对被测零件的测量规划,根据投射面和现场环境的需要选择性投射不同亮度的条纹图像。
步骤六:利用左右相机拍摄投射需要测的各个面的条纹图像。
根据需要调节相机曝光时间拍摄待测面的图像。
步骤七:相位解算与正确性分析。对拍摄的物体横、竖正弦条纹的投影,进行相位解算和相位展开,获得大量匹配点。
相位解算采用四步相移法。对应的解相公式为
其中gi(x,y)为条纹图像,Φ(x,y)为需要解算的相主值,i为相移次数。
利用外差多频相位展开法进行相位展开。
通过投射多级周期相差不大的条纹,根据外差原理,得到一个周期可覆盖整个视场的条纹,利用这个条纹与其中一个原始投射条纹在零相位相同的情况下相位值与周期的反比关系,可进行相位展开。根据实际需要,这里投射三种不同频率的条纹。
经过对物体的横竖投影、拍照,相位解算和展开之后,左右相机的每一个像素对应着一对横竖相位值。这些相位值以二维数组的形式存储,横、竖相位分别存在不同的数组中。
步骤八:利用外差多频相位展开法得到的最后的解相结果进行被测表面的三维复现,得到被测表面的三维点坐标。
本发明与现有技术相比的优点在于:
(1)本发明基于最佳条纹调制度的高动态条纹投射与采集参数自动优化选择方法。在分析投射栅相位法的测量机理、强反光金属结构件的反光特点和现场在位测量环境的基础上,建立强反光表面测量模型,以最佳条纹调制度为准则,分析高动态条纹投射与采集参数优化理论和自动选择方法,提高强反光表面三维形貌在位测量的可行性和实用性。
(2)本发明基于选择性投射的多次反光抑制理论和方法。根据在位测量中的路径规划结果和金属结构件设计模型,自动计算选择性投射的区域,保证多次反光抑制理论和方法的正确性和实用性。
(3)本发明针对强反光金属表面在位测量需求,研制可实现选择性投射的高速、高动态范围数字条纹投射器等专用测量设备,从硬件上满足强反光金属表面在位测量需求,有效推进强反光金属表面在位测量理论和方法的实际应用,从根本上解决该测量难题。
附图说明
图1为本发明方法的流程图。
具体实施方式
为了更好的理解本发明,下面结合附图和实施例对本发明的技术方案作详细的描述。
参照附图1,一种强反光金属结构件的在位测量中多次反光抑制方法,包含如下步骤:
在对强反光金属结构件的在位测量之前,测量系统需要进行测量前的准备工作,包括调试测量软件、投射器等;
步骤一:利用基于平面靶标的相机标定算法,标定左相机、右相机和投射器内参数,标定出左相机与右相机之间的外参数和左相机与投射器之间的外参数。
首先将平面靶标至于投射器和相机前,通过相机拍摄靶标计算圆心点坐标;然后投射横纵条纹,计算圆心点坐标对应的相位,并由相位计算对应的数字投射器坐标,从而建立其靶标圆心三维坐标与投射器像平面上二维图像坐标的关系;最后利用基于平面靶标的相机标定算法,标定投影仪内参数。
利用靶标作为中介标定出投影仪与左相机之间的外参数和左相机与右相机之间的外参数。
步骤二:通过零件的CAD模型进行测量规划。
首先根据高斯球原理,确定测量传感器的测量姿态(方位角);然后根据测量传感器单视场的测量范围,确定每个姿态下测量传感器的空间坐标。再根据坐标转换原理,将传感器的空间坐标转换为机床下的坐标值。最后,就划分出来了需要选择性投射的各个面。
步骤三:计算被测工件和数字投射器之间的坐标关系。
即需要确定旋转矩阵R和平移向量T。数字投射器与左相机坐标系之间的相对位姿可由数字投射器标定给出,左相机坐标系与拼接测量坐标系之间的相对位姿可由拼接测量系统给出,拼接测量坐标系与工件坐标系之间的相对位姿可由坐标对齐系统给出,由坐标系转换关系可得到数字投射器与工件坐标系之间的相对位姿。
步骤四:将工件上需要测量的点投影到数字投射器上,确定出需要投射的区域。
将工件上需要测量的点的坐标转换成数字投射器坐标,由此确定出需要投射的区域,为选择性投射做准备。
步骤五:投射的区域进行腐蚀及膨胀运算后,选择性投射。
根据步骤二对被测零件的测量规划,根据投射面和现场环境的需要选择性投射不同亮度的条纹图像。
步骤六:利用左右相机拍摄投射需要测的各个面的条纹图像。
根据需要调节相机曝光时间拍摄待测面的图像。
步骤七:相位解算与正确性分析。
对拍摄的物体横、竖正弦条纹的投影,进行相位解算和相位展开,获得大量匹配点。
相位解算采用四步相移法。对应的解相公式为
其中gi(x,y)为条纹图像,Φ(x,y)为需要解算的相主值,i为相移次数。
利用外差多频相位展开法进行相位展开。
通过投射多级周期相差不大的条纹,根据外差原理,得到一个周期可覆盖整个视场的条纹,利用这个条纹与其中一个原始投射条纹在零相位相同的情况下相位值与周期的反比关系,可进行相位展开。根据实际需要,这里投射三种不同频率的条纹。
经过对物体的横竖投影、拍照,相位解算和展开之后,左右相机的每一个像素对应着一对横竖相位值。这些相位值以二维数组的形式存储,横、竖相位分别存在不同的数组中。
步骤八:利用外差多频相位展开法得到的最后的解相结果进行被测表面的三维复现,得到被测表面的三维点坐标。
Claims (1)
1.一种强反光金属结构件的在位测量中多次反光抑制方法,其特征在于,它包括以下步骤:
步骤一:利用基于平面靶标的相机标定算法,标定左相机、右相机和投射器内参数,标定出左相机与右相机之间的外参数和左相机与投射器之间的外参数;首先将平面靶标至于投射器和相机前,通过相机拍摄靶标计算圆心点坐标;然后投射横纵条纹,计算圆心点坐标对应的相位,并由相位计算对应的数字投射器坐标,从而建立其靶标圆心三维坐标与投射器像平面上二维图像坐标的关系;最后利用基于平面靶标的相机标定算法,标定投影仪内参数,利用靶标作为中介标定出投影仪与左相机之间的外参数和左相机与右相机之间的外参数;
步骤二:通过零件的CAD模型进行测量规划;首先根据高斯球原理,确定测量传感器的测量姿态即方位角;然后根据测量传感器单视场的测量范围,确定每个姿态下测量传感器的空间坐标;再根据坐标转换原理,将传感器的空间坐标转换为机床下的坐标值;最后,就划分出来了需要选择性投射的各个面;
步骤三:计算被测工件和数字投射器之间的坐标关系;即需要确定旋转矩阵R和平移向量T;数字投射器与左相机坐标系之间的相对位姿由数字投射器标定给出,左相机坐标系与拼接测量坐标系之间的相对位姿由拼接测量系统给出,拼接测量坐标系与工件坐标系之间的相对位姿由坐标对齐系统给出,由坐标系转换关系得到数字投射器与工件坐标系之间的相对位姿;
步骤四:将工件上需要测量的点投影到数字投射器上,确定出需要投射的区域;将工件上需要测量的点的坐标转换成数字投射器坐标,由此确定出需要投射的区域,为选择性投射做准备;
步骤五:投射的区域进行腐蚀及膨胀运算后,选择性投射;根据步骤二对被测零件的测量规划,根据投射面和现场环境的需要选择性投射不同亮度的条纹图像;
步骤六:利用左右相机拍摄投射需要测的各个面的条纹图像;根据需要调节相机曝光时间拍摄待测面的图像;
步骤七:相位解算与正确性分析;对拍摄的物体横、竖正弦条纹的投影,进行相位解算和相位展开,获得大量匹配点;
其中gi(x,y)为条纹图像,Φ(x,y)为需要解算的相主值,i为相移次数;
利用外差多频相位展开法进行相位展开,通过投射多级周期相差不大的条纹,根据外差原理,得到一个周期可覆盖整个视场的条纹,利用这个条纹与其中一个原始投射条纹在零相位相同的情况下相位值与周期的反比关系,进行相位展开;根据实际需要,这里投射三种不同频率的条纹,经过对物体的横竖投影、拍照,相位解算和展开之后,左右相机的每一个像素对应着一对横竖相位值,这些相位值以二维数组的形式存储,横、竖相位分别存在不同的数组中;
步骤八:利用外差多频相位展开法得到的最后的解相结果进行被测表面的三维复现,得到被测表面的三维点坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210210742.8A CN102721375B (zh) | 2012-06-20 | 2012-06-20 | 一种强反光金属结构件的在位测量中多次反光抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210210742.8A CN102721375B (zh) | 2012-06-20 | 2012-06-20 | 一种强反光金属结构件的在位测量中多次反光抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102721375A true CN102721375A (zh) | 2012-10-10 |
CN102721375B CN102721375B (zh) | 2015-02-11 |
Family
ID=46947201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210210742.8A Expired - Fee Related CN102721375B (zh) | 2012-06-20 | 2012-06-20 | 一种强反光金属结构件的在位测量中多次反光抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102721375B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102937425A (zh) * | 2012-10-18 | 2013-02-20 | 北京航空航天大学 | 一种基于高动态条纹投射器的强反射表面三维形貌测量系统 |
CN102967265A (zh) * | 2012-12-15 | 2013-03-13 | 吉林大学 | 基于机器视觉的汽车连杆裂解槽检测方法 |
CN105066962A (zh) * | 2015-07-21 | 2015-11-18 | 中国航空工业集团公司北京长城航空测控技术研究所 | 一种多分辨率大视场角高精度摄影测量装置 |
CN107607040A (zh) * | 2017-08-11 | 2018-01-19 | 天津大学 | 一种适用于强反射表面的三维扫描测量装置及方法 |
CN108507494A (zh) * | 2018-03-07 | 2018-09-07 | 北京航空航天大学 | 一种能够抑制间接光照的三维形貌测量方法 |
CN108519064A (zh) * | 2018-04-20 | 2018-09-11 | 天津工业大学 | 一种应用于多频率三维测量的反光抑制方法 |
CN111415405A (zh) * | 2020-03-16 | 2020-07-14 | 苏州依诺维视智能科技有限公司 | 高反光表面工件的三维高精度视觉测量方法 |
CN113432550A (zh) * | 2021-06-22 | 2021-09-24 | 北京航空航天大学 | 一种基于相位匹配的大尺寸零件三维测量拼接方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005063217A1 (de) * | 2005-12-22 | 2007-07-05 | Pilz Gmbh & Co. Kg | Verfahren zum Konfigurieren einer Überwachungseinrichtung zum Überwachen eines Raumbereichs |
CN102176211A (zh) * | 2010-12-26 | 2011-09-07 | 北京航空航天大学 | 一种基于投影栅相位法测量系统的路径规划方法 |
-
2012
- 2012-06-20 CN CN201210210742.8A patent/CN102721375B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005063217A1 (de) * | 2005-12-22 | 2007-07-05 | Pilz Gmbh & Co. Kg | Verfahren zum Konfigurieren einer Überwachungseinrichtung zum Überwachen eines Raumbereichs |
CN102176211A (zh) * | 2010-12-26 | 2011-09-07 | 北京航空航天大学 | 一种基于投影栅相位法测量系统的路径规划方法 |
Non-Patent Citations (4)
Title |
---|
HONGZHI JIANG: "High dynamic range fringe acquisition: A novel 3-D scanning technique for high-reflective surfaces", 《OPTICS AND LASERS IN ENGINEERING》, vol. 50, no. 10, 12 June 2012 (2012-06-12) * |
姜志宏: "用于强反射表面形貌测量的投影栅相位法", 《光学精密工程》, vol. 18, no. 9, 15 September 2010 (2010-09-15) * |
费键锋: "基于投影二维平面靶标的投影仪现场标定方法研究", 《传感器与微系统》, vol. 27, no. 4, 20 April 2008 (2008-04-20) * |
赵慧洁: "基于相位匹配的大视场视觉检测系统", 《北京航空航天大学学报》, vol. 32, no. 6, 30 June 2006 (2006-06-30) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102937425A (zh) * | 2012-10-18 | 2013-02-20 | 北京航空航天大学 | 一种基于高动态条纹投射器的强反射表面三维形貌测量系统 |
CN102937425B (zh) * | 2012-10-18 | 2015-04-29 | 北京航空航天大学 | 一种基于高动态条纹投射器的强反射表面三维形貌测量系统 |
CN102967265A (zh) * | 2012-12-15 | 2013-03-13 | 吉林大学 | 基于机器视觉的汽车连杆裂解槽检测方法 |
CN102967265B (zh) * | 2012-12-15 | 2015-10-21 | 吉林大学 | 基于机器视觉的汽车连杆裂解槽检测方法 |
CN105066962A (zh) * | 2015-07-21 | 2015-11-18 | 中国航空工业集团公司北京长城航空测控技术研究所 | 一种多分辨率大视场角高精度摄影测量装置 |
CN107607040A (zh) * | 2017-08-11 | 2018-01-19 | 天津大学 | 一种适用于强反射表面的三维扫描测量装置及方法 |
CN107607040B (zh) * | 2017-08-11 | 2020-01-14 | 天津大学 | 一种适用于强反射表面的三维扫描测量装置及方法 |
CN108507494A (zh) * | 2018-03-07 | 2018-09-07 | 北京航空航天大学 | 一种能够抑制间接光照的三维形貌测量方法 |
CN108519064A (zh) * | 2018-04-20 | 2018-09-11 | 天津工业大学 | 一种应用于多频率三维测量的反光抑制方法 |
CN111415405A (zh) * | 2020-03-16 | 2020-07-14 | 苏州依诺维视智能科技有限公司 | 高反光表面工件的三维高精度视觉测量方法 |
CN113432550A (zh) * | 2021-06-22 | 2021-09-24 | 北京航空航天大学 | 一种基于相位匹配的大尺寸零件三维测量拼接方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102721375B (zh) | 2015-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102721375B (zh) | 一种强反光金属结构件的在位测量中多次反光抑制方法 | |
CN102721376B (zh) | 一种大视场三维视觉传感器的标定方法 | |
CN105783775B (zh) | 一种镜面及类镜面物体表面形貌测量装置与方法 | |
CN102183214B (zh) | 一种大口径非球面镜结构光检测方法 | |
CN103335604B (zh) | 一种工作状态下风轮叶片全场三维变形在线监测方法 | |
Shang et al. | Measurement methods of 3D shape of large-scale complex surfaces based on computer vision: A review | |
CN103047944B (zh) | 三维物体测量方法及其测量装置 | |
CN105115560B (zh) | 一种船舱舱容的非接触测量方法 | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
Gai et al. | A flexible multi-view calibration and 3D measurement method based on digital fringe projection | |
CN111322954A (zh) | 一种装配工装位姿测量方法、装置、存储介质及电子设备 | |
JP4743771B2 (ja) | 断面データ取得方法、システム、及び断面検査方法 | |
Malesa et al. | Multi-camera DIC system with spatial data stitching procedure for measurements of engineering objects | |
Cai et al. | Calibration method for the rotating axis in panoramic 3D shape measurement based on a turntable | |
Xiaofei et al. | Three-dimensional reconstruction based on tri-frequency heterodyne principle | |
CN116242277A (zh) | 基于全场三维视觉的电源机箱结构件尺寸自动测量方法 | |
Zheng et al. | Positioning sensor by combining photogrammetry, optical projection and a virtual camera model | |
Jiang et al. | High-precision composite 3D shape measurement of aeroengine blade based on parallel single-pixel imaging and high-dynamic range N-step fringe projection profilometry | |
CN109425308A (zh) | 一种强反光抑制方法 | |
Zhang et al. | Photogrammetric 3D measurement method applying to automobile panel | |
Xu et al. | A geometry and optical property inspection system for automotive glass based on fringe patterns. | |
Ma et al. | Flexible High-Resolution Continuous 3-D Scanning for Large-Scale Industrial Components | |
Knyaz | Automated calibration technique for photogrammetric system based on a multi-media projector and a CCD camera | |
CN104765702A (zh) | 一种船舶外板水火加工过程中的型面数据采集方法 | |
KR100381816B1 (ko) | 삼면반사체를 이용한 6자유도 운동측정장치 |
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: 20150211 Termination date: 20200620 |
|
CF01 | Termination of patent right due to non-payment of annual fee |