基于位相变更的非相干数字全息自适应成像-周宏强.pdf

《基于位相变更的非相干数字全息自适应成像-周宏强.pdf》由会员分享，可在线阅读，更多相关《基于位相变更的非相干数字全息自适应成像-周宏强.pdf（11页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、基于位相变更的非相干数字全息自适应成像周宏强万玉红满天龙Adaptive imaging by incoherent digital holography based on phase changeZhou Hong-Qiang Wan Yu-Hong Man Tian-Long引用信息Citation: Acta Physica Sinica , 67, 044202 (2018) DOI: 10.7498/aps.20172202在线阅读View online: http:/dx.doi.org/10.7498/aps.20172202当期内容View table of contents:

2、 http:/ you may be interested in激光超衍射加工机理与研究进展The mechanisms and research progress of laser fabrication technologies beyond diffraction limit物理学报.2017, 66(14): 144207 http:/dx.doi.org/10.7498/aps.66.144207一种基于全息术的光学系统闭环像差补偿方法A closed-loop aberration compensating method of optics system based on holo

3、graphy物理学报.2015, 64(2): 024206 http:/dx.doi.org/10.7498/aps.64.024206染料掺杂光子晶体荧光带隙边缘的激射研究Lasing in dye-doped photonic crystals at the edge of fluorescence band gaps物理学报.2013, 62(4): 044208 http:/dx.doi.org/10.7498/aps.62.044208大口径多层介质膜光栅衍射效率测量及其在制作工艺中的应用Diffraction efficiency measurement of large ape

4、rture multilayer dielectric grating and its application in thefabrication process物理学报.2012, 61(17): 174203 http:/dx.doi.org/10.7498/aps.61.174203用Stokes参量法实现数字同轴偏振全息的研究Analysis of digital in-line polarization holography by Stokes parameters物理学报.2012, 61(10): 104202 http:/dx.doi.org/10.7498/aps.61.10

5、4202万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202基于位相变更的非相干数字全息自适应成像周宏强万玉红y满天龙(北京工业大学应用数理学院,北京100124)(2017年10月11日收到; 2017年12月1日收到修改稿)菲涅耳非干数字全息作为一种非扫描的三维成像技术具有其独特的优势,但其成像过程中会受到各种像差的影响,导致成像分辨率、再现像的质量降低.为了解决这一问题,可以结合适当的自适应光学技术对波前像差进行探测和校正.位变更是一种基于两幅具有已知位差的强度图像实现波前探测和像差校正的技术.本文发展了基于位变更的非干数字全息自适应成

6、像技术,不需要引入引导星,利用全息记录过程中的两幅移全息图,实现波前像差的探测.本文给出了所发展技术的数值仿真和实验结果,结合位变更算法求解出系统像差的位分布,将像差的共轭位加载到光瞳上,在全息图记录的同时校正像差,从而提高重建像的质量.关键词:数字全息,自适应成像,位变更,像差校正PACS: 42.30.d, 42.40.Lx, 42.30.Rx, 42.15.Fr DOI: 10.7498/aps.67.201722021引言菲涅耳非干全息(Fresnel incoherent corre-lation holography, FINCH)作为一种非扫描三维成像技术得到了广泛的关注和研究1

7、 5.采用准单色的扩展光源(空间非干光)照明物体,物体上每一物点发出的光通过合适的分光方法,如利用在空间光调制器(spatial light modulator, SLM)加载适当的位掩膜进行衍射分光6,分束后的两束光在合适的平进行自干涉产生点源全息图,所有点源全息图的非干叠加便是整个物体的全息图.全息图重建过程在计算机实现,经过数值衍射传播后得到重建像.为充分利用光源的时间干性和系统的空间带宽积, FINCH系统多采用同轴光路配置,结合移技术7实现全息图的记录,在重建过程中去除零级项和孪生像的干扰. FINCH技术作为一种三维显微成像技术8,经过理论和实验研究已经证明其具有超分辨率成像的特性

8、.结合结构照明或者渐变折射率透镜等方法或元件,实现了FINCH的超分辨率成像9;10.文献11对FINCH成像速度的提高方法进行了研究,通过单通道单次曝光的全息图即可再现出物体的三维信息,加快了记录速度.近年来,在FINCH基础上已发展出基于孔径编码的非干全息技术(coded aperturecorrelation holography, COACH)12;13.向实际应用的成像系统不同程度地会受到系统像差的影响,导致系统成像性能下降,成像指标达不到理论预期.基于FINCH技术的显微成像系统中,全息图记录的过程中也会受到某些像差的影响14;15,造成所记录的全息图含有像差以致于重建像质量不高,

9、分辨率降低.为解决这一问题,Kim16;17提出了基于引导星的非干数字全息自适应光学,该方法先引入一个光轴上的点源作为引导星,记录引导星的多幅移全息图,得到复值全息图.利用系统引导星复值全息图(点源全息图)的复共轭与扩展物体的复值全息图进行卷积运算,从而实现系统像差的补偿.该方法需要记录引导星和物体的多幅移全息图,数据量较大,记录速度慢.并且引导星在实际成像系统中的选取往往是非常国家自然科学基金(批准号: 61575009)资助的课题.通信作者. E-mail: 2018中国物理学会Chinese Physical Society http:/044202-1万方数据物理学报Acta Phy

10、s. Sin. Vol. 67, No. 4 (2018) 044202困难和复杂的,其直径大小、光强度都需要根据实际系统优化.位变更(phase diversity, PD)是一种有效的基于两幅具有固定位差异的强度图像计算波前像差的方法18 21. Klapp和Rosen15在FINCH系统中实现了利用PD算法进行波前探测的自适应光学,但其研究工作只给出了简单的模拟结果.本文进一步发展了基于位变更的非干数字全息自适应光学(PD-FINCH),并研究其成像特性.利用全息记录过程中所得到的具有固定位差的两幅全息图作为初始数据,结合位变更算法求解出系统像差的位分布.在本文所提出的方法中,单独提取出

11、光学像差的具体位分布形式,通过将该位的共轭加载在光学系统中的空间光调制器上,在全息图记录的同时实现像差校正,从而提高重建像的质量.与传统PD算法中采取的数值校正方法比,本文提出的方法具有可灵活调制补偿位、便于分析系统像差模式及像差对成像质量的影响等优势.采用标准分辨率板作为物体进行实验,在系统中存在单一像差和多类别像差的情况下,分别实现了波前像差的探测,经过像差校正后,再现像的质量得到明显的提高.2基本原理菲涅耳非干关联全息(FINCH)是一种非干全息成像技术,其全息图记录过程如图1所示.d1LOSLM CCDd2 z1zhz2P1 P2y0 x0yDxDy xz图1 FINCH记录过程示意图

12、Fig. 1. FINCH recording principle diagram.O是扩展型物体在轴上一物点, O到透镜L的距离为d1,透镜L是成像透镜,焦距f1. P1是起偏器,透振方向与x方向成45角.偏振光束传播到SLM平,在SLM加载焦距为fSLM的透镜的包裹位掩膜.根据SLM的偏振敏感特性,只对水平方向偏振的入射光敏感.当入射光偏振方向为45时,只有入射光的水平分量被调制,而垂直分量没有被调制,且水平分量和垂直分量光强等.这样一束光经过SLM就被分成两束光,两束光波传播经过P2 (检偏器,透振方向与P1平行).至探测器CCD(charge coupled device)平干涉形成点

13、源全息图. z1, z2分别是物点经过透镜和SLM调制分光后所成的像到SLM的距离. zh为从SLM到CCD的记录距离.全息图在CCD平的光强分布可表示为理想的物像与系统点源全息图的卷积,物体的全息图表示为i(x;y) = o(x;y) hI(x;y); (1)其中o(x;y)是物体的理想像的强度分布, 代表二维卷积运算, hI(x;y)表示系统的点扩散函数(pointspread function, PSF),这里的PSF即为系统点源全息图.由于FINCH系统为同轴光路配置,在全息图再现过程中,再现像会受到零级项和孪生像的干扰,因此需要结合移全息技术,去除其对再现像的影响(本文中采用四步移技

14、术). FINCH中的PSF系统所记录的点源全息图具体形式为hI =Q( 1d1)Q(1f1)Q( 1d2)Q( 1zh)+ Q( 1d1)Q(1f1)Q( 1d2) ejaQ( 1fSLM)Q( 1zh)ej i2; (2)这里Q(1/a) = expj (x2 + y2)/ /a, 为光源中心波长; i代表移值, i = 1, 2, 3, 4, i分别为0, /2, , 3 /2. a为系统像差在SLM上的位分布(这里加孔径光阑限定SLM平为光瞳),a由泽尼克多项式前N项的线性组合来描述:a =Nn=1CnZn; (3)其中, Zn为泽尼克多项式的第n项数学表达式,它和光学像差有直接对应的

15、关系, n = 1;2; ;N;Cn表示泽尼克多项式对应第n项的系数,为区别仿真与实验中的参数,仿真中为Cn,实验中为Cen.在实验设计中, FINCH系统为同轴共路的系统,为精确使用SLM补偿像差,系统设计时,考虑光瞳限制在SLM平(图1蓝色箭头处).在SLM上加载透镜位和像差位的同时引入一定的移,记录两幅具有固定位差的全息图,分别记为i(x;y)和id(x;y).基于PD算法的核心思想,044202-2万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202为求得系统的像差,定义评价函数E(o;a)22;23:E(o;a) =x2XiohI2

16、 +x2Xid ohId2;(4)其中hI(x, y), hId(x, y)为成像系统的点扩散函数(PSF),这里分别为FINCH系统的两个点源全息图,分别对应于记录全息图i(x, y)和id(x, y)时的情形.根据Parseval定理和傅里叶变换的性质可知,在频域中,E(O;a) =u2UI O OTF2+u2UId O OTFd2; (5)式中I;Id, O;OTF;OTFd分别为i;id;o;hI, hId的傅里叶变换;评价函数E(o;a)表明通过采集到的两幅具有固定移差的图像强度分布,可以估算出波前的位分布(即像差分布).在目标物体分布函数未知情形下,为了在像差估算过程中使评价函数不

17、受物体分布函数的影响,把函数E(O;a)对O求偏导,在偏导数为0时,函数取到最小值,此时评价函数的表达式可以写为E(o) =u2U|I OTFd Id OTF|2|OTF|2 + |OTFd|2 ; (6)(6)式是一个依赖于系统传递函数的一个函数,它并不依赖于物函数.采用遗传算法(genetic algo-rithm, GA)24,在不断改变在光瞳上像差的分布过程中,对上述评价函数的可能解执行大规模并行搜索,找到评价函数的最小值.当某一种像差分布对应于函数取得最小值时,对应的便是最接近实验系统的像差位.根据前述理论以及图2的计算过程流程,便可得到计算的像差位.模拟出 FINCH系统PSF(点

18、源全息图 ) OTFFT评价函数 E(o)改变光瞳面上像差位相系数 CnPDIf min?YES输出位相 a 两幅全息图结束NO评估确定泽尼克像差类型 Zn a=n=1NCnZnE(o)=uUN |I.OTF-Id.OTF|2|OTF|2-|OTF|2图2算法流程图(FT表示二维傅里叶变换)Fig. 2. Flow chart of algorithm (FT: Fourier transform).通过上述算法得到像差的位函数a后,将a的复共轭加载到SLM上,结合移技术记录四幅全息图,合成复值全息图,经衍射传播的一定距离后,可得到经过像差校正的再现像.另外,本文对计算的像差位和剩余残差用峰谷

