Experiment research of accurate wavefront reconstruction of digital holography

随着计算机处理速度的提高以及电荷耦合器件(charge-coupled device，CCD)的问世，数字全息已经成为一个十分活跃的研究领域^[1-10]。数字全息技术使用CCD记录数字全息图，通过衍射的数值计算在计算机的虚拟空间重建物光场的波前。常用的两种波前重建算法是基于菲涅耳衍射的一次傅里叶变换法和基于角谱衍射理论的两次傅立叶变换法^[11-12]。在通常的实验验证过程中，在使用这两种方法进行波前重建时，再现距离通常取物平面到全息面的测量距离，并未考虑物光波在不同介质中传播时对测量距离的影响。一般情况下，这种计算是基本准确的，但在一些特殊的应用场合，如要求高分辨成像，或物光波通过其它介质距离较长时，其影响将不能被忽略。如在显微数字全息中，3维洋葱细胞通常放置在液体介质中进行检测，显然，倘若还是采用测量距离值作为再现距离进行波前重建，并不能得到正确的重建结果。

通过对数字全息系统焦深^[13]的分析得知，只要在系统焦深范围内进行数字全息波前重建，都可得到准确的像面结果，当重建距离超出数字全息系统的焦深范围时，需要分析影响系统准确重建距离的因素。本文中通过对物光波通过非单一介质后对重建距离的影响进行了实验研究。研究表明，物光波通过常用分束镜时引起了重建距离的变化，使得重建距离超出了数字全息重建图像的焦深范围。因而，直接采用测量距离作为再现距离进行波前重建通常都不能得到清晰的再现像。本文中将光波在介质中传播时的波长转换为等效波长，进而计算出等效的重建距离，采用等效距离对全息图波前重建公式进行修正，得到准确的计算结果，并与傍轴近似的柯林斯公式计算结果^[14-17]进行比较。结果表明：使用等效距离修正的重建公式，能简便地实现数字全息波面的准确重建计算。研究结果可为显微数字全息及数字全息检测应用提供有益的参考。

数字全息记录的简图如图 1所示。x₀y₀为物平面，激光自左向右被分束镜S₁分为两束光，其中由S₁透射的光经过全反射镜M₁反射，由准直扩束后投向分束镜S₂，经过S₂反射的光波到形成参考光。通过S₁反射的光波经过反射及扩束以后形成照明光投射向物体，从物体表面散射的光波通过半反半透镜S₂到达CCD形成物光。两束光干涉后由CCD记录得到全息图，计算机上进行实时采集存储，然后在计算机上对数字全息图进行重建。重建时可使用下面两种算法进行。

Figure 1.  Experiment configuration for digital holography recording

设到达CCD的物光波和参考光波的复振幅分别为O和R，则记录的光强I_h可以表示为：

物体到CCD的距离为d，重建光波为单位振幅的平面波，使用重建光波照明全息图后的透射光波经过距离为d的衍射后，在满足傍轴近似的条件下，再现光场复振幅分布为：

式中，λ为波长，波数k=2π/λ，${\rm{j}} = \sqrt { - 1} $，xy和x_oy_o分别为全息面坐标系和重建像面坐标系。设CCD的尺寸为L_x×L_y，像素数为M×N，则取样间距Δx=L_x/M，Δy=L_y/N。(2)式可用快速傅里叶变换(fast Fourier transform，FFT)表示为：

式中，I_h(mΔx, nΔy)为全息面I_h(x, y)的离散形式，m和n为取样点对应的序号; Δx_o和Δy_o为离散傅里叶变换后空域的取样间隔，p和q为取样点对应的序号，根据Shannon采样定理，其取样间隔分别为$\Delta {x_{\rm{o}}} = \frac{{\lambda d}}{{M\Delta x}}, \Delta {y_{\rm{o}}} = \frac{{\lambda d}}{{N\Delta y}}$。

当采用照明光波C(x, y)照射全息图I_h(x, y)，出射光波的傅里叶变换可以表示为：

根据角谱衍射理论，在离全息面距离为d的观察面，再现光场复振幅的离散形式为：

式中，${\mathscr{F}^{ - 1}}$表示快速傅里叶逆变换，频率的取样间隔Δf_x=L_x^－1, Δf_y=L_y^－1。

在实际实验系统中，重建距离通常取记录时物体到CCD的测量距离，下面给出常用系统的实验结果。实验物体为美国空军(United States Air Force, USAF)分辨率测试板，采用图 1的数字全息记录系统。实验研究中，照明光波波长λ=632nm，CCD的像元宽度为0.00465mm，总像素数N×N=1024×1024，对应尺寸为L_x×L_y=4.76mm×4.76mm。在通常的数字全息实验系统中，分束镜是必需的光学元件，故以分束镜为例研究介质空间对波面传播的影响，使用两种不同尺寸的分束镜进行实验研究。采用的分束镜光轴方向的厚度为75mm时，测得物体到CCD的距离d=200mm及d=150mm两组数据；当分束镜光轴方向的厚度为25.4mm时，测得物体到CCD的距离d=65mm，分束镜的折射率n=1.5。采用一次FFT和基于角谱衍射理论两次傅里叶变换法对数字全息图进行重建，图 2和图 3中分别给出再现距离为200mm, 150mm和65mm时的重建图像。观察可知，当光波在不同介质中传播时，对于不同记录距离的数字全息图，采用两种同波前重建算法对其进行再现时，均不能得到清晰的再现像。

