• Photonics Research
  • Vol. 9, Issue 6, 1003 (2021)
Cheng Shen1、*, Mingshu Liang1, An Pan2, and Changhuei Yang1
Author Affiliations
  • 1Department of Electrical Engineering, California Institute of Technology, Pasadena, California 91125, USA
  • 2Xi’an Institute of Optics and Precision Mechanics (XIOPM), Chinese Academy of Sciences (CAS), Xi’an 710119, China
  • show less
    DOI: 10.1364/PRJ.419886 Cite this Article Set citation alerts
    Cheng Shen, Mingshu Liang, An Pan, Changhuei Yang. Non-iterative complex wave-field reconstruction based on Kramers–Kronig relations[J]. Photonics Research, 2021, 9(6): 1003 Copy Citation Text show less

    Abstract

    A non-iterative and non-interferometric computational imaging method to reconstruct a complex wave field called synthetic aperture imaging based on Kramers–Kronig relations (KKSAI) is reported. By collecting images through a modified microscope system with pupil modulation capability, we show that the phase and amplitude profile of the sample at pupil limited resolution can be extracted from as few as two intensity images by using Kramers–Kronig (KK) relations. It is established that as long as each subaperture’s edge crosses the pupil center, the collected raw images are mathematically analogous to off-axis holograms. This in turn allows us to adapt a recently reported KK-relations-based phase recovery framework in off-axis holography for use in KKSAI. KKSAI is non-iterative, free of parameter tuning, and applicable to a wider range of samples. Simulation and experiment results have proved that it has much lower computational burden and achieves the best reconstruction quality when compared with two existing phase imaging methods.
    S(u,v)=F{s}(u,v)·C(u,v),

    View in Article

    Si(u,v)=S(u,v)·D(uui,vvi).

    View in Article

    Ii(x,y)=|F1{Si(u,v)}|2.

    View in Article

    Ii(x,y)=|F1{Si(u,v)}|2=|F1{Si(u,v)+δ(u,v)}|2=|si(x,y)+ej(0·x+0·y)|2,

    View in Article

    I1(x,y)=|F1{S1(u,v)}|2=|F1{S1(u+u1,v+v1)}|2.

    View in Article

    S˜1(u,v)=S1(u+u1,v+v1)

    View in Article

    I1(x,y)=|F1{S˜1(u,v)+δ(u+u1,v+v1)}|2=|s˜1(x,y)+ej(u1·x+v1·y)|2=|s˜1(x,y)|2+1+s˜1*(x,y)·ej(u1·x+v1·y)+s˜1(x,y)·ej(u1·x+v1·y),

    View in Article

    F1{I1}=S˜1(u,v)S˜1(u,v)+δ(u,v)+S˜1*[(u+u1),(v+v1)]+S˜1(uu1,vv1),

    View in Article

    r1(x,y)=F1{δ(u+u1,v+v1)}=ej(u1·x+v1·y).

    View in Article

    H1(u,v)=jsgn(v1)·sgn(v),

    View in Article

    X=ln[F1{S1(u+u1,v+v1)}/r1(x,y)]=ln{[s˜1(x,y)+r1(x,y)]/r1(x,y)}.

    View in Article

    Re{X}=12ln[I1(x,y)/|r1(x,y)|2],Im{X}=F1{F{Re{X}}·H1(u,v)},

    View in Article

    S1(u+u1,v+v1)=F{eRe{X}+jIm{X}·r1(x,y)}.

    View in Article

    Si(u,v),i=1,2,3,4.

    View in Article

    S(u,v)=i=14Si(u,v)/[i=14D(uui,vvi)+ε],

    View in Article

    I1(x,y)I(l)s˜1(x,y)s˜(l)r1(x,y)=ej(u1·x+v1·y)r(l)=ejρr·l,

    View in Article

    I(l)=|s˜(l)+r(l)|2.

    View in Article

    I|r|2=|s˜r+1|2,

    View in Article

    ln|s˜r+1|=12[ln(I)ln(|r|2)].

    View in Article

    eX=s˜r+1=eRe{X}·ejIm{X},

    View in Article

    |s˜r+1|=eRe{X},

    View in Article

    Re{X}=ln|s˜r+1|=12[ln(I)ln(|r|2)].

    View in Article

    α(l)=s˜r=s˜(l)·ejρr·l.

    View in Article

    X=ln(α+1)=n=0(1)nn+1αn+1,

    View in Article

    A(ρ,l)=S˜(ρ|ρr|,l),

    View in Article

    Cheng Shen, Mingshu Liang, An Pan, Changhuei Yang. Non-iterative complex wave-field reconstruction based on Kramers–Kronig relations[J]. Photonics Research, 2021, 9(6): 1003
    Download Citation