Sparse modeling is one of the efficient techniques for imaging that allows recovering lost information. In this paper, we present a novel iterative phase-retrieval algorithm using a sparse representation of the object amplitude and phase. The algorithm is derived in terms of a constrained maximum likelihood, where the wave field reconstruction is performed using a number of noisy intensity-only observations with a zero-mean additive Gaussian noise. The developed algorithm enables the optimal solution for the object wave field reconstruction. Our goal is an improvement of the reconstruction quality with respect to the conventional algorithms. Sparse regularization results in advanced reconstruction accuracy, and numerical simulations demonstrate significant enhancement of imaging.
The conventional sensors detect only the intensity of the light, but the phase is systematically lost in measurements. Phase retrieval is a problem of the phase recovering using a number of intensity observations and some prior on the object wave field. The phase carries important information about an object shape what is necessary for a 3D object imaging and exploited in many areas such as microscopy, astronomy, etc. Moreover, phase-retrieval techniques are often simpler, cheaper and more robust comparing with interferometric ones.
In 1982 Fienup introduced some, for now classical, iterative phase-retrieval algorithms [1]: error-reduction, gradient search and input-output methods. Many phase-retrieval methods are developed based on this pioneer work: the estimated magnitudes at the measurement planes are iteratively replaced by ones obtained from the intensity observations. We are looking for the optimal wave field reconstruction from a number of intensity observations, and the reconstruction problem is formulated in terms of a variational constrained maximum likelihood (ML) approach. The spatial image resolution of the conventional phase-retrieval techniques is limited due to diffraction, what can be one of the main sources of artifacts and image degradations. In order to enhance the imaging quality and recover lost information, in this work we use the novel compressive sensing technique for the variational image reconstruction originated in [2]. The object wave field distribution is assumed to be sparse, and its amplitude and phase are separately analyzed and decomposed using very specific basis functions named BM3D-frames [3]. The proposed phase-retrieval algorithm is derived as a solution of the ML optimization problem using the BM3D-frame based sparse approximation of the object amplitude and phase distributions.
We consider a multi-plane wave field reconstruction scenario: a planar laser beam illuminates an object, and the result of the wave field propagation is detected on a sensor at different distances z r from the object at various measurement (sensor) planes. Here z r = z 1 +(r-1)• z , r=1,…K, where z 1 is the distance from the object to the first measurement plane, z is the distance between the measurement planes, and K is a number of these planes. We assume that the wave field distributions at the object and sensor planes are pixel-wise invariant. In such a discrete-to-discrete model, the forward wave field propagation from the object to the r-th sensor plane is defined in the vector-matrix form as follows: 0 = , 1,… , rr rK u Au (1) where u 0 and u r are n vectors, constructed by columns concatenating 2D discrete complexvalued wave field distributions (N×M matrices) at the object and sensor planes, respectively.
n×n is a discrete forward propagation operator, n=N M. We consider the paraxial approximation of the wave field propagation defined by the Rayleigh-Sommerfield integral.
Depending of the used discretization of this integral, the operators A r in (1) can be, for instance, the angular spectrum decomposition or the discrete diffraction transform in the matrix (M-DDT, [4]) or the Fourier transform domains (F-DDT, [5]). In our numerical experiments, we use F-DDT models enabling the exact pixel-to-pixel mapping u 0 to u r .
According to the used vector-matrix notation, the observation model with the additive zeromean Gaussian noise at the sensor planes, r [k] N(0, r ²), takes the form:
Let us assume that the object amplitude a 0 n and phase 0 n can be separately approximated by small numbers of basic functions with coefficients a for the amplitude and coefficients for the phase. These basic functions are collected in the matrices a and for the amplitude and the phase, respectively. The amplitude and the phase are reconstructed from the noisy intensity data o r .
According to the maximum likelihood approach, the reconstruction of the object wave field is performed by minimization of the criterion:
subject to the following constraints: first of all the forward propagation models (1), and the constraints for the sparse modeling for the object amplitude and phase given as
The quadratic fidelity term in (3) appears due to our assumption that the observation noise is Gaussian. The next two terms define the sparse regularization in the spectral domain, where the positive parameters a and define a balance between the fit of observations, the smoothness of the wave field reconstruction and the complexity of the solution.
The first equations in (4) and (5) give the restrictions for the amplitude and the phase in the synthesis form and the second ones -in the analysis form [2,3]. Note that the priori unknown basic functions for the amplitude and the phase are selected from the synthesis a , and analysis a , matrices defining the redundant sets of the basic functions. The vectors a m are considered as spectra for parametric approximations of the amplitude and the phase. Conventionally, the
This content is AI-processed based on open access ArXiv data.
