CN104111038B - 利用相位融合算法修复饱和产生的相位误差的方法 - Google Patents
利用相位融合算法修复饱和产生的相位误差的方法 Download PDFInfo
- Publication number
- CN104111038B CN104111038B CN201410317007.6A CN201410317007A CN104111038B CN 104111038 B CN104111038 B CN 104111038B CN 201410317007 A CN201410317007 A CN 201410317007A CN 104111038 B CN104111038 B CN 104111038B
- Authority
- CN
- China
- Prior art keywords
- phase place
- modulate intensity
- normalization
- phase
- pixel
- 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
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明为利用相位融合算法修复饱和产生的相位误差的方法,涉及三维测量领域中如何修复饱和引起的相位误差。调节相机使像素点均不饱和,投影机投射纯色灰度图案到物体表面,并使用相机记录扫描后的图像;当相机发生上溢时,投影机投射多组结构光图案扫描物体,相邻两组图案的最大灰度值(幅值)增加,相机记录扫描后的图像;计算扫描后的结构光图像的调制强度及相位;结合捕获的纯色图像归一化调制强度;选取每个像素点的最大归一化调制强度,并将其对应的相位作为最终相位。本发明通过选择最大的归一化调制强度对应的相位进行融合,不受饱和程度的限制,当测量系统存在非线性时,也能有效修复饱和误差,扩大修复饱和误差的适用范围。
Description
技术领域
本发明涉及光学三维数据处理技术领域,具体涉及使用相位融合算法修复饱和产生的相位误差。
背景技术
相位测量轮廓术是一种主动式、非接触式的三维测量技术。它具有高速度、高精度等优点。该技术被广泛地用于产品检测、生物识别以及三维建模等领域。该技术的测量系统是由一台普通计算机、投影机和照相机组成。该技术通过投影机投射一组正弦光栅到被测物体表面,照相机同步地抓取图像数据,然后处理抓取的数据得到物体的三维形貌。但是,当获取的图像中部分像素点的强度值超出照相机的最大量化值,照相机将超出最大量化值的强度量化为其本身的最大量化值。这导致图片中部分像素点的强度值不准确,即光强度饱和。该效应在三维重建过程中对相位的计算有一定影响,降低了相位的精确度。文献(李艳,胡而已.投影条纹相移法中图像饱和误差抑制算法研究[J].光子学报,2010,39(6):1137~1142)提出了饱和系数这一概念,以此来描述照相机抓取的光栅图像中像素点饱和程度,同时也提出了一种相位修复算法,该算法仅使用未饱和像素点的强度值来计算相位。但是当图像中像素点的饱和程度超出饱和系数时,这种方法则无法修复饱和引起的相位误差;当不同的帧图像中的像素点发生饱和时,这种方法修复相位的公式也不相同。并且当测量系统存在非线性畸变时,这种方法也不能修复饱和误差。应用本发明提及的相位融合算法修复饱和引起的相位误差,不仅能获得准确的相位值,而且本发明的方法不受饱和系数的影响,即无论像素点的饱和程度是否超出饱和系数,本发明的方法均能修复。另外,当测量系统存在非线性畸变时,本发明的方法也能有效地修复饱和误差。
发明内容
本发明的目的是针对光强度饱和引起的相位误差问题,提出一种相位融合算法修复该相位误差。这种方法归一化像素点对应的调制强度,然后查找各个像素点的最大的归一化调制强度,并将各个像素点最大的归一化调制强度对应的相位融合。
为实现本发明的目的,本发明提供了以下技术方案:
首先通过投影机投射多组结构光图案到被测物体表面,该结构光图案相邻两组的最大灰度值(振幅)依次增加。其次根据相位测量轮廓术公式计算照相机捕获的各组图像的调制强度以及相位;然后对计算得到的调制强度进行归一化;最后查找各个像素点对应的最大的归一化调制强度,并将各个像素点最大的归一化调制强度对应的相位融合。
与现有技术相比,本发明有如下优点:
本发明采用相位融合算法修复饱和产生的相位误差,是将各个像素点最大的归一化调制强度对应的相位融合。本发明的方法不需要预先逐一判断每个像素点的饱和情况,不仅能获得准确的相位值,而且不受饱和系数的影响,即无论像素点的饱和程度是否超出饱和系数,本发明的方法均能修复。另外,当测量系统存在非线性畸变时,本发明的方法也能有效地修复饱和误差,扩大了修复饱和误差的适用范围。
附图说明
图1是本发明相位融合算法修复饱和误差的流程图。
图2是本发明的实验三维重建前视图和侧视图。图(a)是使用传统的相位测量轮廓术得到的三维重建前视图和侧视图;图(b)是使用本发明的方法得到的三维重建前视图和侧视图。
具体实施方式
下面结合附图和实施方式对本发明作进一步详细说明。
为了使发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
附图1是本发明相位融合算法修复饱和误差的流程图。采用的设备有一台CASIOXJ-A155V型投影机、Prosilica GC650C型照相机以及一台普通计算机组成。其中,投影机的分辨率是800×600,照相机的分辨率是640×480。
附图2是测量系统存在非线性畸变时,被扫面物体墙面的三维重建前视图和侧视图。为了减小噪声对本实验的影响,本实施例中求取25组数据的平均值。图中颜色的深浅表示墙表面的凹凸程度。图2(a)是使用传统的相位测量轮廓术得到的三维重建前视图和侧视图,从图中可以看出,墙表面颜色较深区域由于光强度饱和导致凹凸程度较大,图2(b)是使用本发明中的方法修复后得到的三维重建前视图和侧视图。该区域颜色较浅,凹凸程度明显减小。图2(b)中仍有凹凸的原因是:实验测量系统中存在非线性畸变,该因素也会引起重建墙面的凹凸变化。
本实例具体实施步骤如下:
1. 搭建三维测量系统,使照相机能够捕捉到标靶的全貌,其次,调节投影机的焦距以及照相机的光圈、曝光时间等参数,最后校准三维测量系统获得照相机和投影机的参数矩阵。
2.利用公式(1)生成需要投射的最大灰度值增大的多组正弦图案以及使用公式(2)生成纯色灰度图片。
(1)
(2)
其中,是第k组结构光图案的灰度值,是第k组纯色图案的灰度值,是投影仪的坐标,是本实施例中所述第k的各组最大灰度值。本实施例中第一组图案的最大灰度值是,即相邻两组的最大灰度值递增20,每一组图案的频率f是1和4,相移量为;n是相移指数,其范围从1到4,k是投射的第k组结构光图案,本实施例中投射图案的总组数为5。
3. 调节照相机的曝光时间为7000,像素点均不饱和。通过投影机投射纯色灰度图案到墙表面,并使用照相机抓取扫描后的图像数据。
4. 当相机发生上溢时,通过投影机投射步骤(1)中的图案到墙表面,并且照相机同步抓取图像。然后利用公式(3)和公式(4)分别计算各组数据对应的调制强度及相位。
(3)
(4)
5. 通过公式(5)对各组调制强度进行归一化处理,并且选取同一像素点对应的最大的归一化调制强度,提取该调制强度所在组的相位进行融合,如公式(6)和公式(7)所示。
(5)
(6)
(7)
6. 根据预校准得到的照相机及投影机参数矩阵,以及融合得到的相位,重建物体的三维形貌,计算方法如公式(8):
(8)
其中是物体的三维坐标,照相机及投影仪的参数矩阵分别为:
,
以及其他参数:
,
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种利用相位融合算法修复光强度饱和引起的相位误差的方法,其特征在于以下步骤:首先通过投影机将最大灰度值递增的多组结构光图案投射到物体表面,照相机同时捕捉对应图像,其次根据相位测量轮廓术公式计算出各组图像对应的调制强度及相位,然后分别归一化各组调制强度,然后融合各个像素点的最大的归一化调制强度所对应的相位,从多组归一化后的调制强度中,查找各个像素点的最大的归一化调制强度,并将各个像素点的最大归一化调制强度对应的相位作为每一个像素点的最终相位。
2.按照权利要求1所述的方法,其特征在于,所述最大灰度值递增的多组结构光图案,是指从第一组到最后一组的最大灰度值依次递增的结构光图案。
3.按照权利要求1所述的方法,其特征在于所述的归一化各组调制强度,是指首先根据相位测量轮廓术公式计算出各组数据的调制强度以及相位,且不需要预先判断每一个像素点的饱和程度;当测量系统不存在非线性畸变时,归一化处理调制强度是指将已经计算得到的各组调制强度除以各组对应的最大灰度值,当测量系统存在非线性畸变时,归一化处理调制强度是指将已经计算得到的各组调制强度除以最大灰度值经过非线性畸变后的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317007.6A CN104111038B (zh) | 2014-07-07 | 2014-07-07 | 利用相位融合算法修复饱和产生的相位误差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317007.6A CN104111038B (zh) | 2014-07-07 | 2014-07-07 | 利用相位融合算法修复饱和产生的相位误差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104111038A CN104111038A (zh) | 2014-10-22 |
CN104111038B true CN104111038B (zh) | 2016-08-17 |
Family
ID=51707926
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410317007.6A Active CN104111038B (zh) | 2014-07-07 | 2014-07-07 | 利用相位融合算法修复饱和产生的相位误差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104111038B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108645354A (zh) * | 2018-07-06 | 2018-10-12 | 四川大学 | 高反光对象表面的结构光三维成像方法及系统 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016145582A1 (zh) * | 2015-03-13 | 2016-09-22 | 香港应用科技研究院有限公司 | 相位偏移校准方法、3d形状检测的方法、系统及投影系统 |
CN104713497B (zh) * | 2015-03-13 | 2017-07-04 | 香港应用科技研究院有限公司 | 相位偏移校准方法、3d形状检测的方法、系统及投影系统 |
CN104729666B (zh) * | 2015-03-17 | 2017-10-17 | 东莞市神州视觉科技有限公司 | 一种正弦光栅最佳振幅强度的测量方法及装置 |
CN104697469B (zh) * | 2015-03-30 | 2017-11-10 | 四川大学 | 利用谐波幅值判定饱和的高动态三维测量方法 |
CN108240800B (zh) * | 2016-12-23 | 2020-03-10 | 致茂电子(苏州)有限公司 | 表面形貌的量测方法 |
CN106643559B (zh) * | 2017-03-06 | 2019-04-02 | 中国科学院光电技术研究所 | 一种基于混合干涉条纹的白光显微干涉形貌重建方法 |
CN109307483B (zh) * | 2018-11-20 | 2020-05-19 | 西南石油大学 | 一种基于结构光系统几何约束的相位展开方法 |
CN110057319B (zh) * | 2019-02-20 | 2021-05-25 | 西安理工大学 | 一种大范围反射率变化物体的高速三维测量方法 |
CN111241317B (zh) * | 2020-01-08 | 2022-03-11 | 四川大学 | 一种基于多个二维查找表的相位和调制信息获取方法 |
CN112781522A (zh) * | 2020-12-25 | 2021-05-11 | 复旦大学 | 一种基于彩色相移结构光的去高光轮廓仪 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101694375A (zh) * | 2009-10-23 | 2010-04-14 | 北京航空航天大学 | 一种用于强反射表面三维形貌测量的立体视觉检测方法 |
CN102906536A (zh) * | 2010-05-19 | 2013-01-30 | 株式会社尼康 | 形状测定装置及形状测定方法 |
CN102937425A (zh) * | 2012-10-18 | 2013-02-20 | 北京航空航天大学 | 一种基于高动态条纹投射器的强反射表面三维形貌测量系统 |
CN103615991A (zh) * | 2013-11-27 | 2014-03-05 | 东南大学 | 相位测量轮廓术中相位误差过补偿与欠补偿的解决方法 |
CN103837094A (zh) * | 2014-03-11 | 2014-06-04 | 四川大学 | 一种降低摄像机多路径效应影响的相位融合方法 |
CN103868473A (zh) * | 2014-03-31 | 2014-06-18 | 西北工业大学 | 一种基于递推法的高光物体表面相位快速恢复方法 |
-
2014
- 2014-07-07 CN CN201410317007.6A patent/CN104111038B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101694375A (zh) * | 2009-10-23 | 2010-04-14 | 北京航空航天大学 | 一种用于强反射表面三维形貌测量的立体视觉检测方法 |
CN102906536A (zh) * | 2010-05-19 | 2013-01-30 | 株式会社尼康 | 形状测定装置及形状测定方法 |
CN102937425A (zh) * | 2012-10-18 | 2013-02-20 | 北京航空航天大学 | 一种基于高动态条纹投射器的强反射表面三维形貌测量系统 |
CN103615991A (zh) * | 2013-11-27 | 2014-03-05 | 东南大学 | 相位测量轮廓术中相位误差过补偿与欠补偿的解决方法 |
CN103837094A (zh) * | 2014-03-11 | 2014-06-04 | 四川大学 | 一种降低摄像机多路径效应影响的相位融合方法 |
CN103868473A (zh) * | 2014-03-31 | 2014-06-18 | 西北工业大学 | 一种基于递推法的高光物体表面相位快速恢复方法 |
Non-Patent Citations (3)
Title |
---|
《一种抑制相位测量轮廓术饱和误差的方法》;赵婧 等;《中国激光》;20131010;第40卷(第10期);全文 * |
《多路径模型对相位测量轮廓术的影响》;许晓畅 等;《计算机辅助设计与图形学学》;20140415;第26卷(第4期);全文 * |
《投影条纹相移法中图像饱和误差抑制算法研究》;胡而已 等;《光子学报》;20100615;第39卷(第6期);全文 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108645354A (zh) * | 2018-07-06 | 2018-10-12 | 四川大学 | 高反光对象表面的结构光三维成像方法及系统 |
CN108645354B (zh) * | 2018-07-06 | 2019-10-18 | 四川大学 | 高反光对象表面的结构光三维成像方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104111038A (zh) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104111038B (zh) | 利用相位融合算法修复饱和产生的相位误差的方法 | |
CN110288642B (zh) | 基于相机阵列的三维物体快速重建方法 | |
US9523571B2 (en) | Depth sensing with depth-adaptive illumination | |
US9605961B2 (en) | Information processing apparatus that performs three-dimensional shape measurement, information processing method, and storage medium | |
CN105306922B (zh) | 一种深度相机参考图的获取方法和装置 | |
US20150116460A1 (en) | Method and apparatus for generating depth map of a scene | |
JP6519265B2 (ja) | 画像処理方法 | |
Ben-Hamadou et al. | Flexible calibration of structured-light systems projecting point patterns | |
JP6043735B2 (ja) | 画像評価装置及びパターン形状評価装置 | |
Dufour et al. | Integrated digital image correlation for the evaluation and correction of optical distortions | |
CN109474814A (zh) | 投影仪的二维校准方法、投影仪以及校准系统 | |
Yang et al. | Flexible and accurate implementation of a binocular structured light system | |
KR20170002409A (ko) | 이미지 분석에 의한 신의 3차원 재구성 디바이스 및 방법 | |
JP2015021862A (ja) | 3次元計測装置及び3次元計測方法 | |
IL299069B1 (en) | 3D structure or metrology inspection using deep learning | |
Sun et al. | Lens distortion correction for improving measurement accuracy of digital image correlation | |
JP2019091122A (ja) | デプスマップフィルタ処理装置、デプスマップフィルタ処理方法及びプログラム | |
KR102023087B1 (ko) | 카메라 캘리브레이션 방법 | |
JP6214867B2 (ja) | 計測装置、方法及びプログラム | |
CN109587463A (zh) | 投影仪的校准方法、投影仪及校准系统 | |
JP2011075336A (ja) | 3次元形状計測装置、3次元形状計測方法 | |
JP2021081791A5 (zh) | ||
Gottfried et al. | Time of flight motion compensation revisited | |
JP2020520028A (ja) | 光学測定装置をキャリブレーションするための方法 | |
JP2018151215A (ja) | ひび割れ幅特定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Liu Kai Inventor after: Sun Weiwen Inventor after: Li Na Inventor before: Liu Kai Inventor before: Li Na |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |