高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种两步相移数字全息的设计与实验验证

位文广 韩俊鹤 付玉洲

引用本文:
Citation:

一种两步相移数字全息的设计与实验验证

    作者简介: 位文广(1992-),男,硕士研究生,现主要从事数字全息的重建技术研究.
    通讯作者: 韩俊鹤, junh@henu.edu.cn
  • 中图分类号: O438.1

Design and verification of two-step phase-shifted digital holography

    Corresponding author: HAN Junhe, junh@henu.edu.cn ;
  • CLC number: O438.1

  • 摘要: 为了快速有效地测量微小透明物体表面形貌,设计了一种基于Mach-Zehnder干涉仪的两步相移数字全息系统。该系统使用两个相同的CCD在不同距离同时采集干涉图像,利用光学相移单元在每个CCD记录的干涉图之间形成一个π的相移,然后通过光的空间传递函数及傅里叶变换给出了相位重建的算法。搭建了两步相移干涉光路,以微透镜阵列作为待测物体,验证了该系统的可行性。结果表明,该系统比四步相移法节省一半以上的时间,且能够达到与四步相移法相一致的相位重建结果。该方法对提高相位重建效率具有一定的帮助。
  • Figure 1.  Schematic diagram of two-step phase-shifting digital holography system

    Figure 2.  Simulation of objects to be measured a—amplitude distribution of the simulated object to be measured b—phase distribution of the simulated object to be measured

    Figure 3.  The interference intensity distribution a—D1=61mm and phase shift is 0 b—D2=66mm and phase shift is 0 c—D1=61mm and phase shift is π d—D2=66mm and phase shift is π

    Figure 4.  Phase reconstruction results of simulated interferograms

    Figure 5.  Interference images collected by CCD a, b—residual phase images collected by CCD1 when the phase shift is 0 and π c, d—residual phase images collected by CCD2 when the phase shift is 0 and π e, f—the interferograms collected by CCD1 when the phase shift is 0 and π g, h—the interferograms collected by CCD2 when the phase shift is 0 and π

    Figure 6.  Preliminary phase reconstruction results a—residual phase diagram of the system b—phase distribution of microlens array including residual phase of the system

    Figure 7.  The final result of phase reconstruction a—phase reconstruction result of two-step phase shift b—phase reconstruction result of four-step phase shift

  • [1]

    REN Zh, FAN Z B, LI J Ch, et al. Precision measurement of mini-rotating-angles based on digital holography and its application.Laser Technology, 2012, 36(6):798-801(in Chinese). 
    [2]

    HAO B, SHAN M, ZHONG Z, et al. Common-path interferometer with four simultaneous phase-shifted interferograms using ronchi gra-ting and cube beamsplitter. Optics and Lasers in Engineering, 2013, 51(11):1278-1282. doi: 10.1016/j.optlaseng.2013.05.005
    [3]

    DONG K P, QIAN X F, ZHANG L, et al. Digital holographic microscopy study for cells. Acta Photonica Sinica, 2007, 36(11):2013-2016(in Chinese). 
    [4]

    VESPINI V, COPPOLA S, TODINO M, et al. Forward electrohydrodynamic inkjet printing of optical microlenses on microfluidic devices. Lab on a Chip, 2016, 16(2):326-333. doi: 10.1039/C5LC01386K
    [5]

    CAI H Y, LI G Y, HUANG Zh H. Microsurface profile measurement system based on two-step phase-shifting interference. Laser Technology, 2016, 40(1):20-24(in Chinese). 
    [6]

    ZHANG S Q, ZHOU J Y. A new estimation method for two-step-only quadrature phase-shifting digital holography. Optics Communications, 2015, 335(1):183-188. 
    [7]

    HUANG P S, ZHANG S. Fast three-step phase-shifting algorithm. Applied Optics, 2006, 45(21):5086-5091. doi: 10.1364/AO.45.005086
    [8]

    WANG H J, WANG Zh, ZHAO H, et al. Research on effects of phase error in phase- shifting interferometer.Proceedings of the SPIE, 2007, 6723:67233W. 
    [9]

    CHEN B X, TIAN Y Zh, ZHAO N N, et al. Optimization of two-step phase-shifting digital hologram algorithm and experimental verification. Laser & Optoelectronics Progress, 2015, 52(8):080903(in Chinese). 
    [10]

    DENG L J, YANG Y, SHI B Ch, et al. Two-step phase-shifting digital holography based on extraction of phase shift. Chinese Journal of Lasers, 2014, 41(2):0209014(in Chinese). doi: 10.3788/CJL201441.0209014
    [11]

    MENG X F, CAI L Z, XU X F, et al. Two-step phase-shifting interferometry and its application in image encryption. Optics Le-tters, 2006, 31(10):1414-1416. doi: 10.1364/OL.31.001414
    [12]

    SONG X F, YU M J, WANG H Y, et al. Effect of reference intensity ratio to object on reconstructed image quality in digital holograhpy. Laser Technology, 2014, 38(6):859-862(in Chinese).
    [13]

    LIU J P, POON T C. Two-step-only quadrature phase-shifting digital holography. Optics Letters, 2009, 34(3): 250-252. doi: 10.1364/OL.34.000250
    [14]

    LIU J P, POON T C, JHOU G S, et al. Comparison of two-, three-, and four- exposure quadrature phase- shifting holography. Applied Optics, 2011, 50(16): 2443-2450. doi: 10.1364/AO.50.002443
    [15]

    QIN Y, GONG Q, YANG X Q. A method for accurate phase shift in two-step phase-shifting digital holography. Acta Photonica Sinica, 2011, 40(8): 1282-1286(in Chinese). doi: 10.3788/gzxb20114008.1282
    [16]

    YAMAGUCHI I, ZHANG T. Phase-shifting digital holography. Optics Letters, 1997, 22(16): 1268-1270. doi: 10.1364/OL.22.001268
    [17]

    ABDELSALAM D, YAO G, GAO B L, et al. Single-shot parallel four-step phase shifting using on-axis Fizeau interferometry. A-pplied Optics, 2012, 51(20):4891-4895. 
    [18]

    HAN J H, YAO B L, GAO P, et al. Application of bacteriorhodopsin film for polarization phase-shifting interferometry. Journal of Modern Optics, 2008, 55(14): 2215-2222. doi: 10.1080/09500340802082350
    [19]

    FERRARI J A, GARBUSI E, FRINS E M. Phase modulation by polarization recording in bacteriorhodopsin: Application to phase-shifting interferometry. Optics Letters, 2004, 29(10): 1138-1140. doi: 10.1364/OL.29.001138
    [20]

    ITOH K. Analysis of the phase unwrapping algorithm. Applied Optics, 1982, 21(14): 2470. doi: 10.1364/AO.21.002470
  • [1] 余清婷吕晓旭张庆生刘赣勇 . 基于任意相移量的4步相移数字全息新方法. 激光技术, 2009, 33(4): 409-412. doi: 10.3969/j.issn.1001-3806.2009.04.022
    [2] 张庆生吕晓旭余清婷刘赣勇 . 基于灰度相关的数字全息孔径合成方法. 激光技术, 2009, 33(3): 276-278,282.
    [3] 侯瑞宁 . 希尔伯特变换在同轴数字全息再现中的应用. 激光技术, 2013, 37(3): 362-364. doi: 10.7510/jgjs.issn.1001-3806.2013.03.021
    [4] 裴闯蒋晓瑜王加宗艳桃 . 基于迭代傅里叶变换的3维全息图计算新方法. 激光技术, 2013, 37(3): 347-352. doi: 10.7510/jgjs.issn.1001-3806.2013.03.018
    [5] 魏亚运周昕王梦婷呼有军 . 基于光学扫描全息测量相位物体的相位分布. 激光技术, 2016, 40(6): 820-824. doi: 10.7510/jgjs.issn.1001-3806.2016.06.010
    [6] 姜雨甫周昕王运郭雪帆 . 采用二值化相位光瞳的光学扫描全息方法. 激光技术, 2013, 37(6): 723-727. doi: 10.7510/jgjs.issn.1001-3806.2013.06.004
    [7] 杜源张永安钱晓凡胡振华 . 利用混合编程实现实时数字全息再现系统设计. 激光技术, 2011, 35(6): 841-843,848. doi: 10.3969/j.issn.1001-3806.2011.06.031
    [8] 向东桂进斌刘超郑立婷楼宇丽宋庆和 . 数字全息波前准确重建的实验研究. 激光技术, 2017, 41(3): 406-410. doi: 10.7510/jgjs.issn.1001-3806.2017.03.020
    [9] 陈翠茹王华英赵宝群王学朱巧芬王杰宇王文健雷家良 . 基于UMnet的数字全息相位解包裹. 激光技术, 2023, 47(1): 73-79. doi: 10.7510/jgjs.issn.1001-3806.2023.01.011
    [10] 张中恒李重光李靖哲缪正祥刘营 . 无透镜傅里叶变换数字全息中的相位补偿. 激光技术, 2013, 37(5): 596-600. doi: 10.7510/jgjs.issn.1001-3806.2013.05.008
    [11] 陈希慧焦春妍李俊昌 . 空间载波相移法用于全息CT测量气体温度场. 激光技术, 2006, 30(4): 412-414.
    [12] 秦怡王志鹏李根全 . 部分相干光数字全息消色差相移技术研究. 激光技术, 2012, 36(3): 424-427.
    [13] 桂进斌宋庆和李俊昌楼宇丽 . 彩色数字全息常见波面重建算法的实现与比较. 激光技术, 2015, 39(2): 176-181. doi: 10.7510/jgjs.issn.1001-3806.2015.02.006
    [14] 王华英王大勇谢建军陶世荃 . 数字离轴无透镜傅里叶变换全息重建方法研究. 激光技术, 2007, 31(3): 288-290,294.
    [15] 范美霞张永安李忠芳 . 反射式数字实时全息法对岩石应力场的研究. 激光技术, 2010, 34(4): 502-505. doi: 10.3969/j.issn.1001-3806.2010.04.019
    [16] 张引科王存帅郝劲波 . 改善数字全息再现像质量的背景强度补偿法. 激光技术, 2013, 37(2): 204-206. doi: 10.7510/jgjs.issn.1001-3806.2013.02.016
    [17] 戴欣冉钱晓凡徐天杰 . 同轴菲涅耳全息中提取相位的算法. 激光技术, 2014, 38(2): 172-176. doi: 10.7510/jgjs.issn.1001-3806.2014.02.006
    [18] 盛兆玄孙新利 . 计算全息技术在光学相关识别中的应用. 激光技术, 2008, 32(1): 61-63.
    [19] 浦世亮金其文陈晓锋毛慧吴学成 . 基于二级并行架构的颗粒全息图快速重建方法. 激光技术, 2021, 45(6): 806-810. doi: 10.7510/jgjs.issn.1001-3806.2021.06.022
    [20] 裴瑞景王硕王华英 . 基于深度学习的全息图超分辨重建研究. 激光技术, 2023, 47(4): 485-491. doi: 10.7510/jgjs.issn.1001-3806.2023.04.007
  • 加载中
