CN1302999A - 一种基于移相法的数字投影三维轮廓重构方法 - Google Patents
一种基于移相法的数字投影三维轮廓重构方法 Download PDFInfo
- Publication number
- CN1302999A CN1302999A CN01104163A CN01104163A CN1302999A CN 1302999 A CN1302999 A CN 1302999A CN 01104163 A CN01104163 A CN 01104163A CN 01104163 A CN01104163 A CN 01104163A CN 1302999 A CN1302999 A CN 1302999A
- Authority
- CN
- China
- Prior art keywords
- phase
- projection
- projection pattern
- projector
- reconstructing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000004744 fabric Substances 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 8
- 238000013500 data storage Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 244000309464 bull Species 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 210000001747 pupil Anatomy 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明属于三维传感与计量领域,本发明包括:预先编制投影图案生成程序,用计算机控制数字投影仪将投影图案按一定时间间隔依次投射到被测物体表面上;摄像机依次记录下每一幅投影图案被物体表面调制形成的变形图像;采集视频图像、存储、处理这些图像数据;根据物体上各点的相位信息以及投影仪、摄象机、被测物体之间的几何关系重构出物体的三维轮廓。本发明大大简化了传统测量***中的投影、相移部分,取得了较好的结果。
Description
本发明属于三维传感与计量领域,特别涉及三维物体表面形貌的计算机重构技术。
三维物体表面轮廓的计算机重构技术是目前国内外三维传感与计量研究的热点之一,它被广泛地应用于航空航天、机械制造、在线检测与质量控制、医疗诊断、计算机辅助设计/计算机辅助制造(CAD/CAM)、机器人的视觉***等领域。移相法是所有的三维形貌重构方法中最成熟、最可靠的一种方法,它具有计算量较少、可以用较粗光栅达到很高灵敏度和精度以及抗静态噪声等优点,是条件允许的情况下的首选方案。
移相法的工作原理描述如下:
图1a和图1b所示的是N步移相法。它将投影到物体表面的正弦光栅条纹移动N次,利用离散相移方法获取N幅(N>3)被调制的图像,根据N-PHASE算法计算出位相分布,最后利用几何关系求得物体表面的高度分布。
如图1a所示,正弦光场13以入射角θi入射到参考平面11上,CCD图像传感器12在θv方向上接收被参考平面11调制后的光场,其中z为接收到的光场周期。假设投射正弦光场的周期为x,并且平面移动了Δh,如图1b所示,这时z就会有一个变化量Δz:Δz=Δh[tan(θi)+tan(θv)]cos(θv) (1)这会在调制光场中引起一个相位变化Δφ:
由上式可以看出,参考平面上某处的高度变化Δh,可以由调制光场的相移Δφ求得。当CCD图像传感器垂直于参考平面时,θv为零,于是高度变化为: 只要我们获取到Δφ,就可以重构出物体的三维轮廓。
当一个由正弦光栅调节出来的光强按正弦分布的二维投影图案投射到漫反射物体表面上时,由成像***探测到的被调制的图像可以描述为:Ii(x,y)=a(x,y)[(l+m(x,y)·cos(φ(x,y)+φi)] (4)
其中,a(x,y)代表背景光强度,m(x,y)代表条纹对比度,相位函数φ(x,y)代表了条纹图案的特征,因而包含了物体表面轮廓的信息,φi为初始相位角,N为移相次数。通过调节光栅,可以产生不同初始相位的正弦投影图案。移动光栅N次(N≥3),每次移动P/N(P是正弦光栅的周期)的距离,则产生相位变化是 2π/N的N幅投影图案。用CCD图像传感器记录下每一幅图案调制后的图像,根据N-PHASE算法有: 于是Δφ可求。
传统的移相法测量***利用正弦光栅产生投影光场照到被测物体上,通过步进电机带动光栅移动N次(N≥3),每次移动P/N(P是正弦光栅的周期)的距离,来产生N幅不同初始相位的正弦投影图案,然后用摄像机依次记录下每一幅投影图案被物体表面调制形成的变形图像;计算机通过图像采集卡采集视频图像,并把图像数据存储起来;用N步移相法算法(N-PHASE算法)处理这些图像数据,得到物体上各点的相位信息;根据物体上各点的相位信息以及投影仪、摄象机、被测物体之间的几何关系重构出物体的三维轮廓。由于机械加工和步进电机的失步等各种原因,期望的光栅移动距离常常是不能准确得到的,这给测量结果带来了很大误差。传统的相移装置如图2所示,大都采用精密的微位移传感器来实时的监控光栅移动的距离,将实际移动的距离反馈给计算机,然后再根据偏差通过伺服放大器驱动步进电机带动工作台做相应的补偿。这种方法虽然可以精确的控制相移,但市场上工作稳定的微位移传感器价格昂贵,而且体积较大,使这套反馈式控制***的设计和安装十分不便。
本发明的目的在于为克服传统的移相法测量***相移装置成本高,设计、控制复杂的不足之处,提出了一种数字投影的方法,利用计算机从软件上实现投影光场和相移,大大简化了传统测量***中的投影、相移部分,取得了较好的结果。
本发明提出的数字投影三维轮廓重构方法,其特征在于:预先编制投影图案生成程序,用来生成三幅以上有固定相位差的正弦光栅投影图案;运行所述程序,用计算机控制数字投影仪将投影图案按一定时间间隔依次投射到被测物体表面上;摄像机依次记录下每一幅投影图案被物体表面调制形成的变形图像;计算机通过图像采集卡采集视频图像,并把图像数据存储起来;用N步移相法算法处理这些图像数据,得到物体上各点的相位信息;根据物体上各点的相位信息以及投影仪、摄象机、被测物体之间的几何关系重构出物体的三维轮廓。
上述投影图案的生成可以借助各种编程语言实现。由于生成的投影图案上从左到右光强是按正弦规律周期变化的,所以可以用计算机编程控制投影仪,从左到右依次按正弦灰度变化画一系列直线来实现正弦投影光场图案。
本发明所述的投影图像的采集、处理等步骤的处理均可采用传统的处理方法。
摄像机采集图像对于不同的图像采集卡,其实现方式也可能不同。例如,如果所用的图像采集卡是可编程的,那么可以编程让其在每幅图案投射期间采集几幅图像,选用效果最好的一幅作为这幅投影图案对应的调制图像进入下步的图像处理。如果图像采集卡是不能编程的,可以采用人为控制计算机的方法进行图像采集。
本发明通过采用数字投影的方法,有效的克服了传统的移相法测量***相移装置成本较高,设计、控制复杂的不足,获得了良好的结果。
附图简要说明:
图1为移相法原理示意图。
图2为传统相移装置的光栅移动控制图。
图3为本发明实施例测量装置示意图。
图4为本发明实施例测量装置几何关系示意图。
图5为本发明实施例的投影图案自动生成程序流程图。
图6为本发明实施例的图像处理流程图。
本发明提出的一种计算机控制投影仪来实现投影光场和相移的三维轮廓重构方法的实施例结合附图详细描述如下:
本发明以4幅投影图案为实施例,本实施例的测量***如图3所示,投影仪32采用EPSON7100型LCD投影仪,分辨率1024×768像素、照度650流明;图像传感器35是敏通TMC-56型彩色CCD摄象机,640×480像素、256级灰度(B/W输出);计算机31性能为奔腾233MHz,带有可编程的视频输入卡,33、34、36分别是投影图案、被测物体以及采集到的视频图像。
上述测量***采用的投影和成像***的几何关系如图4所示,P2和I2分别表示投影***和成像***的出瞳和入瞳,O是投影***和成像***光轴的交点,参考平面与成像***的光轴垂直。被测物体应尽量安置在关于O点对称的位置上,这样可以使投影、成像的视场分布均匀。这种位置情况下,计算公式(3)中
,剩下的就是用移相法测出Δφ。
本实施例的数字投影三维轮廓重构方法的具体实现步骤包括:首先,编制投影图案生成程序,其步骤如图5所示,首先定义变量,赋初值;然后定义绘图区域,完成初始化;再采用Visual Basic语言通过三重循环根据各幅图像的初始相位按正弦灰度规律画线。先定义控制投影图案数目、每幅投影图案投影周期数目、每个投影周期内画线数目以及它们各自的控制阈值变量分别为n、i、j、c、b、a,给n、i、j赋零。定义一个延时变量,本实施例中赋值为5秒。定义全屏的0.7倍宽度,全屏高度范围为绘图区域,绘图区域居中设置,完成绘图初始化。本实施例在绘图区域里按每个灰度变化周期画128条线,共画7个周期来绘图,那么只需把图6中的a取为127,b取为6,c取为3,用Visual Basic语言编程实现该三重循环即可。
其次,运行上面编好的程序,用计算机控制LCD投影仪将生成的投影图案按设定的5秒时间间隔依次投射到物体表面上。阅读可编程视频采集卡的编程说明,编制相应的程序,在这段时间间隔内,让CCD摄象机依次记录下每一幅投影图案被物体表面调制形成的变形图像,并把它们存储成BMP文件。
最后,进行图像处理,其过程如图6所示。首先要读入获得的图像文件,把其中的灰度数据存成矩阵,本实施例使用Matlab下的imread()函数直接读入BMP文件,并把其中的灰度数据存成矩阵;接下来对获得的数据进行滤波,去除其中的噪声干扰,本实施例利用了Matlab下功能强大的小波滤波和傅立叶变换工具包来完成这项工作;当依次读入、滤波处理过各幅图像数据后,根据N-PHASE算法公式(6)得到物体上各点的相位信息,通过***的几何关系由式(3)求出各点的高度,从而重构出物体的三维轮廓。
实验表明,本发明的方法大大简化了传统测量***中的投影、相移部分,取得了令人满意的结果。
Claims (1)
1、一种基于移相法的数字投影三维轮廓重构方法,其特征在于包括以下步骤:
(1)预先编制投影图案生成程序,用来生成三幅以上有固定相位差的正弦光栅投影图案;
(2)运行所述程序,用计算机控制数字投影仪将投影图案按一定时间间隔依次投射到被测物体表面上;
(3)摄像机依次记录下每一幅投影图案被物体表面调制形成的变形图像;
(4)计算机通过图像采集卡采集视频图像,并把图像数据存储起来;
(5)用N步移相法算法处理这些图像数据,得到物体上各点的相位信息;
(6)根据物体上各点的相位信息以及投影仪、摄象机、被测物体之间的几何关系重构出物体的三维轮廓。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011041633A CN1149388C (zh) | 2001-02-23 | 2001-02-23 | 一种基于移相法的数字投影三维轮廓重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011041633A CN1149388C (zh) | 2001-02-23 | 2001-02-23 | 一种基于移相法的数字投影三维轮廓重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1302999A true CN1302999A (zh) | 2001-07-11 |
CN1149388C CN1149388C (zh) | 2004-05-12 |
Family
ID=4653721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011041633A Expired - Fee Related CN1149388C (zh) | 2001-02-23 | 2001-02-23 | 一种基于移相法的数字投影三维轮廓重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1149388C (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1302257C (zh) * | 2003-10-09 | 2007-02-28 | 四川大学 | 用双色光栅模板产生π相移实现傅里叶变换轮廓术的方法 |
CN1312461C (zh) * | 2004-06-23 | 2007-04-25 | 香港理工大学 | 薄片型柔性体表面三维重建的***和方法 |
CN100408972C (zh) * | 2006-07-27 | 2008-08-06 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
CN100416466C (zh) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | 单目视觉半实物仿真***及方法 |
CN100417231C (zh) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | 立体视觉半实物仿真***及方法 |
CN101299227B (zh) * | 2008-06-27 | 2010-06-09 | 北京中星微电子有限公司 | 基于三维重构的多人游戏***及方法 |
CN101238488B (zh) * | 2005-08-04 | 2010-08-04 | 皇家飞利浦电子股份有限公司 | 用于产生或重建三维图像的方法和*** |
CN101256672B (zh) * | 2008-03-21 | 2011-10-12 | 北京中星微电子有限公司 | 基于视频摄像设备的物体图像深度重构装置及其投影装置 |
CN102261896A (zh) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | 一种基于相位测量的物体三维形貌测量方法及*** |
CN102589481A (zh) * | 2012-03-13 | 2012-07-18 | 周波 | 最优宽度调制下正弦光栅的设计 |
CN102589477A (zh) * | 2012-02-14 | 2012-07-18 | 乔付 | 光学三维轮廓测量中正弦光栅的设计 |
CN102589480A (zh) * | 2012-03-08 | 2012-07-18 | 乔付 | 结构光三维视觉测量中正弦光栅的设计 |
US8331688B2 (en) | 2009-01-07 | 2012-12-11 | International Business Machines Corporation | Focus-based edge detection |
CN104567721A (zh) * | 2015-01-23 | 2015-04-29 | 清华大学 | 连续剪切干涉测量方法 |
CN105547540A (zh) * | 2015-12-09 | 2016-05-04 | 清华大学 | 实时空间相移的相干梯度敏感干涉方法 |
CN107810384A (zh) * | 2015-06-17 | 2018-03-16 | 惠普德国股份有限公司 | 条纹投影方法、条纹投影装置和计算机程序产品 |
CN108784836A (zh) * | 2018-06-20 | 2018-11-13 | 安徽医科大学第附属医院 | 基于优化镇静管理和区域阻滞骨科麻醉术中图像处理*** |
-
2001
- 2001-02-23 CN CNB011041633A patent/CN1149388C/zh not_active Expired - Fee Related
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1302257C (zh) * | 2003-10-09 | 2007-02-28 | 四川大学 | 用双色光栅模板产生π相移实现傅里叶变换轮廓术的方法 |
CN1312461C (zh) * | 2004-06-23 | 2007-04-25 | 香港理工大学 | 薄片型柔性体表面三维重建的***和方法 |
CN101238488B (zh) * | 2005-08-04 | 2010-08-04 | 皇家飞利浦电子股份有限公司 | 用于产生或重建三维图像的方法和*** |
CN100416466C (zh) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | 单目视觉半实物仿真***及方法 |
CN100417231C (zh) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | 立体视觉半实物仿真***及方法 |
CN100408972C (zh) * | 2006-07-27 | 2008-08-06 | 西安交通大学 | 基于双频彩色条纹投影的三维物体轮廓相位测量方法 |
CN101256672B (zh) * | 2008-03-21 | 2011-10-12 | 北京中星微电子有限公司 | 基于视频摄像设备的物体图像深度重构装置及其投影装置 |
CN101299227B (zh) * | 2008-06-27 | 2010-06-09 | 北京中星微电子有限公司 | 基于三维重构的多人游戏***及方法 |
US8331688B2 (en) | 2009-01-07 | 2012-12-11 | International Business Machines Corporation | Focus-based edge detection |
US8509562B2 (en) | 2009-01-07 | 2013-08-13 | International Business Machines Corporation | Focus-based edge detection |
CN102261896A (zh) * | 2011-04-19 | 2011-11-30 | 长春东瑞科技发展有限公司 | 一种基于相位测量的物体三维形貌测量方法及*** |
CN102589477A (zh) * | 2012-02-14 | 2012-07-18 | 乔付 | 光学三维轮廓测量中正弦光栅的设计 |
CN102589480A (zh) * | 2012-03-08 | 2012-07-18 | 乔付 | 结构光三维视觉测量中正弦光栅的设计 |
CN102589480B (zh) * | 2012-03-08 | 2013-12-11 | 黑龙江科技学院 | 结构光三维视觉测量中正弦光栅的设计 |
CN102589481A (zh) * | 2012-03-13 | 2012-07-18 | 周波 | 最优宽度调制下正弦光栅的设计 |
CN102589481B (zh) * | 2012-03-13 | 2013-12-11 | 黑龙江科技学院 | 最优宽度调制下正弦光栅的设计 |
CN104567721A (zh) * | 2015-01-23 | 2015-04-29 | 清华大学 | 连续剪切干涉测量方法 |
CN107810384A (zh) * | 2015-06-17 | 2018-03-16 | 惠普德国股份有限公司 | 条纹投影方法、条纹投影装置和计算机程序产品 |
US10801834B2 (en) | 2015-06-17 | 2020-10-13 | Hewlett-Packard Development Company, L.P. | Fringe projection for determining topography of a body |
CN107810384B (zh) * | 2015-06-17 | 2020-11-20 | 惠普德国股份有限公司 | 条纹投影方法、条纹投影装置和计算机程序产品 |
CN105547540A (zh) * | 2015-12-09 | 2016-05-04 | 清华大学 | 实时空间相移的相干梯度敏感干涉方法 |
CN108784836A (zh) * | 2018-06-20 | 2018-11-13 | 安徽医科大学第附属医院 | 基于优化镇静管理和区域阻滞骨科麻醉术中图像处理*** |
Also Published As
Publication number | Publication date |
---|---|
CN1149388C (zh) | 2004-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1149388C (zh) | 一种基于移相法的数字投影三维轮廓重构方法 | |
CN102188290B (zh) | 用于获取牙齿的3d表面轮廓图像数据的装置和方法 | |
US7286246B2 (en) | Method and apparatus for non-contact three-dimensional surface measurement | |
EP1994503B1 (en) | Depth-varying light fields for three dimensional sensing | |
Lindner et al. | Lateral and depth calibration of PMD-distance sensors | |
US8050461B2 (en) | Depth-varying light fields for three dimensional sensing | |
JP3525964B2 (ja) | 物体の三次元形状計測方法 | |
US10706570B2 (en) | System and method to acquire the three-dimensional shape of an object using a moving patterned substrate | |
JP7319913B2 (ja) | マシンビジョン検査システムを用いてワークピース表面のz高さ値を測定するための方法 | |
EP3064895A1 (en) | Linear illumination device | |
JPH0338524B2 (zh) | ||
CN111536904A (zh) | 一种基于结构照明的三维测量方法及***、存储介质 | |
Jia et al. | Two-step triangular-pattern phase-shifting method for three-dimensional object-shape measurement | |
Forest Collado | New methods for triangulation-based shape acquisition using laser scanners | |
EP3839599A1 (en) | Metrology system with transparent workpiece surface mode | |
Traxler et al. | Experimental comparison of optical inline 3D measurement and inspection systems | |
WO2004097335A1 (en) | Shape and deformation measurements of large objects by fringe projection | |
CN111982026A (zh) | 一种高反射物体曲面的光栅投影三维测量装置和测量方法 | |
JP2006170744A (ja) | 3次元距離計測装置 | |
JPH07260451A (ja) | 3次元形状測定システム | |
US11928851B2 (en) | Computer-implemented method, computer program, data processing system and device for determining the reflectance behavior of a surface of an object, and storage medium with instructions stored thereon for determining the reflectance behavior of a surface of an object | |
Uyanik et al. | A method for determining 3D surface points of objects by a single camera and rotary stage | |
Kammel | Topography reconstruction of specular surfaces from a series of gray-scale images | |
Gray et al. | Simulation of a tire inspection system | |
CN215114433U (zh) | 一种基于光学位相元件的结构光三维形貌测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |