OSA's Digital Library

Optics Express

Optics Express

  • Editor: Andrew M. Weiner
  • Vol. 22, Iss. 11 — Jun. 2, 2014
  • pp: 13792–13803
« Show journal navigation

Multichannel-Hadamard calibration of high-order adaptive optics systems

Youming Guo, Changhui Rao, Hua Bao, Ang Zhang, Xuejun Zhang, and Kai Wei  »View Author Affiliations


Optics Express, Vol. 22, Issue 11, pp. 13792-13803 (2014)
http://dx.doi.org/10.1364/OE.22.013792


View Full Text Article

Acrobat PDF (919 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

we present a novel technique of calibrating the interaction matrix for high-order adaptive optics systems, called the multichannel-Hadamard method. In this method, the deformable mirror actuators are firstly divided into a series of channels according to their coupling relationship, and then the voltage-oriented Hadamard method is applied to these channels. Taking the 595-element adaptive optics system as an example, the procedure is described in detail. The optimal channel dividing is discussed and tested by numerical simulation. The proposed method is also compared with the voltage-oriented Hadamard only method and the multichannel only method by experiments. Results show that the multichannel-Hadamard method can produce significant improvement on interaction matrix measurement.

© 2014 Optical Society of America

1. Introduction

Adaptive optics (AO) systems can compensate for the turbulence-induced phase distortion in real time so as to improve the image quality of astronomical telescopes [1

1. J. W. Hardy, Adaptive Optics for Astronomical Telescopes (Oxford University, 1998), Chap. 1.

,2

2. H. W. Babcock, “Adaptive optics revisited,” Science 249(4966), 253–257 (1990). [CrossRef] [PubMed]

]. In order to get the proper commands to control the deformable mirror (DM), the gain matrix i.e., interaction matrix (IM) that maps DM commands into wavefront sensor (WFS) measurements has to be calibrated. The traditional zonal method is to drive the different actuators successively and measure the corresponding slope response vectors which compose the columns of the IM [3

3. C. Boyer, V. Michau, and G. Rousset, “Adaptive optics: Interaction matrix measurements and real time control algorithms for the COME-ON project,” Proc. SPIE 1237, 406–421 (1990). [CrossRef]

,4

4. W. Jiang and H. Li, “Hartmann-Shack wavefront sensing and wavefront control algorithm,” Proc. SPIE 1271, 82–93 (1990). [CrossRef]

]. However, with the development of high-order AO systems [5

5. A. Bouchez, R. Dekany, J. Angione, C. Baranec, M. Britton, K. Bui, R. Burruss, J. Cromer, S. Guiwits, J. Henning, J. Hickey, D. Mckenna, A. Moore, J. Roberts, T. Trinh, M. Troy, T. Truong, and V. Velur, “The PALM-3000 high-order adaptive optics system for Palomar Observatory,” Proc. SPIE 7015, 70150Z (2008). [CrossRef]

7

7. S. J. Thomas, L. Poyneer, R. De Rosa, B. Macintosh, D. Dillon, J. K. Wallace, D. Palmer, D. Gavel, B. Bauman, L. Saddlemyer, and S. Goodsell, “Integration and test of the Gemini Planet Imager,” Proc. SPIE 8149, 814903 (2011). [CrossRef]

], the traditional zonal method becomes difficult to implement since the procedure is unacceptably time consuming and sensitive to noise [8

8. S. Oberti, F. Quirós-Pacheco, S. Esposito, R. Muradore, R. Arsenault, E. Fedrigo, M. Kasper, J. Kolb, E. Marchetti, A. Riccardi, C. Soenke, and S. Stroebele, “Large DM AO systems: synthetic IM or calibration on sky?” Proc. SPIE 6272, 627220 (2006). [CrossRef]

]. Considering the saturation of DM actuators, a voltage-oriented Hadamard method was proposed by Kasper to reduce the impact of turbulence and photon noise [9

9. M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, and S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]

]. Nevertheless, in closed-loop systems based on the curvature WFS or the pyramid WFS, the linear range of the WFS is much smaller than the DM so the maximum voltage can hardly be achieved [10

10. S. Oberti, H. Bonnet, E. Fedrigo, L. Ivanescu, M. Kasper, and J. Paufique, “Calibration of a curvature sensor/bimorph mirror AO system: interaction matrix measurement on MACAO systems,” Proc. SPIE 5490, 139–150 (2004). [CrossRef]

,11

11. A. Riccardi, R. Briguglio, E. Pinna, G. Agapito, F. Quiros-Pacheco, and S. Esposito, “Calibration strategy of the pyramid wavefront sensor module of ERIS with the VLT deformable secondary mirror,” Proc. SPIE 8447, 84475M (2012). [CrossRef]

]. Thereby, Meimon suggested a slope-oriented Hadamard method to make the WFS response matrix satisfy a Hadamard matrix pattern. However, an initial estimate of the IM is necessary to calculate the optimal DM commands [12

12. S. Meimon, T. Fusco and C. Petit, “An optimized calibration strategy for high order adaptive optics systems: the Slope-Oriented Hadmard Actuation,” in 1st AO4ELT Conference (2010), paper 07009.

]. Zou has developed an iterative technique for improving the IM calibration to overcome the non-linearity of the DM, in which the static wavefront error is firstly corrected with a generic IM [13

13. W. Zou and S. A. Burns, “High-accuracy wavefront control for retinal imaging with Adaptive-Influence-Matrix Adaptive Optics,” Opt. Express 17(22), 20167–20177 (2009). [CrossRef] [PubMed]

]. Both the slope-oriented method and the iterative method need an initial IM, which influences both the performance and the stability of these methods. Therefore, a good calibration of the initial IM is important for the convergence of these iterative methods. Other calibration methods like those using modulation techniques have also been introduced to reduce the turbulence induced error when calibrating the IM of an adaptive secondary AO system [14

14. F. Wildi and G. Brusa, “Determining the interaction matrix using starlight,” Proc. SPIE 5490, 164–173 (2004). [CrossRef]

,15

15. S. Esposito, R. Tubbs, A. Puglisi, S. Oberti, A. Tozzi, M. Xompero, and D. Zanotti, “High SNR measurement of interaction matrix on-sky and in lab,” Proc. SPIE 6272, 62721C (2006). [CrossRef]

]. The sinusoidal calibration technique has already been tested on sky at the LBT [16

16. E. Pinna, F. Quiros-Pacheco, A. Riccardi, R. Briguglio, A. Puglisi, L. Busoni, C. Arcidiacono, J. Argomedo, M. Xompero, E. Marchetti, and S. Esposito, “First on-sky calibration of a high order adaptive optics system,” Proc. SPIE 8447, 84472B (2012). [CrossRef]

]. These modulation methods are particularly fit for adaptive secondary AO systems because the power spectral density (PSD) of the atmospheric turbulence decreases rapidly with frequency.

