• Photonics Research
  • Vol. 10, Issue 11, 2471 (2022)
Wei Li1, Bingjian Wang2, Tengfei Wu3, Feihu Xu4, and Xiaopeng Shao1、*
Author Affiliations
  • 1School of Optoelectronic Engineering, Xidian University, Xi’an 710071, China
  • 2School of Physics, Xidian University, Xi’an 710071, China
  • 3Laboratoire Kastler Brossel, ENS–Université PSL, CNRS, Sorbonne Université, College de France, F-75005 Paris, France
  • 4Hefei National Laboratory for Physical Sciences at Microscale and School of Physical Science, University of Science and Technology of China, Hefei 230026, China
  • show less
    DOI: 10.1364/PRJ.466065 Cite this Article Set citation alerts
    Wei Li, Bingjian Wang, Tengfei Wu, Feihu Xu, Xiaopeng Shao. Lensless imaging through thin scattering layers under broadband illumination[J]. Photonics Research, 2022, 10(11): 2471 Copy Citation Text show less

    Abstract

    Lensless scattering imaging is a prospective approach to microscopy in which a high-resolution image of an object is reconstructed from one or more measured speckle patterns, thus providing a solution in situations where the use of imaging optics is not possible. However, current lensless scattering imaging methods are typically limited by the need for a light source with a narrowband spectrum. Here, we propose two general approaches that enable single-shot lensless scattering imaging under broadband illumination in both noninvasive [without point spread function (PSF) calibration] and invasive (with PSF calibration) modes. The first noninvasive approach is based on a numerical refinement of the broadband pattern in the cepstrum incorporated with a modified phase retrieval strategy. The latter invasive approach is correlation inspired and generalized within a computational optimization framework. Both approaches are experimentally verified using visible radiation with a full-width-at-half-maximum bandwidth as wide as 280 nm (Δλ/λ=44.8%) and a speckle contrast ratio as low as 0.0823. Because of its generality and ease of implementation, we expect this method to find widespread applications in ultrafast science, passive sensing, and biomedical applications.
    IB=O*PSFB+N,

    View in Article

    IBIB(OO)*(PSFBPSFB),

    View in Article

    PSFB(λ,Δλ)=i=1MωiPSFN(λi,Δλi),

    View in Article

    PSFB(λ,Δλ)PSFB(λ,Δλ)=i=1Mωi2PSFN(λi,Δλi)PSFN(λi,Δλi)+i=1MjiMωiωjPSFN(λi,Δλi)PSFN(λj,Δλj).

    View in Article

    i=1Mωi2PSFN(λi,Δλi)PSFN(λi,Δλi)i=1Mωi2[1+|2J1(πDr/λiz)πDr/λiz|2],

    View in Article

    i=1MjiMωiωjPSFN(λi,Δλi)PSFN(λj,Δλj)i=1MjiMωiωj{1+exp[σh2(n1)|kikj|]×{2J1(πDr/λ¯ijz)πDr/λ¯ijz*F{exp[j(kikj)r22z]}}2},

    View in Article

    Psρ(r)={ρ(r)if  rS0otherwise,

    View in Article

    Pmρ˜(k)=Pm|ρ˜(k)|eiφ(k)=m(k)eiφ(k).

    View in Article

    ρ(n+1)(r)=[PsPm+Ps_(IβPm)]ρ(n)(r)

    View in Article

    ρ(n+1)(r)=[12β(RsRm+I)+(1β)Pm]ρ(n)(r).

    View in Article

    O^=IBPSFB=F1[F(IB)·F(PSFB)¯].

    View in Article

    O^=F1[|F(IB)|α·|F(PSFB)|β·ei(φIBφPSFB)],

    View in Article

    g=H(f)+e,

    View in Article

    f^=argminf012gH(f)22+τfTV+λf1,

    View in Article

    fTV=i(|[Dx(f)]i|+|[Dy(f)]i|),

    View in Article

    f^k+1f^k2/f^k2ε,

    View in Article

    Ek(r,z)=1jλzrEk(r,0)ejk2z(rr)2d2r.(A1)

    View in Article

    Ek(r,0)=tS(r,0)exp[jϕ(r,0)]exp[jk(t^·r^)r],(A2)

    View in Article

    ϕ(r,0)=k(nt^·z^+o^·z^)h(r,0),(A3)

    View in Article

    q=kokt=qt+qzz^.(A4)

    View in Article

    ϕ(r,0)=qzh(r,0).(A5)

    View in Article

    ΓE(r1;r2,z)=Ek1(r1,z)Ek2*(r2,z)¯.(A6)

    View in Article

    ΓE(r1;r2,z)=1λ1λ2z2r1r2ΓE(r1;r2,0)×exp[j2z(k1r12k2r22)]×exp[jz(k1r1r1k2r2r2)]d2r1d2r2.(A7)

    View in Article

    ΓE(r1;r2,0)=Ek1(r1,0)Ek2*(r2,0)¯=|t|2S(r1,0)S*(r2,0)exp{j[qz1h(r1)qz2h(r2)¯]}=|t|2S(r1,0)S*(r2,0)Mh(qz1,qz2)×exp[j(k1t^1·r1^)r1+j(k2t^2·r2^)r2],(A8)

    View in Article

    ΓE(r1;r2,0)κ|t|2|S(r1,0)|2Mh(Δqz)×exp[j(k1t^1·r1^)r1+j(k2t^2·r2^)r2]δ(r1r2).(A9)

    View in Article

    ΓE(r1;r2,z)=κ|t|2λ1λ2z2Mh(Δqz)r|S(r,0)|2exp[j(k1k2)r22z]×exp[j(k1t^1·r^)r+j(k2t^2·r^)r]×exp[j(k1o^1·r^)r+j(k2o^2·r^)r]d2r.(A10)

    View in Article

    qt1=k1t^1·r^+k1o^1·r^qt2=k2t^2·r^+k2o^2·r^.(A11)

    View in Article

    ΓE(r1;r2,z)=κ|t|2λ1λ2z2Mh(Δqz)Ψ(Δqt)r|S(r,0)|2d2r,(A12)

    View in Article

    Ψ(Δqt)=r|S(r,0)|2exp[j(k1k2)r22z]exp[j(Δqt·rt)d2rr|S(r,0)|2d2r,(A13)

    View in Article

    Ψ(Δqt)=F{|S(r,0)|2exp[j(k1k2)r22z]}r|S(r,0)|2d2r.(A14)

    View in Article

    Ψ(Δqt)=2J1(D|Δqt|2)D|Δqt|2*F{exp[j(k1k2)r22z]}.(A15)

    View in Article

    J1(D|Δqt|2)D|Δqt|2J1(πDr/λ¯ijz)πDr/λ¯ijz.(A16)

    View in Article

    |μE(q1,q2)|2=|Mh(Δqz)|2|Ψ(Δqt)|2.(A17)

    View in Article

    |Mh(Δqz)|2=exp(σh2Δqz2).(A18)

    View in Article

    Δqz=(n1)|kikj|.(A19)

    View in Article

    ΓI(Δr)=I¯2[1+|μE(Δr)|2].(A20)

    View in Article

    i=1MjiMωiωjPSFN(λi,Δλi)PSFN(λj,Δλj)i=1MjiMωiωj{1+exp[σh2(n1)|kikj|]×{2J1(πDr/λ¯ijz)πDr/λ¯ijz*F{exp[j(kikj)r22z]}}2}.(A21)

    View in Article

    R[z,r0,r1]{Ek(r0)}=1jλzr0Ek(r0,0)ejk2z(r1r0)2d2r0.(B1)

    View in Article

    PSFB=|kR[v,r1,r2]{R[u,r0,r1]{Ek(r0,0)}tk(r1)p(r1)}dk|2.(B2)

    View in Article

    Wei Li, Bingjian Wang, Tengfei Wu, Feihu Xu, Xiaopeng Shao. Lensless imaging through thin scattering layers under broadband illumination[J]. Photonics Research, 2022, 10(11): 2471
    Download Citation