# Two-Step Phase Retrieval Algorithm Using Single-Intensity Measurement

Authors: Cheng Zhang; Meiqin Wang; Qianwen Chen; Dong Wang; Sui Wei
Journal: International Journal of Optics
 (2018)

---

## Abstract

Aiming at the problem that the single-intensity phase retrieval method has poor reconstruction quality and low probability of successful recovery, an improved method is proposed in this paper. Our method divides the phase retrieval into two steps: firstly, the GS algorithm is used to recover the amplitude in the spatial domain from the single-spread Fourier spectrum, and then the classical GS algorithm using two intensity measurements (one is recorded and the other is estimated from the first step) measurements is used to recover the phase. Finally, the effectiveness of the proposed method is verified by numerical experiments.


---

## Body

## 1. Introduction

Most of the information in the optical field is contained in the phase, such as the depth and shape of the object and so on. The oscillation frequency of the optical field is as high as 1015 Hz. However, existing detection devices can only record intensity, directly [1]. Phase retrieval method recovers lost phase from recorded intensity measurements with some prior knowledge, which plays an important role in the imaging field. Phase retrieval problem arises in numerous of areas, such as crystallography [2–4], optics [5–9], astronomical imaging [10], microscopy [11, 12], biomedical [13], and holographic imaging[14, 15].In 1972, Gerchberg and Saxton first proposed an alternate projection-based phase retrieval algorithm, the Gerchberg-Saxton (GS) algorithm [16]. The main idea of GS algorithm is to use the intensity of the spatial domain and the Fourier domain to recover the phase of the optical field. Subsequently, Fienup proved that the GS algorithm had obvious error-decreasing properties, Error Reduction (ER) algorithm, and Hybrid Input-Output (HIO) algorithm were proposed [17]. At present, ER algorithm and HIO algorithm are considered to be the most effective methods in the field of phase retrieval [18]. Since the above algorithms are only suitable for a positive linear transformation system, nobody cares about any linear transformation system. Therefore, Yang and Gu proposed the amplitude-phase detection theory in arbitrary linear transformation systems, namely, the Yang-Gu (Y-G) algorithm [19]. In 2015, Guo et al. optimized the iterative algorithm and proposed two improved GS iterative phase retrieval algorithms, which were the spatial phase perturbation GS algorithm and the combined GS HIO algorithm. For the two improved algorithms, the squared value of the squared error value rapidly drops to an acceptable value, and the lost phase can be successfully recovered in the spatial domain and Fourier domain, which means that both algorithms can jump out of the local minimum value and convergence to the global minimum [9].The GS algorithm was originally proposed in connection with the problem of reconstructing phase given only two-intensity measurements in the spatial domain and Fourier domain, i.e., given two intensity measurements (one is recorded in the spatial domain and the other is estimated in the Fourier domain) measurements. Unfortunately, such two intensity measurements at two different planes cannot be measured in some cases. Therefore, a GS algorithm using single-intensity measurement with some prior knowledge combined was proposed to recover the phase. However the single-intensity phase retrieval method has the drawbacks of poor reconstruction quality and low probability of successful recovery. This paper proposes an improved method based on single-intensity measurement, Two-Step Phase Retrieval (TSPR) from single-intensity measurement algorithm. The main idea of the TSPR algorithm is to first recover the amplitude (i.e., the square root of the intensity) in the spatial domain using the Single-Intensity Phase Retrieval (SIPR) algorithm with some prior knowledge and then recover the phase in the spatial domain using the Two-Intensity Phase Retrieval (TIPR) algorithm. Experimental results demonstrate that our proposed TSPR algorithm can effectively enhance the quality of phase recovered and the probability of successful recovery.


## 2. Basic Principles

The complex amplitudeg(x,y) is(1)gx,y=gx,yexp⁡iφx,yIts Fourier transform form is(2)Gξ,η=Gξ,ηexp⁡iψξ,ηwhere g(x,y) and Gξ,η denote the complex amplitude in the spatial domain and Fourier domain, respectively, gx,y and Gξ,η denote the amplitude in the spatial domain and Fourier domain, respectively, φx,y and ψξ,η denote the phase in the spatial domain and Fourier domain, respectively. Phase retrieval refers to recovering the losing phase φx,y in the spatial domain from two amplitudes gx,y and Gξ,η.The flow chart of the classical GS algorithm is shown in Figure1. This algorithm was first proposed and was easily used to solve the phase retrieval problem using two intensity measurements. exp⁡(iϕ0(x,y)) denotes the initialized phase distribution function. The algorithm repeatedly iterates between the spatial domain and the Fourier domain until the error ε satisfies the termination condition [20].Figure 1
Flow chart of classical GS algorithms using two intensity measurements, i.e., TIPR.


## 3. Two-Step Phase Retrieval Using Single-Intensity Measurement

The TIPR algorithm uses the two intensity measurements recorded in the spatial domain and the Fourier domain, respectively, to recover the amplitude and phase in the spatial domain with high quality. Although the TIPR algorithm can achieve good results, in some case, it is difficult to obtain two intensity measurements simultaneously in the both spatial and Fourier domains.Therefore, the SIPR algorithm is proposed, and it can recover the phase from single-intensity measurement in the Fourier domain with some prior knowledge. However, the SIPR method has low probability of successful recovery and poor quality. As shown in Figure2, the purpose of the experiment is to test the SIPR method’s performance to recover the amplitude and phase. The coded aperture is generated by a uniform distribution. The sampling rate of the coded aperture is used as a variable parameter in the experiment. The sampling rate gradually increased from 0.1 to 0.7; the step length is 0.05. The experiment runs 500 times independently to compute the probability of successful recovery. If the SNR of the phase reconstructed is greater than 25 dB, the reconstruction is deemed as successful. The probability of successful recovery versus sampling rate for both amplitude and phase is shown in Figure 2.Figure 2
Probability of successful recovery of amplitude and phase versus sampling rate.As depicted in Figure2, it can be found that the SIPR method can achieve high probability of successful recovery of the amplitude, but the probability of successful recovery of the phase is very low and unstable. Hence, our TSPR method is proposed to solve the phase retrieval problem from single-intensity measurement. Firstly, SIPR method is used to recover the amplitude in the spatial domain from the intensity in the Fourier domain and the support set in the spatial domain. Then, TIPR method using two intensity measurements (one of the two intensity measurements is recorded in the Fourier domain and the other is estimated in the spatial domain from the first step) is used to recover the phase in the spatial domain.The flow chart of TSPR method is shown in Figure3, which includes the following steps: (1) initialize the amplitude A1 and phase φ0(x,y) in the spatial domain, i.e., use all-one amplitude and randomized phase. (2) Perform Fourier transform on the complex amplitudes and replace the amplitude after the Fourier transform in the Fourier domain with the recorded Fourier amplitudes A2. (3) Perform the inverse Fourier transform on the synthesized complex amplitude in the Fourier domain and multiply it with the support in the spatial domain. (4) Repeat steps (2)-(3) until the amplitude converges in the spatial domain. (5) The classical TIPR algorithm using the recorded Fourier amplitude A2 and the estimated amplitude gest(x,y) from the first step is used to iteratively recover the phase in the spatial domain.Figure 3
Flow diagram of the TSPR algorithm.The complex amplitude of the initialization is as follows:(3)gx,y=A1.∗expφ0x,yThe optical setup of imaging process of the TSPR method is shown in Figure4. The coded aperture and CCD are placed on the front and back focal planes of the lens, respectively.Figure 4
The optical setup of imaging process.Firstly, the complex optical fieldux,y is filtered through the coded aperture M to obtain a complex optical field uMx,y and then imaged by a Fourier lens; finally the intensity is recorded in the CCD plane. The intensity recorded at the CCD plane is expressed as follows: (4)ICCD=FuMx,y2The main purpose of this paper is to recover the phase of the complex optical fielduMx,y from the recorded intensity measurement ICCD.


## 4. Numerical Experiments

In order to verify the effectiveness and superiority of our proposed method, three experiments are presented in this section. The first simulation experiment is the single reconstruction experiment. The second simulation experiment is to test the performance of TSPR algorithm under different coded apertures. The final simulation experiment is to compare the performance of the different phase retrieval algorithms.


### 4.1. Experiment 1: Single Reconstruction Experiment

The purpose of the experiment is to verify the feasibility of our proposed TSPR method. Two grayscale images (“Lena” and “Cameraman” with256×256 pixels) are chosen in the experiment. The two images are multiplied by the coded aperture (256×256 pixels) to obtain the amplitude and phase, respectively. Then the complex amplitude is synthesized by amplitude and phase. The coded aperture is 0/1 random distribution. The sampling rate in this experiment is 0.4, and the phase retrieval results are shown in Figure 5.Figure 5
Comparison of SIPR algorithm and TSPR algorithm. (a) A 0/1 randomly distributed coded aperture. (b) and (c) are the amplitude and phase of the complex optical field behind the coded aperture. (d) The intensity recorded in the Fourier domain. (e) and (f) are the amplitude and phase recovered using SIPR algorithm. (g) and (h) are the amplitude and phase recovered using TSPR algorithm.
(a) (b) (c) (d) (e) (f) (g) (h)Figure5(a) is the coded aperture with 0/1 randomly distributed, and the sampling rate is 0.4. Figures 5(b) and 5(c) are the amplitude and phase of the complex optical field behind the coded aperture to be retrieved. Figure 5(d) is the intensity recorded in the Fourier domain. Figures 5(e) and 5(f) are the amplitude and phase recovered using the SIPR algorithm, and the Signal-to-Noise Ratio (SNR) is 343.04 dB and 19.85 dB, respectively. Figures 5(g) and 5(h) are the amplitude and phase recovered using the TSPR algorithm, and their SNR is 343.04dB and 45.58dB, respectively. Comparing Figures 5(f)–5(h), it is easy to find that the quality of phase recovered is greatly improved by our TSPR method.The trend of the amplitude and phase with the number of iterations in our TSPR method is shown in Figure6. In the first step (i.e., SIPR) of TSPR, the amplitude can be well recovered, but the phase recovered is poor. As the number of iterations increases, the SNR of the phase recovered can only reach about 19dB. However, in the second step (i.e., TIPR) of the TSPR, the SNR of the phase recovered increases as high as 45 dB. Compared with SIPR method, TSPR method can greatly improve the SNR of phase reconstructed.Figure 6
Amplitude and phase versus the number of iterations.


### 4.2. Experiment 2: Reconstruction Experiments under Different Coded Apertures

The coded aperture in Experiment 1 is a uniform random sampling pattern. The purpose of Experiment 2 is to test the performance of the TSPR method under different coded apertures.In Figure7, three different types of the coded apertures are selected, such as Uniform Random (UR) sampling, Radial Line (RL) sampling, and Variable Density (VD) sampling. The sampling rate of the coded aperture is used as a variable parameter in the experiment, and the sampling rate is set from 0.1 to 0.7. Each set of experiments with the same parameters runs 500 times independently. Compute the recovered Average Signal-to-Noise Ratio (ASNR) and probability of successful recovery. The ASNR and probability of successful recovery versus sampling rate under different coded apertures are shown in Figures 8 and 9, respectively.Figure 7
Three types of coded apertures with sampling rate of 0.3. (a) Uniform random sampling pattern. (b) Radial line sampling pattern. (c) Variable density sampling pattern.
(a) (b) (c)Figure 8
ASNR versus sampling rate.Figure 9
Probability of successful recovery versus sampling rate.As shown in Figure8, as the sampling rate increases, the TSPR algorithm achieves higher ASNR for phase recovered than the SIPR algorithm under three different coded apertures. In Figure 9, the best choice for coded aperture is uniform random sampling pattern, which achieve the optimal performance. The probability of successful recovery has also been significantly improved and has a good stability. The reason is that the support set of the uniform random sampling pattern is relatively uniform distributed in the spatial domain, while the support set of the other two coded apertures is more concentrated to the center, so the phase cannot be retrieved well.


### 4.3. Experiment 3: Reconstruction Experiment under Different Sampling Rates

The purpose of Experiment 3 is to compare the performance of phase retrieval between our TSPR method and the other two algorithms (TIPR algorithm and SIPR algorithm) under different sampling rates. The coded aperture uses a uniform random sampling pattern with a sampling rate varies from 0.1 to 0.7 and the step size is 0.05. Under the same parameters, experiments runs 500 times independently. The ASNR and probability of successful recovery versus sampling rate for different algorithms are presented in Figures10 and 11, respectively.Figure 10
ASNR versus sampling rate.Figure 11
Probability of success versus sampling rate.As shown in Figures10 and 11, the TIPR method is the best, the TSPR method is the second, and the SIPR method is the worst. The reason is that TIPR method has the most known information, i.e., the two intensity measurements at two different planes. The difference between the TSPR method and the SIPR method is that the TSPR method draws on the idea of the TIPR method. Firstly, use the SIPR method to estimate the amplitude in the spatial domain, then use the TIPR method to recover the phase. In summary, the TSPR method is superior to the SIPR method.


## 4.1. Experiment 1: Single Reconstruction Experiment

