CN104132628A - 用相位计算实现线结构光三维测量的方法 - Google Patents
用相位计算实现线结构光三维测量的方法 Download PDFInfo
- Publication number
- CN104132628A CN104132628A CN201410358743.6A CN201410358743A CN104132628A CN 104132628 A CN104132628 A CN 104132628A CN 201410358743 A CN201410358743 A CN 201410358743A CN 104132628 A CN104132628 A CN 104132628A
- Authority
- CN
- China
- Prior art keywords
- phase
- structured light
- phase place
- line
- testee
- 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)
Abstract
本发明是三维传感技术中一种用相位计算实现线结构光三维测量的方法。用投影设备对被测物体进行线扫描,用摄像装置记录下线结构光图像序列,采用相位测量轮廓术的相位计算公式分别计算出基频相位和多个高频截断相位,根据时间相位展开方法得到绝对相位,再根据亮度调制滤除相位图中信噪比低的区域,最后用相位和系统标定参数计算出被测物体表面的三维坐标。本发明可用于包括激光线扫描的线结构光三维测量技术。本发明方案具有计算效率高、抗干扰能力强和测量精度高的优点。
Description
技术领域
本发明涉及光学三维传感技术,特别是涉及通过投影线结构光实现对目标物体表面的三维测量。
背景技术
基于结构光投影的三维测量技术因其非接触、测量精度高的优点在工业生产和科学研究领域有广泛的应用,如工业检测、逆向工程、三维生物特征识别、辅助医疗诊断等。结构光技术包括以下基本步骤:用投影装置向目标物体投射结构光编码图案,再用摄像装置同步拍摄,通过解码得到摄像空间与投影空间的对应关系,进而利用三角原理获得目标物体的三维坐标。相位测量轮廓术(Phase Measuring Profilometry,简称PMP)是一种经典的结构光三维测量技术,具有测量精度高、对环境光和被测物体纹理等干扰因素不敏感的优点。相位测量轮廓术采用具有相移的正弦结构光图案序列,利用捕捉到的图像序列计算相位,并通过相位信息获取被测物体表面的三维坐标。线结构光扫描,特别是激光线扫描,被广泛应用于工业三维测量。现有线结构光图案的解码方案主要是基于对摄像头获取的激光条纹的峰值进行检测。峰值检测的精度直接影响三维测量结果的精度,论文“Naidu, D. K., and Robert B. Fisher. "A comparative analysis of algorithms for determining the peak position of a stripe to sub-pixel accuracy." BMVC91. Springer London, 1991. 217-225”对多种峰值检测算法进行了综述。现有基于峰值检测的线结构光测量方法有以下缺点:(1)由于峰值检测采用邻域操作,如差分、查找最大值等操作,检测过程受环境光、被测物体反射率不一致等干扰因素影响较大,抗干扰能力较弱,在复杂测量环境下精度难以保证。(2)为了保证峰值检测算法的稳定性和测量精度,一般需要额外使用复杂的图像处理方法,如曲线拟合、异常点检测与间断点插值等,计算过程繁琐,计算效率不高。如何克服现有方法计算效率低、抗干扰能力弱、精度有限的缺点,在保证测量精度的同时提高计算效率以及增强抗干扰能力,应用本发明提及的方案就可以解决这一关键技术问题。
发明内容
本发明的目的是针对目前基于峰值检测的线结构光三维测量方法计算效率低、抗干扰能力弱、精度有限的缺陷,提供一种基于相位计算的方法,该方法在保证高测量精度的同时具有较高的计算效率和很强的抗干扰能力。
本发明的目的是采用下述技术方案来实现的:
用投影设备对被测物体进行线扫描,用摄像装置记录下线结构光图像序列,用相位测量轮廓术的相位计算公式分别计算出基频相位和多个高频截断相位,根据时间相位展开方法得到绝对相位,再利用相位测量轮廓术的计算公式得到图像序列的亮度调制,根据亮度调制信息滤除相位图中信噪比低的区域,用相位和系统标定参数计算出被测物体表面的三维坐标。
本发明与现有技术相比有如下优点:
与现有基于峰值检测的线结构光三维测量方法相比,因为本发明无需繁琐的图像处理,直接计算相位,所以算法实现简单,计算效率更高;由于本发明在对图像序列的相位计算过程中不涉及邻域操作,即对每个像素点的操作是独立进行的,且对环境光和被测物体纹理等干扰因素不敏感,所以抗干扰能力更强,系统运行更稳定;因为本发明通过对高频截断相位展开得到最终相位,所以能效抑制随机噪声,且计算出的相位具有子像素精度,从而能保证较高的三维测量精确度。
附图说明
图1为本发明三维测量系统结构图。
图2为本发明三维测量方法的工作流程图。
具体实施方式
下面结合附图、工作原理对本发明作进一步详细说明。
附图1为本发明三维测量系统结构图,图中101为摄像头,102为投影机,103为计算机,104为被测样品。采用的装置有1台CASIO XJ-M140投影机,投影机缓存帧大小为 像素,灰度量化等级为8bit,投影机最大输出频率为150帧/s;1个Prosilica GC650工业摄像头,分辨率为像素,灰度量化等级为8bit,摄像头最大捕捉频率为62帧/s。1台具有Core i3 3530 CPU,4GB内存的计算机。由计算机对结构光投影和拍摄过程进行控制。附图2为本实施例线结构光三维测量方法流程图。本实例具体实施步骤如下:
(1)对摄像头和投影机进行标定,分别得到摄像头与投影机大小为的投影矩阵、。
(2)生成线结构光图案。垂直方向扫描的线结构光图案序列可以表示为:
(1)
其中,表示投影机空间坐标;表示线结构光图案在像素点处的灰度值;为投影线条纹的幅值;表示单位冲激函数,即满足,且当时;表示投影机空间的高度;表示扫描线的位移系数,当分别取时,分别对应扫描线处于第行的图案。这里线结构光图案参数取值为:,。
(3)对目标物体进行结构光扫描。用投影机把上述结构光图案序列依次投射到目标物体表面,并用一个摄像头对目标场景进行同步拍摄。拍摄到的线结构光图像序列表示为,其中表示摄像头空间坐标,表示位移系数为的摄像头图像在像素点处的灰度值。
(4)用拍摄到的线结构光图像解码得到相位。计算相位的公式为:
(2)
其中,表示摄像头图像在像素点处的相位值;表示该相位对应的空间频率。当取时,计算出的是基频相位;而当取高频时,上式计算出的是出未展开的高频相位。分别取基频和高频,这里具体取值为,用上式分别计算出基频相位和未展开的高频相位,然后通过时间相位展开法用基频相位依次对高频相位进行展开,最终得到展开后的高频相位。多频相位展开的具体方法在技术报告“Yalla V G, Hassebrook L G. Very high resolution 3D surface scanning using multi-frequency phase measuring profilometry[C]//Defense and Security. International Society for Optics and Photonics, 2005: 44-53”中有详细介绍。设展开后高频相位的取值范围为,摄像头与投影机的坐标对应关系可以表示为:。需要指出的是,当时,上述相位计算公式,即式(2),是相位测量轮廓术中标准的相位计算公式;当取高频时,式(2)为双频相位测量轮廓术中高频相位的计算公式,其中双频相位测量轮廓术在论文“Liu, Kai, et al. "Dual-frequency pattern scheme for high-speed 3-D shape measurement." Optics express 18.5 (2010): 5229-5244”中有详细介绍。
(5)根据亮度调制滤除相位图中信噪比低的区域。由计算亮度调制的公式为:
(3)
其中,表示摄像头图像在像素点处的亮度调制。摄像头图像中某像素点值的大小能反映该点相位是否有效。值较大区域对应有效测量区域,而值接近于0的区域则对应信噪比较低的区域或背景阴影区域。故用作为滤波器过滤掉相位图中的无效区域,即设定一个阈值,这里取值为,把摄像头坐标区域判定为有效测量区域,保留该区域内的相位,而把区域判定为无效区域,舍弃该区域内的相位。
(6)根据相位以及系统标定参数和可以计算出被测物体表面的三维坐标。具体计算方法在论文“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所述的方法,其特征在于所说的用相位测量轮廓术的相位计算公式得到图像序列的相位,是通过用相位测量轮廓术的相位计算公式分别计算出基频相位和多个高频截断相位,根据时间相位展开方法得到绝对相位,这些操作的目的是利用展开高频相位抑制随机噪声,并获得高精确度的相位信息。
4.按照权利要求1所述的方法,其特征在于所说的滤除相位图中信噪比低的区域,是利用相位测量轮廓术的计算公式得到图像序列的亮度调制,设定一个阈值,把相位图中亮度调制小于该阈值的像素区域滤除,目的是去除信噪比低的无效相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410358743.6A CN104132628B (zh) | 2014-07-27 | 2014-07-27 | 用相位计算实现线结构光三维测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410358743.6A CN104132628B (zh) | 2014-07-27 | 2014-07-27 | 用相位计算实现线结构光三维测量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104132628A true CN104132628A (zh) | 2014-11-05 |
CN104132628B CN104132628B (zh) | 2017-06-23 |
Family
ID=51805403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410358743.6A Active CN104132628B (zh) | 2014-07-27 | 2014-07-27 | 用相位计算实现线结构光三维测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104132628B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697469A (zh) * | 2015-03-30 | 2015-06-10 | 四川大学 | 利用谐波幅值判定饱和的高动态三维测量方法 |
CN105698708A (zh) * | 2016-01-29 | 2016-06-22 | 四川大学 | 一种三维视觉重建方法 |
CN107018398A (zh) * | 2017-03-08 | 2017-08-04 | 清华大学 | 一种用于光场三维显示量化标定的方法 |
KR101819141B1 (ko) | 2015-12-30 | 2018-01-16 | 인천대학교 산학협력단 | 푸리에 해석을 이용한 라인 구조화된 광 패턴 디코딩 방법 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008031648A1 (de) * | 2006-09-13 | 2008-03-20 | Phoenix Conveyor Belt Systems Gmbh | Einrichtung zur überwachung einer förderanlage |
CN101762243A (zh) * | 2010-01-04 | 2010-06-30 | 北京航空航天大学 | 一种受限空间三维形貌结构光视觉测量方法 |
CN101806587A (zh) * | 2010-04-29 | 2010-08-18 | 浙江师范大学 | 一种采用绝对相位测量的光学三维测量方法 |
CN101813462A (zh) * | 2010-04-16 | 2010-08-25 | 天津理工大学 | 单处理器控制的三维形貌光学测量系统及测量方法 |
CN103729251A (zh) * | 2013-11-06 | 2014-04-16 | 中国科学院上海光学精密机械研究所 | 并行计算光学条纹图相位提取方法 |
US20140132730A1 (en) * | 2007-05-17 | 2014-05-15 | Ilya Blayvas | Compact 3D Scanner with Fixed Pattern Projector and Dual Band Image Sensor |
CN103913128A (zh) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | 一种基于质量图区域分割的快速相位展开方法 |
CN103940371A (zh) * | 2014-05-12 | 2014-07-23 | 电子科技大学 | 一种用于跃变物体的高精度三维面形测量的方法 |
-
2014
- 2014-07-27 CN CN201410358743.6A patent/CN104132628B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008031648A1 (de) * | 2006-09-13 | 2008-03-20 | Phoenix Conveyor Belt Systems Gmbh | Einrichtung zur überwachung einer förderanlage |
US20140132730A1 (en) * | 2007-05-17 | 2014-05-15 | Ilya Blayvas | Compact 3D Scanner with Fixed Pattern Projector and Dual Band Image Sensor |
CN101762243A (zh) * | 2010-01-04 | 2010-06-30 | 北京航空航天大学 | 一种受限空间三维形貌结构光视觉测量方法 |
CN101813462A (zh) * | 2010-04-16 | 2010-08-25 | 天津理工大学 | 单处理器控制的三维形貌光学测量系统及测量方法 |
CN101806587A (zh) * | 2010-04-29 | 2010-08-18 | 浙江师范大学 | 一种采用绝对相位测量的光学三维测量方法 |
CN103729251A (zh) * | 2013-11-06 | 2014-04-16 | 中国科学院上海光学精密机械研究所 | 并行计算光学条纹图相位提取方法 |
CN103913128A (zh) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | 一种基于质量图区域分割的快速相位展开方法 |
CN103940371A (zh) * | 2014-05-12 | 2014-07-23 | 电子科技大学 | 一种用于跃变物体的高精度三维面形测量的方法 |
Non-Patent Citations (3)
Title |
---|
吴庆阳等: "线结构光双传感器测量系统的标定方法", 《中国激光》 * |
国家自然科学基金委员会工程与材料科学部等组编: "《国家自然科学基金机械学科自助项目结题成果汇编(2005年)》", 31 July 2010 * |
张建保等编: "《中国生物医学工程进展(上册)》", 30 April 2007 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697469A (zh) * | 2015-03-30 | 2015-06-10 | 四川大学 | 利用谐波幅值判定饱和的高动态三维测量方法 |
CN104697469B (zh) * | 2015-03-30 | 2017-11-10 | 四川大学 | 利用谐波幅值判定饱和的高动态三维测量方法 |
KR101819141B1 (ko) | 2015-12-30 | 2018-01-16 | 인천대학교 산학협력단 | 푸리에 해석을 이용한 라인 구조화된 광 패턴 디코딩 방법 |
CN105698708A (zh) * | 2016-01-29 | 2016-06-22 | 四川大学 | 一种三维视觉重建方法 |
CN105698708B (zh) * | 2016-01-29 | 2018-06-01 | 四川大学 | 一种三维视觉重建方法 |
CN107018398A (zh) * | 2017-03-08 | 2017-08-04 | 清华大学 | 一种用于光场三维显示量化标定的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104132628B (zh) | 2017-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9965870B2 (en) | Camera calibration method using a calibration target | |
JP5317169B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
US9633439B2 (en) | Image processing system, and image processing method | |
KR20190065432A (ko) | 뷰로부터 거리 정보를 획득하는 장치 및 방법 | |
JP2012058076A (ja) | 3次元計測装置及び3次元計測方法 | |
JP5633058B1 (ja) | 3次元計測装置及び3次元計測方法 | |
CN110337674B (zh) | 三维重建方法、装置、设备及存储介质 | |
CN102721376A (zh) | 一种大视场三维视觉传感器的标定方法 | |
CN112184811B (zh) | 单目空间结构光系统结构校准方法及装置 | |
CN104132628A (zh) | 用相位计算实现线结构光三维测量的方法 | |
CN111256628A (zh) | 墙面平整度检测方法、装置、计算机设备和存储介质 | |
CN104200456B (zh) | 一种用于线结构光三维测量的解码方法 | |
KR20230065978A (ko) | 구조화된 광을 사용하여 장면에서 평면 표면들을 직접 복구하기 위한 시스템, 방법 및 매체 | |
CN113074667A (zh) | 基于标志点的全局绝对相位对齐方法、存储介质和系统 | |
JP5656018B2 (ja) | 球体の検出方法 | |
CN113551617B (zh) | 基于条纹投影的双目双频互补三维面型测量方法 | |
Hongsheng et al. | Three-dimensional reconstruction of complex spatial surface based on line structured light | |
CN113450460A (zh) | 一种基于人脸面形空间分布的无相位展开三维人脸重建方法及系统 | |
KR101770866B1 (ko) | 뎁스 센서와 촬영 카메라 간의 관계에 기초한 뎁스 영상 보정 장치 및 방법 | |
EP2953096B1 (en) | Information processing device, information processing method, system and carrier means | |
US20230090275A1 (en) | Systems and Methods for Stereo Depth Sensing | |
Liu et al. | Vision measurement method for impact point in large planar region | |
JP2016061622A (ja) | 情報処理システム、及び情報処理方法 | |
CN117808858A (zh) | 一种管线数据采集系统、方法、装置、设备及存储介质 | |
Zhao et al. | A review on 3D measurement of highly reflective objects using structured light projection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |