OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 17 — Aug. 17, 2009
  • pp: 15257–15263
« Show journal navigation

Linear phase retrieval with a single far-field image based on Zernike polynomials

Min Li and Xin-Yang Li  »View Author Affiliations


Optics Express, Vol. 17, Issue 17, pp. 15257-15263 (2009)
http://dx.doi.org/10.1364/OE.17.015257


View Full Text Article

Acrobat PDF (168 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Wavefront aberrations can be represented accurately by a number of Zernike polynomials. We develop a method to retrieve small-phase aberrations from a single far-field image with a Zernike modal-based approach. The difference between a single measured image with aberration and a calibrated image with inherent aberration is used in the calculation process. In this paper, the principle of linear phase retrieval is introduced in a vector–matrix format, which is a kind of linear calculation and is suitable for real-time calculation. The results of numerical simulations on atmosphere-disturbed phase aberrations show that the proposed Zernike modal-based linear phase retrieval method works well when the rms of phase error is less than 1 rad, and it is valid in a noise condition when the signal-to-noise ratio (SNR)>3.

© 2009 Optical Society of America

1. Introduction

The purpose of phase retrieval is to obtain the unknown phase from far-field image data. It has been proposed for a variety of applications and has received much attention over many years. Solutions to the problem include measuring the modulus in one or two conjugate domains, using a single intensity measurement plus a non-negativity constraint, using two defocused images, and so on [1

1. J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Eng. 21, 2758–2769 (1982).

6

6. R. A. Gonsalves, “Small-phase solution to the phase-retrieval problem,” Opt. Lett. 26, 684–685 (2001).

].

R. A. Gonsalves [6

6. R. A. Gonsalves, “Small-phase solution to the phase-retrieval problem,” Opt. Lett. 26, 684–685 (2001).

] proposed an odd–even function decomposition method that can be used to provide evidence for the solution to the phase retrieval problem. Based on this idea, Min Li et al. found that in a small-phase condition, the odd and even parts of a phase aberration can be obtained with a simple linear calculation method. The difference between a single measured image with aberration and a calibrated image with inherent aberration was used in the calculation process. But it was impractical due to the point-to-point calculations. Its calculative quantity was huge and it was too sensitive to noise [7

7. M. Li, X.-Y. Li, and W.-H. Jiang. “Small-phase retrieval with a single far-field image,” Opt. Express 16, 8190–8197 (2008). [PubMed]

].

Generally, a phase can be accurately represented by a number of orthogonal bases such as Zernike polynomials. With the Hartmann–Shack wavefront sensor [8

8. L. Xinyang and J. Wenhan, “Zernike modal wavefront reconstruction error of Hartmann wavefront sensor,” Acta. Optica Sinica 22, 1236–1240 (2002).

], there is a linear relationship between the Zernike polynomials and the measure parameters; that is, a=R·m·a is the vector of Zernike polynomials, R is the linear matrix, and m is the vector of the measure parameter. This linear phase reconstruction is favorable for actual application in wavefront sensing.

2. Theory

The principle of linear phase retrieval is shown in Fig. 1 [7

7. M. Li, X.-Y. Li, and W.-H. Jiang. “Small-phase retrieval with a single far-field image,” Opt. Express 16, 8190–8197 (2008). [PubMed]

]. It is performed on a conventional imaging system. An ideal beam source is used first to calibrate inherent aberration B of the imaging system itself. This calibrated far-field image intensity P serves as the standard.

Fig. 1. Principal diagram of the linear phase retrieval method.

After getting the calibrated image, the same imaging system is used to measure the beam with small disturbed phase W. We break W into two parts:

W=We+Wo,
(1)

where We and Wo are the even and odd parts of W, respectively. It is based on the theory of odd–even function decompositions. Any function f(x,y) can be decomposed uniquely into the sum of an even and an odd function:

f(x,y)=fo(x,y)+fe(x,y),
(2)

where fo(x,y)=f(x,y)f(x,y)2,fe(x,y)=f(x,y)+f(x,y)2..

In the same way, the vector–matrix description of an image f(x,y) can be decomposed as

f=fJxyf2+f+Jxyf2=fo+fe,
(3)

in which f, fo, and fe are the vector formats of functions f(x, y), fo(x, y), and fe(x, y); and Jxy is the transform matrix to make the x-direction and the y-direction flip over a matrix.

We also break inherent phase aberration B into two parts: B=B e+Bo. It is proved in [7

7. M. Li, X.-Y. Li, and W.-H. Jiang. “Small-phase retrieval with a single far-field image,” Opt. Express 16, 8190–8197 (2008). [PubMed]

] that if BeBo or Bo=0 and W and B are small, the odd part of disturbed phase W is proportional to the difference of odd part of intensities PBo and Po:

ŴoPBoPo.
(4)

The even part of an aberration is approximately proportional to the difference of the even part of intensities PBe and Pe,

ŴePBePe,
(5)

in which P is the calibrated far-field image intensity, P=Pe+Po, and PB is the far-field intensity with inherent system aberration and measuring aberrations PB=PBe+PBo.

The relationships can be described in a vector–matrix format:

ŵo=Ro·Δo,
(6)
ŴeRe·Δe.
(7)

ŵo and ŵe are the vector formats of Ŵo and Ŵe. R o is the response matrix between phase aberration ŵo and intensity Δ o o, which is the vector format of ΔI o; ΔI o=PBo-Po; and Re is the response matrix between phase aberration ŵe and intensity Δe, which is the vector format of ΔIe, ΔIe=PBe-Pe.

Then the phase aberration can be retrieved as

ŵ=ŵo+ŵeRoΔo+ReΔe=(RoRoJxy2+Re+ReJxy2)Δ=R·Δ,
(8)

in which R is the linear retrieval matrix between phase aberration ŵ and image difference Δ, Δ=Δo+Δe. Therefore, there is an inverse linear relationship between them:

Δ=R+·ŵ,
(9)

in which R + is the pseudo inverse of matrix R and describes the linear matrix between them as well.

We use the Zernike circle polynomials, which were defined by Noll [9

9. X.-Y. Li, M. Li, B. Chen, and W.-H. Jiang “A kind of novel linear phase retrieval wavefront sensor and its application in close-loop adaptive optics system,” in Proc. Sixth International Workshop on Adaptive Optics for Industry and Medicine, C. Dainty, ed. (Imperial College Press, 2007), pp. 212–218.

], as the orthonormal basis functions. So, the phase can be described as the sum of several Zernike modes :

W(x,y)=Σi=1pciMi(x,y),
(10)

in which P is the Zernike number, ci is the coefficient, and Mi(x, y) is a description of the Zernike polynomials. Or, in matrix format,

w=D·c,
(11)

in which D is the response matrix between w and Zernike polynomial c, and c is the vector format of ci. Equation (9) becomes

Δ=R+D·c=Z·c.
(12)

It is proved that there is a linear relationship between the change of Zernike coefficients and the changes of intensities, too. The matrix Z represents the relationship between them. Matrix Z can be obtained formally according to the parameters of the imaging system. Then we can reconstruct the Zernike coefficients by

c=Z+·Δ,
(13)

in which Z + is the pseudo inverse of matrix Z. The wavefront also can be reconstructed.

From Eqs. (12) and (13), a linear relationship has been established between the intensities and the Zernike coefficients of the aberration phase, and the aberration phase can be retrieved from intensities by using a linear calculation algorithm. Equations (12) and (13) are the principle of linear phase retrieval in the vector–matrix format, which is a kind of linear calculation and is suitable for real-time calculation.

3. Numerical Simulations

Numerical simulations are employed to analyze the performance of this linear phase retrieval algorithm. First, the inherent system phase aberration characteristic is analyzed. Second, the dynamic range of this method is tested without noise to judge how small the aberration must be to satisfy the method. Last, the performance of the method is tested with different levels of measurement noise added to the far-field images.

A series of phase aberrations are generated to test the performance of this linear phase retrieval method. The method proposed by N. Roddier is used to generate the phase screen with the Kolmogorov atmosphere spectrum, and a total of 65 Zernike polynomials are used [10

10. N. Roddier, “Atmospheric wavefront simulation using Zernike polynomials,” Opt. Eng. 29, 1174–1180 (1990).

]. The entrance pupil is a 64×64 grid and the far-field image plane is in 128×128 pixels. The far-field image is calculated using a fast Fourier transform (FFT) program. During the FFT process, the far field is expanded 7 times experientially, and only 128×128 pixels in the center area are used.

In order to compare the difference between the unknown phase W(x, y) and the estimated phase Ŵ(x, y), which are expressed by 65 Zernike polynomials, we define the error wavefront as

E(x,y)=W(x,y)Ŵ(x,y).
(14)

An error coefficient η, which is the ratio between the rms of the error wavefront E(x, y) and the unknown wavefront W(x, y), is used as one criterion to determine the validity of the phase retrieval method:

η=RMS[E(x,y)]RMS[W(x,y)].
(15)

If η<1, the retrieval effect is valid.

Another metric is the residual Strehl ratio (SRe), which is the ratio between the peak intensity of the far-field image produced by the error wavefront E(x, y) and the maximum intensity of the Airy spot. If the SRe is closer to 1, the performance of this method is better.

The performance of this linear phase retrieval algorithm depends on the inherent phase aberration of the imaging system. It is proved in [7

7. M. Li, X.-Y. Li, and W.-H. Jiang. “Small-phase retrieval with a single far-field image,” Opt. Express 16, 8190–8197 (2008). [PubMed]

] that most of the inherent phase aberration of the imaging system must be an even type, such as defocus, to keep the method working well. Here we have done some simulations to test the tolerance of the inherent system phase aberration. The results are shown in Fig. 2. The atmosphere-disturbed phase aberration is set as σ=0.4 rad and the main type of system aberration is defocus.

Fig. 2. Relationship between the error coefficient η and the proportion of coma to the total system aberration; the defocus is set as σ=4 rad.

From Fig. 2, if there is some other odd aberration such as coma in the system, the conclusion is that when the proportion of coma to the total system aberration is smaller than 0.1, it has faint effect on the retrieval result.

Fig. 3. Numerical simulation results of the linear phase retrieval method proposed in this paper without noise. (a) Relationship between the average Strehl ratio SRˉ of aberration and the average RMS of atmosphere disturbed aberration σ̄ ; (b) relationship between the average error coefficient η̄ and the average rms of atmosphere-disturbed aberration σ̄.

It can be seen from the results that with the increase of phase aberration, the retrieval effects become worse. When σ̄≥1 rad, the average residual Strehl ratio SReˉ<0.85 and the average error coefficient η̄>0.6. So we can conclude under the conditions of this paper that the valid dynamic range of this method is approximately σ<1 rad.

To examine the applicability of this method in a noise condition, we investigate the sensitivity of the method to noise. Random noise, which satisfies the Gaussian distribution, is added to the imaging plane. All the values of noise are positive. The SNR is defined as

SNR=P1σn,
(16)

where PI is the peak value of the far-field image without noise, and σn is the rms value of noise.

Because the mean value of noise is nonzero, a threshold has to be subtracted from the noisy image during data processing. In this paper, the threshold is denoted as

Threshold=3×σn.
(17)

After subtraction, if the intensity value of a pixel is negative, it is set to zero.

For an average wavefront aberration σ̄=0.4 rad whose initial average Strehl ratio is 0.84, different levels of noise are added to the far-field image, and 100 frame simulations are performed. The system aberration is defocus and its rms is σ=4 rad. The retrieval results are presented in Table 1.

Table 1. Results of phase retrieval with different SNRs

table-icon
View This Table

4. Conclusions

In the future we will provide evidence of the performance of this method in detailed experiments. Because of its limited measuring range and linear calculation, this linear phase retrieval method will be more suitable for closed-loop adaptive optics systems.

Acknowledgments

This project is funded by the National Natural Science Foundation of China with contract 60408005. The authors thank W. Jiang, C. Rao, H. Xian, and S. Feng for their helpful comments and suggestions.

References and links

1.

J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Eng. 21, 2758–2769 (1982).

2.

J. M. Wood, M. A. Fiddy, and R. E. Burge, “Phase retrieval using two intensity measurements in the complex plane,” Opt. Lett. 6, 514–516 (1981). [PubMed]

3.

J. R. Fienup, “Phase retrieval using a support constraint,” presented at the Institute of Electrical and Electronics Engineers ASSP Workshop on Multidimensional Digital Signal Processing, Leesburg, Va., October 28–30 (1985).

4.

R. A. Gonsalves, “Phase retrieval from modulus data,” J. Opt. Soc. Am. 66, 961–964 (1976).

5.

B. Ellerbroek and D. Morrison, “Linear methods in phase retrieval,” Proc. SPIE 351, 90–95 (1983).

6.

R. A. Gonsalves, “Small-phase solution to the phase-retrieval problem,” Opt. Lett. 26, 684–685 (2001).

7.

M. Li, X.-Y. Li, and W.-H. Jiang. “Small-phase retrieval with a single far-field image,” Opt. Express 16, 8190–8197 (2008). [PubMed]

8.

L. Xinyang and J. Wenhan, “Zernike modal wavefront reconstruction error of Hartmann wavefront sensor,” Acta. Optica Sinica 22, 1236–1240 (2002).

9.

X.-Y. Li, M. Li, B. Chen, and W.-H. Jiang “A kind of novel linear phase retrieval wavefront sensor and its application in close-loop adaptive optics system,” in Proc. Sixth International Workshop on Adaptive Optics for Industry and Medicine, C. Dainty, ed. (Imperial College Press, 2007), pp. 212–218.

10.

N. Roddier, “Atmospheric wavefront simulation using Zernike polynomials,” Opt. Eng. 29, 1174–1180 (1990).

OCIS Codes
(010.1080) Atmospheric and oceanic optics : Active or adaptive optics
(100.5070) Image processing : Phase retrieval

ToC Category:
Image Processing

History
Original Manuscript: April 16, 2009
Revised Manuscript: July 20, 2009
Manuscript Accepted: July 23, 2009
Published: August 14, 2009

Citation
Min Li and Xin-Yang Li, "Linear phase retrieval with a single far-field image based on Zernike polynomials," Opt. Express 17, 15257-15263 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-17-15257


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. R. Fienup, "Phase retrieval algorithms: a comparison," Appl. Eng. 21, 2758-2769 (1982).
  2. J. M. Wood, M. A. Fiddy, and R. E. Burge, "Phase retrieval using two intensity measurements in the complex plane," Opt. Lett. 6, 514-516 (1981). [PubMed]
  3. J. R. Fienup, "Phase retrieval using a support constraint," presented at the Institute of Electrical and Electronics Engineers ASSP Workshop on Multidimensional Digital Signal Processing, Leesburg, Va., October 28-30 (1985).
  4. R. A. Gonsalves, "Phase retrieval from modulus data," J. Opt. Soc. Am. 66, 961-964 (1976).
  5. B. Ellerbroek and D. Morrison, "Linear methods in phase retrieval," Proc. SPIE 351, 90-95 (1983).
  6. R. A. Gonsalves, "Small-phase solution to the phase-retrieval problem," Opt. Lett. 26, 684-685 (2001).
  7. M. Li, X.-Y. Li, and W.-H. Jiang. "Small-phase retrieval with a single far-field image," Opt. Express 16, 8190-8197 (2008). [PubMed]
  8. L. Xinyang and J. Wenhan, "Zernike modal wavefront reconstruction error of Hartmann wavefront sensor," Acta.Optica Sinica 22,1236-1240 (2002).
  9. X.-Y. Li, M. Li, B. Chen, and W.-H. Jiang, "A kind of novel linear phase retrieval wavefront sensor and its application in close-loop adaptive optics system," in Proc. Sixth International Workshop on Adaptive Optics for Industry and Medicine, C. Dainty, ed. (Imperial College Press, 2007), pp. 212-218.
  10. N. Roddier, "Atmospheric wavefront simulation using Zernike polynomials," Opt. Eng. 29, 1174-1180 (1990).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited