• Photonics Research
  • Vol. 10, Issue 3, 758 (2022)
Li Song and Edmund Y. Lam*
Author Affiliations
  • Department of Electrical and Electronic Engineering, The University of Hong Kong, Pokfulam, Hong Kong, China
  • show less
    DOI: 10.1364/PRJ.447862 Cite this Article Set citation alerts
    Li Song, Edmund Y. Lam. Fast and robust phase retrieval for masked coherent diffractive imaging[J]. Photonics Research, 2022, 10(3): 758 Copy Citation Text show less

    Abstract

    Conventional phase retrieval algorithms for coherent diffractive imaging (CDI) require many iterations to deliver reasonable results, even using a known mask as a strong constraint in the imaging setup, an approach known as masked CDI. This paper proposes a fast and robust phase retrieval method for masked CDI based on the alternating direction method of multipliers (ADMM). We propose a plug-and-play ADMM to incorporate the prior knowledge of the mask, but note that commonly used denoisers are not suitable as regularizers for complex-valued latent images directly. Therefore, we develop a regularizer based on the structure tensor and Harris corner detector. Compared with conventional phase retrieval methods, our technique can achieve comparable reconstruction results with less time for the masked CDI. Moreover, validation experiments on real in situ CDI data for both intensity and phase objects show that our approach is more than 100 times faster than the baseline method to reconstruct one complex-valued image, making it possible to be used in challenging situations, such as imaging dynamic objects. Furthermore, phase retrieval results for single diffraction patterns show the robustness of the proposed ADMM.
    li(x,y)=ejkzijλzi+exp{jk2zi[(xξ)2+(yη)2]}dξdη,

    View in Article

    u^(x,y)=a(x,y)ejθ(x,y),

    View in Article

    u^(x,y)=mb(x,y)r(x,y)+mo(x,y)u(x,y),

    View in Article

    le(x,y)=li(x,y)u^(x,y).

    View in Article

    lo(x,y)=F{le(x,y)},

    View in Article

    o(x,y)=|lo(x,y)|.

    View in Article

    o=|FL(Wbr+Wou)|,

    View in Article

    minimizeR(u)+C(ϕ)subjecttoFLu^=oϕu^=Wbr+Wou,

    View in Article

    C(ϕ)={0,|ϕ|=1,otherwise.

    View in Article

    u=Wo(u^Wbr).

    View in Article

    minimizeR(u)+C(ϕ)subjecttoFLu^=oϕ,u=Wo(u^Wbr).

    View in Article

    L(u^,ϕ,u,λ,μ)=R(u)+C(ϕ)+ρ2FLu^oϕ+λ22+τ2uWo(u^Wbr)+μ22,

    View in Article

    u^(k+1)=argminu^L(u^,ϕ(k),u(k),λ(k),μ(k)).

    View in Article

    u^L(u^,ϕ(k),u(k),λ(k),μ(k))|u^=u^(k+1)=0.

    View in Article

    u^L=ρ(FL)H(FLu^oϕ+λ)+τWoH(Wou^uWoWbrμ).

    View in Article

    u^(k+1)=(ρLHFHFL+τWoHWo)1[ρ(FL)H(oϕ(k)λ(k))+τWoH(u(k)+WoWbr+μ(k))].

    View in Article

    ϕ(k+1)=Po(λ(k)+FLu^(k+1)),

    View in Article

    u(k+1)=proxRτ[Wo(u^(k+1)Wbrμ(k))],

    View in Article

    λ(k+1)=λ(k)+FLu^(k+1)oϕ(k+1),

    View in Article

    μ(k+1)=μ(k)+u(k+1)Wo(u^(k+1)Wbr),

    View in Article

    Po(s)(i)={s(i)/o(i)s(i)/o(i),o(i)0,s(i)01,otherwise.

    View in Article

    proxRτ(v)=argminx(R(x)+τ2xv).

    View in Article

    R(x)=xH[xD(x)],

    View in Article

    R(x)=(xr+jxi)H{[xrD(xr)]+j[xiD(xi)]}.

    View in Article

    R(x)=γxHQx,

    View in Article

    S=G*(uuH),

    View in Article

    R(u^,o)=|FLu^|o2.

    View in Article

    ϵrel(k+1)=u(k+1)u(k)2u(k+1)2,

    View in Article

    FL[Wou+Wbr+(1WoWb)n]=yϕ,

    View in Article

    Li Song, Edmund Y. Lam. Fast and robust phase retrieval for masked coherent diffractive imaging[J]. Photonics Research, 2022, 10(3): 758
    Download Citation