基于光强传输方程离散余弦变换求解的数字全息相位解调方法

本文链接：https://blog.csdn.net/weixin_43608857/article/details/146873187

在众多干涉测量应用中，例如数字全息［1- 5］、合成孔径雷达［6- 8］、磁共振成像［9］和条纹投影轮廓仪［10- 11］等，由于检索到的相位被包裹到区间［−π，π］之间，必须通过相位解包裹（PU）才能获得真实的连续相位分布。在过去的几十年中，已经提出了很多 PU 方法，可将其分为两类：时间［12- 16］和空间方法［17- 29］。时间方法，如多频方法［13- 14］和外差方法［15- 16］，PU 相位图是随时间变化来获取，并且不需要相邻像素信息，因此可以避免缠绕相位缺陷的影响，但是这些方法需要较多包裹相位或多频条纹图。空间方法利用空间信息展开相位，与时间方法相比约束条件更少。其中路径跟踪方法通过解决路径积分问题来展开相位图［17- 20］，例如分支方法［19- 22］、质量引导路径跟踪（QGPF）方法［21- 22］和掩码切割算法［23］。路径跟踪方法采用的是逐点比较法，对于理想条件，该方法是有效的，它的计算量小、效率高。但是，如果存在噪声或相位跳跃，这些方法可能会在展开的相位图上带来条纹噪声；此外，由于 PU 过程依赖于积分路径，结果并不唯一［18］。最小范数方法是另一种类型的空间解调方法［21- 27］，它擅长处理带有噪声或相位跳跃的展开相位图［29］。局部噪声可以通过最小化定义的误差度量来过滤。最小二乘方法（LS）是最常用的最小范数 PU 方法之一。它通过泊松方程将包裹相位的梯度与绝对相的相位关联，然后利用展开结果的导数与包裹相位数据之间的差异来获得展开相位图。研究学者提出了诸多方法来提高 LS 方法的效率，例如基于傅里叶变换的 LS 方法［24］、基于离散余弦变换的 LS（DCTLS）方法［23］，以及多重网格方法［24］。这些方法通常具有良好的鲁棒性，但有时倾向于过度平滑相位图。为了克服这个缺陷，加权最小二乘 PU 算法［26，29］被提出。加权系数通常首先从质量图中生成并用于相位展开过程，以限制来自平滑效应带来的误差，该算法可以有效屏蔽无效相位数据对邻近区域的不利影响，具有较好的稳定性。然而，获得质量图和权重系数可能很耗时，这将限制 PU 方法的实时应用。

近年来，为提高 PU 工艺的效率，新的技术已被引用。CruzSantos 等［30］引入了一种基于 redblack tree 数据结构的新剪枝策略来加速 PU 过程，但应首先实施一个额外的过程来获得数据结构。Kulkarni 等［31］提出了一种基于多项式相位逼近和线性卡尔曼滤波器的噪声鲁棒 PU 算法。然而该方法使用了线扫描方法或像素选择的质量引导策略。张祥朝等［32］将非下采样 Contourlet 变换应用于 PU 算法以提高鲁棒性，然而区域分割必须在展开过程之前进行，必然会增加处理时间

最近，学术界还提出了基于光强传输方程（TIE）的 PU 方法［33- 35］。它从强度信息中获得绝对相位图。首先计算强度的轴向导数。然后将结果用作泊松方程的输入，并且可以通过使用快速傅里叶变换（FFT）求解方程来获得展开的相位图。基于 FFTTIE 的 PU 方法被证明是准确和快速的，但效率仍有提高空间。文献［33- 35］提出将泊松方程用 Volkov 方案求解［36］，该方法将相位图对称扩展，但处理时间也随之增加。 Zuo 等［37- 38］在他们的出色工作中证明，用离散余弦变换（DCT）并且利用硬边界条件（边界外光强为零）可以检索 TIE 真实相位，并且不需要 Volkov 方案。本文将这一研究成果引入到数字全息显微 PU 算法中，通过模拟与实验数据测试，证明基于 DCTTIE 的数字全息 PU 方法与传统的 FFT 算法效率更高。通过将边缘部分反射（RMP）来扩展相位图，表明采用 RMP 方案与 Volkov 方案相比可以更高的速度展开相位图，且不损失精度。同时通过与传统的 DCTLS 全息相位解调相比，表明该方法在处理相位剧烈变化的物体时具有更强的鲁棒性。

2 基本原理

2. 1 TIE 原理简介

TIE 是一种非干涉测量光学相位检索的方法。根据文献［39- 40］，在近菲涅耳衍射区，光波的强度（I）和相位（ϕ）之间的关系可以描述为

式中：k 是波数；∇ 是空间坐标（x，y）上的横向梯度算子；∂I/∂z 是强度的轴向导数，一般来说，在式（1）中选择 z 为 0。Neumann 边界条件通常应用于式（1）的解。则有

式中：Ω 是由 TIE 控制的区域，∂Ω 是边界；g 是一个光滑函数。在此边界条件下，TIE 的解在加性常数之前是唯一的。这里可直接给出绝对相位，无需 PU 处理。如果给出包裹相位图 φW，则可以通过数值方式生成波函数

然后可以计算强度的轴向导数。最后，可以基于 TIE 获得展开的相位图。当 TIE 用于 PU 时，式（1）可以进一步简化。当波函数 U 的振幅设置为常数时，强度 I 对于所有坐标（x，y）都是常数。为方便计算，此处将振幅设置为 1，对于用 TIE 方法进行相位 PU 而言，该设置不影响最终结果。然后式（1）可以写为