In this paper we describe a novel IM calibration method called the multichannel-Hadamard (MH) method for traditional high-order AO systems with the advantages of both the voltage-oriented Hadamard only (HO) method and a kind of multichannel only (MO) method. The IM acquired by this method can be used for other iterative or closed-loop calibration methods as the initial IM. The IM calibration metric is defined in section 2. The principle of the multichannel-Hadamard method is introduced in detail in section 3. Section 4 shows the experimental results.

2. Calibration error

In order to quantify the performance of different IM calibration methods, the definition of the calibration error is necessary. The calibration error used here is a minor modification of the one used by S. Oberti et. al. [8

8. S. Oberti, F. Quirós-Pacheco, S. Esposito, R. Muradore, R. Arsenault, E. Fedrigo, M. Kasper, J. Kolb, E. Marchetti, A. Riccardi, C. Soenke, and S. Stroebele, “Large DM AO systems: synthetic IM or calibration on sky?” Proc. SPIE 6272, 627220 (2006). [CrossRef]

]. To give the calibration error more meaning, the normalization is used here. Thus, the calibration error is defined as the normalization of the mean-square error between the exact IM and the measured IM. If the exact IM is D0 and the estimated IM is Dm, the calibration error here is defined as

J=trace[(D0Dm)T(D0Dm)]/trace(D0TD0).
(1)

Obviously, the smaller the IM error is, the smaller the calibration error will be. Owing to the normalization, if the estimated IM is a null matrix, J equals 1, while if the estimated IM equals the exact IM, i.e., the perfect calibration performance achieved, J equals 0. Although this metric is different to the one used in [9

9. M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, and S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]

] which is the geometric mean of the standard deviation of the IM elements, both of them reflect the same physical meaning that is the deviation of the measured IM from the exact IM.

3. Principle of the multichannel-Hadamard calibration method

3.1. Basic principle

As each actuator of the DM can only influence local areas of the WFS, slopes measured in subapertures farther from the pushed actuator will have lower signal-to-noise ratios (SNR) so it is better to replace some of the slope responses with extremely low SNRs by zeros. Therefore, the IM of a large AO system is a sparse matrix and actuators far enough from each other can be driven at the same time with little coupling. These actuators are treated as a channel and the whole actuators of a DM can be divided into a series of channels. Then, the voltage-oriented Hadamard actuation is applied to these channels to reject the calibration noise. Finally, the IM mapping the multichannel pattern commands into the WFS measurements is translated into the one that maps the DM commands into the WFS measurements. This calibration method is called as the multichannel-Hadamard method. Taking the 595-element AO system as an example (The layout of the DM actuators and the WFS subapertures is shown in Fig. 1
Fig. 1 Multichannel mode of the 595-actuator DM. The circles denote the DM actuators and the numbers in the circles represent which channel the actuators belong to. The squares mean the subapertures.
), the coupling of the DM is about 9.5%. The detailed description of this method is below. The whole procedure can be summarized as three steps:

1. The DM actuators are firstly divided into 19 channels by the rule that the distance between any of the two actuators in the same channel is not closer than 3.4dact, where dact is the spacing of two adjacent actuators in the same line. This rule makes actuators in the same channel almost uncoupled with each other. The detailed description of the channel dividing result is also shown in Fig. 1. The channel dividing result can also be expressed by a channel dividing matrix (Fig. 2
Fig. 2 Channel dividing matrix (C) of the 595-actuator DM. The x axis represents the channel number and the y axis denotes the actuator number. Elements that equal 1 are black and those equal 0 are white.
) which is defined as
C(i,j)={=1whiletheithactuatorbelongstothejthchannel=0otherwise.
(2)
Here r equaling 1.7dact is defined as the effective influence radius (EIR), and the circle with center located at the actuator’s position and radius equaling r is defined to be the effective influence area (EIA) of the actuator. The EIA can be described by an effective influence matrix (EIM) defined as
Seffect(2p1,q)2nsub×nact=Seffect(2p,q)2nsub×nact={1;whilethecenterofthepthsubapertureisintheEIAoftheqthactuator0;otherwise,
(3)
where nsub is the subaperture number and nact is the actuator number. In Eq. (3), the row number of Seffect is 2nsub because both x-slope and y-slope responses are contained in the IM. The EIM of the 595-element AO system is shown in Fig. 3
Fig. 3 The effective influence matrix Seffect. The x axis represents the actuator number and the y axis denotes the slope number. Elements that equal 1 are black and those equal 0 are white.
, which indicates that this matrix is a sparse matrix. It is certain that the 9.5% coupling relationship in our system is a bit low, but one can easily adjust the EIR with different actuator coupling as required. If the coupling factor of the AO system increases, the EIR should increase too.

2. The 19 channels are treated as 19 “big” actuators and then a new DM with 19 “big” actuators is another form of the 595-actuator DM. Afterwards, if the MO method is used, it will do a push-pull on the channels and gets the average responses one by one. In that case, an identity matrix pattern is used, and the voltage matrix would be
VMO=vMOInMO,
(4)
where vMO is the amplitude of the driving voltage and nMO is the channel number of the MO method. Different from that, our MH method employs a Hadamard matrix pattern for the channels to benefit from the noise rejection effect of the HO method [9

9. M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, and S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]

]. Thus, the voltage matrix for the MH method can be described as
VMH=vMHHnMH,
(5)
where vMH is the amplitude of the driving voltage and nMH is channel number of the MH method. HnMHis a Hadamard matrix of dimension nMH. Depending on the value of nMO, several channels without any actuators may be necessary to make the channel number nMH satisfy the dimension requirement of a Hadamard matrix. In that case,
CMH=[C0],
(6)
where the dimension of the block matrix 0 is nact×(nMHnMO). When driving the channels, the WFS measurement matrix
G=DcV+N,
(7)
where V is the voltage matrix, N is the measurement noise matrix and Dc is the relationship matrix that maps the channels’ commands into the WFS measurements which can be estimated by

Dcm=GV1.
(8)

3. In order to translate Dcm into the IM mapping DM actuators’ commands into the WFS measurements, two steps have to be done. Firstly, the slope response of each channel is copied to the columns of the corresponding slope response of the actuators that belong to that channel by
Dtemp=DcmCMHT.
(9)
Secondly, as the slope response of a channel is composed of the slope response of each actuator in that channel, elements in the slope response that belong to other actuators in the same channel should be set to zeros. This process depends on the EIA of the actuators. As the centre of a subaperture can only be in the EIA of one actuator in the same channel, one can easily determine which subaperture slope response is caused by which actuator in that channel. This procedure can be expressed in mathematics as
DMH=Dtemp×Seffect.
(10)
Thus, the matrix DMH is the IM that we need which maps the DM actuators’ commands into the WFS measurements.

3.2. Error analysis

To compare the performance of the MO, HO and MH methods, the calibration errors are analyzed in theory. For the MH method, if the EIR is large enough to make the elements that are going to be set to zeros small enough, the IM error
ΔDMH=DMHD01vMHnMHNMHHnMHTCMHT×Seffect,
(11)
where NMH is the noise matrix of the MH method whose dimension is 2nsub×nMH. Similarly, the IM error of the MO method
ΔDMO=DMOD01vMONMOCT×Seffect,
(12)
where NMO is the noise matrix of the MO method whose dimension is 2nsub×nMO. According to [9

9. M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, and S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]

], the IM error of the HO method here is
ΔDHO=DHOD0=1vHOnHONHOHnHOT,
(13)
where NHO is the noise matrix of the HO method whose dimension is 2nsub×nHO, vHO is the driving voltage amplitude and HnHO is a Hadamard matrix of dimension nHO. As the slope measurement noise of each subaperture can be considered as a Gaussian white noise, the noise matrix meets
NMHTNMH=2nsubσn2InMH,
(14)
where InMH is the identity matrix of dimension nMH and σn2 is the average variance . If one push-pull is applied when driving the actuators, this variance would be
σn2=[noise(j,k)noise(j,k1)2]2j,k
(15)
where noise(j,k) is the jth element of the slope noise vector at the kth time of driving the DM, and <>j,k means the ensemble average. Besides, the Hadamard matrix meets
HnMHHnMHT=nMHInMH,
(16)
and the trace of the channel dividing matrix
trace(CMHCMHT)=nact.
(17)
Because of the EIM Seffect, lots of elements in the DMH with small SNRs are set to zeros. Thus, with Eqs. (1), (11), (14), (16) and (17), the calibration error of the MH method can be approximated by
JMH2(nDn0)nactnsubσn2nDnMHvMH2trace(D0TD0),
(18)
where nD is the total number of elements in Seffect and n0 is the number of the elements equaling 0. Similarly, it can be inferred that the calibration error of the MO method
JMO2(nDn0)nactnsubσn2nDvMO2trace(D0TD0)
(19)
and the calibration error of the HO method
JHO2nactnsubσn2nHOvHO2trace(D0TD0).
(20)
With Eqs. (18), (19) and (20)
JMH:JMO:JHOnDn0nDnMHvMH2:nDn0nDvMO2:1nHOvHO2.
(21)
For a zonal IM, the factor (nD-n0)/nD is from the EIA definition which makes the elements with low SNRs set to zeros. Obviously, this factor increases with the EIR which is partly determined by the DM coupling. If the EIR is set large enough, the factor will become 1 and nMH will equal nHO. In that case, the MH method is the same as the HO method. Therefore, the HO method is a special case of the MH method which especially fits the AO system with a large coupling DM. What’s more, the factors 1/ nMH and 1/ nHO are brought by the Hadamard matrix pattern of the driving voltage matrix because with a Hadamard matrix pattern, each actuator is driven up and down repeatedly to average the measurement noise. If the channel number nMH becomes larger, the calibration error will be smaller. However, the time used to finish the calibration will be longer like the principle of increasing the push-pull times to average the measurement noise [17

17. X. Zhang, C. Arcidiacono, A. R. Conrad, T. M. Herbst, W. Gaessler, T. Bertram, R. Ragazzoni, L. Schreiber, E. Diolaiti, M. Kuerster, P. Bizenberger, D. Meschke, H. W. Rix, C. Rao, L. Mohr, F. Briegel, F. Kittmann, J. Berwein, and J. Trowitzsch, “Calibrating the interaction matrix for the LINC-NIRVANA high layer wavefront sensor,” Opt. Express 20(7), 8078–8092 (2012). [PubMed]

].