19、值(peak-valleyvalue, P-V)和位标准差(standard deviation, )描述和对比校正前后的光瞳位的变化程度.P V = max min; (7)max和min分别表示光瞳位的最大值和最小值;=vuut 1NNi=1(xi )2; (8)N表示有效光瞳平区域内像素数, xi表示各个像素值, i = 1;2;N, 表示所有像素值的平均.3仿真与实验根据图1中的光路结构设计仿真,在FINCH系统中引入彗差,其数学表达式为25Z7 = 8(3r3 2r)sin : (9)根据(3)和(9)式,使用彗差位a = C7 Z7,其系数为C7 = 0:37作为像差.极坐标参数r

20、 = (x2 + y2),非干光源的中心波长为627 nm,物体到透镜的距离d1 = 260 mm,透镜L的焦距为f = 150 mm,透镜到SLM的距044202-3万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202离d2 = 85 mm; SLM上加载的透镜的焦距为fSLM = 4000 mm,大小为10241024像素,每个像素的大小为6.4 m6.4 m. CCD到SLM的距离zh是279 mm.物光经过起偏后形成与水平方向成45的线偏振光,衍射传播通过透镜L后照射到SLM上,经过SLM加载的掩膜的调制分光.最后经过检偏器P2使两束

21、光同一偏振态(与P1透振方向平行)后在CCD上干涉产生全息图,改变SLM上的移间隔,分别记录两幅移全息图作为PD计算波前像差的强度图.模拟结果如图3所示,图3(a)为校正前全息图,图3(c)是四步移全息图的重建像. PD计算后得到的泽尼克彗差系数为C7 = 0:371的像差位.将得到的共轭位加载到SLM上,补偿系统像差.重新记录四步移全息图,合成复值全息图,经过角谱传播后得到再现像.图3(b)为补偿像差后记录的全息图.校正后的再现像如图3(d)所示.模拟过程中引入系统的像差位如图3(e)所示,校正后剩余的像差位如图3(f)所示.abcd1005-5-1015001500100010005005

22、00001005-5-1015001500100010005005000086420-2-4-6-886420-2-4-6-8(f)(e)(d)(c)(b)(a)图3仿真结果(a)记录的全息图(含有彗差); (b)加载共轭位补偿像差后记录的全息图; (c)校正前的全息图重建像; (d)校正后全息图重建像(红框中为细节部分放大); (e)模拟中引入的彗差位(C7 = 0:37); (f)像差校正后的剩余位Fig. 3. Simulated results: (a) Recorded hologram with coma aberration; (b) recorded hologram afte

23、r correction coma;(c) reconstructed image with coma; (d) reconstructed image after correction(enlarge the details in the red block); (e) comaaberration phase (C7 = 0:37); (f) residual phase after correction.044202-4万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202由图3(c)和图3(d)可以看出,经过PD校正后的重建像质量更高,

24、分辨率也提高了.像差位的P-V由校正前(图3(e) 13.06 rad降低到(图3(f) 0.035 rad;位的标准差 由校正前(图3(e) 0.9859 rad降低到(图3(f) 0.0027 rad.为进一步研究所发展的方法的像差校正特性,在系统中引入复杂像差的模式,像散(Z3)、彗差(Z9)、球差(Z12),像差系数为C3 = 0:4, C9 = 0:2,C12 = 0.3,泽尼克多项式数学表达式如下25:Z3 = 2rsin ; (10)Z9 = 8r3 sin3 , (11)Z12 = 10(4r4 3r2)cos2 , (12)即加入像差位a = C3Z3+C9Z9+C12Z12

25、.记录的全息图如图4(a),经过PD计算之后的解为C3 = 0:457, C9 = 0:18, C12 = 0:284.将其共轭位加载到SLM上补偿系统像差后记录全息图如图4(b).全息图校正前和校正后的再现像分别如图4(c)和图4(d)所示.模拟过程中,引入系统的像差如图4(e),校正后剩余的像差位如图4(f).(d)(c)(b)(a)100012000 0400 80015001500100010005005000 0(f)(e)20100-10-2020100-10-2020151050-5-10-15-2020151050-5-10-15-202000图4模拟结果(a)记录的全息图(含

26、有像差); (b)加载共轭位补偿像差后记录的全息图; (c)校正前的全息图重建像; (d)校正后全息图重建像(红框中为细节部分放大); (e)模拟的像差位(像散C3 = 0:4,彗差C9 = 0:2,高阶球差C12 = -0:3); (f)校正后的剩余位Fig.4. Simulationresults: (a)Hologramwithaberrations; (b)hologramafteraberrationscorrection; (c)reconstructedimagewith aberrations; (d) reconstructed image after correction

27、(enlarge the details in the red block); (e) simulated aberrationphase (astigmatism C3 = 0:4, coma aberration C9 = 0:2, spherical aberration C12 = -0:3); (f) the residual phase.044202-5万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202综上模拟结果分析,经过PD校正后,重建像的分辨率得到明显的提高,能看到第6组第6线对(图4(d).像差位的峰谷值P-V由校正前(

28、图4(e) 26.93 rad降低到(图4(f) 2.12 rad;位的标准差(standard deviation, )由校正前(图4(e) 2.230 rad降低到(图4(f) 0.146 rad.根据基本理论,设计实验光路如图5所示.实验系统中,照明光源采用LED,中心波长为627 nm,带宽是16 nm. CL是准直透镜,焦距为75 mm.光束经过准直后起偏器P1后变成线偏振光.将P1的透振方向调节到与水平方向成45角的方向. A为可变光阑.采用标准分辨率板(USAF1951)作为物体O,分辨率板到透镜L的距离大约CCDLEDCL F A O LBSSLMP1 P2图5实验光路配置Fi

29、g. 5. Experimental set-up.cd(d)(c)ab(b)(a)1005-5-101500150010001000500500001005-5-1015001500100010005005000086420-2-4-6-886420-2-4-6-8(f)(e)图6实验结果(a)记录的全息图(含有彗差); (b)补偿像差后记录的全息图; (c)校正前的重建像; (d)校正后的重建像(红框中为细节部分放大); (e)光瞳彗差位(Ce7 = 0:37); (f)校正后像差的剩余位Fig. 6. Experimental results: (a) recorded hologram

30、 with coma aberration; (b) recorded hologram after correction coma;(c) reconstructed image with coma (enlarge the details in the red block); (d) reconstructed image after correction; (e) comaaberration phase (Ce7 = 0:37); (f) residual phase after correction.044202-6万方数据物理学报Acta Phys. Sin. Vol. 67, N

31、o. 4 (2018) 044202cd0.70.60.50.40.30.20.100 10像素值归一化灰度值20 30 40 50 60 70 80 9020100-10-201500150010001000500 5000020100-10-2015001500100010005005000020151050-5-10-15-2020151050-5-10-15-20(d)(c)(f)(g)(e)(b)(a)图7实验结果(a)记录的全息图含有像差; (b)加载共轭位补偿像差后记录的全息图; (c)校正前的全息图重建像; (d)校正后全息图重建像; (e)光瞳像差位(像散Ce3 = -0:3

32、,彗差Ce9 = -0:4,高阶球差Ce12 = 0:2); (f)校正后的剩余位; (g)归一化校正前(c)和校正后(d)条纹对比度曲线(红线(d);蓝线(c)Fig.7. Experimentalresults: (a)Hologramwithaberrations; (b)hologramaftercorrectingaberrations; (c)reconstructedimagewithaberrations; (d) reconstructed image after correction (bottom right: normalized fringe contrast cur

33、ve, the red is after correction);(e) simulated aberration phase (astigmatism Ce3 = -0:3, coma aberration Ce9 = -0:4, spherical aberration Ce12 = 0:2); (f) phasedistribution of the residual aberration; (g) the normalized contrast curve before and after correction (red (d); blue (c).044202-7万方数据物理学报Ac

34、ta Phys. Sin. Vol. 67, No. 4 (2018) 044202为262 mm.透镜L的焦距为f = 150 mm.经过分光棱镜(beam splitter, BS)后照射到空间光调制器(SLM, Holoeye LETO 3.0)上. SLM像素数为1920 1080,像元尺寸6.4 m 6.4 m. SLM是LCOS纯位型的反射式空间光调制器.由于该SLM只对水平偏振的入射光敏感,将P1的透振方向旋转至与水平方向呈45角,此时入射到SLM上的光波的为线偏振,并为与水平方向夹角为45,光强比近似为1 : 1只有水平分量的光波前被调制,垂直分量没有被调制.物光经过透镜和SL

35、M调制分光后在CCD平自干涉产生全息图. CCD与SLM间距为279 mm.其中CCD(Thorlabs DCU224 M,1280 1024)像元尺寸4.65 m4.65 m.实验过程中,在SLM平加载直径为800像素的圆形掩膜作为孔径光阑将系统的光瞳限制在SLM平.并将透镜位掩膜和像差(彗差系数Ce7 = 0:37)加载到SLM上,记录四幅移全息图,组合成复值全息图.经过角谱传播一段距离得到再现像,实验中再现距离为4.5 mm.实验结果如图6所示,其中图6(a)和图6(c)为校正前全息图和重建像.利用记录的前两幅全息图和移值,经过PD算法计算后,得到的彗差系数为Ce7 = 0:36:为补偿

36、系统像差,将得到的共轭位加载到SLM上,重新记录四幅移全息图,合成复值全息图,经过角谱传播到指定距离,便可得到清晰的重建像.校正后的全息图和再现像分别如图6(b)和图6(d)所示.同时可以看到系统的彗差像差位以及补偿后的剩余位如图6(e)和图6(f).利用PD算法求解得到系统的泽尼克系数为Ce7 = 0:36.比较图6(c)和图6(d)中校正前后的重建像,红色方框中细节部分可以清楚地看到,校正后的图像质量得到明显的改善,分辨率也有明显的提高.像差位的峰谷值P-V由校正前(图6(e) 13.06 rad降低到(图6(f) 0.353 rad;位的标准差 由校正前(图6(e) 0.9859 rad

37、降低到(图6(f) 0.0266 rad.在系统中引入复杂像差模式,以研究本文所提出的自适应成像方法对不同像差类型的校正性能.该方法对于泽尼克多项式描述的前12项像差均有校正作用,本文针对荧光显微成像中可能存在对成像质量影响较大的低阶像差(前12项如像散,彗差,球差),模拟和实验都选取了不同的像差模式和不同系数大小进行实验验证,都得到了良好的校正效果.这里展示的是在系统中引入像散(Z3),彗差(Z9),球差(Z12),其中像差系数分别为Ce3 = 0:3, Ce9 = 0:4, Ce12 = 0:2:经过实验并计算求解之后的解为Ce3 = 0:29, Ce9 = 0:48,Ce12 = 0:2

38、8.实验结果如图7所示.其中图7(a)和图7(c)为像差校正前的全息图和重建像,校正后的全息图和再现像如图7(b)和图7(d)所示.同时可以看到系统所引入的像差位以及补偿后的剩余像差位如图7(e)和图6(f).图7(g)给出了图7(c)和图7(d)中沿红线标示方向的再现像图像强度分布.从实验结果上看,校正效果明显,分辨率提高,校正前后的条纹对比度都有明显的提高.像差位的峰谷值P-V由校正前(图7(e) 34.36 rad降低到(图7(f) 3.73 rad;位的标准差 由校正前(图7(e) 2.34 rad降低到(图7(f) 0.30 rad.4结论与讨论本文结合位变更算法发展了一种基于非干数

39、字全息的自适应光学成像技术(PD-FINCH).所提出的方法仅需利用全息记录过程中所得到的具有固定位差的两幅全息图即可实现系统像差的探测.并且光学像差的具体位分布形式被单独提取出来,校正过程采用加载共轭位到SLM进行灵活补偿,校正过程可分离出系统像差分析,校正效果可定量评价.与基于引导星的非干数字全息自适应光学的方法比,该方法记录速度更快,校正补偿方便灵活并且校正效果可以定量地评价.本文通过模拟和实验验证了该方法的可行性,并研究了PD-FINCH系统的自适应成像特性.研究表明,利用所发展的方法, FINCH系统内存在的像差得到了有效校正,成像性能提升,分辨率和校正前后的条纹对比度都有明显的提高

40、. PD-FINCH系统成像性能的提升,将有力推进其在荧光成像系统,三维显微成像领域的应用进展.参考文献1 Rosen J, Brooker G 2007 Opt. Lett. 32 9122 Man T L, Wan Y H, Wu F, Wang D Y 2017 Appl. Opt.56 F913 WanYH,ManTL,WuF,WangDY2014 Opt. Express22 8565044202-8万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 0442024 Shi X, Yuan B, Zhu W F, Liang E J, Liu

41、X M, Ma F Y2014 Proceedings of the Holography, Diractive Optics,and Applications VI 9271 92710W5 Siegel N, Rosen J, Brooker G 2013 Opt. Lett. 38 39226 Katz B, Rosen J, Kelner R 2012 Opt. Express 20 91097 Yamaguchi I, Zhang T 1997 Opt. Lett. 22 12688 Rosen J, Brooker G 2008 Nature Photonics 2 1909 Ka

42、shter Y, Vijayakumar A, Miyamoto Y, Rosen J 2016Opt. Lett. 41 155810 Brooker G, Siegel N, Rosen J, Hashimoto N, KuriharaM, Tanabe A 2013 Opt. Lett. 38 526411 Kelner R, Rosen J 2012 Opt. Lett. 37 372312 Vijayakumar A, Ashter Y, R Kelner, Rosen J 2016 Opt.Express 24 1243013 Vijayakumar A, Rosen J 2017

43、 Opt. Express 25 1388314 Man T L,Wan Y H, Wang D Y 2015 Proceedings ofthe Holography: Advances and Modern Trends IV 950895080415 Klapp I, Rosen J 2013 Proceeding of Imaging and Ap-plied Optics. CTH3C16 Kim M K 2012 Opt. Lett. 37 269417 Kim M K 2013 Appl. Opt. 52 A11718 Gonsalves R A, Chidlaw R 1979

44、Proceedings of the Ap-plications of Digital Image Processing 207 3219 Gonsalves R A 1982 Opt. Engineering 21 82920 Li F 2012 Acta Phys. Sin. 61 230203 (in Chinese) 李斐2012物理学报61 23020321 Li F, Rao C H 2012 Acta Phys. Sin. 61 029502 (in Chi-nese) 李斐,饶长辉2012物理学报61 02950222 Paxman R G, Fienup J R 1988 J

45、. Opt. Soc. Am. A 591423 Paxman R G, Schulz T J, Fienup J R 1992 J. Opt. Soc.Am. A 9 107224 Lei Y J, Zhang S W 2014 Matlab Genetic AlgorithmToolbox and Applications (2nd Ed.) (Xian: Xidian Uni-versity Press), pp143182 (in Chinese) 雷英杰,张善文2014 MATLAB遗传算法工具箱及应用(第二版) (西安:西安电子科技大学出版社第143182页25 Noll R J

46、1976 J. Opt. Soc. Am. 66 207044202-9万方数据物理学报Acta Phys. Sin. Vol. 67, No. 4 (2018) 044202Adaptive imaging by incoherent digital holographybased on phase changeZhou Hong-Qiang Wan Yu-Hong Man Tian-Long(College of Applied Sciences, Beijing University of Technology, Beijing 100124, China)( Received 11 O

47、ctober 2017; revised manuscript received 1 December 2017 )AbstractFresnel incoherent correlation holography (FINCH) is a unique scanning-free three-dimensional imaging techniquewhich enables holograms to be created from incoherent light illumination. However, the image quality is inclined tobe destr

48、oyed by various optical aberrations, in the practical application of microscopic imaging. In order to solve thisproblem, some kinds of adaptive optics are combined with imaging technologies to detect the distorted wavefront andcompensate for the aberrations. Phase diversity is an image-based adaptive optics method where two intensity imageswith a certain phase dierence are used