基于强度传输方程法相位恢复技术实例演示

简单光学

已于 2023-07-07 15:48:06 修改

阅读量4.2k

点赞数 6

分类专栏： # 2.1 相位恢复技术

于 2022-03-15 23:22:42 首次发布

本文链接：https://blog.csdn.net/qq_36584460/article/details/123515029

版权

相位恢复强度传输方程傅里叶变换三维重建光学显微成像

关键词由CSDN通过智能技术生成

2.1 相位恢复技术专栏收录该内容

5 篇文章 15 订阅

订阅专栏

一、引言

CCD相机无法直接记录相位分布。因此，如何通过测量的强度信息来恢复相位信息，即对相位恢复技术的研究尤为重要。

在这里插入图片描述

图 1 相位恢复技术的重要性示意图[1]

目前，基于强度传输方程(Transport of Intensity Equation, TIE)的相位恢复技术已广泛应用于三维重建、光学显微成像、全息技术、衍射光学和光学测量等众多领域，解决了这些领域的许多研究问题。

二、强度传输方程的推导

假设一束单色光沿着z轴自由传播，其复振幅的表示形式可以写为

在这里插入图片描述
其中， $U\left( x,y,z \right)$ 为光场复振幅， $I\left( x,y,z \right)$ 为光场强度， $\varphi \left( x,y,z \right)$ 为光场相位。那么在傍轴近似的条件下，光场的传播满足下面抛物线方程：

在这里插入图片描述

其中 ${{\partial }_{z}}={\partial }/{\partial z}\;$ 代表沿z方向的强度分布的变化量，

在这里插入图片描述
表示梯度算子

表示二维拉普拉斯算子，k为波数，与波长λ之间满足k=2π/λ。分别对 $U\left( x,y,z \right)$ 在(x,y,z)方向上求一维偏导数;

在这里插入图片描述
分别对 $U\left( x,y,z \right)$ 在x,y方向上求二维偏导数;

在这里插入图片描述
先求 $U\left( x,y,z \right)$ 的共轭函数 ${{U}^{*}}\left( x,y,z \right)$

在这里插入图片描述
用 ${{U}^{*}}\left( x,y,z \right)$ 左乘公式(2)得

在这里插入图片描述
用 $U\left( x,y,z \right)$ 左乘公式(4)

在这里插入图片描述
用式(6)减去式(7)得：

化简可得：

在这里插入图片描述
将式(9)进一步化简整理得：

在这里插入图片描述
公式(10)即是推导出的强度传输方程，该方程指出，明确了强度分布以及光轴方向上的强度微分两项重要信息以后，接下来就可以对于二阶偏微分方程进行求解，最终可以得到位置z具体的平面相位分布信息。一般来说，可以采用两个或者两个以上的散焦图像，通过有限差分的方式来对于轴向微分进行预测。

三、强度传输方程的求解

关于强度传输方程的求解，较常用到有傅里叶解法、多重网格方法、格林函数方法以及把扩展成多项式方法等。我们主要介绍的是如何采用傅里叶来进行推导，这种方法可以运用在两种情况下，一种是光波场强度分布均匀的情况，还有一种是非均匀情况。

当强度分布均匀的时候，即I等于常数时，可以通过傅里叶变换有效求解强度传输方程，公式(10)简化为：
在这里插入图片描述
公式(11)是一个泊松方程，可以利用傅里叶变换求解两次泊松方程得到相位[2]:

在这里插入图片描述
当待测光波场强度分布不均匀时，Paganin 等将傅里叶求解方法进行了拓展，其通过将辅助函数引入方程求解过程中，这样再求解公式（11）得到的相位如下：

其中

${{f}_{x}}$ 和 ${{f}_{y}}$ 表示频率域坐标。

公式（12）和式（13）中都有强度轴向微分 ${\partial I}/{\partial z}\;$ 这一项，通常情况下强度轴向微分都是用有限差分法近似得到：

在这里插入图片描述
式中z0为传播距离， $I\left( x,y,{{z}_{0}}+\Delta z \right)$ 和 $I\left( x,y,{{z}_{0}}-\Delta z \right)$ 分别为记录面 ${{z}_{0}}+\Delta z$ 和 ${{z}_{0}}-\Delta z$ 距离处的强度分布， $\Delta z$ 为距离步长。