3.3. Optimal EIR

When dividing the DM actuators, the EIR is one of the most important factors that influence the performance of the MH method. For different AO systems, the actuator couplings and the WFS measurement error may be different. During the IM calibration, the slope response’s SNR of a subaperture is in direct proportion to the actuator coupling while inversely proportional to the measurement noise. Therefore, the optimal EIR indeed exists.

The influence of the actuator coupling and the measurement noise on the EIR was analyzed by Monte Carlo simulations. Based on the influence function of the DM and the geometry relationship between the actuators and the subapertures, the ideal IM was firstly computed without adding any measurement noise, i.e., the slope response’s measurement SNR was infinite. Then, the IM errors were computed by simulating the IM calibration with different EIRs in different SNR conditions. During the simulation, the calibration SNR was increased from 5 to 15, and the actuator couplings were 9.5%, 20% and 30%, respectively. The calibration SNR here is defined as
snr=vmax(|D0|)σnnnMH
(22)
where v is the amplitude of the driving voltage, max(|D0|) is the maximum of the absolute value of the IM, n is the number of push-pulls added during the calibration, nMH is that of the MH method

During the simulation, in spite of the different channel numbers (all less than 64), an order 64 Hadamard matrix was used as the voltage matrix pattern to exclude the influence of different numbers of push-pulls required by different Hadamard matrices, so both n and nMH equaled 64 here. The results (Fig. 4
Fig. 4 The relationship between optimal EIR, coupling and the measurement SNR. The markers are the average of 10 times repeated results and the dashed lines are the fitted results.
) show that the optimal EIRs increase almost linearly with snr. Besides, the larger the coupling is, the larger the optimal EIR will be if snr is fixed. For the 595-element AO system we used (9.5% coupling), the optimal EIR is about 1.7~1.8dact when snr increases from 5 to 15. That’s why the EIR used in our experiment is 1.7dact.

4. Experiment

In order to test the performance of the MO, HO and MH calibration methods, experiments have been done on the 595-element AO system. The 595-element AO system (located in the Institute of Optics and Electronics, Chinese Academy of Sciences) mainly contains a 595-actuator piezoelectric DM (9.5% actuator coupling), a 600-subaperture Shack-Hartmann WFS (28x28 subapertures). The optical schematic of this AO system is shown in Fig. 5
Fig. 5 Optical schematic of the 595-element AO system.
. The source is a 600 nm semiconductor laser.

In different methods, the DM actuators were divided into different number of groups, so different number of push-pulls were needed to drive the DM actuators during the whole calibration period. The first task was to define the calibration condition to exclude the impact of the different number of push-pulls. The calibration SNR defined in Eq. (22) was taken as the calibration condition. Generally, The RMS of a measured signal is inversely proportional to the square root of the measurements times when the measurement noise is white. If the numbers of push-pulls of driving the DM by using the MO, HO and MH methods in our system were denoted by nMO, nHO and nMH, the driving voltages of these methods should satisfy
vMH:vMO:vHO=1:nMH/nMO:nMH/nHO
(23)
to keep these methods working in the same calibration condition. With Eqs. (21) and (23),
JMH:JMO:JHO=1:nMO:nDnDn0
(24)
In our experiment, nMO is 19, nD is 714000 and n0 is 699552 so
JMH:JMO:JHO=1:19:49.4
(25)
In addition, the real IM was also required to calculate the calibration error J defined in Eq. (1). Although it couldn’t be really measured, we could use an IM calibrated in a high SNR condition to replace. In our experiment, an IM calibrated with snr equal to 35 was used as the real IM. Sets of ten IMs were taken by using each of the three methods mentioned above with different SNRs. The corresponding RMs were then calculated. The parameter σn during the experiment was about 1.99x10−5 rad. The simulation was also done on the computer using the real IM and the slope noise collected on the 595-element AO system with the calibration condition the same as the experiment. The experimental and simulation results are shown in Fig. 6
Fig. 6 Calibration errors of the MO, HO and MH methods with calibration SNR from 5 to 15. These errors are calculated by Eq. (1). The dash lines are simulation results by using the three methods and the markers are the experimental results of the average calibration error of ten IMs.
. Because of the relatively small coupling value of the DM in our experiment, a great number of elements approaching zeros exist in the IM, which resulted in the calibration error of the MO method smaller than the HO method. However, the MO method isn’t always better than the HO method, especially if a larger coupling DM is used which makes less elements of the IM close to zeros. With the advantages of both the MO and HO methods, the MH method had the smallest error.

Since the IMs were used to calculate the RM and make the AO closed-loop system work, the most important way to judge the quality of them is through closed-loop performance. As the pseudo inverse of the IMs measured above, RMs were calculated and put to use in a kind of “one-shot” close, in which the compensation voltage vector calculated by the RM was directly applied to the DM. To simplify the closed-loop procedure, none iteration was used. It was a test of the wavefront reconstruction capability of different RMs. After “one-shot” close, the slope residual variance caused by the IM error can be given by [12

12. S. Meimon, T. Fusco and C. Petit, “An optimized calibration strategy for high order adaptive optics systems: the Slope-Oriented Hadmard Actuation,” in 1st AO4ELT Conference (2010), paper 07009.

]
σgIM2(D+)=σg2(D+)σg2(D0+),
(26)
where σg2(D+),σg2(D0+) were the slope residual variances obtained by using D+ and D0+. Since D0 was measured with a high SNR about 35 by using the MH method, it was considered as the real IM in our experiment. Therefore, σg2(D0+) was considered as the slope error variance which was only caused by the fitting error of the DM. σgIM2(D+) was the slope residual variance caused by the IM error. The initial slope error in the experiment was caused by the static aberration of the AO system which scarcely varied. The initial slope error was then compensated by the ‘one-shot’ close. The slope variances induced by the IM error after the ‘one-shot’ close are shown in Fig. 7
Fig. 7 Slope variances induced by the IM error after the “one-shot” close with calibration SNR from 5 to 15.
. The similar trend of Fig. 6 and Fig. 7 demonstrates that the calibration error defined in Eq. (1) can indeed represent the slope residual variance caused by the IM error.

Although the slope error variance could indicate the quality of a zonal IM, the phase error variance would be a more convincing index. The phase error variances induced by the IM error after the “one-shot” close are shown in Fig. 8
Fig. 8 Phase variances induced by the IM error after the “one-shot” close with calibration SNR from 5 to 15. These variances were calculated by subtracting the phase variance caused by the fitting error (this variance equaled the phase residual variance when snr equaled 35 in our experiments) from the phase residual variances. The unit λ denotes wavelength.
. Generally speaking, the phase error variances induced by the IM error decrease with the calibration SNR and the variances of the MH method are the smallest among the three methods. Even if the calibration SNRs of the other two methods reach 15, their variances are still larger than the variances of the MH method with SNR 5. These results agree with the slope variances in Fig. 7 very well.

In fact, almost every AO system works in real on-going closed-loop scene. Obviously, when correcting the same static aberration, a better IM should compensate the aberration faster and if the same closed-loop time is used, a larger Strehl ratio (SR) of the focal plane image should be achieved. In our experiment, the IMs and the corresponding RMs measured above were used to make the AO system closed-loop work and compensate the static aberration of the AO system. The sampling frequency was set at a small value about 10Hz for good compare of the closed-loop convergence process. Each close-loop procedure made 200 closed-loop iterations before the focal plane image was acquired. The controller was based on a PID control algorithm whose parameters were fixed all through. Some of the focal plane intensities and corresponding SRs are shown in Fig. 9
Fig. 9 Focal plane intensities of the imaging system by using MO, HO and MH methods for calibration.
. As it shows, the SRs acquired by using the MH method when the calibration SNRs equal 35 and 9 are approaching. The reason is that when snr is high, fitting error is the dominant factor that influences the focal plane intensity and the IM error becomes trivial. More importantly, although the SR acquired by using the MH method decreases with snr decreasing, the one acquired by the MH method with snr equivalent to 5 is still larger than those acquired by the MO and HO methods with snr equal to 10. These results can match the corresponding phase error variances shown in Fig. 8, which also indicts the superiority of the MH method. It means that if the aberration is static, the MH method can make a least phase residual variance in the three methods with the same closed-loop iterations and if the aberration is dynamic, the MH method will track the change of the turbulence best, too. The PSDs of the residual phases with snr 10 are shown in Fig. 10
Fig. 10 PSD of the closed-loop residual phase error by using the IM calibrated with snr 10.
. As it shows, the MH method performs very well not only in high spatial frequencies but also in low spatial frequencies. Considering the calibration efficiency requirement of high-order AO systems, the MH method is particular useful for high-order AO systems.

5. Conclusion

In conclusion, a new multichannel-Hadamard method for IM calibration has been presented. The method is based on the sparse characteristic of the IM and the noise rejection effect of the voltage-oriented Hadamard method. The IM calibration error has been defined and validated by the experiments. The excellent slope error variances, phase error variances and focal plane intensity performance are also demonstrated by experiments. The MH calibration method is useful for high-order AO systems thanks to its strong anti-noise ability. It may either work independently or be used to provide an initial IM for other closed-loop or iteration calibration methods

Acknowledgments

This work has been supported by the National Natural Science Joint Foundation (NO. 11178004), the Hi-tech project of China and the Graduate Student Innovation Foundation of the Institute of Optics and Electronics, Chinese Academy of Sciences (2013.2).

References and links

1.

J. W. Hardy, Adaptive Optics for Astronomical Telescopes (Oxford University, 1998), Chap. 1.

2.

H. W. Babcock, “Adaptive optics revisited,” Science 249(4966), 253–257 (1990). [CrossRef] [PubMed]

3.

C. Boyer, V. Michau, and G. Rousset, “Adaptive optics: Interaction matrix measurements and real time control algorithms for the COME-ON project,” Proc. SPIE 1237, 406–421 (1990). [CrossRef]

4.

W. Jiang and H. Li, “Hartmann-Shack wavefront sensing and wavefront control algorithm,” Proc. SPIE 1271, 82–93 (1990). [CrossRef]

5.

A. Bouchez, R. Dekany, J. Angione, C. Baranec, M. Britton, K. Bui, R. Burruss, J. Cromer, S. Guiwits, J. Henning, J. Hickey, D. Mckenna, A. Moore, J. Roberts, T. Trinh, M. Troy, T. Truong, and V. Velur, “The PALM-3000 high-order adaptive optics system for Palomar Observatory,” Proc. SPIE 7015, 70150Z (2008). [CrossRef]

6.

