OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 5 — Mar. 2, 2009
  • pp: 3052–3062
« Show journal navigation

Restoration of turbulence-degraded extended object using the stochastic parallel gradient descent algorithm: numerical simulation

Huizhen Yang, Xinyang Li, Chenglong Gong, and Wenhan Jiang  »View Author Affiliations


Optics Express, Vol. 17, Issue 5, pp. 3052-3062 (2009)
http://dx.doi.org/10.1364/OE.17.003052


View Full Text Article

Acrobat PDF (542 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

An adaptive optics (AO) system with Stochastic Parallel Gradient Descent (SPGD) algorithm and a 61-element deformable mirror is simulated to restore the image of a turbulence-degraded extended object. SPGD is used to search the optimum voltages for the actuators of the deformable mirror. We try to find a convenient image performance metric, which is needed by SPGD, merely from a gray level distorted image and without any additional optics elements. Simulation results show the gray level variance function acts more promising than other metrics, such as metrics based on the gray level gradient of each pixel. The restoration capability of the AO system is investigated with different images and different turbulence strength wave-front aberrations using SPGD with the above resultant image quality criterion. Numerical simulation results verify the performance metric is effective and the AO system can restore those images degraded by different turbulence strengths successfully.

© 2009 Optical Society of America

1. Introduction

2. Model of space object imaging

The symbol (x,y) denotes coordinates of a point in the image plane and f(x, y) is gray level of the point (x, y). For incoherent imaging, the image f(x, y) is given by the convolution of the object function I(x, y) and the intensity point spread function (PSF), h(x, y) , of the system[7

7. J. W. Goodman. Introduction to Fourier Optics (Publishing House of Electronics Industry, 2006).

]:

f(x,y)=I(x,y)*h(x,y)
(1)

For clarity, magnification factors are omitted. All aberration effects are manifested in the PSF which can be obtained through Fast Fourier Transform of the wave-front in simulation. The original images used are given in Figs. 1(a) and 1(b), and are called Image A and Image B respectively.

Image A has a clear background and the object has clear edges, while the object in Image B almost fills the image plane. The image size is 128 × 128 pixels and the gray level is from 0 to 255. The diffraction limit angle (λ/D) is supposed to be 5 pixels, so the full Field of View (FOV) is about 24 times the size of diffraction limited angle. The ideal imaging results of Image A and Image B are given in Figs. 1(c) and 1 (d) under the telescope parameters in this paper.

We use the method proposed by N. Roddier, which makes use of a Zernike expansion of randomly weighted Karhunen-Loeve functions, to simulate atmospherically distorted wavefronts [8

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

]. Considering that the low-order aberrations (tilts, defocus, astigmatism, etc) have the most significant impact on image quality, we use the first 104 Zernike polynomial orders. Different phase screens generated according to this method are not correlated to each other and represent the Kolmogorov spectrum. The phase screens ϕ(r) are defined over 128 × 128 pixels which is also the grid of the wave-front corrector and don’t include the tip/tilt aberrations. The tip/tilt aberrations are usually controlled by another control loop and are considered as being removed completely in our simulation. Atmospheric turbulence strength for a receiver system with aperture size D can be characterized by the following two parameters: the ratio D/r 0 and the averaged Strehl Ratio (SR) of phase fluctuations. Phase screens of different atmospheric turbulence strength can be obtained through changing r 0 in the simulation program.

The correction capability of the AO imaging system is analyzed through three different turbulence strength wave-front aberrations where D/r 0 is set at 5, 10 and 20 respectively.

Fig. 1. Original image A (a), with clear edges and Background; Original image B (b), more complicated and almost fills the image plane; ideal imaging result (c) of Image A and that (d) of Image B when the FOV is about 24 times the size of diffraction limit angle.

3. Definition of several image sharpness functions

Four candidate image-quality metrics are examined in this paper. The first metric F 1 is defined according to the gray level variance of image, which expresses the discrete degree of gray levels distribution. Gray levels are distributed in a bigger range when the image is much clearer. Here

F1=xy[f(x,y)A]2
(2)

and A=(1/MN)xyf(x,y) is the mean of all pixel gray level, and M × N is the image size.

A high resolution image has clear edges and the gray level gradient expresses the edge information of an image, i.e. a bigger gray level gradient can offer us a better image. The second metric F 2 and the third F 3 are defined according to the gray level gradient information of each pixel.

F2=xy{[f(x,y)f(x+1,y)]2+[f(x,y)f(x,y+1)]2}
(3)

We define F 3 by using the Laplacian operator, which implements a second derivative operation on the image.

F3=xy{4f(x,y)f(x,y+1)f(x+1,y)f(x1,y)f(x,y1)}2
(4)

The square operation in F 2 and F 3 gives prominence to the contribution of some points, which have a large gray level gradient, to image quality evaluation function. F 2 and F 3 are also regarded as modification of sharpness function S 1 in Ref. 2

2. R. A. Muller and A. Buffington. “Real-time correction of atmospherically degraded telescope images through image sharpening,” J. Opt. Am . A. 64, 1200–1210 (1974). [CrossRef]

.

F4=uv(u2+v2)P(u,v)
(5)

To give emphasis to high frequency content, we weight the energy spectrum, P(u,v), with the weighting coefficient (u 2 + v 2), so F 4 can also be called as the mean square radius of the image spectrum.

4. Description of AO system

The high resolution imaging model with an AO system is shown in Fig. 2. The model includes an imaging system that records the image plane gray level f(x,y), an image quality analyzer that calculates the image quality metric and the SPGD algorithm that produces control signals u = {u 1,u 2,…u 61} for a 61-element Deformable Mirror (DM) according to changes of metric. The phase compensation m(x,y) , introduced by the DM, can be combined linearly with influence functions of actuators:

m(x,y)=j=161ujSj(x,y)
(6)

where Sj(x,y) is the influence function and uj is the control signal of the jth actuator. On the basis of experimental measurements, we know the actuator influence function of a 61-element DM actuators is approximately Gaussian [9

9. W. H. Jiang, N. Ling, X. J. Rao, and F. shi. “Fitting capability of deformable mirror,” SPIE 1542130–137 (1991). [CrossRef]

]:

Sj(r)=Sj(x,y)=exp{Inω[(xxj)2+(yyj)2/d]α}
(7)

where (xj,yj) is the location of the jth actuator, α is the coupling value between actuators and is set to 0.08, and a is the Gaussian index and is set to 2. The distance between actuators is d which is about 0.1367 normalized in a unit circle. Fig. 3 gives actuators location distribution of DM. The circled line in the figure denotes the effective aperture and the layout of all actuators is hexagonal. We suppose the stroke of the DM is enough to correct the wave-front aberrations in the simulation.

The SPGD algorithm [5

5. M. A. Vorontsov and G. W. Carhart. “Adaptive optics based on analog parallel stochastic optimization: analysis and experimental demonstration,” J. Opt. Soc. Am . A. 17, 1440–1453 (2000) [CrossRef]

] requires small random perturbations ∆u = {∆u 1,∆u 2,…∆u 61} with fixed amplitude |∆uj|=σ and random signs with equal probabilities for Pr(∆uj = ±σ) = 0.5 [10

10. J. C. Spall, “Multivariate stochastic approximation using a simultaneous perturbation gradient approximation,” IEEE Trans. Autom. Control . 37, 332–341 (1992). [CrossRef]

], to be applied to all 61 DM control channels simultaneously. Then for a given random ∆u , the control signals are updated with the rule:

u(k+1)=u(k)+γΔu(k)ΔF(k),k=0,1,
(8)

Where the scalar ∆F = F(u + ∆u) - F(u - ∆u) is the corresponding perturbation of the image quality metric and γ is a gain coefficient which scales the size of the control parameter corrections (positive for the case of metric maximization, negative for minimization and positive in this paper).

Fig. 2. Schematic diagram of a high resolution imaging model.
Fig. 3. Actuator distribution of 61-element DM.

5. Numerical results and analysis

5.1 Adaptation process

We perform the adaptation process over 50 phase realizations. The averaged evolution curves and the standard deviation evolution curves of the metric are the recorded simulation results. For comparing correction ability of the four image-quality metrics, we also give corresponding averaged SR evolution curves.

First, correction performance of Image A is investigated. Fig. 4 to Fig. 7 show simulation results when D/r 0 is 5, 10 and 20 respectively. Averaged evolution curves of four image-quality metrics are given in Fig. 4(a) to Fig. 7(a), in which the averaged evolution curves are normalized to be 1 in the optimal case. Corresponding standard deviation curves and Strehl Ratio curves during the control algorithm’s 1500 iterations are presented in Fig. 4(b) to Fig. 7(b) and Fig. 4(c) to Fig. 7(c).

All image-quality metric curves have converged after 1500 iterations in Fig. 4(a) to Fig. 7(a). The iteration number needed for convergence increases and correction performance decreases with the growth of turbulence strength. From Fig. 4(b) to Fig. 7(b), we can see that F 1 and F 4 have relatively smaller standard deviations than F 2 and F 3 , which shows that F 1 and F 4 have stronger adaptability to different turbulence realizations than F 2 and F 3. Fig. 4(c) to Fig. 7(c) indicate that all the four image-quality metrics have strong correction ability for D/r 0 = 5 , F 1 and F 4 are much better than F 2 and F 3 for D/r 0 = 10 , and F 1 is the best for D/r 0 = 20 .

The same adaptation process is carried out for Image B. Results are very similar to those of Image A and aren’t given in the section.

Fig. 4. Averaged curve of F 1 (a), the corresponding standard deviation curves (b) and SR curves(c) during 1500 iterations.
Fig. 5. Averaged curve of F 2 (a), the corresponding standard deviation curves (b) and SR curves(c) during 1500 iterations.
Fig. 6. Averaged curve of F 3 (a), the corresponding standard deviation curves (b) and SR curves(c) during 1500 iterations.
Fig. 7. Averaged curve of F 4 (a), the corresponding standard deviation curves (b) and SR curves(c) during 1500 iterations.

5.2 Zernike order and PSF of single frame phase screen

Figure 8 gives Zernike coefficients 3-104 decomposed from the same phase screen (D/r 0 = 20) when F 1, F 2, F 3 and F 4 are used as performance metrics optimized respectively. Corresponding PSFs are shown in Fig. 9. We also fit the DM figure to the phase screen using least squares to obtain the best correction achievable with the given 61-element DM. The fitting results are also shown in Fig. 8 and Fig. 9.

Fig. 8. Comparison of Zernike coefficients 3-104 before correction (a) and after correction with F 1 (c), F 2 (d), F 3 (e) and F 4 (f) ; (b) is Zernike coefficients of the residual wave-front with the least squares fitting.
Fig. 9. Comparison of PSF before correction (a) and after correction with F 1 (c), F 2 (d), F 3 (e) and F 4 (f); (b) is PSF of the residual wave-front with the least squares fitting.

From Fig. 8, we can see relatively low order content occupies the main part of the wave-front aberrations. Fig. 8 and Fig. 9 all show F 1 and F 4 have much higher correction ability than F 2 and F 3. Compared with the least squares fitting, F 1 almost obtains the best correction achievable for the 61-element DM.

5.3 Cormparison of imaging results

Figure 10 presents averaged imaging results of Image A after different turbulence strength wave- front aberrations are corrected by using F 1 , F 2 , F 3 and F 4 respectively.

Fig. 10. Imaging results comparison of Image A before correction (A) and after correction (B), (C), (D) and (E) with F 1, F 2 F 3 ,and F 4 respectively.

Figure 11 shows imaging results of Image B after different turbulence strength wave-front aberrations are corrected by using F 1 , F 2 , F 3 and F 4 respectively.

Fig. 11. Imaging results comparison of Image B before correction (A) and after correction (B), (C), (D) and (E) with F 1, F 2 F3 ,and F 4 respectively.

5.4 Correction ability analysis of performance metrics

For comparing correction performance of four image-quality metrics, we examine the SR before and after correction by these metrics on Image A. Results are shown in Fig. 12. The result of the least squares fitting is also given in Fig. 12. From Fig. 12, we see that four image-quality metrics have strong correction ability when the turbulence strength is relatively weak. F 1 is similar in correction ability to F 4 , F 2 is similar to F 3 , and F 1 and F 4 are better than F 2 and F 3 when D/r 0 is 10. F 1 is the best when D/r 0 is 20. The correction capability of F 1 is very close to the least squares fitting.

Fig. 12. Comparison of four image-quality metrics and the least squares fitting on correction ability

Fig. 13. Different image metrics versus D/r 0.

6. Conclusion

We have simulated a high resolution imaging AO system with a 61-element deformable mirror and SPGD control algorithms. Numerical simulation results verify the performance metric offered is effective, which can provide us with guidance on how to choose appropriate performance metrics in extended object imaging: image metrics based on image global information are better than those based on image local information and the grey level variance function is the most promising image-quality metric. The compensation effect is very close to the optimal correction on combining the grey level variance function and SPGD control algorithm for a given DM. The AO system can restore those images degraded by different turbulence strength successfully. Note that the iteration number needed for convergence increases and correction performance decreases with the growth of turbulence strength. The system can be applied to static or slowly changing wave-front aberrations, while the convergence rate must be improved for rapidly changing aberrations. Further research should be focused on improving the convergence rate so that the system can be used for the dynamic environment.

Acknowledgments

We express our thanks for the support of the National High Technology Project of China.

References and links

1.

R. K. Tyson, Principle of Adaptive Optics (Academic Press, 1991).

2.

R. A. Muller and A. Buffington. “Real-time correction of atmospherically degraded telescope images through image sharpening,” J. Opt. Am . A. 64, 1200–1210 (1974). [CrossRef]

3.

S. Zommer, E. N. Ribak, S. G. Lipson, and J. Adler. “Simulated annealing in ocular adaptive optics,” Opt. Lett . 31, 1–3 (2000).

4.

P. Yang, M. W. Ao, Y. Li, B. Xu, and W. H. Jiang, “Intracavity transverse modes controlled by a genetic algorithm based on Zernike mode coefficients,” Opt. Express 15, 17051–17062 (2007). [CrossRef] [PubMed]

5.

M. A. Vorontsov and G. W. Carhart. “Adaptive optics based on analog parallel stochastic optimization: analysis and experimental demonstration,” J. Opt. Soc. Am . A. 17, 1440–1453 (2000) [CrossRef]

6.

M. S. Zakynthinaki and Y. G. Saridakis, “Stochastic optimization for a tip-tilt adaptive correcting system,” Comput. Phys. Commun . 150, 274–292 (2003) [CrossRef]

7.

J. W. Goodman. Introduction to Fourier Optics (Publishing House of Electronics Industry, 2006).

8.

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

9.

W. H. Jiang, N. Ling, X. J. Rao, and F. shi. “Fitting capability of deformable mirror,” SPIE 1542130–137 (1991). [CrossRef]

10.

J. C. Spall, “Multivariate stochastic approximation using a simultaneous perturbation gradient approximation,” IEEE Trans. Autom. Control . 37, 332–341 (1992). [CrossRef]

OCIS Codes
(010.0010) Atmospheric and oceanic optics : Atmospheric and oceanic optics
(010.1080) Atmospheric and oceanic optics : Active or adaptive optics
(110.0110) Imaging systems : Imaging systems
(350.4600) Other areas of optics : Optical engineering

ToC Category:
Atmospheric and Oceanic Optics

History
Original Manuscript: November 3, 2008
Revised Manuscript: January 4, 2009
Manuscript Accepted: January 21, 2009
Published: February 17, 2009

Citation
Huizhen Yang, Xinyang Li, Chenglong Gong, and Wenhan Jiang, "Restoration of turbulence-degraded extended object using the stochastic parallel gradient descent algorithm: numerical simulation," Opt. Express 17, 3052-3062 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-5-3052


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. R. K. Tyson, Principle of Adaptive Optics (Academic Press, 1991).
  2. R. A. Muller and A. Buffington, "Real-time correction of atmospherically degraded telescope images through image sharpening," J. Opt. Am. A. 64, 1200 -1210 (1974). [CrossRef]
  3. S. Zommer, E. N. Ribak, S. G. Lipson and J. Adler, "Simulated annealing in ocular adaptive optics," Opt. Lett. 31, 1-3 (2000).
  4. P. Yang, M. W. Ao, Y. Li, B. Xu, and W. H. Jiang, "Intracavity transverse modes controlled by a genetic algorithm based on Zernike mode coefficients," Opt. Express 15, 17051-17062 (2007). [CrossRef] [PubMed]
  5. M. A. Vorontsov and G. W. Carhart, "Adaptive optics based on analog parallel stochastic optimization: analysis and experimental demonstration," J. Opt. Soc. Am. A. 17, 1440 -1453 (2000) [CrossRef]
  6. M. S. Zakynthinaki and Y. G. Saridakis, "Stochastic optimization for a tip-tilt adaptive correcting system," Comput. Phys. Commun. 150, 274 -292 (2003) [CrossRef]
  7. J. W. Goodman. Introduction to Fourier Optics (Publishing House of Electronics Industry, 2006).
  8. N. Roddier, "Atmospheric wavefront simulation using Zernike polynomials," Opt. Eng. 29, 1174 -1180 (1990). [CrossRef]
  9. W. H. Jiang, N. Ling, X. J. Rao, and F. Shi. "Fitting capability of deformable mirror," Proc. SPIE 1542, 130 -137 (1991). [CrossRef]
  10. J. C. Spall, "Multivariate stochastic approximation using a simultaneous perturbation gradient approximation," IEEE Trans. Autom. Control. 37, 332 -341 (1992). [CrossRef]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited