光学学报, 2018, 38 (9): 0905002, 网络出版: 2019-05-09  

基于鬼衍射方法的相位定量重建 下载: 754次

Quantitative Phase Retrieval Based on Ghost Diffraction Method
作者单位
1 华南理工大学电子与信息学院, 广东 广州 510641
2 华南理工大学物理系, 广东 广州 510641
摘要
在鬼成像装置的参考光路上安放特定的透射屏,提出了一种新的鬼衍射方案来实现物体相位的定量重建,并给出了对应的理论解释与实验仿真。仿真结果与理论基本吻合,验证了实验方案的正确性与可行性。
Abstract
A novel ghost diffraction scheme for the quantitative phase retrieval of objects is proposed, in which a special transmission screen is inserted in the reference light path of the ghost imaging setup. The corresponding theoretical analysis and the experimental simulation are shown. The simulation results are basically consistent with the theoretical ones, which confirms the correctness and the feasibility of this proposed experimental scheme.

1 引言

鬼成像是一种新型的非局域成像技术[1-7],可以通过非局域的光场强度关联测量,在不包含物体信息的光路中恢复出物体的像,且在待测物体的光路上只需安放一个单像素的探测器,故该技术引起了研究者们极大的研究兴趣。鬼成像技术已在显微成像[8]、**侦查[9-13]、对地观测[14]等领域表现出了巨大的应用潜力。鬼成像技术主要用于获取物体的透射分布,近年来也有致力于利用鬼成像技术恢复物体相位的研究[15-21]。Borghi等[15]通过相干干涉实现了对物体相位的恢复。Ying等[16-17]通过两步法同样实现了对物体相位的恢复。

物体的相位通常需要通过计算量较大的数值迭代或较复杂的光路干涉测量得到,本文改进了常规的鬼衍射方案来实现对物体相位的直接恢复,并给出了相应的理论解释与实验仿真。通过在参考光路上放置4个不同的参考屏幕以获得不同的鬼衍射测量结果,对得到的参考光路与测试光路之间的4个关联结果进行处理,便可以定量恢复出未知物体的相对相位分布。仿真结果与理论基本吻合,验证了实验方案的正确性与可行性。

2 理论解释

图(1)所示为改进的鬼衍射方案的光路示意图,源光场通过分束器(BS)变成了两束光,探测光路包含了一个未知物体t(y),z1z0分别为光源到物体的距离及物体到测试光路探测器Dt的距离。参考光路包含了一个透射率可调的参考屏幕r(u),z2z3分别为光源到参考屏幕的距离及参考屏幕到参考光路探测器Dr的距离。xyuxtxr分别为光源平面、物体平面、参考屏幕平面、测试光路探测器平

图 1. 改进的鬼衍射方案的示意图

Fig. 1. Schematic of improved ghost imaging scheme

下载图片 查看所有图片

面与参考光路探测器平面的坐标。

描述鬼成像图像的强度涨落的关联公式[2]

G(xt,xr)=<δI(xtI(xr)>=<It(xt)Ir(xr)>-<It(xt)><Ir(xr)>,(1)

式中δIt,r=It,r-<It,r>为强度的涨落。未知物体的鬼成像或鬼衍射的图像样式为

I(xr)=DG(xt,xr)dxt,(2)

式中D为探测光路探测器的面积。

可以直接得到

G(xt,xr)=1λ4z0z1z2z3dx1dx2dyduΓ(x1,x2)hz1(x1,y)t(y)hz0(y,xt)hz3(x2,u)r(u)hz4(u,xr)*2,(3)

式中λ为所使用光场的波长,x1x2均为光源面上的位置坐标,Γ(x1,x2)为光源的一阶相干函数,hz(x,y)=exp λz(x-y)2为菲涅耳衍射公式的积分核。

若光源是完全非相干的,则Γ(x1,x2)=Is(x1)δ(x1-x2),其中Is(x1)为光源的空间强度分布,在这种条件下,若z3=z0=l,z1=z2=d,则有

G(xt,xr)=1λ4z0z1z2z3Is(x)t(y)r*(u)expλdy2-2xy-u2-2xu×expλdy2-2yxt-u2-2uxrdxdydu2(4)

若光源足够大而且均匀分布,则Is(x)=I0(I0为均匀光场强度),且其关于x的积分将会得到δ(y-u)项,故最终的结果为

G(xt,xr)=I0λ4l2d22t(y)r*(y)exp-i2πy(xr-xt)λldy2(5)

假设在探测光路坐标xt=0处放置一个点探测器,得到鬼衍射的图像为

I(xr)=I0λ4l2d22TxrλlRxrλl2,(6)

式中T(k)与R(k)分别为t(y)跟r(u)的傅里叶变换,⊗表示卷积。至此得到了图1所示实验框架中最终的鬼成像关联结果。对于一个完全非相干的均匀分布的热光源,当z3=z0=l,z1=z2=d时,鬼衍射图像I(xr)是t(u) ×r(u)的傅里叶变换,即I(xr)∝ TxrλlRxrλl2。若未放置参考屏幕r(u),即r(u)=1,可以得到熟知的鬼衍射图像结果[2]:I(xr)∝ Txrλl2

选取四个完全不同的参考屏幕[22]:

r1(u)=1+exp(-2iπqu)r2(u)=1+iexp(-2iπqu)r3(u)=1-exp(-2iπqu)r4(u)=1-iexp(-2iπqu),(7)

式中q为反映参考屏的空间调制频率的一个参数。

得到的鬼衍射图像分别为

I1(xr)=I0λ4l2d22Txrλl+Txrλl-q2I2(xr)=I0λ4l2d22Txrλl-iTxrλl-q2I3(xr)=I0λ4l2d22Txrλl-Txrλl-q2I4(xr)=I0λ4l2d22Txrλl+iTxrλl-q2,(8)

对这四个鬼衍射图像进行数据处理可以得到

I1(xr)-I3(xr)-i I2(xr)-I4(xr)=

4 I0λ2l2d22TxrλlT*xrλl-q。(9)

S=I1(xr)-I3(xr)-i I2(xr)-I4(xr),故

ΦSxrλl=ΦTxrλl-ΦTxrλl-q,(10)

式中ΦS(k)为Sk点处的相位,ΦT(k)为T(k)在k点处的相位。不失一般性,令ΦT(0)=0,则可以通过(10)式得到ΦT(±q),ΦT(±2q),…的相位值。

至此得到了幅值|T(k)|与抽样相位ΦT(k),可以恢复出物体t(y)的相位值。不过由于ΦT(0)=0,通过这种方式得到的相位与真实相位之间存在一个常数误差,且这个误差值大小正是真实的ΦT(0)值。有意义的测量结果通常只与相位的相对值分布有关,故这个常数的相位偏差不影响实际的相位复原效果。

3 数值模拟

前面的理论均基于光源完全非相干的假设,但是自然界中完全非相干的光源并不存在。为了明确光源的相干性对相位恢复的影响,在仿真实验中,通过高斯谢尔模型光源模拟光源的部分相干性。高斯谢尔模型光源的一阶相干函数[23-25]

<E(x1)E*(x2)>=1bexp -x12+x22a2-(x1-x2)2b2,(11)

式中a为光源的横向尺寸,b为光源的横向相干长度。

将(11)式代入(3)式,得到鬼成像I(xr)为

I(xr)=π2λ8d4l4|1ba1c1-b12exp-π2λ2d2(a1c1-b12)(a1u2+c1y2+2b1uy)t(y)r*(u)·×expλ1l+1d(u2-y2)expλl(-2uxr)dydu|2,(12)

式中参数a1=λl+1a2+1b2,b1=-1b2,c1=-λl+1a2+1b2。当a趋于无穷大而b=0,即当光源完全非相干时,有

I(xr)=π2λ8d4l4dydua2exp-π2a22λ2d2(u-y)2t(y)r*(u)expλ1l+1d(u2-y2)×expλl(-2uxr)2=π2λ8d4l4dut(u)r*(u)expλl(-2uxr)2(13)

可以看出,(13)式与(6)式吻合,验证了理论的正确性。

实验仿真使用一个有固定相位差的双缝物体和一个具有高斯相位分布的透射板。双缝物体的数学表达式为

t(y)=al×recty+w/2w'expiθ1+ar×recty-w/2w'expiθ2,(14)

式中w为双缝的缝宽,w'为双缝的缝间距,alar为双缝的透射率,rect(·)为矩形函数,θ1θ2分别为两条缝所对应的相位延迟。在实验方案中,al=0.5,ar=1,w=80 μm,w'=40 μm,θ1=0.5π,θ2=0.2π。

具有高斯相位分布的透射板的数学表达式为

t(y)=rectycexp-i0.4y2d12,(15)

式中c为透射板的长度,d1为相位分布的宽度参数,且c=160 μm,d1=40 μm。

为了了解光源的横向相干长度b对物体相位恢复的影响,选取了8组不同的数据并对其进行了仿真模拟,取l=20 mm,d=200 mm,双缝物体与高斯相位物体的相位恢复结果如图2、3所示,对应的恢复相位与原始相位的差值如图4所示。这里主要研究的是光源的横向相干长度对相位的影响,故默认光源的横向尺寸a无穷大。

图 2. 不同b下的双缝物体相位的恢复结果。(a) b=0 μm; (b) b=2 μm; (c) b=4 μm; (d) b=6 μm; (e) b=8 μm; (f) b=12 μm; (g) b=16 μm; (h) b=20 μm

Fig. 2. Phase retrieval results of double slits under different b. (a) b=0 μm; (b) b=2 μm; (c) b=4 μm; (d) b=6 μm; (e) b=8 μm; (f) b=12 μm; (g) b=16 μm; (h) b=20 μm

下载图片 查看所有图片

图 3. 不同b下高斯相位物体的相位恢复结果。(a) b=0 μm; (b) b=2 μm; (c) b=4 μm; (d) b=6 μm; (e) b=8 μm; (f) b=12 μm; (g) b=16 μm; (h) b=20 μm

Fig. 3. Phase revival results of Gaussian phase objects under different b. (a) b=0 μm; (b) b=2 μm; (c) b=4 μm; (d) b=6 μm; (e) b=8 μm; (f) b=12 μm; (g) b=16 μm; (h) b=20 μm

下载图片 查看所有图片

图 4. 恢复相位与原始相位的差值。(a)双缝物体;(b)高斯相位透射板

Fig. 4. Difference between reconstructed phase and original phase. (a) Double slits; (b) Gaussian phase transmission plate

下载图片 查看所有图片

图2(a)、图3(a)及图(4)可以看出,b=0即光源完全非相干时,在幅值不为零处成功地得到了与原始物体一致的相对相位分布,恢复出来的物体相位与真实物体相位的差值是一个恒定的常数,这进一步验证了理论的正确性。从图2~4可以看出,当b的取值大于零且小于4 μm时,恢复相位与原始相位的差值基本也是一个恒定的常数,只有边缘处极个别点不符合要求。随着b的取值增大,恢复相位与原始相位之间的差值已经不再是一个恒定的常数,即恢复出来的相位偏离了真实相位,相位恢复结果越来越差。

在具体实验中,可能会存在参考屏幕不理想、参考屏幕位置与待测物体不对应重合(z3z0)等问题。通过具体计算,分别分析了这两种不理想状态对最终相位重建结果的影响。在参考屏幕不理想即混有噪声信号的情况下,参考屏rs(u)=r(u)+ au×n(u),其中au为噪声系数,n(u)为在[0,1]均匀随机分布的变量函数,对au选取了3组数值进行数值模拟。当参考屏幕与待测物体不对应重合(z3z0)时,也就是z2=z1=d,z0=l,z3=l-Δ,其中Δ为距离误差,对Δ选取了4组数值进行仿真模拟,结果如图5所示。

图 5. 双缝物体在不同条件下的相位恢复结果。(a)不同au;(b)不同Δ

Fig. 5. Phase revival results of double slits under different conditions. (a) Different au ; (b) different Δ

下载图片 查看所有图片

图5可以看出,当噪声系数(距离误差)较小时,能较好地恢复物体的相位,但当噪声系数(距离误差)增大到一定范围时,相位恢复效果越来越差。

故具体的实验操作最好通过单光路计算鬼成像[26-28]的方法来实现,通过空间光调制器产生赝热光源,并通过计算准确得到四种不同参考屏幕对应的鬼衍射结果,这样可以避免双光路鬼成像过程中可能出现的参考屏幕不理想、参考屏幕位置与待测物体不对应重合(z3z0)等问题。

4 结论

为了解决鬼衍射技术中的相位恢复问题,提出了一种改进的鬼衍射实验方案。通过在参考光路依次放入四个不同的参考屏幕,得到了四个鬼衍射图像,对其进行简单的计算处理便可以得到物体的相对相位分布。在部分相干的高斯谢尔光源模型下,通过数值计算,证实了理论方案在相位恢复方面的可行性。光源的相干性较差时,鬼衍射方案仍能得到较好的相位恢复结果,可以正确重构物体的相对相位分布。

参考文献

[1] Bennink R S, Bentley S J, Boyd R W. “Two-photon” coincidence imaging with a classical source[J]. Physical Review Letters, 2002, 89(11): 113601.

[2] Cheng J, Han S S. Incoherent coincidence imaging and its applicability in X-ray diffraction[J]. Physical Review Letters, 2004, 92(9): 093903.

[3] Zhang D, Zhai Y H, Wu L A, et al. Correlated two-photon imaging with true thermal light[J]. Optics Letters, 2005, 30(18): 2354-2356.

[4] Scarcelli G, Vincenzo B, Shih Y. Can two-photon correlation of chaotic light be considered as correlation of intensity fluctuations[J]. Physical Review Letters, 2006, 96(6): 063602.

[5] Pelliccia D, Rack A, Scheel M, et al. Experimental X-ray ghost imaging[J]. Physical Review Letters, 2016, 117(11): 119902.

[6] Liu X F, Chen X H, Yao X R, et al. Lensless ghost imaging with sunlight[J]. Optics Letters, 2014, 39(8): 2314-2317.

[7] Padgett M J, Boyd R W. An introduction to ghost imaging: Quantum and classical[J]. Philosophical Transactions of the Royal Society A, 2017, 375(2099): 20160223.

[8] Shen Q, Bai Y F, Shi X H, et al. Ghost microscope imaging system from the perspective of coherent-mode representation[J]. Laser Physics Letters, 2018, 15(3): 035207.

[9] Magaña-Loaiza O S, Howland G A, Malik M, et al. . Compressive object tracking using entangled photons[J]. Applied Physics Letters, 2013, 102(23): 231104.

[10] Zhao C, Gong W, Chen M, et al. Ghost imaging lidar via sparsity constraints[J]. Applied Physics Letters, 2012, 101(14): 141123.

[11] Meyers R E, Deacon K S, Shih Y. Turbulence-free ghost imaging[J]. Applied Physics Letters, 2011, 98(11): 111115.

[12] 梅笑冬, 龚文林, 严毅, 等. 可预置强度关联激光三维成像雷达实验研究[J]. 中国激光, 2016, 43(7): 0710003.

    Mei X D, Gong W L, Yan Y, et al. Experimental research on prebuilt three-dimensional imaging lidar[J]. Chinese Journal of Lasers, 2016, 43(7): 0710003.

[13] 王靖茹, 山卓玉, 张勇. 构造雷达关联成像随机辐射场的方法分析[J]. 光学学报, 2017, 37(8): 0811004.

    Wang J R, Shan Z Y, Zhang Y. Methodology analysis on stochastic radiation field of radar correlated imaging[J]. Acta Optica Sinica, 2017, 37(8): 0811004.

[14] 韩申生, 龚文林, 陈明亮, 等. 基于稀疏和冗余表面的鬼成像雷达研究进展[J]. 红外与激光工程, 2015, 44(9): 2548-2555.

    Han S S, Gong W L, Chen M L, et al. Research process of GISC lidar[J]. Infrared and Laser Engineering, 2015, 44(9): 2548-2555.

[15] Borghi R, Gori F, Santarsiero M. Phase and amplitude retrieval in ghost diffraction from field-correlation measurements[J]. Physical Review Letters, 2006, 96(18): 183901.

[16] Ying G R, Wei Q, Shen X, et al. A two-step phase-retrieval method in Fourier-transform ghost imaging[J]. Optics Communications, 2008, 281(20): 5130-5132.

[17] Yu H, Lu R H, Han S S, et al. Fourier-transform ghost imaging with hard X rays[J]. Physical Review Letters, 2016, 117(11): 113901.

[18] Zhang M H, Wei Q, Shen X, et al. Lensless Fourier-transform ghost imaging with classical incoherent light[J]. Physical Review A, 2007, 75(2): 021803.

[19] Gong W L, Han S S. Phase-retrieval ghost imaging of complex-valued objects[J]. Physical Review A, 2010, 82(2): 023828.

[20] Zhang M H, Wei Q, Shen X, et al. Sub-wavelength Fourier-transform imaging of a pure-phase object with thermal light[J]. Physics Letters A, 2007, 366(6): 569-574.

[21] Shirai T, Setala T, Friberg A T. Ghost imaging of phase objects with classical incoherent light[J]. Physical Review A, 2011, 84(4): 041801.

[22] Kashter Y, Vijayakumar A, Miyamoto Y, et al. Enhanced super resolution using Fresnel incoherent correlation holography with structured illumination[J]. Optics Letters, 2016, 41(7): 1558-1561.

[23] Lin J, Cheng J. Lensless ghost diffraction with partially coherent sources: Effects of the source size, transverse coherence, detector size and defocusing length[J]. Chinese Physics Letters, 2011, 28(9): 094203.

[24] Foley J T, Zubairy M S. The directionality of Gaussian Schell-model beams[J]. Optics Communications, 1978, 26(3): 297-300.

[25] Cai Y J, Zhu S Y. Ghost imaging with incoherent and partially coherent light radiation[J]. Physical Review E, 2005, 71(5): 056607.

[26] Bromberg Y, Katz O, Silberberg Y. Ghost imaging with a single detector[J]. Physical Review A, 2009, 79(5): 053840.

[27] Shapiro H J. Computational ghost imaging[J]. Physical Review A, 2008, 78(6): 061802.

[28] Kuhn S, Hartmann S, Elsäβer W. Photo-statistics-based classical ghost imaging with one single detector[J]. Optics Letters, 2016, 41(12): 2863-2866.

陈紫艳, 徐海红, 程静. 基于鬼衍射方法的相位定量重建[J]. 光学学报, 2018, 38(9): 0905002. Ziyan Chen, Haihong Xu, Jing Cheng. Quantitative Phase Retrieval Based on Ghost Diffraction Method[J]. Acta Optica Sinica, 2018, 38(9): 0905002.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!