Photonics Research, Volume. 9, Issue 6, 1003(2021)

Non-iterative complex wave-field reconstruction based on Kramers–Kronig relations

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
    Figures & Tables(16)
    Principle of KKSAI. (a) Schematic of experimental setup, where pupil modulation is achieved by an SLM-based module. (b) Simplified 4f system corresponding to (a). (c) Simulated complex-valued sample. (d) Amplitude pupil modulation indicated by the green circle, whose center is (ui,vi). (e) Measured images corresponding to (d). (f) KKSAI reconstruction algorithm flowchart. It finally recovers the pupil-limited sample spectrum. MMF, multimode fiber; CL, collimating lens; RL, relay lens; P, polarizer; BS, beam splitter; SLM, spatial light modulator; TL, tube lens.
    Analogy between KKSAI measurement and off-axis hologram. (a) Measurement I1 between its corresponding complex-valued spectrum subregion and the amplitude of its FT. (b) Shifted spectrum subregion still brings in the same measurement due to the frequency shifting property of FT and phase loss of the square-law detector. (c) The shifted spectrum subregion can be hypothetically decomposed into a Dirac delta function and the shifted scattered complex-valued function.
    Titchmarsh theorem applied to a band-limited signal. (a) Amplitude and (b) phase of s˜(l→) with bandwidth of ρNAape. (c) Logarithm of its 2D Fourier amplitude spectrum. (d) Logarithm of its 1D Fourier amplitude spectrum along the l∥ axis and its shifted copies by (e) |ρ→r|<ρNAape and (f) |ρ→r|=ρNAape.
    Scanning scheme examples to cover the entire pupil. (a) Four circular apertures; (b) two rectangular apertures. The circled numbers are used to label the measurement sequence.
    Reconstruction of phase-only samples by two existing methods and our proposed method. (a) Weak phase sample; (b) strong phase sample.
    Reconstruction of complex-valued sample by two existing methods and our proposed method. (a) Phase; (b) amplitude.
    Effect of distance between aperture edge and pupil center on the final reconstruction accuracy.
    KKSAI based on the scanning scheme with only two measurements. (a) Measurements; (b) reconstruction.
    Experimental results for a microlens array. (a) Reconstructions of a single microlens by PM-DPC, PM-FPM, and KKSAI using four measurements. (b) A close-up view of the SEM image of the microlens array (adapted from Thorlabs website). (c) Radial average profile of three phase recoveries compared with the ground truth (GT).
    Experimental results for a thyroid carcinoma pap smear slide. (a) Two out of four measurements acquired by KKSAI and their Fourier amplitude spectrum. (b) Amplitude reconstruction by PM-FPM and KKSAI compared with ground truth. (c) Phase reconstruction by PM-DPC, PM-FPM, and KKSAI compared with ground truth.
    Chromatic aberration correction by digital refocusing ability of KKSAI. (a) Reconstructed amplitudes by KKSAI of three channels. (b) Reconstructed phases by KKSAI of three channels. (c) Digitally refocused amplitudes with the corresponding refocusing distance labeled at the bottom. (d) Color composite of three channels before and after digital refocusing with the enlargements showing improved image quality. R, red (638 nm); G, green (532 nm); B, blue (405 nm).
    Complex wave-field reconstruction by KKSAI based on only two measurements. (a) Scanning scheme, raw measurements, and their spectrum amplitude; (b) and (c) reconstructed amplitude and phase, respectively, by KKSAI from two measurements, four measurements, and PM-FPM with 47 measurements. Here PM-FPM reconstruction is taken as the reference to calculate the FSIM metric for KKSAI reconstruction.
    • Table 1. Quantitative Evaluation of Reconstructions in Fig. 5

      View table
      View in Article

      Table 1. Quantitative Evaluation of Reconstructions in Fig. 5

      MetricPM-DPCPM-FPMKKSAI
      (a)MSE4.80×1094.53×1097.84×1010
      FSIM0.99991.00000.9998
      Time (s)2.0128.432.12
      (b)MSE0.17110.06400.0136
      FSIM0.98941.00000.9973
      Time (s)1.6777.332.03
    • Table 2. Quantitative Evaluation of Reconstructions in Fig. 6

      View table
      View in Article

      Table 2. Quantitative Evaluation of Reconstructions in Fig. 6

      MetricPM-DPCPM-FPMKKSAI
      PhaseMSE0.05310.01200.0037
      FSIM0.99340.99970.9976
      AmplitudeMSE/1.55×1083.78×104
      FSIM/1.00000.9965
      Time (s)3.83109.372.82
    • Table 3. Similarity Evaluation of Overlapping Spectrum Regions in Fig. 6

      View table
      View in Article

      Table 3. Similarity Evaluation of Overlapping Spectrum Regions in Fig. 6

      Overlapping Region
      Phase0.99400.99670.99900.9926
      Amplitude0.99970.99980.99990.9998
    • Table 4. Similarity Evaluation of Overlapping Spectrum Regions in Fig. 10

      View table
      View in Article

      Table 4. Similarity Evaluation of Overlapping Spectrum Regions in Fig. 10

      Overlapping Region
      Phase0.98720.99270.98640.9691
      Amplitude0.99910.99980.99860.9991
    Tools

    Get Citation

    Copy Citation Text

    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

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

    Category: Imaging Systems, Microscopy, and Displays

    Received: Jan. 19, 2021

    Accepted: Mar. 23, 2021

    Published Online: May. 27, 2021

    The Author Email: Cheng Shen (cshen3@caltech.edu)

    DOI:10.1364/PRJ.419886

    Topics