Figure 2.  Wavefront reconstruction based on fast Fourier transform

Figure 3.  Wavefront reconstruction based on angular spectrum diffraction theory

当光在不同的介质中传播时，其改变的是光波的波长。当物光波在厚度为d_b分束镜中传播，相较于光在真空中传播，其相当于在距离为d_b的长度内填充折射率为n的介质。设光在真空、介质中传播的波长分别为λ和λ₀，传播速率分别为c和v，则在介质中波长为：

对于菲涅耳积分表达式，将(6)式带入(2)式中，令物光波传播到分束镜入射面的复振幅为U_b(x_b, y_b)，出射面复振幅为U(x_{o, b}, y_{o, b})，则光波介质中传播时的菲涅耳衍射积分为：

对于图 1中数字全息记录系统，考虑到光波在折射率为n中传播时，将(7)式与(2)式对比可以知道，在分束镜中光波实际走过的距离为d_b/n，因而实际的距离d_e可以等效为：

即物光波在实际的介质中传播时的等效距离可以通过(8)式计算得到，当采用厚度为75mm和25.4mm的分束镜时，其计算得到对应的等效记录距离分别为d_e=175mm，d_e=125mm和d_e=56mm。采用一次FFT和基于角谱衍射理论的两次傅里叶变换法对数字全息图进行重建，结果如图 4和图 5所示。观察可知，采用等效距离重建数字全息图时，得到了比较清晰的再现像。与第1.3节中的重建结果比较知，由于分束镜厚度对重建距离的影响，使得实际的记录距离的光线的光程发生了变化，因而直接采用测量距离作为再现距离进行重建，并不能获得准确的重建像。

Figure 4.  Wavefront reconstruction based on fast Fourier transform

Figure 5.  Wavefront reconstruction based on angular spectrum diffraction theory

对于一次傅里叶变换系统的焦深范围，参考文献[13]中给出了详细的推导证明。根据文献中的研究，数字全息重建图像的系统焦深定义如下：

式中，z_a和z_i分别是实际的再现距离和理想的成像距离。不难看出，随着理想像重建距离z_i的增加，重建图像的系统焦深会增加；而随着全息图宽度的增加，系统焦深会随之下降。根据第2节中的结论，当再现距离为d_e=175mm，L取CCD的面尺寸，代入(9)式，|z_a－z_i|=1.3mm，即重建时实际能成清晰像的系统焦深范围为(175±1.3)mm，显然，若采用200mm进行数字全息波前重建，已经超出了焦深范围，因而不能得到清晰的再现像。

对于采用等效距离计算得到的重建图像, 如图 6所示。在一系列对图像清晰度的评价方法中，有一种方法是通过图像的高频与低频部分的比值来衡量重建像的清晰度，比值的最大值为图像最清晰处^[18]。采用等效距离修改正后的公式对图像进行重建，对重建图像做傅里叶变换，将零级视为低频，其余部分视为高频，在系统的焦深范围内，采用自动聚焦的方式，对于不同记录距离的重建图像，其高频与低频的比值曲线如图 7所示。由分析可知，重建图像在175.3mm处有最大值，即此处成像是最清晰的位置，这与之前计算得到等效距离的值是基本一致的，因而采用等效距离修正后的公式能较好地成像。

Figure 6.  Wavefront reconstruction based on equivalent distance

Figure 7.  Ratio of high frequency and low frequency

物光波通过实际的光学系统时，通常采用柯林斯公式进行研究。光波通过介质时，轴对称傍轴光学系统可以由2×2的矩阵$\left[ {{array}{*{20}{c}} A&B\\ C&D {array}} \right]$描述^[14]。设入射平面为全息面，出射平面为观察面，坐标分别由xy和x_oy_o表示。根据柯林斯建立的入射平面光波场I_h(x, y)与出射平面光波场U(x_o, y_o)的关系^[15]，观察面复振幅的离散形式可以表示为：

光传输矩阵中元素ABCD的值分别是A=1，$ - B = d - {d_{\rm{b}}} + \frac{{{d_{\rm{b}}}}}{n}$，C=0，D=1，L_axis为ABCD光学系统的轴上光程，N为CCD的像素数，L为CCD面尺寸，L_o为出射面光波场的空域宽度(即出射平面的尺寸)，L_o可以表示为：

为便于讨论，选取记录距离d=200mm时的数字全息图进行重建, 图像如图 8所示。采用柯林斯公式重建像与采用等效距离修正后的公式进行波前重建像，取第30行的复振幅的值做差, 如图 9所示。从图 9中可知，光波通过不同尺寸的分束镜时，采用柯林斯公式和采用基于等效距离对数字全息图重建时有基本一致的结果。

Figure 8.  Reconstruction based on Collins formula

Figure 9.  Complex amplitude subtract of two reconstruction images in different method

当光波在非单一介质中传播时，给出了将记录距离等效为介质中传播距离的方法。结果表明，计算光波在介质中传播的等效距离，并采用等效距离对数字全息图进行波前重建的结果与采用柯林斯公式重建的结果是一致的。此方法相较于采用柯林斯公式对数字全息图进行波前重建更简单方便，期望此方法能对数字全息图的波前重建提供有益的参考。