CN105158921B - 一种基于互补随机抽样的无透镜衍射成像方法 - Google Patents
一种基于互补随机抽样的无透镜衍射成像方法 Download PDFInfo
- Publication number
- CN105158921B CN105158921B CN201510698056.3A CN201510698056A CN105158921B CN 105158921 B CN105158921 B CN 105158921B CN 201510698056 A CN201510698056 A CN 201510698056A CN 105158921 B CN105158921 B CN 105158921B
- Authority
- CN
- China
- Prior art keywords
- sampling
- complex amplitude
- complementary
- amplitude
- screen
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005070 sampling Methods 0.000 title claims abstract description 92
- 230000000295 complement effect Effects 0.000 title claims abstract description 29
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 49
- 238000004422 calculation algorithm Methods 0.000 claims description 24
- 238000009826 distribution Methods 0.000 claims description 17
- 238000002834 transmittance Methods 0.000 claims description 5
- 230000000149 penetrating effect Effects 0.000 claims description 2
- 238000011084 recovery Methods 0.000 abstract description 11
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000007796 conventional method Methods 0.000 abstract 1
- 230000001427 coherent effect Effects 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000010363 phase shift Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 239000002585 base Substances 0.000 description 1
- 238000012984 biological imaging Methods 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000002425 crystallisation Methods 0.000 description 1
- 230000008025 crystallization Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- CPBQJMYROZQQJC-UHFFFAOYSA-N helium neon Chemical compound [He].[Ne] CPBQJMYROZQQJC-UHFFFAOYSA-N 0.000 description 1
- 238000005210 holographic interferometry Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/42—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
- G02B27/4233—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect having a diffractive element [DOE] contributing to a non-imaging application
- G02B27/4238—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect having a diffractive element [DOE] contributing to a non-imaging application in optical recording or readout devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Holo Graphy (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种基于互补随机抽样的无透镜衍射成像方法,该方法是引入一组具有互补特性的二值随机抽样屏;将这组互补抽样屏依次置于被测物体和记录平面之间,并依次记录透过抽样屏的物波在记录平面处的强度图样;从所记录的强度图样中恢复出被测物波的复振幅分布,利用所恢复出的复振幅分布再经过数字衍射运算,即得到被测物体在空间任意位置的衍射成像。该方法通过在物体和记录面之间***互补随机抽样屏,大大提高了相位恢复的精确度和迭代效率;同时,可以有效消除传统方法中由于二值随机抽样引起的波前抽样损失。另外,由于本方法中迭代运算过程仅应用于抽样面和记录面之间,可适用于更一般的复值物体,并且不需要对物体进行特殊的限制。
Description
技术领域
本发明涉及一种无透镜衍射成像方法,属于衍射成像技术领域。
背景技术
无透镜相干衍射成像(CDI)(参见文献1)是一种利用物波的衍射强度图样计算物波波前复振幅分布的技术,该技术在自适应光学、光学精密测量、高分辨显微成像等领域有重要实际应用。CDI技术的关键问题是怎样从物波的衍射强度图样中恢复出物波波前的复振幅信息(包括振幅信息和相位信息,其中主要是相位信息的恢复)。为解决波前复振幅的恢复问题,人们已经提出了很多方法。其中,全息干涉测量技术是最常用的一类方法,比如同轴或离轴全息(参见文献2-文献10),但是全息测量需要引入合适的参考光或者使用相移算法。另一种是使用迭代算法(参见文献11-文献30),这种方法的好处是不需要引入参考光,但待测波前的复振幅信息需要从其衍射强度图样中通过一种迭代运算恢复出来。传统的迭代算法通常需要大量的迭代次数,并且存在迭代结果的非唯一性。为了提高迭代效率和波前复振幅恢复的准确性,人们也已经提出了一些改进的方法,如采用超抽样法(参见文献14)、采用多幅衍射图样(参见文献15-文献22)、采用结构光照明技术(参见文献23和文献24)、采用相位板波前调制技术等(参见文献25-文献30)。但其中还有很多问题亟待解决。例如,在已有得方法中,为了满足迭代收敛条件,需要采用特殊支架,物体横向尺寸必须远小于图像传感器有效记录平面的大小等。
最近,Horisaki等人提出了一种基于二值随机抽样的CDI方法(参见文献31和文献32)。该方法通过理论分析和计算机模拟证明,通过在被测物体与图像传感器之间***一幅二值随机抽样(RS)屏可以解决波前复振幅恢复中存在的非唯一性和非收敛性问题。为简单起见,在以下叙述中将这种方法称为RS-CDI法。RS-CDI法较好地克服了传统CDI方法中对物体尺寸的限制。但由于在物体和记录面之间引入了一个二值RS屏,该方法只能得到物波波前复振幅的一组稀疏抽样数据。尽管用压缩感知算法(参见文献33)也可以基于这种稀疏抽样数据实现衍射成像,但对待测物物波的冗余性提出了较高的要求。
上述所述及的文献是指:
文献1:J.Miao,T.Ishikawa,I.K.Robinson,and M.M.Murnane,"Beyondcrystallography:Diffractive imaging using coherent x-ray light sources,"Science 348,530-535(2015)
文献2:A.Rosenhahn,R.Barth,F.Staier,T.Simpson,S.Mittler,S.Eisebitt,andM.Grunze,"Digital in-line soft x-ray holography with element contrast,"J.Opt.Soc.Am.A 25,416-422(2008).
文献3:M.Guizar-Sicairos and J.R.Fienup,"Holography with extendedreference by autocorrelation linear differential operation,"Opt.Express 15,17592-17612(2007).
文献4:S.Marchesini,S.Boutet,A.E.Sakdinawat,M.J.Bogan,S.Bajt,A.Barty,H.N.Chapman,M.Frank,S.P.Hau-Riege,A.Szoeke,C.Cui,and D.A.Shapiro,"Massivelyparallel X-ray holography,"Nat.Photonics 2,560-563(2008).
文献5:M.Guizar-Sicairos and J.R.Fienup,"Direct image reconstructionfrom a Fourier intensity pattern using HERALDO,"Opt.Lett.33,2668-2670(2008).
文献6:R.L.Sandberg,D.A.Raymondson,C.La-o-vorakiat,A.Paul,K.S.Raines,J.Miao,M.M.Murnane,H.C.Kapteyn,and W.F.Schlotter,"Tabletop soft-x-ray Fouriertransform holography with 50nm resolution,"Opt.Lett.34,1618-1620(2009).
文献7:J.Geilhufe,C.Tieg,B.Pfau,C.M.Günther,E.Guehrs,S.Schaffert,andS.Eisebitt,“Extracting depth information of 3-dimensional structures from asingle-view X-ray Fouriertransform hologram,"Opt.Express 22,24959(2014).
文献8:X.F.Xu,L.Z.Cai,Y.R.Wang,X.F.Meng,W.J.Sun,H.Zhang,X.C.Cheng,G.Y.Dong,and X.X.Shen,”Simple direct extraction of unknown phase shift andwavefront reconstruction in generalized phase-shifting interferometry:algorithm and experiments,”Opt.Lett.33,776-778(2008).
文献9:J.Deng,H.Wang,D.Zhang,L.Zhong,J.Fan,and X.Lu,"Phase shiftextraction algorithm based on Euclidean matrix norm,"Opt.Lett.38,1506-1508(2013)
文献10:N.T.Shaked,Y.Zhu,M.T.Rinehart,and A.Wax,"Two-step-only phase-shifting interferometry with optimized detector bandwidth for microscopy oflive cells,"Opt.Express 17,15585-15591(2009).
文献11:J.R.Fienup,"Phase retrieval algorithms:a comparison,"Appl.Opt.21,2758-2769(1982).
文献12:U.Weierstall,Q.Chen,J.C.H.Spence,M.R.Howells,M.Isaacson,andR.R.Panepucci,"Image reconstruction from electron and X-ray diffractionpatterns using iterative algorithms,"Ultramicroscopy 90,171-195(2002).
文献13:D.Shapiro,P.Thibault,T.Beetz,V.Elser,M.Howells,C.Jacobsen,J.Kirz,E.Lima,H.Miao,A.M.Neiman,and D.Sayre,"Biological imaging by soft x-raydiffraction microscopy,"Proc.Natl.Acad.Sci.U.S.A.102,15343-15346(2005).
文献14:J.Miao,D.Sayre,and H.N.Chapman,"Phase retrieval from themagnitude of the Fourier transforms of nonperiodic objects,"J.Opt.Soc.Am.A15,1662-1669(1998)。
文献15:M.Agour,P.F.Almoro,and C.Falldorf,“Investigation of smoothwave fronts using SLM-based phase retrieval and a phase diffuser,”J.Europ.Opt.Soc.Rap.Public.7,12046(2012).
文献16:J.Rodrigo,H.Duadi,T.Alieva,and Z.Zalevsky,"Multi-stage phaseretrieval algorithm based upon the gyrator transform,"Opt.Express 18,1510-1520(2010).
文献17:H.I.Campbell,S.Zhang,A.H.Greenaway,and S.Restaino,"Generalizedphase diversity for wave-front sensing,"Opt.Lett.29,2707-2709(2004).
文献18:H.M.L.Faulkner and J.M.Rodenburg,"Movable aperture lenslesstransmission microscopy:A novel phase retrieval rlgorithm,"Phys.Rev.Lett.93,023903(2004).
文献19:J.M.Rodenburg,A.C.Hurst,A.G.Cullis,B.R.Dobson,F.Pfeiffer,O.Bunk,C.David,K.Jefimovs,and I.Johnson,"Hard-X-ray lensless imaging ofextended objects,"Phys.Rev.Lett.98,034801(2007).
文献20:P.Thibault,M.Dierolf,A.Menzel,O.Bunk,C.David,and F.Pfeiffer,"High-resolution scanning X-ray diffraction microscopy,"Science 321,379-382(2008).
文献21:J.M.Rodenburg,A.C.Hurst,and A.G.Cullis,"Transmissionmicroscopy without lenses for objects of unlimited size,"Ultramicroscopy 107,227-231(2007).文献22:A.Anand,G.Pedrini,W.Osten,and P Almoro,"Wavefrontsensing with random amplitude mask and phase retrieval,"Opt.Lett.32,1584-1586(2007).
文献23:A.Fannjiang,"Absolute uniqueness of phase retrieval withrandom illumination,"Inverse Problems 28,75008(2012).
文献24:E.J.Candès,Y.C.Eldar,T.Strohmer,and V.Voroninski,"Phaseretrieval via matrix completion,"SIAM J.Imaging Sci.6,199-225(2013).
文献25:I.Johnson,K.Jefimovs,O.Bunk,C.David,M.Dierolf,J.Gray,D.Renker,and F.Pfeiffer,"Coherent diffractive imaging using phase frontmodifications,"Phys.Rev.Lett.100,155503(2008).
文献26:F.Zhang and J.M.Rodenburg,“Phase retrieval based on wave-frontrelay and modulation,”Phys.Rev.B 82,121104(2010)
文献27:C.Falldorf,M.Agour,C.von.Kopylow,and R.B.Bergmann,"Phaseretrieval by means of a spatial light modulator in the Fourier domain of animaging system,"Appl.Opt.49,1826-1830(2010)
文献28:A.Fannjiang,and W.Liao,"Phase retrieval with random phaseillumination,"J.Opt.Soc.Am.A 29,1847-1859(2012).
文献29:H.Tao,S.P.Veetil,J.Cheng,X.Pan,H.Wang,C.Liu,and J.-Q.Zhu,"Measurement of the complex transmittance of large optical elements withmodulation coherent imaging,"Appl.Opt.54,1776-1781(2015).
文献30:H.-Y.Wang,C.Liu,S.P Veetil,X.-C.Pan,and J.-Q.Zhu,"Measurementof the complex transmittance of large optical elements with PtychographicalIterative Engine,"Opt.Express 22,2159-2166(2014).
文献31:R.Horisaki,Y.Ogura,M.Aino,and J.Tanida,"Single-shot phaseimaging with a coded aperture,"Opt.Lett.39,6466-6469(2014).
文献32:R.Horisaki and J.Tanida,"Multidimensional object acquisitionby single-shot phase imaging with a coded aperture,"Opt.Express 23,9696-9704(2015).
文献33:Y.Rivenson,A.Stern,and B.Javidi,"Compressive FresnelHolography,"J.Disp.Technol.6,506-509(2010).
发明内容
本发明针对现有无透镜相干衍射成像技术存在的问题,提供一种基于互补随机抽样的无透镜衍射成像方法,该方法能够提高物波波前复振幅的恢复效率和质量,而且还能有效避免由物波波前的二值随机抽样造成的波前抽样损失。
本发明的基于互补随机抽样的无透镜衍射成像方法,是:
引入一组具有互补特性的二值随机抽样屏(简称互补随机抽样(CRS)屏);将这组互补抽样屏依次置于被测物体和记录平面之间,并依次记录透过抽样屏的物波在记录平面处的强度图样;从所记录的强度图样中恢复出被测物波的复振幅(即振幅和相位)分布,利用所恢复出的复振幅分布再经过数字衍射运算,即得到被测物体在空间任意位置的衍射成像。
随机抽样屏的幅数大于等于2,最佳值为4。
从所记录的强度图样中恢复出被测物波波前的复振幅的过程为:
(1)设定记录面上的初始复振幅的振幅和相位分别为和
(2)由记录面的初始复振幅通过逆衍射算法计算出该复振幅逆衍射回抽样面处得到的抽样面复振幅;
(3)将得到的抽样面复振幅乘以对应的随机抽样屏函数Sn得到修正抽样面复振幅然后通过正向衍射算法计算该修正抽样面复振幅再衍射到记录面处的记录面复振幅;
(4)保持记录面复振幅的相位而用代替其振幅部分,得到修正记录面复振幅,再通过逆衍射算法计算出该修正记录面复振幅逆衍射回抽样面处的抽样面复振幅;
(5)重复步骤(3)和(4),直到计算出的修正抽样面复振幅收敛;
(6)将所有的相加,得到最终恢复出的物波波前复振幅分布。
本发明通过在物体和记录面之间***互补随机抽样屏,大大提高了相位恢复的精确度和迭代效率;同时,可以有效消除传统方法中由于二值随机抽样引起的波前抽样损失。另外,由于本方法中迭代运算过程仅应用于抽样面和记录面之间,可适用于更一般的复值物体,并且不需要对物体进行特殊的限制。
附图说明
图1给出了实现本发明方法的光路原理示意图。
图2是本发明方法的实验结果示意图。(a)-(d)是设计的一组二值CRS图样(其中插图为相应CRS屏的局部放大图),设计参数N=4。(e)-(h)给出了将图(a)-d)所示的CRS图样分别显示到SLM上时由图像传感器得到的记录平面RP处的衍射光强度图样。(i)和(j)是基于本发明方法从所记录的强度图样中回复出来的抽样面SP处物波波前复振幅的振幅和相位分布。(k)和(l)则为通过对抽样面上恢复出的复振幅进行逆衍射运算得到的在原物体所在位置处的物波衍射成像的振幅和相位分布。
具体实施方式
如图1所示,实现本发明基于互补随机抽样的无透镜衍射成像方法的光路,主要由照明光S、被测物体Obj,置于抽样面SP处的随机抽样屏和置于记录面RP处的图像传感器构成。假设透过物体Obj的光经过一个菲涅耳衍射在抽样面SP处形成的待测物波波前的复振幅分布为o(x,y),置于抽样面SP处的随机抽样屏的透过率分布为S(x,y),抽样面SP到记录面RP的距离为z,则物波波前o(x,y)经随机抽样屏抽样后继续衍射到达记录面RP处的衍射光场的强度分布I(u,v)可用以下公式表示:
I(u,v)=|Frz{S(x,y)o(x,y)}|2, (1)
其中,Frz{}表示衍射距离为z的衍射积分算符。
文献31中已通过理论分析和计算机模拟证明,透过抽样屏的抽样物波复振幅S(x,y)o(x,y)可通过常规的迭代算法如输入-输出迭代算(参见文献13)从公式(1)所示的强度图样I(u,v)中高效率地恢复出来,该复振幅恢复过程可用公式表示为:
S(x,y)o(x,y)≈Iter{I(u,v),S(x,y)}, (2)
式中,Iter{}表示以S(x,y)为限制条件从强度信息I(u,v)中恢复抽样波前复振幅S(x,y)o(x,y)的迭代运算。但在RS-CDI方法中,由于只采用了单幅二值随机抽样,该方法只能得到被测物波波前的稀疏抽样信息。由二值随机抽样屏的引入导致的抽样损失和噪声虽然通过进一步采用压缩感知算法(参见文献33)可以在一定程度上得到抑制,但代价是需要在物面和抽样面之间引入二次迭代运算和物体限制条件。
本发明提出了一种从原理上消除由二值随机抽样屏导致的抽样损失,从而进一步提高波前复振幅恢复效率和质量的新方法。与已有方法所不同的是采用了一组具有互补特性的二值随机抽样屏(以下简称为CRS屏);这组CRS屏的透过率Sn(p,q)是由以下公式确定的:
其中,N为该组CRS屏的幅数,其取值为大于或等于2的整数;(p,q)是二维离散坐标;R(p,q)是一个二维随机整数阵列,其取值范围为整数1到整数N。CRS屏的互补特性是指它们的透过率之和满足以下条件:
公式(4)称为互补条件。
将一组CRS屏中的一幅放入图1的抽样面SP上,在记录面RP上,就可以用图像传感器记录下对应的衍射强度,其中第n幅CRS屏所对应的强度图样可以表示为:
In(u,v)=|Frz{Sn(p,q)o(p,q)}|2,n=1,2,…,N. (5)
利用该强度分布通过常规的迭代算法可以将透过第n幅CRS屏的物波波前复振幅信息恢复出来,即:
Sn(p,q)o(p,q)exp(iδn)≈Iter{In(u,v),Sn(p,q)},n=1,2,…,N, (6)
其中,δn是由迭代算法引入的一个附加相位。进一步的模拟和实验结果表明,如果在对所有强度图样In(u,v)进行迭代运算时都采用相同的初始迭代相位(可为任意常数),则附加相位δn将保持不变,即δn≈δ0,其中δ0是一个与初始迭代相位有关的常数。利用CRS屏的互补特性,将(6)式所示的所有恢复出的物波抽样复振幅相加并取附加相位δn为零就可以得到没有抽样损失的物波波前o(p,q),即:
从公式(7),可以看到,利用CRS屏的互补特点可以很好地消除单个抽样屏对被测物波波前影响。利用公式(7)得到的物波复振幅(包含相位和振幅)就可在计算机中通过简单的数字逆菲涅耳衍射算法实现物体的无透镜衍射成像。利用这种方法可以实现复数物体或三维物体的波前再现和数字成像。
上述基于互补抽样的衍射成像方法的主要实施步骤可归纳如下:
(a)根据公式(3)设计一组CRS屏,幅数大于等于2,其最佳值为4。
(b)将CRS屏依次置于图1所示的记录光路的抽样面SP处,并用图像传感器依次获取记录面RP处的衍射强度图样In(u,v)并输入计算机中。
(c)利用常规的迭代算法从所记录的强度图样中恢复出被测物波波前的复振幅,其运算过程为:
(1)设定记录面RP上的初始复振幅的振幅和相位分别为和
(2)由记录面的初始复振幅通过常规的逆衍射算法计算出该复振幅逆衍射回抽样面SP处得到的抽样面复振幅;
(3)将得到的抽样面复振幅乘以对应的CRS屏函数Sn得到修正抽样面复振幅然后通过常规的正向衍射算法计算该修正抽样面复振幅再衍射到记录面RP处的记录面复振幅;
(4)保持记录面复振幅的相位而用代替的其振幅部分,得到修正记录面复振幅,再通过常规的逆衍射算法计算出该修正记录面复振幅逆衍射回抽样面SP处的抽样面复振幅;
(5)反复重复步骤(3)和(4),直到计算出的修正抽样面复振幅收敛;(6)将所有的相加,得到最终恢复出的物波波前复振幅分布。
(d)利用最终恢复出的物波波前复振幅信息,通过简单的数字衍射计算就可以得到物体在空间任意位置的衍射成像。
原理上,本发明的方法可适用于任何波长,比如红外线、可见光、X射线以及电子束。为了实验验证本发明方法的可行性,建立了一个无透镜相干衍射成像实验装置,采用了波长λ=632.8nm的氦氖激光作为照明光光源;物体是一个透射型的标准USAF分辨率板,物体到抽样屏的距离设置为185mm;抽样屏到记录面的距离为201mm。将一个透射式的空间光调制器(SLM)置于如图1所示的抽样面SP处,并将SLM设为强度调制模式。实验所需的CRS是通过将根据公式(3)设计好的CRS图样分别显示到SLM上来实现的。图2(a)-2(d)是基于公式(3)设计的一组二值CRS图样(其中插图为相应CRS屏的局部放大图),该组CRS图样的设计参数N=4。透过CRS屏到达记录平面RP处的衍射光强度图样则用一个商用的图像传感器记录。
图2(e)-(h)给出了将图2(a)-2(d)所示的CRS图样分别显示到SLM上时,由图像传感器得到的记录平面RP处的衍射光强度图样。图2(i)和2(j)则是基于本发明方法从所记录的强度图样中回复出来的抽样平面SP处物波波前复振幅的振幅和相位分布,其中所用的迭代运算次数为30次。图2(k)和2(l)则进一步给出了通过对抽样面上恢复出的复振幅进行逆衍射运算得到的在原物体所在位置处的物波衍射成像的振幅和相位分布。由图2可见,USAF分辨率板的振幅和相位得到了成功恢复和成像。
总之,本发明通过理论分析和实验证明,物波的全部复振幅信息可从一组通过二值CRS的衍射强度图样中准确恢复出来,而物体的衍射成像则可通过直接计算复振幅的数字衍射来实现。该方法中引入了一组具有(4)式所示的互补特性的二值随机抽样屏(CRS);将这组互补抽样屏依次置于被测物体和记录平面之间,并用图像传感器依次记录透过抽样屏的物波在记录平面RP处的强度图样。从所记录的强度图样中利用常规的迭代算法可准确恢复出被测物波的振幅和相位分布,利用所恢复出的振幅和相位再经过数字衍射运算,即得到被测物体在空间任意位置的衍射成像。
实验结果表明,通过在物体和记录面之间***CRS屏可以大大提高相位恢复的精确度和迭代效率;同时,可以有效消除传统方法中由于二值随机抽样引起的抽样损失。另外,由于本方法中迭代运算过程仅应用于抽样面SP和记录面RP之间,这种方法可适用于更一般的复值物体,并且不需要对物体进行特殊的限制。
Claims (3)
1.一种基于互补随机抽样的无透镜衍射成像方法,其特征是:
引入一组具有互补特性的二值随机抽样屏;将这组互补抽样屏依次置于被测物体和记录平面之间,并依次记录透过抽样屏的物波在记录平面处的强度图样;从所记录的强度图样中恢复出被测物波的复振幅分布,利用所恢复出的复振幅分布再经过数字衍射运算,即得到被测物体在空间任意位置的衍射成像;
具有互补特性的二值随机抽样屏的互补特性是指它们的透过率之和满足以下条件:其中上述公式称为互补条件,互补特性的二值随机抽样屏的透过率Sn(p,q)是由以下公式确定的:
其中,N为该组CRS屏的幅数,其取值为大于或等于2的整数;(p,q)是二维离散坐标;R(p,q)是一个二维随机整数阵列,其取值范围为整数1到整数N。
2.根据权利要求1所述的基于互补随机抽样的无透镜衍射成像方法,其特征是:所述具有互补特性的二值随机抽样屏的幅数为4。
3.根据权利要求1所述的基于互补随机抽样的无透镜衍射成像方法,其特征是:所述从所记录的强度图样中恢复出被测物波波前的复振幅的过程为:
(1)设定记录面上的初始复振幅的振幅和相位分别为和
(2)由记录面的初始复振幅通过逆衍射算法计算出该复振幅逆衍射回抽样面处得到的抽样面复振幅;
(3)将得到的抽样面复振幅乘以对应的随机抽样屏函数Sn得到修正抽样面复振幅然后通过正向衍射算法计算该修正抽样面复振幅再衍射到记录面处的记录面复振幅;
(4)保持记录面复振幅的相位而用代替的其振幅部分,得到修正记录面复振幅,再通过逆衍射算法计算出该修正记录面复振幅逆衍射回抽样面处的抽样面复振幅;
(5)重复步骤(3)和(4),直到计算出的修正抽样面复振幅收敛;
(6)将所有的相加,得到最终恢复出的物波波前复振幅分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698056.3A CN105158921B (zh) | 2015-10-26 | 2015-10-26 | 一种基于互补随机抽样的无透镜衍射成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510698056.3A CN105158921B (zh) | 2015-10-26 | 2015-10-26 | 一种基于互补随机抽样的无透镜衍射成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105158921A CN105158921A (zh) | 2015-12-16 |
CN105158921B true CN105158921B (zh) | 2017-08-08 |
Family
ID=54799828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510698056.3A Active CN105158921B (zh) | 2015-10-26 | 2015-10-26 | 一种基于互补随机抽样的无透镜衍射成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105158921B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11175627B2 (en) * | 2016-11-04 | 2021-11-16 | miDiagnostics NV | System and method for reconstruction of holographic lens-free images by multi-depth sparse phase recovery |
CN107828654B (zh) * | 2017-10-20 | 2019-12-03 | 江苏大学 | 基于无透镜衍射成像的细胞活性无标记监测装置与方法 |
CN109143607B (zh) * | 2018-09-17 | 2020-09-18 | 深圳奥比中光科技有限公司 | 补偿显示屏、屏下光学***及电子设备 |
CN111474188B (zh) * | 2019-01-24 | 2021-07-27 | 中国科学院上海光学精密机械研究所 | 基于动态调制的单次曝光波前重建和相位成像装置和方法 |
CN112683794A (zh) * | 2020-12-11 | 2021-04-20 | 中国科学院上海光学精密机械研究所 | 基于波前调制的相位成像及元件检测的装置和方法 |
CN117631246A (zh) * | 2023-12-25 | 2024-03-01 | 山东理工大学 | 一种定量相位成像方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4778662B2 (ja) * | 2000-06-29 | 2011-09-21 | パナソニック株式会社 | 光ヘッド装置および光情報処理装置 |
WO2012119094A2 (en) * | 2011-03-03 | 2012-09-07 | California Institute Of Technology | Light guided pixel |
CN102660457A (zh) * | 2012-04-17 | 2012-09-12 | 南昌航空大学 | 无透镜全息衍射成像血细胞分析计数的装置及方法 |
CN103984979B (zh) * | 2014-06-06 | 2017-11-03 | 南昌航空大学 | 无透镜衍射成像的藻类自动检测计数装置及方法 |
-
2015
- 2015-10-26 CN CN201510698056.3A patent/CN105158921B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105158921A (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105158921B (zh) | 一种基于互补随机抽样的无透镜衍射成像方法 | |
Shechtman et al. | Phase retrieval with application to optical imaging: a contemporary overview | |
Poon et al. | Introduction to modern digital holography: with MATLAB | |
CN106990694B (zh) | 一种部分相干光照明下的非迭代相位恢复装置及方法 | |
CN106094487B (zh) | 基于多个记录距离的太赫兹同轴全息成像方法 | |
CN105404128B (zh) | 多帧相移数字全息方法及装置 | |
US20140016183A1 (en) | Reconstruction of nonlinear wave propagation | |
CN1260052A (zh) | 直接一数字全息术,全息干涉测量术,和全息影象 | |
EP1230576B1 (en) | System and method for recovering phase information of a wave front | |
CN103322940B (zh) | 一种获取三维形貌显微图像的方法 | |
CN103323938B (zh) | 一种获取立体显微图像的方法 | |
CN103149827A (zh) | 消除单光束同轴数字全息直流项和共轭像的方法 | |
CN109283821B (zh) | 基于涡旋透镜的相移数字全息单次曝光成像装置及方法 | |
CN111736442A (zh) | 基于各向同性双焦超颖透镜的单次曝光产生全息图的方法 | |
Özcan et al. | Digital holography image reconstruction methods | |
Hoffmann et al. | Holographic method for site-resolved detection of a 2D array of ultracold atoms | |
Gautam et al. | Reconstruction of complex-object using edge point referencing | |
Tsang et al. | Non-diffractive optical scanning holography for hologram acquisition | |
He et al. | Resolution improvement in off-axis digital holography using an iterative scheme | |
Mansoor et al. | A versatile digital holography software for three-dimensional investigations | |
CN114764220B (zh) | 一种基于离轴数字全息改善散斑自相关重建效果的方法 | |
Hamed | Scanning holography using a modulated linear pupil: simulation | |
Chaab et al. | An adaptive anisotropic magnification algorithm of numerical reconstructed holographic image via a cascaded Fresnel transform | |
JPWO2012115042A1 (ja) | デジタルホログラフィ装置、及びデジタルホログラフィによる3次元像再生方法 | |
Chen et al. | High-quality phase-only holographic projection based on bandwidth constraint optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |