摘抄自博士论文第3.2.2节
虽然数字相移干涉技术基于其准确的相位恢复能力获得了大量的应用,然而其需要复杂的相移装置。
希尔伯特变换是在时域中发展起来的一种特殊的正交变换。已被广泛的应用于相位提取等方面。实函数
式中,*表示卷积运算。对式(3.8)进行傅里叶变换得到希尔伯特变换的频域表示为:
式中,
后来,时变场信号的复解析形式在光学领域也获得了广泛应用。特别的,复杂分析信号的实部与虚部之间的希尔伯特转换关系,通常被用于对单个空间干涉图进行相移。这种定量希尔伯特相位提取技术,允许从单个空间干涉图中提取光场相位信息。此外,该技术可直接对相位实现解包裹运算,这有助于研究远大于一个光波长的相位目标。
下面详细介绍基于希尔伯特变换的相位获取方法:
如图3.11所示,假设在成像平面上的待恢复散斑光场表示为
强度均匀的参考波为
之后让参考光与散斑光场进行干涉,则CCD上的干涉强度表示为
同时若能在上述参考光中引入π/2的相移,则可在CCD上得到光强分布:
运用上述记录的4幅图像,即可解调得出CCD平面上物光复振幅为
以上方法实验原理比较简单,但是式(3.12)的获取需要对参考光进行相移,这通常需要额外引入相移装置,进而加大了装置复杂度,减慢了成像时间。不过由上面式(3.9)的分析可知,希尔伯特变换能直接实现π/2相移,并滤除直流分量。
为了减小振幅的波动对相位恢复效果的影响,这里首先依据式(3.11)单独提取出其相位项:
之后将式(3.14)进行一次希尔伯特变换即可得到
H[.]表示希尔伯特变换运算。进一步可获得散斑光场相位信息:
一般来说,此处的
基于以上分析可知希尔伯特变换相位提取技术具有简单高效的优点,不过其也存在着局限性,特别是对于相位随机变化的散斑光场。这是因为希尔伯特变换的有效性严重依赖于相邻像素相位之间的相关性,因此原则上讲只对空间变化比较缓慢的信号有效。为此,进行了如下实验以验证其有效性及鲁棒性:
(1)相位变化缓慢的散斑光场
如图3.12所示,用两幅随机散斑图案来模拟散斑光场的振幅和相位分布,其中考虑到希尔伯特变换技术需要依赖于相邻像素相位之间的相关性,故作为相位的散斑图案其散斑颗粒相对更大。
参考光的振幅如图3.13(a)所示为常值,相位如图3.13(b)所示自左向右逐渐减小。
物光与参考光的干涉强度图如图3.13(c)所示,将其代入式(3.14),再结合式(3.15)和式(3.16)就能提取出散斑光场的初始相位,其结果如图3.14(a)所示。最后再去除由于参考光倾斜照明而叠加的额外相位,就能恢复出散斑光场的最终相位信息,结果如图3.14(b)所示。比较3.12(b)和图3.14(b)可知两者具有高度的相似性,进而验证了希尔伯特变换获取缓变散斑相位的有效性。
最后,为了定量地评估重建相位的准确性,运用式(3.7)计算了重建相位与原目标散斑相位的平均测量误差为0.0257弧度。由此可见,虽然希尔伯特变换技术对于相位缓变的散斑光场具有比较好的重建效果,但是不同于数字相移干涉技术在不存在光学噪声的情况下,具有接近0的相位重建误差。希尔伯特变换技术即使在不存在光学噪声的情况下,其重建误差也不为0。
(2)相位变化剧烈的散斑光场
更进一步地,当光场发生严重畸变时,散斑光场相邻像素之间的相位可能就不存在相关性,这时从理论上来说希尔伯特变换就会失效,为此这里也对相位变化更剧烈的散斑光场进行了仿真验证。当原始光场相位分布为如图3.15(a)所示的完全随机的散斑时,相位的重建结果如图3.15(b)所示,比较图3.15(a)和图3.15(b)可见两者的确存在较大的差异,之后通过式(3.7)计算表明其平均测量误差高达1.1395弧度。由于对于未解包裹的相位而言,其有效取值仅为0~2π,因此平均1.1395弧度的相位测量误差是非常大的。