CN104315996A - 用二进制编码策略实现傅里叶变换轮廓术的方法 - Google Patents
用二进制编码策略实现傅里叶变换轮廓术的方法 Download PDFInfo
- Publication number
- CN104315996A CN104315996A CN201410554298.0A CN201410554298A CN104315996A CN 104315996 A CN104315996 A CN 104315996A CN 201410554298 A CN201410554298 A CN 201410554298A CN 104315996 A CN104315996 A CN 104315996A
- Authority
- CN
- China
- Prior art keywords
- binary
- pattern
- fourier transform
- image
- phase
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种用二进制编码策略实现傅里叶变换轮廓术的方法。把一幅高频正弦结构光图案按位平面分解为一系列二进制编码图案，使用投影装置向三维目标投射该二进制图案序列，并用摄像装置进行同步拍摄，对拍摄的二进制图像进行二值化处理，再根据位平面与灰度图像的对应关系合成一幅高频正弦结构光图像，利用傅里叶变换轮廓术处理合成的正弦结构光图像，得到截断相位，对截断相位进行相位展开得到绝对相位，用绝对相位和***标定参数计算出物体表面的三维坐标。本发明可用于测量具有漫反射特性物体的三维面形。本发明能有效避免传统傅里叶变换轮廓术存在的频谱混叠问题和背景光干扰，具有较高的测量精度。
Description
技术领域
本发明涉及光学三维传感技术，特别是涉及基于傅里叶变换轮廓术对空间物体的三维面形测量。
背景技术
傅里叶变换轮廓术（Fourier Transform Profilometry，简称FTP）是一种应用广泛的非接触式光学三维测量技术，论文“Su, X., & Chen, W. (2001). Fourier transform profilometry: a review. Optics and lasers in Engineering, 35(5), 263-284”详细介绍了傅里叶变换轮廓术的原理。因为傅里叶变换轮廓术一般只需要一帧高频结构光图案，所以它的显著优点是能实时对运动物体进行三维测量。傅里叶变换轮廓术的核心步骤是在频域用带通滤波器滤除直流分量和高次谐波，提取主频分量。然而，背景光亮度和物体反射率在空间上的变化、随机噪音以及摄像头-投影机***的非线性均会产生谐波，而这些谐波与主频分量在频谱上混叠，从而使傅里叶变换轮廓术无法准确提取主频分量，导致严重的测量误差。此外，传统傅里叶变换轮廓术对背景光也很敏感。上述原因致使傅里叶变换轮廓术虽然实时性较好，但测量精确度较低。二进制编码技术，作为另一种典型的结构光三维测量技术，由于其使用的投影图案只具有两级灰度，所以解码的鲁棒性较强，相对于连续灰度编码具有诸多优势：一方面，二进制编码技术能有效抑制随机噪音；另一方面，它不受测量***的非线性、背景光、物体表面反射率等因素影响。此外，现有投影设备扫描二进制图案的帧率大大高于扫描灰度图案的帧率，这为利用二进制编码技术实施高速三维扫描提供了保证，例如论文“Li, B., Wang, Y., Dai, J., Lohry, W., & Zhang, S. (2014). Some recent advances on superfast 3D shape measurement with digital binary defocusing techniques. Optics and Lasers in Engineering, 54, 236-246”报道了采用二值图案投影技术实现了KHz的三维测量速度。如何把二进制编码技术与傅里叶变换轮廓术有机结合起来，在保证傅里叶变换轮廓术测量快速性的基础上，解决背景光亮度和物体反射率在空间上的变化、随机噪音以及摄像头-投影机***的非线性导致的频谱混叠问题，消除背景光的对傅里叶变换轮廓术的影响，从而提高傅里叶变换轮廓术的测量精确度，应用本发明提及的方案就可以解决这一关键技术问题。
发明内容
本发明的目的是针对传统傅里叶变换轮廓术由于背景光亮度和物体反射率在空间上的变化、随机噪音以及摄像头-投影机***的非线性导致的频谱混叠问题，提出一种三维传感技术中采用二进制编码策略的傅里叶变换轮廓术，这种方法能解决频谱混叠问题，且不受背景光的影响，在保证快速测量的基础上，具有较高的测量精度。
本发明的目的是采用下述技术方案来实现的：
把一幅高频正弦结构光图案按位平面分解为一系列二进制编码图案，使用投影装置向三维目标投射该二进制图案序列，并用摄像装置进行同步拍摄，对拍摄的二进制图像进行二值化处理，提取出二进制编码信息，根据位平面与灰度图像的对应关系合成一幅高频正弦结构光图像，利用傅里叶变换轮廓术处理合成的正弦结构光图像，从而计算出截断相位，对截断相位进行相位展开得到绝对相位，用绝对相位和测量***标定参数计算出物体表面的三维坐标。
本发明与现有技术相比有如下优点：
因为本发明在扫描过程中使用二进制编码图案，所以本发明方法不受测量***的非线性、背景光、物体表面反射率等因素影响，且能有效抑制随机噪声，所以本发明能克服传统傅里叶变换轮廓术的频谱混叠问题和背景光干扰，提高测量精确度；相对于传统傅里叶变换轮廓术，本发明由于采用了二进制编码图案，对摄像头亮度饱和具有更强的鲁棒性。
附图说明
图1为结构光三维测量***结构图
图2为本发明采用二进制编码策略的傅里叶变换轮廓术的工作流程图。
具体实施方式
下面结合附图、工作原理对本发明作进一步详细说明。
附图1为本发明三维测量***结构图。采用的装置有1台CASIO XJ-M140投影机，投影机缓存帧大小为 像素，灰度量化等级为8bit，投影机最大输出频率为150帧/s；1个Prosilica GC650工业摄像头，分辨率为像素，灰度量化等级为8bit，摄像头最大捕捉频率为62帧/s；1台具有Core i3 3530 CPU，4GB内存的计算机，由计算机对结构光投影和拍摄过程进行控制。附图2为本实施例二进制编码的傅里叶变换轮廓术流程图。本实例具体实施步骤如下：
（1）对摄像头和投影机进行标定，分别得到摄像头与投影机大小为的投影矩阵、。
（2）生成一幅高频正弦结构光图案，该结构光图案可以表示为：
（1）
其中，为投影机空间坐标，表示结构光图案在像素点处的灰度值，是投影机空间的高度，是结构光图案的空间频率，为结构光图案的直流分量，为结构光图案的交流振幅，满足。本实例正弦结构光图案的参数为，空间频率取f=32。
（3）把正弦结构光图案按位平面分解为二进制结构光图案序列。正弦结构光图案与其分解产生的二进制结构光图案的关系可以表示为：
（2）
其中，表示第个比特平面的二进制图案在像素点处的灰度值，表示连续灰度结构光图案分解产生的二进制图案的数目。由于格雷码比直接二进制编码误码率更低，所以根据二进制数与格雷码的对应关系把直接二进制图案序列转换为对应的格雷码图案序列。本实例中。
（4）用二进制图案序列对目标物体进行结构光扫描。向目标物体投射格雷码图案序列。中只有0和1两级灰度，在实际投影过程中分别对应投影机的最低亮度级和最高亮度级。用摄像头同步捕捉每一幅格雷码图像，并用阈值法对捕捉图像进行二值化，二值化处理后的图像表示为，其中表示摄像头图像坐标，表示在像素点处的灰度值。
（5）根据位平面与灰度图像的对应关系，用二进制图像序列合成高频正弦结构光图像。由二进制数与格雷码的对应关系把格雷码图像序列转换为对应的直接二进制图像序列。具有连续灰度级的正弦结构光图像可以用合成，合成公式为：
（3）
其中表示在摄像头空间合成的余弦图像在像素点处的灰度值。也可以表示为：
（4）
其中表示摄像头空间像素点对应的相位值。
（6）用傅里叶变换轮廓术从合成的正弦结构光图像中提取相位。为了去除结构光图像中的直流分量，我们通过减去已知常量对进行规范化处理。规范化后的可以表示为：
（5）
式（5）的复数形式可以写作：
（6）
利用傅里叶变换轮廓术从中提取相位，傅里叶变换轮廓术的相位计算过程在论文“Lohry, W., & Zhang, S. (2012). Fourier transform profilometry using a binary area modulation technique. Optical Engineering, 51(11), 113602-113602.”中有详细介绍，现简述如下：首先，固定，沿方向对进行一维傅里叶变换；第二，在频域用带通滤波器提取的主频分量；第三，进行傅里叶逆变换，滤波后的信号可以表示为：
（7）
用可以计算出截断相位：
， （8）
其中和分别表示取复数的虚部和实部。最后，对截断相位进行空间相位展开得到绝对相位。
（7）用绝对相位以及***标定参数和可以计算出被测物体表面的三维坐标。具体计算方法在论文“Li J, Hassebrook L G, Guan C. Optimized two-frequency phase-measuring-profilometry light-sensor temporal-noise sensitivity[J]. JOSA A, 2003, 20(1): 106-115.”中有详细介绍。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种用二进制编码策略实现傅里叶变换轮廓术的方法，其特征在于把一幅高频正弦结构光图案按位平面分解为一系列二进制编码图案，使用投影装置向三维目标投射该二进制图案序列，并用摄像装置进行同步拍摄，对拍摄的二进制编码图像进行处理，提取出二进制编码信息，根据位平面与灰度图像的对应关系合成一幅高频正弦结构光图像，利用傅里叶变换轮廓术处理合成的正弦结构光图像，从而计算出截断相位，对截断相位进行相位展开得到绝对相位，用绝对相位和测量***标定参数计算出物体表面的三维坐标。
2.按照权利要求1所述的方法，其特征在于所说的二进制编码图案，是指与高频结构光图案具有唯一对应关系的二进制图案序列，其包括直接二进制编码图案和格雷码图案。
3.按照权利要求1所述的方法，其特征在于所说的二进制编码图案，是指灰度图案或颜色编码图案，其中颜色编码图案是在RGB颜色通道中存放多幅二进制图案，使用颜色编码图案的目的是减少投影图案数量，提高扫描速度。
4.按照权利要求1所述的方法，其特征在于所说的对拍摄的二进制编码图像进行处理，是从拍摄的颜色编码图像中分离出各颜色通道的二进制图像，以及对二进制图像进行二值化处理，从而得到每个像素点上的二进制编码信息（即0或1）。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410554298.0A CN104315996B (zh) | 2014-10-20 | 2014-10-20 | 用二进制编码策略实现傅里叶变换轮廓术的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410554298.0A CN104315996B (zh) | 2014-10-20 | 2014-10-20 | 用二进制编码策略实现傅里叶变换轮廓术的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104315996A true CN104315996A (zh) | 2015-01-28 |
| CN104315996B CN104315996B (zh) | 2018-04-13 |
Family
ID=52371262
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410554298.0A Active CN104315996B (zh) | 2014-10-20 | 2014-10-20 | 用二进制编码策略实现傅里叶变换轮廓术的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104315996B (zh) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104596444A (zh) * | 2015-02-15 | 2015-05-06 | 四川川大智胜软件股份有限公司 | 一种基于编码图案投影的三维照相***及方法 |
| CN104634323A (zh) * | 2015-02-15 | 2015-05-20 | 四川川大智胜软件股份有限公司 | 一种多级摄像三维照相***及方法 |
| CN104634277A (zh) * | 2015-02-12 | 2015-05-20 | 北京唯创视界科技有限公司 | 拍摄设备和方法、三维测量***、深度计算方法和设备 |
| CN105066904A (zh) * | 2015-07-16 | 2015-11-18 | 太原科技大学 | 基于相位梯度阈值的流水线产品三维面型检测方法 |
| CN105806259A (zh) * | 2016-04-29 | 2016-07-27 | 东南大学 | 一种基于二值光栅离焦投影的三维测量方法 |
| CN106979757A (zh) * | 2017-03-03 | 2017-07-25 | 浙江华睿科技有限公司 | 一种三维测量方法及装置 |
| WO2018171385A1 (zh) * | 2017-03-24 | 2018-09-27 | 南京理工大学 | 一种基于改进的傅立叶变换轮廓技术的超快三维形貌测量方法及其*** |
| CN109000587A (zh) * | 2018-08-30 | 2018-12-14 | 无锡信捷电气股份有限公司 | 获取精确高密度点云的方法 |
| CN109799928A (zh) * | 2017-11-16 | 2019-05-24 | 清华大学深圳研究生院 | 投影触摸板中用户手指参数的获取方法和*** |
| CN110057317A (zh) * | 2018-11-16 | 2019-07-26 | 中山大学 | 一种基于奇异值分解和查找表构造的投影噪声消除方法 |
| CN110440714A (zh) * | 2019-09-05 | 2019-11-12 | 南昌航空大学 | 一种基于多频和二进制条纹的相位解包裹方法 |
| CN115670392A (zh) * | 2023-01-04 | 2023-02-03 | 北京智拓视界科技有限责任公司 | 用于对被扫描对象的脸部表情进行采集的三维扫描装置 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61153507A (ja) * | 1984-12-26 | 1986-07-12 | Toyota Central Res & Dev Lab Inc | 三次元形状認識装置 |
| CN101608906A (zh) * | 2009-07-21 | 2009-12-23 | 西安交通大学 | 一种空间编码的光学三维轮廓测量方法 |
| CN101655358A (zh) * | 2009-07-01 | 2010-02-24 | 四川大学 | 采用颜色编码提高正交复合光栅位相测量轮廓术的动态特性 |
| CN101697233A (zh) * | 2009-10-16 | 2010-04-21 | 长春理工大学 | 一种基于结构光的三维物体表面重建方法 |
| CN101986098A (zh) * | 2010-09-21 | 2011-03-16 | 东南大学 | 基于三色光栅投影的傅里叶变换三维测量法 |
| CN102798354A (zh) * | 2012-08-28 | 2012-11-28 | 西北工业大学 | 基于二进制条纹叠加的正弦光栅生成方法 |
| JP2013190406A (ja) * | 2012-03-12 | 2013-09-26 | Sekisui Chemical Co Ltd | 三次元形状測定装置 |
| JP2014178124A (ja) * | 2013-03-13 | 2014-09-25 | Ricoh Co Ltd | ３次元計測システム、プログラム及び方法。 |
-
2014
- 2014-10-20 CN CN201410554298.0A patent/CN104315996B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61153507A (ja) * | 1984-12-26 | 1986-07-12 | Toyota Central Res & Dev Lab Inc | 三次元形状認識装置 |
| CN101655358A (zh) * | 2009-07-01 | 2010-02-24 | 四川大学 | 采用颜色编码提高正交复合光栅位相测量轮廓术的动态特性 |
| CN101608906A (zh) * | 2009-07-21 | 2009-12-23 | 西安交通大学 | 一种空间编码的光学三维轮廓测量方法 |
| CN101697233A (zh) * | 2009-10-16 | 2010-04-21 | 长春理工大学 | 一种基于结构光的三维物体表面重建方法 |
| CN101986098A (zh) * | 2010-09-21 | 2011-03-16 | 东南大学 | 基于三色光栅投影的傅里叶变换三维测量法 |
| JP2013190406A (ja) * | 2012-03-12 | 2013-09-26 | Sekisui Chemical Co Ltd | 三次元形状測定装置 |
| CN102798354A (zh) * | 2012-08-28 | 2012-11-28 | 西北工业大学 | 基于二进制条纹叠加的正弦光栅生成方法 |
| JP2014178124A (ja) * | 2013-03-13 | 2014-09-25 | Ricoh Co Ltd | ３次元計測システム、プログラム及び方法。 |
Non-Patent Citations (2)
| Title |
|---|
| YANJUN FU等: "Fourier transform profiometry based on defocusing", 《OPTICS & LASER TECHNOLOGY》 * |
| 宋倩等: "基于格雷编码投影的三维测量技术研究", 《激光与光电子学进展》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104634277A (zh) * | 2015-02-12 | 2015-05-20 | 北京唯创视界科技有限公司 | 拍摄设备和方法、三维测量***、深度计算方法和设备 |
| CN104634323A (zh) * | 2015-02-15 | 2015-05-20 | 四川川大智胜软件股份有限公司 | 一种多级摄像三维照相***及方法 |
| CN104596444B (zh) * | 2015-02-15 | 2017-03-22 | 四川川大智胜软件股份有限公司 | 一种基于编码图案投影的三维照相***及方法 |
| CN104596444A (zh) * | 2015-02-15 | 2015-05-06 | 四川川大智胜软件股份有限公司 | 一种基于编码图案投影的三维照相***及方法 |
| CN105066904A (zh) * | 2015-07-16 | 2015-11-18 | 太原科技大学 | 基于相位梯度阈值的流水线产品三维面型检测方法 |
| CN105806259A (zh) * | 2016-04-29 | 2016-07-27 | 东南大学 | 一种基于二值光栅离焦投影的三维测量方法 |
| CN105806259B (zh) * | 2016-04-29 | 2018-08-10 | 东南大学 | 一种基于二值光栅离焦投影的三维测量方法 |
| CN106979757B (zh) * | 2017-03-03 | 2019-03-26 | 浙江华睿科技有限公司 | 一种三维测量方法及装置 |
| CN106979757A (zh) * | 2017-03-03 | 2017-07-25 | 浙江华睿科技有限公司 | 一种三维测量方法及装置 |
| WO2018171385A1 (zh) * | 2017-03-24 | 2018-09-27 | 南京理工大学 | 一种基于改进的傅立叶变换轮廓技术的超快三维形貌测量方法及其*** |
| US11029144B2 (en) | 2017-03-24 | 2021-06-08 | Nanjing University Of Science And Technology | Super-rapid three-dimensional topography measurement method and system based on improved fourier transform contour technique |
| CN109799928A (zh) * | 2017-11-16 | 2019-05-24 | 清华大学深圳研究生院 | 投影触摸板中用户手指参数的获取方法和*** |
| CN109799928B (zh) * | 2017-11-16 | 2022-06-17 | 清华大学深圳研究生院 | 投影触摸板中用户手指参数的获取方法和*** |
| CN109000587A (zh) * | 2018-08-30 | 2018-12-14 | 无锡信捷电气股份有限公司 | 获取精确高密度点云的方法 |
| CN110057317A (zh) * | 2018-11-16 | 2019-07-26 | 中山大学 | 一种基于奇异值分解和查找表构造的投影噪声消除方法 |
| CN110440714A (zh) * | 2019-09-05 | 2019-11-12 | 南昌航空大学 | 一种基于多频和二进制条纹的相位解包裹方法 |
| CN110440714B (zh) * | 2019-09-05 | 2021-01-08 | 南昌航空大学 | 一种基于多频和二进制条纹的相位解包裹方法 |
| CN115670392A (zh) * | 2023-01-04 | 2023-02-03 | 北京智拓视界科技有限责任公司 | 用于对被扫描对象的脸部表情进行采集的三维扫描装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104315996B (zh) | 2018-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104315996A (zh) | 用二进制编码策略实现傅里叶变换轮廓术的方法 | |
| Zhang | Recent progresses on real-time 3D shape measurement using digital fringe projection techniques | |
| KR101475382B1 (ko) | 광학적 3차원 측량의 자기 적응 윈도우 푸리에 위상추출방법 | |
| CN107578464A (zh) | 一种基于线激光扫描的传送带工件三维轮廓测量方法 | |
| CN105069789B (zh) | 基于编码网格模板的结构光动态场景深度获取方法 | |
| CN106091986B (zh) | 一种适用于光亮表面的三维测量方法 | |
| EP2881702A1 (en) | Image processing system, and image processing method | |
| CN105844633B (zh) | 基于De序列和相位编码的单帧结构光深度获取方法 | |
| CN101788274B (zh) | 一种彩色复合光栅的三维面形测量方法 | |
| CN111028295A (zh) | 一种基于编码结构光和双目的3d成像方法 | |
| CN104236482A (zh) | 结合几何标定的相位测量轮廓术***非线性校正方法 | |
| CN105157614B (zh) | 基于二值相移图案的三维测量方法 | |
| CN108596008B (zh) | 针对三维人脸测量的面部抖动补偿方法 | |
| Li et al. | A three dimensional on-line measurement method based on five unequal steps phase shifting | |
| CN111222472A (zh) | 一种基于结构光频域特征的人脸识别方法 | |
| Chen et al. | High-speed 3D surface profilometry employing trapezoidal phase-shifting method with multi-band calibration for colour surface reconstruction | |
| Je et al. | Color-phase analysis for sinusoidal structured light in rapid range imaging | |
| Deng et al. | High-speed and high-accuracy fringe projection profilometry without phase unwrapping | |
| Wang et al. | Highly reflective surface measurement based on dual stereo monocular structured light system fusion | |
| CN104132628B (zh) | 用相位计算实现线结构光三维测量的方法 | |
| Xu et al. | Realtime 3D profile measurement by using the composite pattern based on the binary stripe pattern | |
| CN114136236A (zh) | 一种结合极线约束的相位轮廓自适应投影测量方法 | |
| CN111462331A (zh) | 扩展对极几何并实时计算三维点云的若干方法 | |
| WO2019238583A1 (en) | Deflectometric techniques | |
| KR20060095173A (ko) | 조명-반사율 모델에 기반한 영상의 왜곡 보정 시스템 및 방법 |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liu Kai Inventor after: Long Yunfei Inventor after: Liu Yiguang Inventor after: Wu Wei Inventor after: Yang Xiaomin Inventor after: Zheng Xiaojun Inventor before: Liu Kai Inventor before: Long Yunfei Inventor before: Zheng Xiaojun Inventor before: Wu Wei Inventor before: Yang Xiaomin |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |