▒▒本文目录▒▒
一、引言
相移技术是目前使用最多的相位提取技术之一。
几种常用的相移法有:五步相移法[1]、四步相移法[2]、三步相移法[3]和两步相移法[4]等。因其具有非常高的测量精度高而广泛地受到学者的关注。该方法的基本原理是为了求解相位信息,在不同时刻或在同一时刻的不同空间位置,通过对其中一支光路引入不同的相移来得到足够数量的散斑干涉条纹图。根据基本是三角函数变换公式对这几幅散斑干涉图进行相应计算就可得到散斑干涉图中所包含的与位移或变形等有关的相位信息。比其他相位提取方法相比相移法的测量精度更高,一般测量精度可以高达激光波长的百分之一。相移法通常可以分为时间相移法,空间相移法和空间载波相移法。在研究静态问题时,一般都首选时间相移法。该方法是用CCD相机在不同时刻对同一空间位置采集三幅或三幅以上的干涉条纹图进行相位求解。因而该方法在动态或振动过程中不能加以采用。空间相移法与时间相移法不同,是用CCD相机在同一时刻的不同空间位置同时采集三幅或三幅以上的散斑干涉条纹图进行相位求解。因为这些散斑干涉条纹图可以在同一时刻进行采集,因而该方法可以用于动态相位测量。但是其缺点是,要求CCD相机的性能要一致,并且需要匹配在不同空间的位置得到的散斑干涉条纹,这大大影响了测量精度。
二、时间相移法(Temporal Phase-shifting Method)
时间相移技术是目前应用在散斑干涉领域中较常见也发展的较成熟的一种相移技术。此方法由 Carre 于 1966 年最早提出方法构想,并由 Bruning 等人由 1974年付诸实践。它的思路是在光路中增加一个动态的相移器,对一路激光引入固定步长相移。在物体的单一变形状态保持稳定的条件下,拍摄多幅不同相移值的干涉图,并通过前述的相移公式进行解算。时间相移法有着理论精度较高、测量可全程自动化实现等优势。目前存在压电陶瓷法、光栅相移法、倾斜玻璃法等。其中压电陶瓷法因为其使用简单、工作稳定、重复性较好等优势被广泛使用。
2.1 三步相移算法原理与实例分析
假设三步法中三次相移量分别为-α、0、和α,有:
{ I 0 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) − α ] I 1 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) ] I 2 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) + α ] \left\{ \begin{matrix} {{I}_{0}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right)-\alpha \right] \\ {{I}_{1}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right) \right] \\ {{I}_{2}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right)+\alpha \right] \\ \end{matrix} \right. ⎩ ⎨ ⎧I0(x,y)=A(x,y)+B(x,y)cos[φ(x,y)−α]I1(x,y)=A(x,y)+B(x,y)cos[φ(x,y)]I2(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+α]
可改写为以下形式:
{ I 0 ( x , y ) = A ( x , y ) + B ( x , y ) [ cos φ ( x , y ) cos α + sin φ ( x , y ) sin α ] I 1 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) ] I 2 ( x , y ) = A ( x , y ) + B ( x , y ) [ cos φ ( x , y ) cos α − sin φ ( x , y ) sin α ] \left\{ \begin{matrix} {{I}_{0}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\left[ \cos \varphi \left( x,y \right)\cos \alpha +\sin \varphi \left( x,y \right)\sin \alpha \right] \\ {{I}_{1}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right) \right] \\ {{I}_{2}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\left[ \cos \varphi \left( x,y \right)\cos \alpha -\sin \varphi \left( x,y \right)\sin \alpha \right] \\ \end{matrix} \right. ⎩ ⎨ ⎧I0(x,y)=A(x,y)+B(x,y)[cosφ(x,y)cosα+sinφ(x,y)sinα]I1(x,y)=A(x,y)+B(x,y)cos[φ(x,y)]I2(x,y)=A(x,y)+B(x,y)[cosφ(x,y)cosα−sinφ(x,y)sinα]
通过求解方程组得:
当每步相移量为π/2时有:
通过实验获取三幅相移干涉图,根据三步相移公式,解得的包裹相位如下图所示。
2.2 四步相移算法原理与实例分析
四步相移算法是最经典的相移算法之一,取相移量分别为 0 π 2 π 3 π 2 0\frac{\pi }{2}\pi \frac{3\pi }{2} 02ππ23π
则对应的强度表达式为:
{ I 0 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) ] I 1 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) + π / 2 ] I 2 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) + π ] I 3 ( x , y ) = A ( x , y ) + B ( x , y ) cos [ φ ( x , y ) + 3 π / 2 ] \left\{ \begin{matrix} {{I}_{0}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right) \right] \\ {{I}_{1}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right)\text{+}{\pi }/{2}\; \right] \\ {{I}_{2}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right)\text{+}\pi \right] \\ {{I}_{3}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right)\text{+}{3\pi }/{2}\; \right] \\ \end{matrix} \right. ⎩ ⎨ ⎧I0(x,y)=A(x,y)+B(x,y)cos[φ(x,y)]I1(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+π/2]I2(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+π]I3(x,y)=A(x,y)+B(x,y)cos[φ(x,y)+3π/2]
将其改写为:
{
I
0
(
x
,
y
)
=
A
(
x
,
y
)
+
B
(
x
,
y
)
cos
[
φ
(
x
,
y
)
]
I
1
(
x
,
y
)
=
A
(
x
,
y
)
−
B
(
x
,
y
)
sin
[
φ
(
x
,
y
)
]
I
2
(
x
,
y
)
=
A
(
x
,
y
)
+
B
(
x
,
y
)
cos
[
φ
(
x
,
y
)
]
I
3
(
x
,
y
)
=
A
(
x
,
y
)
+
B
(
x
,
y
)
sin
[
φ
(
x
,
y
)
]
\left\{ \begin{matrix} {{I}_{0}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right) \right] \\ {{I}_{1}}\left( x,y \right)=A\left( x,y \right)-B\left( x,y \right)\sin \left[ \varphi \left( x,y \right) \right] \\ {{I}_{2}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\cos \left[ \varphi \left( x,y \right) \right] \\ {{I}_{3}}\left( x,y \right)=A\left( x,y \right)+B\left( x,y \right)\sin \left[ \varphi \left( x,y \right) \right] \\ \end{matrix} \right.
⎩
⎨
⎧I0(x,y)=A(x,y)+B(x,y)cos[φ(x,y)]I1(x,y)=A(x,y)−B(x,y)sin[φ(x,y)]I2(x,y)=A(x,y)+B(x,y)cos[φ(x,y)]I3(x,y)=A(x,y)+B(x,y)sin[φ(x,y)]
则包裹的相位信息可根据下式求解:
通过实验获取四幅相移干涉图,根据四步相移公式,解得的包裹相位如下图所示。
2.3 两步广义相移算法实例分析
通过实验获取两幅相移干涉图,通过两步广义算法[5],求解得到的包裹相位如下图所示。
2.4 对比分析
两步、三步、四步相移算法求解得到的包裹相位如下图所示。
三、其他相移算法
3.1 五步相移算法理论及其验证
五步相移算法理论推导及其验证,可参考博主这篇博文:一种数字全息散斑干涉测量仿真模拟系统|五步相移算法理论及其验证|点击链接即可跳转
参考文献
[1] 罗志勇, 陈朝晖, 顾英姿, et al. 基于数值模拟的高准确度五步相移算法研究 [J]. 光学学报, 2006, 26(11): 1687-90.
[2] 杨国标, 章彰. 电子散斑干涉法在铝合金材料结构内部柱状缺陷检测中的应用研究 [J]. 计量学报, 2008, 29(b09): 191-4.
[3] Peisen S Huang, Qingying J Hu, Chiang Fu-Pen. Double three-step phase-shifting algorithm [J]. Appl Opt, 2002, 41(22): 4503-9.
[4] Yang Fujun, He Xiaoyuan. Two-step phase-shifting fringe projection profilometry: intensity derivative approach [J]. Appl Opt, 2007, 46(29): 7172.
[5] 李校宇. 基于数字全息干涉光学图像加密新算法的研究 [D]; 天津大学, 2015.
⭐️◎⭐️◎⭐️◎⭐️ · · · **博 主 简 介** · · · ⭐️◎⭐️◎⭐️◎⭐️ ♪▁▂▃▅▆▇ 博士研究生生 ,研究方向主要涉及定量相位成像领域,具体包括干涉相位成像技术(如**全息干涉☑**、散斑干涉☑等)、非干涉法相位成像技术(如波前传感技术☑,相位恢复技术☑)、此外,还对各种相位解包裹算法☑,相干噪声去除算法☑等开展过深入的研究。