图(7)
计量
  • 文章访问数:  5190
  • HTML全文浏览量:  3915
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-01
  • 录用日期:  2019-04-16
  • 刊出日期:  2020-01-25

一种两步相移数字全息的设计与实验验证

    通讯作者: 韩俊鹤, junh@henu.edu.cn
    作者简介: 位文广(1992-),男,硕士研究生,现主要从事数字全息的重建技术研究
  • 河南大学 物理与电子学院,开封 475004

摘要: 为了快速有效地测量微小透明物体表面形貌,设计了一种基于Mach-Zehnder干涉仪的两步相移数字全息系统。该系统使用两个相同的CCD在不同距离同时采集干涉图像,利用光学相移单元在每个CCD记录的干涉图之间形成一个π的相移,然后通过光的空间传递函数及傅里叶变换给出了相位重建的算法。搭建了两步相移干涉光路,以微透镜阵列作为待测物体,验证了该系统的可行性。结果表明,该系统比四步相移法节省一半以上的时间,且能够达到与四步相移法相一致的相位重建结果。该方法对提高相位重建效率具有一定的帮助。

English Abstract

    • 相移干涉数字全息术主要利用光电探测元件(如CCD、CMOS等)代替全息记录干板来记录全息图像,然后在计算机上模拟光的干涉过程对数字全息图像进行重建,从而实现全息技术从数据采集到重建全过程的数字化。相移干涉数字全息技术以光的波长为计量单位且具有全域、无损的特性,在物体表面形变测量、细胞分析、无损检测及光学元件缺陷检测等领域得到广泛应用[1-5]。光学相移干涉的基本原理是在两束相干光,即物光和参考光之间引入有序可控的相移,当物光或参考光的相位发生变化时,干涉条纹的强弱分布也随之发生相应的变化;在此过程中用光电探测器如电荷耦合器件(charge coupled device, CCD)逐步采集对应的干涉图,并将其转化为数字信号后按帧存储在计算机内;然后由计算机按照特定的相移干涉的波前重建算法由采集的数字干涉图样解出物光的复振幅分布。根据所需相移次数不同,可将相移法分为两步相移法[6]、三步相移法[7]和四步相移法[8]等。显然,随着相移次数的增加,所采集的干涉图像和耗费的时间也随之增加,这也是2步相移法较之其它相移法优势所在[9-10]。虽然传统的两步相移法与四步相移法相比,节省了近一半时间,但其算法和光路结构也相对复杂。2006年,MENG等人提出了一种两步相移干涉法,只需要两幅干涉图和一个单独的参考光强度分布来进行相位测量[11]。该方法要求参考波的平均强度比物光波最高强度高两倍,且随着物光和参考光相对强度的变化会对全息重建结果产生不同影响[12]。LIU和POON证明了相位成像可以用两个干涉图和参考波的估计而非测量来完成[13],其后又对该算法进行了优化[14],此算法虽然可以很好地去除零级像和孪生像,但计算原理较为复杂而且计算量较大,且使用压电陶瓷器实现相移操作不便[15]

      本文中所述的两步相移干涉数字全息算法,主要利用两个相同的CCD放置于距离像平面后的不同位置,同时采集不同位置的两幅干涉图以缩短采集干涉图像所需时间,然后利用光学相移单元在参考光中引入π的相移后再次采集两幅干涉图,并根据光学传递函数的性质得到相移前后干涉图间的数学关系,从而解出物光的波前分布,这样只需一次相移就可实现相位重建的目的。以微透镜阵列作为待测物体,验证了本文中提出的两步相移算法的可行性。

    • 图 1为两步相移干涉数字全息系统示意图。系统的主体为一个Mach-Zehnder干涉仪,其中,光源为一个波长为632.8nm的He-Ne激光器。He-Ne激光器发出的光束经过透振方向沿竖直方向的偏振器P后入射扩束器(beam expander, BE),经过扩束器BE扩束后的光束由消偏振分光棱镜(non-polarizating beam splitter, NPBS)NPBS1分为参考波和物光波。在参考波光路中,放置一个光学相移单元(phase shifting unit, PSU),用于调节参考光的相位。在物光波光路中,一个微透镜阵列(microlens array,MA)作为待检测物体放置于4倍显微镜物镜(microscopic objective, MO)的前焦平面上。通过待测物体的物光波经一个由显微镜物镜MO和消色差透镜L组成的显微放大系统放大。Im为像平面, 放大后的物光波和参考波合束于消偏振分光棱镜NPBS2。型号相同的两个CCD,即CCD1和CCD2分别放置于距像平面为D1D2的位置,用于采集对应位置的干涉图样。CCD型号为DMK 23G274, 分辨率为1600pixel×1200pixel,像素尺寸为4.4μm。

      Figure 1.  Schematic diagram of two-step phase-shifting digital holography system

      CCD1与像平面(x0, y0)的距离为d1。设Od1为物光波O(x0, y0)在自由空间中传播d1距离后的复振幅分布,可表示为:

      $ {\mathit{\boldsymbol{O}}_{{d_1}}}\left( {{x_1}, {y_1}, {z_1}} \right) = {{\mathscr{F}}^{ - 1}}\left\{ {{\mathscr{F}}\left[ {\mathit{\boldsymbol{O}}\left( {{x_0}, {y_0}} \right)} \right] \cdot {\mathit{\boldsymbol{G}}_{{d_1}}}} \right\} $

      (1)

      式中,${\mathscr{F}}\left\{ {} \right\}$和${{\mathscr{F}}^{ - 1}}\left\{ {} \right\}$分别表示傅里叶变换和逆傅里叶变换。Gd1为自由空间传递函数,表示为:

      $ {\mathit{\boldsymbol{G}}_{{d_1}}}\left( {\xi , \eta , {\mathit{\boldsymbol{d}}_1}} \right) = \exp \left[ {{\rm{i}}k{\mathit{\boldsymbol{d}}_1}\sqrt {1 - {{(\lambda \xi )}^2} - {{(\lambda \eta )}^2}} } \right] $

      (2)

      式中,ζη为频域空间坐标,λ为激光波长。设平面参考波为R,则在d1位置处的干涉强度分布A11可表示为:

      $ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{A}}_{11}}\left( {{x_1}, {y_1}, {\mathit{\boldsymbol{d}}_1}} \right) = }\\ {{{\left| {{\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}} \right|}^2} + |\mathit{\boldsymbol{R}}{|^2} + {\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}{\mathit{\boldsymbol{R}}^*} + {\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}} \cdot \mathit{\boldsymbol{R}}} \end{array} $

      (3)

      在参考光中引入π的相移后,d1处干涉图样的强度分布A12可写为:

      $ \begin{array}{*{20}{c}} {{\mathit{\boldsymbol{A}}_{12}}\left( {{x_1}, {y_1}, {\mathit{\boldsymbol{d}}_1}} \right) = }\\ {{{\left| {{\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}} \right|}^2} + |\mathit{\boldsymbol{R}}{|^2} - {\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}{\mathit{\boldsymbol{R}}^*} - {\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}} \cdot \mathit{\boldsymbol{R}}} \end{array} $

      (4)

      则由(3)式和(4)式可得:

      $ {\mathit{\boldsymbol{A}}_{11}} - {\mathit{\boldsymbol{A}}_{12}} = 2{\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}\mathit{\boldsymbol{R}}_1^* + 2\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}^*\mathit{\boldsymbol{R}} $

      (5)

      同时对(5)式两边进行傅里叶变换可得:

      $ {\mathscr{F}}\left\{ {{\mathit{\boldsymbol{A}}_{11}} - {\mathit{\boldsymbol{A}}_{12}}} \right\} = {\mathscr{F}}\left\{ {2{\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}}{\mathit{\boldsymbol{R}}^*} + 2\mathit{\boldsymbol{O}}_{{\mathit{\boldsymbol{d}}_1}}^*\mathit{\boldsymbol{R}}} \right\} $

      (6)

      为简化计算,假设平面参考波R的复振幅为1。根据(2)式,可知Gd1的复共轭Gd1*G-d1相等,再将(1)式和(2)式代入(6)式可得:

      $ \begin{array}{*{20}{c}} {{\mathscr{F}}\left\{ {{\mathit{\boldsymbol{A}}_{11}} - {\mathit{\boldsymbol{A}}_{12}}} \right\} = }\\ {2{\cal T}\left\{ \mathit{\boldsymbol{O}} \right\} \cdot {\mathit{\boldsymbol{G}}_{{\mathit{\boldsymbol{d}}_\mathit{\boldsymbol{1}}}}} + 2{\cal T}\left\{ {{\mathit{\boldsymbol{O}}^*}} \right\} \cdot {\mathit{\boldsymbol{G}}_{ - {\mathit{\boldsymbol{d}}_1}}}} \end{array} $

      (7)

      同理,通过位于距像平面(x0, y0)为d2的CCD2上所记录的两幅干涉图也可以得到如下式所示关系:

      $ \begin{array}{*{20}{c}} {{\mathscr{F}}\left\{ {{\mathit{\boldsymbol{A}}_{21}} - {\mathit{\boldsymbol{A}}_{22}}} \right\} = }\\ {2{\mathscr{F}}\{ \mathit{\boldsymbol{O}}\} \cdot {G_{{\mathit{\boldsymbol{d}}_2}}} + 2{\mathscr{F}}\left\{ {{\mathit{\boldsymbol{O}}^*}} \right\} \cdot {G_{ - {\mathit{\boldsymbol{d}}_2}}}} \end{array} $

      (8)

      由(7)式×G-d1-(6)式×G-d2可解出待测物体物光波O(x0, y0)的频谱分布, 如下式所示:

      $ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\mathscr{F}}\left\{ \mathit{\boldsymbol{O}} \right\} = \\ \frac{{{\mathscr{F}}\left\{ {{\mathit{\boldsymbol{A}}_{21}} - {\mathit{\boldsymbol{A}}_{22}}} \right\} \cdot {\mathit{\boldsymbol{G}}_{ - {\mathit{\boldsymbol{d}}_1}}} - {\mathscr{F}}\left\{ {{\mathit{\boldsymbol{A}}_{11}} - {\mathit{\boldsymbol{A}}_{12}}} \right\} \cdot {\mathit{\boldsymbol{G}}_{ - {\mathit{\boldsymbol{d}}_2}}}}}{{2\left( {{\mathit{\boldsymbol{G}}_{{\mathit{\boldsymbol{d}}_2} - {\mathit{\boldsymbol{d}}_1}}} - {\mathit{\boldsymbol{G}}_{{\mathit{\boldsymbol{d}}_1} - {\mathit{\boldsymbol{d}}_2}}}} \right)}} \end{array} $

      (9)

      式中, G-d1, G-d2, Gd1-d2Gd2-d1分别为在距离-d1, -d2, d1-d2d2-d1的自由空间传递函数。对(9)式进行逆傅里叶变换, 可以解出待测物体物光波O(x0, y0)的复振幅。这种方法不仅有着和传统四步相移相同的相位重建效果,而且只用了一半的时间。与传统相移算法类似[16-19],通过上述方法解出的复振幅分布仍然是包裹的[20],作者使用最小二乘法对包裹的复振幅分布进行解包裹,就可得到物光波的真实相位分布。

    • 为了验证上述理论的可行性,首先使用MATLAB对其进行模拟验证。图 2为模拟的待测物体的振幅分布和相位分布。图 3I11, I12, I21I22为按照(3)式和(4)式等生成的数字全息干涉强度分布图。其中I11I12D1=61mm时生成的干涉强度分布图,且I11I12之间有π的相移;I21I22D2=66mm时生成的干涉强度分布图,I21I22之间相移同样为π。

      Figure 2.  Simulation of objects to be measured a—amplitude distribution of the simulated object to be measured b—phase distribution of the simulated object to be measured

      Figure 3.  The interference intensity distribution a—D1=61mm and phase shift is 0 b—D2=66mm and phase shift is 0 c—D1=61mm and phase shift is π d—D2=66mm and phase shift is π

      按照上面所述方法对模拟出的干涉强度分布图进行相位重建,可得到图 4。对比图 4图 2b中相位分布可知, 两步相移干涉全息算法确实能够实现相位重建的目的。

      Figure 4.  Phase reconstruction results of simulated interferograms

      为了进一步验证两步相移干涉全息术的有效性,以微透镜阵列作为待测物体重建其相位分布。微透镜阵列为THORLABS公司生产的刻在紫外融石英底板上的熔融石英凸透镜阵列,型号为MLA300-14AR-M。

      在实验过程中,由于光的衍射,干涉图应靠近物体像平面采集,因此,将光学相移单元放置于参考光路,且将两个CCD靠近合束棱镜NPBS2放置。作者将CCD1和CCD2垂直于光路放置,其与物体经显微放大

      系统所成像的距离分别为47mm和52mm。由于实验所用的激光束并非理想的平面波,所用的光学器件也会引入附加残余相位,因此, 需要在不放置待测物体时采集干涉图重建系统的残余相位分布,实验结果如图 5所示。

      Figure 5.  Interference images collected by CCD a, b—residual phase images collected by CCD1 when the phase shift is 0 and π c, d—residual phase images collected by CCD2 when the phase shift is 0 and π e, f—the interferograms collected by CCD1 when the phase shift is 0 and π g, h—the interferograms collected by CCD2 when the phase shift is 0 and π

      图 5a图 5b分别为相移为0和π时CCD1所采集干涉图,图 5c图 5d分别为相移为0和π时CCD2所采集干涉图。把微透镜阵列放置于显微放大系统的物平面,继续按要求采集4幅干涉图,结果如图 5e图 5f图 5g图 5h所示。图 5e图 5f为相移为0和π时CCD1所采集干涉图,图 5g图 5h分别为相移为0和π时CCD2所采集干涉图。按照本文中所述的相位重建算法,经过拟合重建,当D1=47.2mm和D2=52.15mm(即ΔD=D2-D1=4.95mm)时对图 5a~图 5d的相位重建结果最好,如图 6a所示。在相同位置对图 5e~图 5h进行相位重建,重建结果如图 6b所示。图 6a是系统的残余相位图,图 6b是包含系统残余相位时微透镜阵列的相位分布图。

      Figure 6.  Preliminary phase reconstruction results a—residual phase diagram of the system b—phase distribution of microlens array including residual phase of the system

      在包含系统残余相位的微透镜阵列相位分布图中除去系统的残余相位可以得到微透镜阵列的真实相位分布,结果如图 7a所示。为验证相位重建结果的有效性,同时使用四步相移法对微透镜阵列的相位分布进行重建,重建结果如图 7b所示。对比图 7a图 7b可知,两种方法所得相位重建结果相一致,由此可知本文中所用的两步相移算法可以有效的实现微透镜阵列的相位重建。

      Figure 7.  The final result of phase reconstruction a—phase reconstruction result of two-step phase shift b—phase reconstruction result of four-step phase shift

      实验中系统测量误差来源主要有:光学系统误差和相移误差。其中光学系统误差主要是透镜引入的像差可能引起干涉图样的变形,此外透镜对焦的不准确也会造成图像变形,使用高质量的傅里叶变换透镜可以有效减少此类误差。相移误差主要在于相移单元相移的准确度以及实验环境引起的干涉条纹晃动,提高相移单元相移的准确度以及在CCD采集干涉图时使用降低噪声的功能能够减少此类误差。

    • 从理论模拟和实验结果可知,通过将两个相同CCD放置于不同位置同时采集物光波在不同衍射距离的干涉图,同时使用光学相移单元在参考光中引入特定的相移π,并根据光学传递函数的性质可得干涉图之间的函数关系,最终实现相位重建的方法是可行的。这种方法不仅结构相对简单、操作方便,而且比四步相移节约了一半以上时间,可以快速有效地实现透明物体的相位重建,理论模拟和实验结果都表明了这种方法的可靠性和有效性。

参考文献 (20)

目录

    /

    返回文章
    返回