J.-F. Sauvage, T. Fusco, C. Petit, S. Meimon, E. Fedrigo, M. Suarez Valles, M. Kasper, N. Hubin, J.-L. Beuzit, J. Charton, A. Costille, P. Rabou, D. Mouillet, P. Baudoz, T. Buey, A. Sevin, F. Wildi, and K. Dohlen, “SAXO, the eXtreme Adaptive Optics System of SPHERE. Overview and calibration procedure,” Proc. SPIE 7736, 77360F (2010). [CrossRef]

7.

S. J. Thomas, L. Poyneer, R. De Rosa, B. Macintosh, D. Dillon, J. K. Wallace, D. Palmer, D. Gavel, B. Bauman, L. Saddlemyer, and S. Goodsell, “Integration and test of the Gemini Planet Imager,” Proc. SPIE 8149, 814903 (2011). [CrossRef]

8.

S. Oberti, F. Quirós-Pacheco, S. Esposito, R. Muradore, R. Arsenault, E. Fedrigo, M. Kasper, J. Kolb, E. Marchetti, A. Riccardi, C. Soenke, and S. Stroebele, “Large DM AO systems: synthetic IM or calibration on sky?” Proc. SPIE 6272, 627220 (2006). [CrossRef]

9.

M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, and S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]

10.

S. Oberti, H. Bonnet, E. Fedrigo, L. Ivanescu, M. Kasper, and J. Paufique, “Calibration of a curvature sensor/bimorph mirror AO system: interaction matrix measurement on MACAO systems,” Proc. SPIE 5490, 139–150 (2004). [CrossRef]

11.

A. Riccardi, R. Briguglio, E. Pinna, G. Agapito, F. Quiros-Pacheco, and S. Esposito, “Calibration strategy of the pyramid wavefront sensor module of ERIS with the VLT deformable secondary mirror,” Proc. SPIE 8447, 84475M (2012). [CrossRef]

12.

S. Meimon, T. Fusco and C. Petit, “An optimized calibration strategy for high order adaptive optics systems: the Slope-Oriented Hadmard Actuation,” in 1st AO4ELT Conference (2010), paper 07009.

13.

W. Zou and S. A. Burns, “High-accuracy wavefront control for retinal imaging with Adaptive-Influence-Matrix Adaptive Optics,” Opt. Express 17(22), 20167–20177 (2009). [CrossRef] [PubMed]

14.

F. Wildi and G. Brusa, “Determining the interaction matrix using starlight,” Proc. SPIE 5490, 164–173 (2004). [CrossRef]

15.

S. Esposito, R. Tubbs, A. Puglisi, S. Oberti, A. Tozzi, M. Xompero, and D. Zanotti, “High SNR measurement of interaction matrix on-sky and in lab,” Proc. SPIE 6272, 62721C (2006). [CrossRef]

16.

E. Pinna, F. Quiros-Pacheco, A. Riccardi, R. Briguglio, A. Puglisi, L. Busoni, C. Arcidiacono, J. Argomedo, M. Xompero, E. Marchetti, and S. Esposito, “First on-sky calibration of a high order adaptive optics system,” Proc. SPIE 8447, 84472B (2012). [CrossRef]

17.

X. Zhang, C. Arcidiacono, A. R. Conrad, T. M. Herbst, W. Gaessler, T. Bertram, R. Ragazzoni, L. Schreiber, E. Diolaiti, M. Kuerster, P. Bizenberger, D. Meschke, H. W. Rix, C. Rao, L. Mohr, F. Briegel, F. Kittmann, J. Berwein, and J. Trowitzsch, “Calibrating the interaction matrix for the LINC-NIRVANA high layer wavefront sensor,” Opt. Express 20(7), 8078–8092 (2012). [PubMed]

OCIS Codes
(010.1080) Atmospheric and oceanic optics : Active or adaptive optics
(090.1000) Holography : Aberration compensation
(010.1285) Atmospheric and oceanic optics : Atmospheric correction

ToC Category:
Adaptive Optics

History
Original Manuscript: March 6, 2014
Revised Manuscript: May 19, 2014
Manuscript Accepted: May 21, 2014
Published: May 30, 2014

Citation
Youming Guo, Changhui Rao, Hua Bao, Ang Zhang, Xuejun Zhang, and Kai Wei, "Multichannel-Hadamard calibration of high-order adaptive optics systems," Opt. Express 22, 13792-13803 (2014)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-22-11-13792


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. W. Hardy, Adaptive Optics for Astronomical Telescopes (Oxford University, 1998), Chap. 1.
  2. H. W. Babcock, “Adaptive optics revisited,” Science 249(4966), 253–257 (1990). [CrossRef] [PubMed]
  3. C. Boyer, V. Michau, G. Rousset, “Adaptive optics: Interaction matrix measurements and real time control algorithms for the COME-ON project,” Proc. SPIE 1237, 406–421 (1990). [CrossRef]
  4. W. Jiang, H. Li, “Hartmann-Shack wavefront sensing and wavefront control algorithm,” Proc. SPIE 1271, 82–93 (1990). [CrossRef]
  5. A. Bouchez, R. Dekany, J. Angione, C. Baranec, M. Britton, K. Bui, R. Burruss, J. Cromer, S. Guiwits, J. Henning, J. Hickey, D. Mckenna, A. Moore, J. Roberts, T. Trinh, M. Troy, T. Truong, V. Velur, “The PALM-3000 high-order adaptive optics system for Palomar Observatory,” Proc. SPIE 7015, 70150Z (2008). [CrossRef]
  6. J.-F. Sauvage, T. Fusco, C. Petit, S. Meimon, E. Fedrigo, M. Suarez Valles, M. Kasper, N. Hubin, J.-L. Beuzit, J. Charton, A. Costille, P. Rabou, D. Mouillet, P. Baudoz, T. Buey, A. Sevin, F. Wildi, K. Dohlen, “SAXO, the eXtreme Adaptive Optics System of SPHERE. Overview and calibration procedure,” Proc. SPIE 7736, 77360F (2010). [CrossRef]
  7. S. J. Thomas, L. Poyneer, R. De Rosa, B. Macintosh, D. Dillon, J. K. Wallace, D. Palmer, D. Gavel, B. Bauman, L. Saddlemyer, S. Goodsell, “Integration and test of the Gemini Planet Imager,” Proc. SPIE 8149, 814903 (2011). [CrossRef]
  8. S. Oberti, F. Quirós-Pacheco, S. Esposito, R. Muradore, R. Arsenault, E. Fedrigo, M. Kasper, J. Kolb, E. Marchetti, A. Riccardi, C. Soenke, S. Stroebele, “Large DM AO systems: synthetic IM or calibration on sky?” Proc. SPIE 6272, 627220 (2006). [CrossRef]
  9. M. Kasper, E. Fedrigo, D. P. Looze, H. Bonnet, L. Ivanescu, S. Oberti, “Fast calibration of high-order adaptive optics systems,” J. Opt. Soc. Am. A 21(6), 1004–1008 (2004). [CrossRef] [PubMed]
  10. S. Oberti, H. Bonnet, E. Fedrigo, L. Ivanescu, M. Kasper, J. Paufique, “Calibration of a curvature sensor/bimorph mirror AO system: interaction matrix measurement on MACAO systems,” Proc. SPIE 5490, 139–150 (2004). [CrossRef]
  11. A. Riccardi, R. Briguglio, E. Pinna, G. Agapito, F. Quiros-Pacheco, S. Esposito, “Calibration strategy of the pyramid wavefront sensor module of ERIS with the VLT deformable secondary mirror,” Proc. SPIE 8447, 84475M (2012). [CrossRef]
  12. S. Meimon, T. Fusco and C. Petit, “An optimized calibration strategy for high order adaptive optics systems: the Slope-Oriented Hadmard Actuation,” in 1st AO4ELT Conference (2010), paper 07009.
  13. W. Zou, S. A. Burns, “High-accuracy wavefront control for retinal imaging with Adaptive-Influence-Matrix Adaptive Optics,” Opt. Express 17(22), 20167–20177 (2009). [CrossRef] [PubMed]
  14. F. Wildi, G. Brusa, “Determining the interaction matrix using starlight,” Proc. SPIE 5490, 164–173 (2004). [CrossRef]
  15. S. Esposito, R. Tubbs, A. Puglisi, S. Oberti, A. Tozzi, M. Xompero, D. Zanotti, “High SNR measurement of interaction matrix on-sky and in lab,” Proc. SPIE 6272, 62721C (2006). [CrossRef]
  16. E. Pinna, F. Quiros-Pacheco, A. Riccardi, R. Briguglio, A. Puglisi, L. Busoni, C. Arcidiacono, J. Argomedo, M. Xompero, E. Marchetti, S. Esposito, “First on-sky calibration of a high order adaptive optics system,” Proc. SPIE 8447, 84472B (2012). [CrossRef]
  17. X. Zhang, C. Arcidiacono, A. R. Conrad, T. M. Herbst, W. Gaessler, T. Bertram, R. Ragazzoni, L. Schreiber, E. Diolaiti, M. Kuerster, P. Bizenberger, D. Meschke, H. W. Rix, C. Rao, L. Mohr, F. Briegel, F. Kittmann, J. Berwein, J. Trowitzsch, “Calibrating the interaction matrix for the LINC-NIRVANA high layer wavefront sensor,” Opt. Express 20(7), 8078–8092 (2012). [PubMed]

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