Laser & Optoelectronics Progress, Volume. 58, Issue 18, 1811003(2021)

Optical-Field Coherence Measurement and Its Applications in Computational Imaging

Runnan Zhang1,2,3, Zewei Cai1,2,3、**, Jiasong Sun1,2,3, Linpeng Lu1,2,3, Haitao Guan1,2,3, Yan Hu1,2,3, Bowen Wang1,2,3, Ning Zhou1,2,3, Qian Chen3、***, and Chao Zuo1,2,3、*
Author Affiliations
  • 1Smart Computational Imaging Laboratory, School of Electronic and Optical Engineering, Nanjing University of Science & Technology, Nanjing, Jiangsu 210094, China;
  • 2Smart Computational Imaging Research Institute, Nanjing University of Science & Technology, Nanjing, Jiangsu 210019, China;
  • 3Jiangsu Key Laboratory of Spectral Imaging & Intelligent Sense, Nanjing, Jiangsu 210094, China;
  • show less
    Figures & Tables(61)
    Several typical examples of light sources with different degrees of temporal coherence and spatial coherence
    Representation of optical signal in time and space. (a) Optical signal is a function of time at a certain point in space; (b) optical signal is a function of space at a certain point in time
    Superposition of light waves with different frequencies. (a) Waves with different frequencies are coherently superimposed into one pulse wave packet; (b) waves with different frequencies are incoherently superimposed into one continuous wave (non-periodic and infinite width), and its phase and amplitude vary randomly
    Basic principle of Fourier transform spectrometer
    Michelson stellar interferometer measures the spatial coherence of a quasi-monochromatic wave field by interferometry to infer the size of the light source. (a) Optical configuration; (b) photograph of a real system
    Relationship between classical coherence theory and phase space optics
    Characterization of common signal transformations in phase space. (a) Fresnel propagation; (b) Chirp modulation (lens); (c) Fourier transform; (d) fractional Fourier transform; (e) magnifier
    Wigner distribution function of special signals. (a) Point source; (b) plane wave; (c) spherical wave; (d) phase slow-varying wave; (e) Gaussian signal
    Parameterization of the light field. (a) Seven-dimensional plenoptic function; (b) two planes parameterization of four-dimensional light field; (c) position-angular parameterization of four-dimensional light field
    Relationship between Wigner distribution function and light field of a smooth coherent wavefront. Phase is represented as the localized spatial frequency (instantaneous frequency) in the Wigner distribution function. Rays travel perpendicularly to the wavefront (phase normal) [128]. (a) Wavefront in real space; (b) Wigner distribution function in phase space; (c) light field in position-angle space
    Classification of phase imaging techniques
    Classification of the coherence measurement techniques
    Young’s interferometry with two holes[22]
    Reversed-wavefront Young’ interferometry[77]
    Distributions of nonredundant array method and experimental scheme[78]. (a1) Superior in points on axes; (a2) central distribution of points; (a3) superior in points out of axes; (b) experimental scheme of nonredundant array method
    Experimental scheme of self-referencing interferometry[80]
    Correspondence between Wigner distribution function and ambiguity function
    Basic principle of phase space tomography. (a) Vertical projection; (b) quarter rotation projection; (c) rotation by 90° projection; (d) superposition of all projections from different angles
    Two different transformations of WDF for phase space tomography. (a) WDF of complex signal; (b) WDF after Fresnel diffraction; (c) phase space WDF after fractional Fourier transform; (d) correspondence between Fresnel diffraction and fractional Fourier transform
    Optical path structure of phase space tomography[85]
    Experimental device for measuring spatial coherence based on edge diffraction[84]
    Direct phase space measurement. (a) Direct phase space measurement based on pinhole scanning; (b) direct phase space measurement based on microlens array
    Schematic of a simplistic view of coherent field and partially (spatially) coherent field. (a) A coherent field requires a 2D complex amplitude representation, the surface of the constant phase is interpreted as wavefronts with geometric light rays traveling normal to them; (b) a partially coherent field requires a 4D coherence function to accurately represent its properties such as propagation and diffraction. The “phase” (generalized phase) of a partially coherent light field is the statistical average of phases (spatial frequency, direction of propagation) at each position in space
    Principle of the Shack-Hartmann sensor and light field camera. (a) For coherent field, the Shack-Hartmann sensor forms a focus spot array sensor signal; (b) for partially coherent field, the Shack-Hartmann sensor forms an extended source array sensor signal; (c) for incoherent imaging, the light field camera produces a 2D sub-aperture image array
    Classification of light field imaging techniques
    Light field capture based on camera arrays. (a) Light field gantry[126]; (b) large camera arrays[110]; (c) micro light field acquisition acquired by the 5×5 camera array system[189]
    Various light field cameras based on microlens array
    Computational light field imaging based on coded mask. (a) Light field acquisition of mask enhanced camera[191]; (b) light field acquisition of compressive photography[192]
    Light field imaging based on programmable aperture. (a) Programmable aperture light field camera[104]; (b) programmable aperture light field microscope[106]
    Light field representation of a slowly varying object under spatially stationary illumination[128]
    Light field microscope model. (a) Traditional bright field microscope; (b) light field microscope[111]; (c) light field microscopic model based on wave optics theory[112]; (d) Fourier light field microscope[113]
    Comparison between TIE and WOTF[207]. (a1)(b1) Physical implications of TIE and WOTF; (a2)(b2) geometric illustrations for deriving the PGTF and WOTF; (a3)(b3) PGTF and WOTF for phase imaging under different s
    Direct visualization of coherent images reconstructed from coherent holograms[209]
    Photon correlation holography[213]. (a) Concept diagram of photon-dependent holography; (b) schematic diagram of intensity interferometer
    Representative optical setup for incoherent holography. (a) Optical path of modified triangular interferometer[217]; (b) optical path of FINCH[214]; (c) Michelson interferometer[119]; (d) Sagnac interferometer[215]; (e) Mach-Zehnder interferometer[216,218]
    Typical imaging optical path for COACH. (a) Structure of COACH[115]; (b) structure of I-COACH[117]; (c) structure of LI-COACH[116]
    Schematic of non-invasive scattering imaging through strong scattering layers[237]
    Single frame imaging based on speckle autocorrelation through strong scattering layer[238]. (a) Experimental setup; (b) raw camera image; (c) autocorrelation; (d) image reconstructed by an iterative phase-retrieval algorithm; (e) photograph of the experiment; (f) raw camera image; (g)--(k) Left column is calculated autocorrelation, middle column is reconstructed object; right column is image of the real object
    Conventional incoherent synthetic aperture structure. (a) Michelson interferometer; (b) common secondary structure; (c) multiple telescopes structure
    Design model of the initial generation of SPIDER imaging conceptual system. (a) Explosive view of SPIDER; (b) PIC schematics of the two physical baselines and three spectral bands; (c) arrangement of SPIDER microlens; (d) corresponding frequency-spectrum coverage
    Incoherent synthetic aperture technology based on FINCH[243]
    RSI of visible cone-beam tomography[79]
    Application of light field microscopy in bioscience.(a) Mouse with a head-mounted MiniLFM[252]; (b) imaging Golgi-derived membrane vesicles in living COS-7 cells using HR-LFM [74]; (c) dynamics during neutrophil migration in mouse liver using DAOSLIMIT[73]; (d) hunting activity of zebrafish and the neural activity of mouse brain observed by confocal light field microscope[75]
    Microscopy imaging based on FINCH. (a) FINCHSCOPE schematic; (b) FINCHSCOPE fluorescence sections of pollen grains[118]; (c) wide-field image and reconstructed FINCH image of pollen grains captured using a 20×(0.75 NA) objective[253]; (d) comparative imaging of three different Golgi apparatus proteins in HeLa cells using wide-field (left) and FINCH(right)[255]
    Light field imaging in computational photography. (a) Light field refocusing[101]; (b) synthetic aperture imaging based on light field[259]
    Computational photography refocusing based on FINCH. (a) Digital refocusing based on FINCH[214]; (b) colorful digital holography refocusing[260]; (c) full color holographic digital refocusing under natural light illumination[119]
    X-ray characterization via phase space tomography. (a) Measured intensity distribution of the X rays as a function of lateral position and along the direction of propagation[263]; (b) phase space density reconstructed from the data in Fig.47(a); (c)(d) measured complex degree of coherence for the beams in the two conditions [263]
    Optical beam characterization via phase space tomography. (a) 1D signal[265]; (b) optical beams separable in Cartesian coordinates[264]; (c) rotationally symmetric beams[266]; (d) intensity distributions of the test beams with different degrees of coherence (first row), the Wigner distribution function of the beams exhibits hidden differences associated with their coherence state (second row)[267]
    Schematic diagram of phase retrieval and factor M2 calculation[268]. (a)(b) Axial intensity images at two different longitudinal positions; (c) phase retrieval by TIE; (d) reconstructed intensity distribution at any selected plane; (e) performing a hyperbolic fit to the beam widths and calculating the M2
    Under different numerical apertures, the phase is recovered directly through the gravity of the light field[128]. (a) 0.05; (b) 0.15; (c) 0.2; (d) 0.25
    Reconstructed phases with and without mode decomposition method under partially coherent illumination[269]
    Stack imaging based on coherent mode decomposition. (a) Decoherence in scattering imaging[270];(b) experimental scheme of Fourier stack imaging with single-mode and multi-mode multiplexing[271]
    Synthetic aperture technique based on FINCH. (a)--(c) Three phase functions loaded on SLM; (d) single aperture reconstruction result; (e) synthetic multi aperture reconstruction result[243]; (f) image obtained by the conventional imaging system; (g) reconstructed image corresponding to the hologram produced by the 360×360 FINCH system; (h) reconstructed image corresponding to the hologram produced by synthetic aperture of double lens FINCH; (i) reconstructed image corresponding to the hologram produced by the 1080×1080 FINCH system[275]
    Experimental results of SPIDER imaging[277]. (a) PIC image experimental platform; (b) iterative image reconstruction result of Fig.54(g); (c)(g) two images of the target; (d)(h) corresponding principle simulation results of target image; (e)(i) imaging results obtained by inverse Fourier transform reconstruction; (f)(j) after correcting the swing error of the turntable, the imaging results are reconstructed by inverse Fourier transform
    Lensless noninterference coded aperture dependent holography[116]. (a) Two LEDs; (b) two one-dime coins
    Incoherent lensless imaging based on Fresnel region aperture. (a) Real time image capture and reconstruction of lens less camera[279]; (b) binary, gray, and color images are reconstructed by Fresnel region aperture single frame lensless camera[282]
    • Table 1. Coherence measurement: classical coherence theory and phase space optics theory

      View table

      Table 1. Coherence measurement: classical coherence theory and phase space optics theory

      TheoryFunctionDefinitionTemporal/Spatial coherence
      Classical coherence theoryMutual coherence functionΓ(x1,x2,τ)=<U(x1,t)U*(x2,t+τ)>Temporal and spatial
      Complex degree of coherenceγ(x1,x2,τ)=Γ(x1,x2,τ)Γ(x1,x1,0)Γ(x2,x2,0)1/2
      Cross-spectral density functionW(x1,x2,ω)=∫Γ(x1,x2,τ)exp(2πiντ)dτ
      Self-coherence functionΓ(x,τ)=<U(x,t)U*(x,t+τ)>Note: I(x)=Γ(x,0)Temporal
      Self complex degree of coherence functionγ(x,τ)=Γ(x,τ)Γ(x,0)
      Mutual intensityJ(x1,x2)≡Γ(x1,x2,0)=<U(x1,t)U*(x2,t)>Spatial quasi-monochromatic
      Complex coherence factorj(x1,x2)≡γ(x1,x2,0)=J(x1,x2)J(x1,x1)J(x2,x2)1/2
      Phase space optics theoryWigner distribution functionW(x,u)=∫Wx+x'2,x-x'2exp(-j2πux')dx'=∫Γu+u'2,u-u'2exp(j2πxu')du'Spatial
      Ambiguity functionA(u',x')=∫Wx+x'2,x-x'2exp(-j2πux')dx=∫Γu+u'2,u-u'2exp(j2πux')du
    • Table 2. Properties of Wigner distribution function

      View table

      Table 2. Properties of Wigner distribution function

      PropertyRepresentationExplanation
      RealnessW(x,u)∈ℝW is always a real function
      Spatial marginal propertyI(x)=∫W(x,u)duI(x) is the intensity
      Spatial frequency marginal propertyS(u)=∫W(x,u)dxS(u) is the power spectrum
      Convolution propertyU(x)=U1(x)U2(x) W(x,u)=W1(x,u)uW2(x,u)U(x)=U1(x)xU2(x) W(x,u)=W1(x,u)xW2(x,u)x is the convolution over xu is the convolution over u
      Instantaneous frequencyuW(x,u)duW(x,u)du=12πÑϕ(x)ϕ(x) is the phase component Ñϕ(x) is the instantaneous frequency
    • Table 3. Common optical transformation of Wigner distribution function

      View table

      Table 3. Common optical transformation of Wigner distribution function

      Optical transformationRepresentationExplanation
      Fresnel diffractionWz(x,u)=W0(x-λzu,u)λ is the wavelengthz is diffraction distance
      Chirp modulation (lens)W(x,u)=W0x,u+xλfλ is the wavelengthf is the focal length of lens
      Fourier transform(Fraunhofer diffraction)WÛ(x,u)=WU(-u,x)Û is the Fourier transform of signal
      Fractional Fourier transformWÛθ(x,u)=WU(xcos θ-usin θ,ucos θ+xsin θ)Ûθ is the fractional Fourier transform, θ is the rotation angle
      Beam amplifier (compressor)W(x,u)=W0(x,u/M)M is the magnification
      First order optical systemx'u'=ABCDxuA,B,C,D corresponding to first order optical system
    • Table 4. Spatial and phase space characterization of common optical signals

      View table

      Table 4. Spatial and phase space characterization of common optical signals

      Optical signalSpatial representationPhase space representationExplanation
      Point sourceU(x)=δ(x-x0)W(x,u)=δ(x-x0)Line perpendicular to the x-axis in phase space
      Plane waveU(x)=exp(i2πu0x)W(x,u)=δ(u-u0)Line perpendicular to the u-axis in phase space
      Spherical waveU(x)=exp(i2πax2)W(x,u)=δ(u-ax)Straight line across the origin of phase space
      Slow-varying waveU(x)=A(x)expiϕ(x)W(x,u)≈I(x)δu-12πϕCurve in phase space
      Gaussian signalU(x)=exp-πσ2(x-x0)2W(x,u)=exp-πσ2(x-x0)2+σ22πu22D Gaussian function in phase space
      Spatially incoherent fieldWx+x'2,x-x'2=I(x)δ(x')W(x,u)=cI(x)c is a constant only related to x
      Spatially stationary fieldWx+x'2,x-x'2=I0μ(x')W(x,u)=cμ̅(u)μ̅(u) is the Fourier transform of μ(x')
      Quasi-homogeneous fieldWx+x'2,x-x'2I(x)μ(x')W(x,u)≈I(x)μ̅(u)I is a relatively slow-varying signal compared with μ
    • Table 5. Optical field transmission: from coherent to partially coherent

      View table

      Table 5. Optical field transmission: from coherent to partially coherent

      CoherenceCoherentPartially coherent
      RepresentationU(x,z)W(x1,x2)W(x,u)
      Wave equation(2+k2)U(x,z)=0x12W(x1,x2)+k2W(x1,x2)=0x22W(x1,x2)+k2W(x1,x2)=0k2-4π2|u|2kW(x,u)z=-λuxW(x,u)
      Spatial convolutionUz(x)=exp(jkz)jλz∫U0(x0expλzx-x02dx0Wz(x1,x2)=W0(x1,x2)x1,x2hz(x1,x2)Wz(x,u)=W0(x,u)xWhz(x,u)
      Angular spectrumHFz(ux,uy)=exp(jkzexp-λz(ux2+uy2)Ŵz(u1,u2)=Ŵ0(u1,u2)Hz(u1,u2)
      Transport of intensity equation-kI(x)z=Ñ·I(x,z)φ(x)zWx+x'2,x-x'2=-jkxx'Wx+x'2,x-x'2I(x)z=-λx·uWω(x,u)du
    Tools

    Get Citation

    Copy Citation Text

    Runnan Zhang, Zewei Cai, Jiasong Sun, Linpeng Lu, Haitao Guan, Yan Hu, Bowen Wang, Ning Zhou, Qian Chen, Chao Zuo. Optical-Field Coherence Measurement and Its Applications in Computational Imaging[J]. Laser & Optoelectronics Progress, 2021, 58(18): 1811003

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Imaging Systems

    Received: Jul. 14, 2021

    Accepted: Aug. 10, 2021

    Published Online: Sep. 2, 2021

    The Author Email: Cai Zewei (zeweicai@njust.edu.cn), Chen Qian (chenqian@njust.edu.cn), Zuo Chao (zuochao@njust.edu.cn)

    DOI:10.3788/LOP202158.1811003

    Topics