The purpose of the experiment is to verify the feasibility of our proposed TSPR method. Two grayscale images (“Lena” and “Cameraman” with256×256 pixels) are chosen in the experiment. The two images are multiplied by the coded aperture (256×256 pixels) to obtain the amplitude and phase, respectively. Then the complex amplitude is synthesized by amplitude and phase. The coded aperture is 0/1 random distribution. The sampling rate in this experiment is 0.4, and the phase retrieval results are shown in Figure 5.Figure 5
Comparison of SIPR algorithm and TSPR algorithm. (a) A 0/1 randomly distributed coded aperture. (b) and (c) are the amplitude and phase of the complex optical field behind the coded aperture. (d) The intensity recorded in the Fourier domain. (e) and (f) are the amplitude and phase recovered using SIPR algorithm. (g) and (h) are the amplitude and phase recovered using TSPR algorithm.
(a) (b) (c) (d) (e) (f) (g) (h)Figure5(a) is the coded aperture with 0/1 randomly distributed, and the sampling rate is 0.4. Figures 5(b) and 5(c) are the amplitude and phase of the complex optical field behind the coded aperture to be retrieved. Figure 5(d) is the intensity recorded in the Fourier domain. Figures 5(e) and 5(f) are the amplitude and phase recovered using the SIPR algorithm, and the Signal-to-Noise Ratio (SNR) is 343.04 dB and 19.85 dB, respectively. Figures 5(g) and 5(h) are the amplitude and phase recovered using the TSPR algorithm, and their SNR is 343.04dB and 45.58dB, respectively. Comparing Figures 5(f)–5(h), it is easy to find that the quality of phase recovered is greatly improved by our TSPR method.The trend of the amplitude and phase with the number of iterations in our TSPR method is shown in Figure6. In the first step (i.e., SIPR) of TSPR, the amplitude can be well recovered, but the phase recovered is poor. As the number of iterations increases, the SNR of the phase recovered can only reach about 19dB. However, in the second step (i.e., TIPR) of the TSPR, the SNR of the phase recovered increases as high as 45 dB. Compared with SIPR method, TSPR method can greatly improve the SNR of phase reconstructed.Figure 6
Amplitude and phase versus the number of iterations.


## 4.2. Experiment 2: Reconstruction Experiments under Different Coded Apertures

The coded aperture in Experiment 1 is a uniform random sampling pattern. The purpose of Experiment 2 is to test the performance of the TSPR method under different coded apertures.In Figure7, three different types of the coded apertures are selected, such as Uniform Random (UR) sampling, Radial Line (RL) sampling, and Variable Density (VD) sampling. The sampling rate of the coded aperture is used as a variable parameter in the experiment, and the sampling rate is set from 0.1 to 0.7. Each set of experiments with the same parameters runs 500 times independently. Compute the recovered Average Signal-to-Noise Ratio (ASNR) and probability of successful recovery. The ASNR and probability of successful recovery versus sampling rate under different coded apertures are shown in Figures 8 and 9, respectively.Figure 7
Three types of coded apertures with sampling rate of 0.3. (a) Uniform random sampling pattern. (b) Radial line sampling pattern. (c) Variable density sampling pattern.
(a) (b) (c)Figure 8
ASNR versus sampling rate.Figure 9
Probability of successful recovery versus sampling rate.As shown in Figure8, as the sampling rate increases, the TSPR algorithm achieves higher ASNR for phase recovered than the SIPR algorithm under three different coded apertures. In Figure 9, the best choice for coded aperture is uniform random sampling pattern, which achieve the optimal performance. The probability of successful recovery has also been significantly improved and has a good stability. The reason is that the support set of the uniform random sampling pattern is relatively uniform distributed in the spatial domain, while the support set of the other two coded apertures is more concentrated to the center, so the phase cannot be retrieved well.


## 4.3. Experiment 3: Reconstruction Experiment under Different Sampling Rates

The purpose of Experiment 3 is to compare the performance of phase retrieval between our TSPR method and the other two algorithms (TIPR algorithm and SIPR algorithm) under different sampling rates. The coded aperture uses a uniform random sampling pattern with a sampling rate varies from 0.1 to 0.7 and the step size is 0.05. Under the same parameters, experiments runs 500 times independently. The ASNR and probability of successful recovery versus sampling rate for different algorithms are presented in Figures10 and 11, respectively.Figure 10
ASNR versus sampling rate.Figure 11
Probability of success versus sampling rate.As shown in Figures10 and 11, the TIPR method is the best, the TSPR method is the second, and the SIPR method is the worst. The reason is that TIPR method has the most known information, i.e., the two intensity measurements at two different planes. The difference between the TSPR method and the SIPR method is that the TSPR method draws on the idea of the TIPR method. Firstly, use the SIPR method to estimate the amplitude in the spatial domain, then use the TIPR method to recover the phase. In summary, the TSPR method is superior to the SIPR method.


## 5. Conclusion

This paper proposes an improved TSPR algorithm to solve the problem of poor quality and low success rate of SIPR method. TSPR algorithm is completed in two steps: firstly, SIPR algorithm is used to recover the amplitude in the spatial domain, then TIPR algorithm is used to recover the phase in the spatial domain. Finally, the effectiveness and superiority of the TSPR method are verified. The experiment comparison results demonstrate that the proposed method can effectively improve the ASNR and probability of successful recovery under the same parameters. TSPR algorithm can recover the lost phase from known amplitude in the Fourier domain and support set in the spatial domain.



---

*Source: 8643819-2018-10-01.xml*