CN101655359A - 用非等步相移正交光栅实现实时三维面形测量的方法 - Google Patents
用非等步相移正交光栅实现实时三维面形测量的方法 Download PDFInfo
- Publication number
- CN101655359A CN101655359A CN200910059847A CN200910059847A CN101655359A CN 101655359 A CN101655359 A CN 101655359A CN 200910059847 A CN200910059847 A CN 200910059847A CN 200910059847 A CN200910059847 A CN 200910059847A CN 101655359 A CN101655359 A CN 101655359A
- Authority
- CN
- China
- Prior art keywords
- grating
- phase
- step phase
- fourier transform
- phase shift
- 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.)
- Pending
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明是三维传感技术中用非等步相移正交光栅对物体的三维面形和变形进行实时测量的方法。运用计算机进行设计编码非等步相移正交光栅,采用光学投影手段将光栅的面结构调制模板投影成像在物体表面上,再用摄像装置记录下物体的变形条纹图,随后用傅里叶变换频谱滤波方法和非等步相移位相测量轮廓术处理这些变形条纹,能够精确地恢复出物体的三维面形,进一步分析数据处理结果可以获得诸如物体变形等一些数字化信息。本发明可以用于具有漫反射特性的物体表面特征的测量研究。本发明具有测量精度高,响应时间快,能够实时获取物体的三维面形数据等优点。
Description
一、技术领域
本发明涉及光学三维传感技术,特别是涉及基于位相测量轮廓术方法的用非等步相移正交光栅实现对物体的实时三维面形的测量。
二、技术背景
三维面形测量,在机器视觉、生物医学、工业检测、快速成型、影视特技、产品质量控制等领域具有重要意义。光学三维传感技术,由于其具有非接触、精度高、大面积测量、易于自动控制等优点获得广泛的研究和应用。现有的光学三维传感方法主要包括:三角测量法、莫尔条纹法(MoiréTopography,简称MT)、傅里叶变换轮廓术(Fourier Transform Profilometry,简称FTP)、空间相位检测术(Spatial Phase Detection,简称SPD)、位相测量轮廓术(Phase MeasuringProfilometry,简称PMP)等,这些方法都是通过对受三维物体面形调制的空间结构光场进行解调制,来获得物体的三维面形信息。其中最常用的空间结构光场三维传感方法是傅立叶变换轮廓术和位相测量轮廓术。傅里叶变换轮廓术是通过对变形条纹图像进行傅里叶变换、频域滤波和逆傅里叶变换等步骤实现的。傅里叶变换轮廓术只需要用一帧条纹图来重建三维面形,实时性较好,可以用于动态过程的三维传感;但由于其涉及到滤波操作,频谱混叠会降低测量精度,该方法对环境光也比较敏感。位相测量轮廓术则需要从多帧相移条纹图形来重建三维面形,具有很高的精度,但由于采用多次相移,实时性较差。如何综合位相测量轮廓术和傅里叶变换轮廓术的优点而尽量减少其缺点,应用本发明提及的方法就可以解决这一关键技术难题。
三、发明内容
本发明的目的则是针对对一帧条纹进行解码三维重建精度较低的缺陷,提出一种在三维传感技术测量中用非等步相移正交光栅实现实时的三维面形测量方法。这种方法能非常好地做到实时获得物体表面变形条纹分布的图像信息,具有较高的测量精度,能真正地实现动态和瞬态测量。
本发明的目的是采用下述技术方案来实现的:
采用计算机设计编码产生所需要的正交光栅图案,即对透射光场进行面结构编码,然后运用光学手段将数字的光栅图案成像在胶片上,即将数字的透过率分布转化成模拟的透过率分布,脱离了数字投影仪等电子器件本身的特性对系统的干扰。
本发明与现有技术相比有如下优点:
1.本发明使用非等步相移正交光栅,使三维传感光学系统结构能够实时测量物体的三维面形。
2.同样是获取一帧变形条纹图案,本发明相比傅立叶变换轮廓术具有更高的测量精度;而相比位相测量轮廓术则不需要进行相移,以及不需要采集多帧变形条纹图案。
3.本发明通过计算机可以灵活设计所需要的光栅编码,通过光学手段产生胶片的过程也非常简单,因此可以在很短的时间内得到所需要的光栅。
四、附图说明
图1位相测量轮廓术的光路示意图。
图2待测心形物体。
图3通过计算机设计编码的非等步相移正交光栅。
图4受到物体面形调制的非等步相移正交光栅变形条纹。
图5用上述技术方案恢复的心形物体。
五、具体实施方式
下面结合附图、工作原理对本发明作进一步详细说明。
用非等步相移正交光栅实现实时三维面形测量的光路与传统的位相测量轮廓术的光路相似。图1是PMP方法的投影光路,P1和P2是投影系统的入瞳和出瞳,I2和I1是成像系统的入瞳和出瞳。成像光轴垂直于参考平面,并与投影光轴的夹角为θ,它们相交于参考平面上的O点,d为探测器光心到投影设备光心之间的距离,l为探测器光心到参考平面之间的距离。
用计算机设计一组非等步相移光栅,其投影强度像表示为:
a表示直流偏置,b/a表示条纹对比度,表示每次移动的相位,将这组非等步相移光栅分别调制在一组载频正弦光栅上,载频光栅的栅线方向和相移光栅的栅线方向垂直,将各个调制成分叠加,得到的非等步正交光栅的投影强度像为:
频率ξn为载频正弦光栅的频率,c为正交光栅的直流偏置,d/c为条纹对比度。当一块非等步正交光栅被投影到物体表面时,不考虑物体的反射率和环境光强度,CCD接收到的变形条纹可以表示为:
其中φ为物体的高度引起的相移光栅条纹的相位变化。对(3)进行二维快速傅里叶变换(FFT),得到的频谱IFFT(ξ,η)如下:
式中ξ,η为空间频率,从(4)式可知I(x,y)经过FFT后得到的频谱图中一共有6N+1个的频谱成分,其中cδ(ξ,η)为携带背景信息的零频;含有ηy的才是携带相移信息的频谱。为了能够有效提取出相应的相移光栅条纹,要求尽可能的使各个频谱成分相互独立,没有混叠,因此可通过选择合适的载频光栅的频率来使各个频谱在频谱图中尽量独立。当满足上述条件时,选择适当的滤波器,可从整个频谱中提取第n个相移光栅变形条纹的频谱成分IFFTn(x,y):
对(5)进行快速逆傅里叶变换(IFFT)得IRn,即:
IRn是一虚数,对IRn取模即可解调出一相移光栅变形条纹IIMAGn:
可以看出,从非等步相移正交光栅中解调出的相移光栅变形条纹IIMAGn与所编码非等步相移正交光栅中设计的非等步相移光栅变形条纹In(x,y)只是在强度上有一个倍数关系,可以从经解调后的一组相移光栅变形条纹IIMAGn用进行相位提取。
当非等步正交光栅投影到一个三维漫反射物体表面上时,通过成像系统可以获得受物体表面面形调制的变形条纹,根据离散相移技术获取3幅变形光栅的数学表达式为(8)、(9)及(10)式:
I0=R[A+Bcos(2πx/P+φ)] (8)
I1=R[A+Bcos(2πx/P+φ+α)] (9)
I2=R[A+Bcos(2πx/P+φ+β)] (10)
其中R(x,y)是表示物体表面的反射率,A(x,y)表示背景光强,B(x,y)为受物体面形反射率影响的光场调制强度,P为光栅周期,φ(x,y)为调制光场的相位分布函数,它表征了条纹的变形,并且与物体的三维面形有关。为了叙述方便,式子中的各量均省略(x,y),α与β表示,此时α≠β、并且都在2π范围内。整理(8)、(9)、(10)式得:
I1-I0=RB[cos(2πx/P+φ+α)-cos(2πx/P+φ)](11)
I2-I0=RB[cos(2πx/P+φ+β)-cos(2πx/P+φ)](12)
令M=2πx/P+φ,当I2-I0≠0、cosM≠0时,由(11)、(12)式联立求解,得到(12)式:
即:
φ=M-2πx/P (14)
(14)式即为物体截断相位的表达式,但此时计算的结果都是在分母不为零的情况下获得的,下面考虑当分母为零时的情况:
当I2-I0=0时:
φ=-β/2-2πx/P (15)
同样当cosM=0时:
φ=π/2-2πx/P (16)
综合(14)、(15)、(16)三式得到物体的截断相位为:
通过(17)式的反正切计算,相位被截断在-π→π之间,因此此时的相位是离散的,而实际物体表面的相位是连续的,故需要对截断相位进行展开,设展开后的相位分布函数为ψ(x,y),从投影光路的几何关系可以看出,物体的高度信息被编码在三维相位分布ψ(x,y)中,而Δψ(x,y)=ψ(x,y)-ψ0(x,y)对应着待测物体的真实高度分布h(x,y),式中ψ0(x,y)为参考平面的展开位相,因此只需要求出Δψ(x,y)后再利用相位和高度的对应关系式:
即可恢复出物体的三维物体表面面形高度分布。
图2是用计算机设计生成的非等步相移正交光栅。
图3是受到物体面形调制的非等步相移正交光栅变形条纹。
图4是用上述技术方案恢复的物体。
Claims (5)
1.一种适用于三维重建中使用非等步相移正交光栅实现对物体三维面形实时测量的方法,其特征在于使用照明光源,将非等步相移正交光栅的模板图案投影到物体表面,用摄像装置记录下变形的光场图像,采用二维傅里叶变换获取光场图像的空间频谱,通过空间滤波和二维傅里叶逆变换提取非等步相移光栅变形条纹图,用非等步相移解相和高度映射手段得到物体瞬时的三维面形分布。
2.按照权利要求1所述的方法,其特征在于所说的非等步相移正交光栅的模板,是指运用计算机设计编码、采用光学手段实现的投影面结构调制模板,包括非等步相移正交正弦光栅、非等步相移正交罗奇光栅编码模板。
3.按照权利要求1所述的方法,其特征在于所说的对获得的变形光场图像进行处理,是利用快速傅里叶变换的方法对变形条纹进行傅里叶变换、滤波、逆傅里叶变换、相位提取和展开等处理和操作,目的是通过解调相位的方法获取物体的面形分布,也包括利用空间位相检测方法对变形光场进行处理,获取物体的面形分布。
4.按照权利要求1所述的方法,其特征在于所说的对获得的变形光场图像的非等步解相算法解相位,目的是获取更加准确的解调相位,提高物体的三维面形分布重构精度。
5.按照权利要求1所述的方法,其特征在于运用非等步相移正交光栅,研究物体的瞬间状态,进行模式识别或提取物体的某些特征参数,如面形、姿态、变形量、特征点距离等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910059847A CN101655359A (zh) | 2009-07-01 | 2009-07-01 | 用非等步相移正交光栅实现实时三维面形测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910059847A CN101655359A (zh) | 2009-07-01 | 2009-07-01 | 用非等步相移正交光栅实现实时三维面形测量的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101655359A true CN101655359A (zh) | 2010-02-24 |
Family
ID=41709741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910059847A Pending CN101655359A (zh) | 2009-07-01 | 2009-07-01 | 用非等步相移正交光栅实现实时三维面形测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101655359A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003948A (zh) * | 2010-09-21 | 2011-04-06 | 暨南大学 | 一种高精度光学条纹位相提取方法 |
CN102288589A (zh) * | 2011-06-23 | 2011-12-21 | 深圳大学 | 单个荧光探针的定位方法、装置及超分辨成像方法、系统 |
CN102519393A (zh) * | 2011-11-15 | 2012-06-27 | 四川大学 | 用两个正交正弦光栅实现快速调制度测量轮廓术的方法 |
CN102589479A (zh) * | 2012-03-06 | 2012-07-18 | 天津大学 | 三维形貌中心摄动复合光栅投影测量方法及装置 |
CN102589478A (zh) * | 2012-02-29 | 2012-07-18 | 天津工业大学 | 一种应用于多频率三维测量的全局相位解相方法 |
CN104270591A (zh) * | 2014-09-12 | 2015-01-07 | 深圳市索恩达电子有限公司 | 一种可调节成像图像的光栅投影器及其光学测量系统 |
CN104330052A (zh) * | 2014-11-21 | 2015-02-04 | 天津工业大学 | 外差式三频不等步相移解相位方法 |
CN105066907A (zh) * | 2015-07-27 | 2015-11-18 | 南京理工大学 | 一种针对傅里叶轮廓术的载频相位去除方法 |
CN108613637A (zh) * | 2018-04-13 | 2018-10-02 | 深度创新科技(深圳)有限公司 | 一种基于参考图像的结构光系统解相方法及系统 |
CN109903377A (zh) * | 2019-02-28 | 2019-06-18 | 四川川大智胜软件股份有限公司 | 一种无需相位展开的三维人脸建模方法及系统 |
CN111238398A (zh) * | 2020-03-02 | 2020-06-05 | 四川大学 | 一种基于概率分布函数的相移误差检测方法 |
CN113029042A (zh) * | 2021-05-25 | 2021-06-25 | 四川大学 | 一种高温熔融态金属表面形貌动态测量装置及方法 |
CN113514009A (zh) * | 2021-08-06 | 2021-10-19 | 哈尔滨理工大学 | 一种移位阶梯相位码与相移条纹非对称组合三维测量方法 |
-
2009
- 2009-07-01 CN CN200910059847A patent/CN101655359A/zh active Pending
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003948A (zh) * | 2010-09-21 | 2011-04-06 | 暨南大学 | 一种高精度光学条纹位相提取方法 |
CN102003948B (zh) * | 2010-09-21 | 2012-02-29 | 暨南大学 | 一种高精度光学条纹位相提取方法 |
CN102288589A (zh) * | 2011-06-23 | 2011-12-21 | 深圳大学 | 单个荧光探针的定位方法、装置及超分辨成像方法、系统 |
CN102288589B (zh) * | 2011-06-23 | 2013-05-08 | 深圳大学 | 单个荧光探针的定位方法、装置及超分辨成像方法、系统 |
CN102519393A (zh) * | 2011-11-15 | 2012-06-27 | 四川大学 | 用两个正交正弦光栅实现快速调制度测量轮廓术的方法 |
CN102589478A (zh) * | 2012-02-29 | 2012-07-18 | 天津工业大学 | 一种应用于多频率三维测量的全局相位解相方法 |
CN102589479A (zh) * | 2012-03-06 | 2012-07-18 | 天津大学 | 三维形貌中心摄动复合光栅投影测量方法及装置 |
CN104270591A (zh) * | 2014-09-12 | 2015-01-07 | 深圳市索恩达电子有限公司 | 一种可调节成像图像的光栅投影器及其光学测量系统 |
CN104330052B (zh) * | 2014-11-21 | 2016-09-28 | 天津工业大学 | 外差式三频不等步相移解相位方法 |
CN104330052A (zh) * | 2014-11-21 | 2015-02-04 | 天津工业大学 | 外差式三频不等步相移解相位方法 |
CN105066907A (zh) * | 2015-07-27 | 2015-11-18 | 南京理工大学 | 一种针对傅里叶轮廓术的载频相位去除方法 |
CN105066907B (zh) * | 2015-07-27 | 2017-09-01 | 南京理工大学 | 一种针对傅里叶轮廓术的载频相位去除方法 |
CN108613637A (zh) * | 2018-04-13 | 2018-10-02 | 深度创新科技(深圳)有限公司 | 一种基于参考图像的结构光系统解相方法及系统 |
CN108613637B (zh) * | 2018-04-13 | 2020-04-07 | 深度创新科技(深圳)有限公司 | 一种基于参考图像的结构光系统解相方法及系统 |
CN109903377A (zh) * | 2019-02-28 | 2019-06-18 | 四川川大智胜软件股份有限公司 | 一种无需相位展开的三维人脸建模方法及系统 |
CN109903377B (zh) * | 2019-02-28 | 2022-08-09 | 四川川大智胜软件股份有限公司 | 一种无需相位展开的三维人脸建模方法及系统 |
CN111238398A (zh) * | 2020-03-02 | 2020-06-05 | 四川大学 | 一种基于概率分布函数的相移误差检测方法 |
CN111238398B (zh) * | 2020-03-02 | 2021-06-04 | 四川大学 | 一种基于概率分布函数的相移误差检测方法 |
CN113029042A (zh) * | 2021-05-25 | 2021-06-25 | 四川大学 | 一种高温熔融态金属表面形貌动态测量装置及方法 |
CN113029042B (zh) * | 2021-05-25 | 2021-08-03 | 四川大学 | 一种高温熔融态金属表面形貌动态测量装置及方法 |
CN113514009A (zh) * | 2021-08-06 | 2021-10-19 | 哈尔滨理工大学 | 一种移位阶梯相位码与相移条纹非对称组合三维测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101655359A (zh) | 用非等步相移正交光栅实现实时三维面形测量的方法 | |
CN101655360A (zh) | 一种可绝对位相展开的十六步双频光栅相移轮廓术 | |
CN105783775B (zh) | 一种镜面及类镜面物体表面形貌测量装置与方法 | |
CN103940371A (zh) | 一种用于跃变物体的高精度三维面形测量的方法 | |
CN104330052B (zh) | 外差式三频不等步相移解相位方法 | |
CN102589479B (zh) | 三维形貌中心摄动复合光栅投影测量方法及装置 | |
CN101655358A (zh) | 采用颜色编码提高正交复合光栅位相测量轮廓术的动态特性 | |
CN105066906B (zh) | 一种快速高动态范围三维测量方法 | |
CN106032976B (zh) | 基于波长选择的三条纹投影相位展开方法 | |
EP3012576B1 (en) | Method for measuring a contour of an object | |
CN102607466A (zh) | 高反射自由曲面光栅投影快速非接触测量方法及装置 | |
CN102261896A (zh) | 一种基于相位测量的物体三维形貌测量方法及系统 | |
CN104006765A (zh) | 单幅载频干涉条纹相位提取方法及检测装置 | |
CN103791856A (zh) | 基于四幅光栅条纹图像的相位求解与去包裹方法 | |
CN102410819B (zh) | 一种测量膜基反射镜三维面形的方法 | |
CN108955571A (zh) | 双频外差与相移编码相结合的三维测量方法 | |
CN109631798A (zh) | 一种基于π相移方法的三维面形垂直测量方法 | |
CN105953749A (zh) | 一种光学三维形貌测量方法 | |
CN104680534B (zh) | 基于单帧复合模板的物体深度信息获取方法 | |
CN109141290A (zh) | 一种大偏离量光学自由曲面面形的检测方法 | |
CN101298981A (zh) | 平面倾斜角度的光栅衍射条纹测量法 | |
CN101718533B (zh) | 镜像莫尔测量装置及方法 | |
CN100529656C (zh) | 利用电子散斑干涉载频调制技术测量物体形貌的方法 | |
CN111947600B (zh) | 基于相位级次代价滤波的鲁棒立体相位展开方法 | |
CN106840039B (zh) | 一种基于复合光栅投影的快速三维测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20100224 |