Variability of wavefront aberration measurements in small pupil sizes using a clinical Shack-Hartmann aberrometer.

BACKGROUND: Recently, instruments for the measurement of wavefront aberration in the living human eye have been widely available for clinical applications. Despite the extensive background experience on wavefront sensing for research purposes, the information derived from such instrumentation in a clinical setting should not be considered a priori precise. We report on the variability of such an instrument at two different pupil sizes.
METHODS: A clinical aberrometer (COAS Wavefront Scienses, Ltd) based on the Shack-Hartmann principle was employed in this study. Fifty consecutive measurements were performed on each right eye of four subjects. We compared the variance of individual Zernike expansion coefficients as determined by the aberrometer with the variance of coefficients calculated using a mathematical method for scaling the expansion coefficients to reconstruct wavefront aberration for a reduced-size pupil.
RESULTS: Wavefront aberration exhibits a marked variance of the order of 0.45 microns near the edge of the pupil whereas the central part appears to be measured more consistently. Dispersion of Zernike expansion coefficients was lower when calculated by the scaling method for a pupil diameter of 3 mm as compared to the one introduced when only the central 3 mm of the Shack - Hartmann image was evaluated. Signal-to-noise ratio was lower for higher order aberrations than for low order coefficients corresponding to the sphero-cylindrical error. For each subject a number of Zernike expansion coefficients was below noise level and should not be considered trustworthy.
CONCLUSION: Wavefront aberration data used in clinical care should not be extracted from a single measurement, which represents only a static snapshot of a dynamically changing aberration pattern. This observation must be taken into account in order to prevent ambiguous conclusions in clinical practice and especially in refractive surgery.

Background

Since the application of excimer laser technology for the correction of eye's simple refractive errors (i.e. <span class="Disease">defocus and astigmatism), there has been a considerable debate concerning the visual impact of correcting the higher order monochromatic aberrations of the eye (e.g. spherical aberration, <span class="Disease">coma and secondary astigmatism), which also degrade retinal image quality [1-4]. Advances in the measurement of the eye's wave aberration have led to the emergence of sophisticated instrumentation for the clinical evaluation of the ocular higher-order aberrations. In general, these devices typically represent the aberrations as a wavefront-error map at the corneal or pupil plane.

Among other subjective and objective techniques which are now available for measuring <span class="Disease">ocular aberrations (e.g. the <span class="Chemical">Spatially Resolved Refractometer (SRR)[5,6], the Tscherning aberrometer[7] and the Retinal ray-tracing[8], Shack-Hartmann based instruments [9-11] have become the most widely adopted. Aberrometry is rapidly making its way into the clinic and has already been applied in measuring aberrations of normal[12,13] or clinically abnormal eyes (eg dry, keratoconic) [14-16], eyes undergone refractive surgery [17-20], as well as in situ aberration structures of soft and RGP contact lenses [21,22] and intraocular lenses (IOLs)[23,24].

Recently, the clinical measurement of higher-order aberrations has become important for <span class="Species">patient care. There is currently a major ongoing effort to refine laser refractive surgery, with the aim to eliminate higher order aberrations. In principle, wavefront aberration is measured using devices such as the Hartmann-Shack wavefront sensor. This information is then fed to a computer that generates the excimer laser's scanning pattern to allow simple refractive errors as well as higher order aberrations to be corrected. Preliminary results are still tentative as no clinically significant difference between conventional and wavefront-guided ablations has been demonstrated[17,25-27].

Moreover, there have been studies using state-of-the-art aberrometers to evaluate the refractive state[12,13,28] and the accommodative response[29,30] of the <span class="Species">human eye. Therefore, an obvious requirement of each of these devices is accuracy and repeatability of the measurement of the low order (sphero-cylindrical error) as well as the higher-order aberrations of the eye. Several studies have addressed the accuracy and the repeatability of static measurements of wavefront aberration[10,11,31,32]. These studies have shown that, although there is some variation arising from a combination of misalignment errors and small drifts in the measuring equipment, these are well beyond the clinicians' normal operation range. However, the use of a single measurement of the wavefront error in the planning of a custom correction is not recommended[32].

In this study, we used a clinical aberrometer to evaluate the variability of low and higher order aberrations at different pupil sizes, since no standard pupil has been established for reporting <span class="Disease">ocular aberrations. Morever, we compared the variance of individual Zernike coefficients as determined by the aberrometer with the variance of coefficients calculated with a matrix method that reconstructs a new set of expansion coefficients appropriate for any reduced-size pupil.

Methods

Subjects

Four right eyes of four subjects aged between 23 and 33 years (mean age: 29.2 years) were tested. One subject was emmetropic (AP: plano), two were low myopes (SP: -2.00 / -0.25 × 30, HG: -1.75 / -1.25 × 86), and one was an intermediate myope (OL: -4.75 / -0.25 × 10). None of the subjects had any ocular pathology or had undergone any kind of refractive surgery. Subjects were familiarised with the measurements. Prior to data collection, institutional research board approval was obtained.

Instrumentation

The monochromatic aberration function of the eye was measured using the Complete Ophthalmic Analysis System (COAS, Wavefront Sciences Ltd), which is based on the Shack-Hartmann principle as described elsewhere[9]. COAS also provides a real time display of the pupil image, which is used to objectively measure pupil size to the nearest 0.1 mm. COAS uses an 840 nm infrared super-luminescent diode as the light source and utilises a square lenslet array of 33 × 44 (a total of 1452 lenslets). The diameter of each lenslet is 144 μm. According to the manufacturers the pupil magnification factor is about 0.685, which means that the lenslet array samples the exiting wavefront every 210 μm in the pupil plane. This allows approximately 600 sample points within a 6.0 mm diameter pupil (150 sample point within a 3.0 mm diameter pupil), providing very high resolution sampling of the aberration. The software allows continuous recording of Shack – Hartmann images and pupil size with an exposure time of about 130 ms for each frame capture, i.e. a temporal frequency of 7.7 Hz. The data extracted from COAS, consist of a set of Zernike coefficients (up to 4th order) in Malacara format, that quantify the type and the magnitude of aberrations present. The ordering of Zernike coefficients was transposed to the format recommended by the Optical Society of America (OSA)[33].

Procedure

All measurements were performed on natural pupils without the use of any dilating or cycloplegic drug. Room illuminance was set at mesopic light levels. Large pupil analysis was based on the full pupil, which, at these lighting conditions, varied for each subject (ranged between 4.5 and 7.1 mm). The subject positioned his head on the chin rest and fixated on the centre of a circular grid. The operator manually aligned the subject's pupil centre with the optical axis of the device by means of six dots (that lie on a circle concentric with the pupil) displayed on a video monitor. This ensured that subject's line-of-sight was coaxial with the instrument's optical axis. A series of fifty consecutive measurements (total recording time 6.5 sec) for each eye were taken for the full pupil without re-alignment. Subjects were asked to blink prior to the measurement. In addition to subjects, we measured the wave aberration of an artificial eye supplied by the manufacturer as a test object.

Data analysis

Data analysis was performed using MATLAB (V 5.2, The Mathworks, Inc Natick, MA) mathematical software. The Zernike expansion coefficients derived from the wave inclination data for the full pupil, were initially transposed to the OSA format and then corrected for chromatic aberration (from 840 to 550 nm) (see Appendix A).

The corrected coefficients were "scaled" to a smaller pupil diameter (3 mm) using two different techniques: (i) the "direct" technique (the standard employed by COAS software), which re-calculated Zernike expansion coefficients (up to 4th order) after discarding the Shack-Hartmann image outside the 3 mm zone (ii) the "scaling" technique, which uses a matrix method to reconstruct a new set of Zernike coefficients that describe a wavefront aberration corresponding to the central 3 mm of the pupil using all available raw data. To achieve this, we used formulas developed by Schwiegerling[34] implemented in a MATLAB file, as previously described by Campbell[35].

Figure 1 depicts wavefront difference map of the central 3 mm of the full-size pupil and the calculated central 3 mm as obtained by the "direct" and the "scaling" methods. It is evident that the difference map for the "scaling" method is practically zero for every point of the entrance pupil. In contrast, the "direct" method produces a map that does not correspond closely to the initial data. Apparently, the calculated Zernike polynomial coefficients are different when peripheral data points are discarded.

Figure 1

Difference maps. Wavefront aberration difference maps between the central 3 mm pupil of the full-size pupil and the re-calculated central 3 mm as obtained by the "scaling" (left) and "direct" (right) methods. The vertical colour bar on the right shows corresponding wavefront aberration error in micrometers.

Difference maps. Wavefront aberration difference maps between the central 3 mm pupil of the full-size pupil and the re-calculated central 3 mm as obtained by the "scaling" (left) and "direct" (right) methods. The vertical colour bar on the right shows corresponding wavefront <span class="Disease">aberration error in micrometers.

Results

Variance of the wave aberration data for the full-size pupil

Prior to the presentation of the results for the scaled pupils it is of interest to report on the variance of the raw wavefront aberration data (50 measurements) for the full-size pupil of each subject. Figure 2 illustrates the spatial distribution of the standard deviation of the measured wavefront aberration as a function of pupil position.

Figure 2

Colour patterns of wavefront aberration at full-size pupils. Colour patterns of the standard deviation of the wavefront error as a function of horizontal (x) and vertical (y) pupil position for the full-size pupils of the four subjects tested. Map size is 150 × 150 pixels.

It is obvious that wavefront aberration variance at peripheral points of the pupil and especially near the edge of the pupil, is increased. Reasonably, this effect is more pronounced for subjects OL and SP having larger pupil diameters.

Comparison of the two scaling methods

Figure 3 shows the variation with time of Zernike expansion coefficients C20 and <span class="Gene">C40, determined for a pupil diameter of 3 mm. The upper graphs correspn>ond to fifty consecutive measurements made on an artificial eye. There is little doubt that discarding peripheral pupil data points in the "direct" method does not affect the variance of the low order coefficient, C20. However, it does introduce noise to the value determined for the higher order coefficient, <span class="Gene">C40. On the other hand, the "scaling" method produces minimal noise in both coefficients. It is noteworthy that noise introduced by the direct method is not observed in the coefficients corresponding to the large pupil.

Figure 3

Variation of aberration coefficients CVariation of wavefront aberration coefficient C20 (left) and C40 (right) with time for a 3 mm pupil as calculated by the direct (filled symbols) and scaling (open symbols) methods. Data from an artificial eye (upper graphs) and two subjects are shown. The dotted lines are least-square regression coefficients. Note, that the scaling of y axis is different for the artificial eye.

Variation of aberration coefficients CVariation of wavefront aberration coefficient C20 (left) and <span class="Gene">C40 (right) with time for a 3 mm pupil as calculated by the direct (filled symbols) and scaling (open symbols) methods. Data from an artificial eye (upper graphs) and two subjects are shown. The dotted lines are least-square regression coefficients. Note, that the scaling of y axis is different for the artificial eye.

The two lower graphs depict similar plots for two of the subjects tested. The variance for the C20 is much higher in this case, as would be expected due to the dynamic nature of the <span class="Species">human eye. Moreover, note that there is no difference between the two methods, which also supports the above statement. For the 4th order term, the "scaling" method improves substantially the dispn>ersion of the coefficient. Similar improvement in the variability of all higher order aberration coefficients is observed when the scaling method is applied for the calculation at 3 mm pupil.

Another point worth mentioning is that a drift in the value of the <span class="Gene">C40 with time indicated by the direct method (see subject SP) may be ambiguous, as this is not the case when the scaling method is used. Moreover, there is no evidence of change of any aberration coefficient with time, as no statistically significant correlation can be established. This implies that during time period of the recordings (6.5 sec) needed to capture the fifty measurements, subjects maintained fixation and wavefront changes related to the tear film quality as well as accommodation state did not interfere significantly with the estimated dispn>ersion.

Figures 4a and 4b show frequency histograms of the spherical aberration coefficient (<span class="Gene">C40) and the higher-order RMS error for all subjects tested. It is clear from these data that the dispn>ersion of the values is wider when the direct method is used, where as the scaling method calculates coefficients with higher variability. Furthermore, there is an obvious difference in the mean value, with a trend of the direct method to overestimate the amount of the wavefront error.

Figure 4

Frequency histograms of C(a) Frequency histograms of the spherical aberration coefficient, C40 for all subjects tested. Comparison between results derived from the direct (upper) and scaling (lower) methods for a 3 mm pupil. Bin width is 0.001. (b) Frequency histograms of the higher-order RMS error for all subjects tested. Comparison between results derived from the direct (upper) and scaling (lower) methods for a 3 mm pupil. Bin width is 0.0025.

Frequency histograms of C(a) Frequency histograms of the spherical aberration coefficient, <span class="Gene">C40 for all subjects tested. Comparison between results derived from the direct (upper) and scaling (lower) methods for a 3 mm pupil. Bin width is 0.001. (b) Frequency histograms of the higher-order RMS error for all subjects tested. Comparison between results derived from the direct (upper) and scaling (lower) methods for a 3 mm pupil. Bin width is 0.0025.

Impact of error on different aberration terms and radial orders

Figure 5 presents signal-to-noise (S/N) ratio charts (mean / SD) for different radial orders of the wavefront aberration for all subjects tested. High S/N values imply high variability of the measured order of aberration. For example, subject's HG 2nd order terms are measured to have a value about thirty times higher than their standard deviation (noise level), whereas the magnitude of the 4th order terms is only five times higher than the standard deviation. There are three points to note from these data. First, the signal-to-noise ratios are much higher for the 2nd order compared to the higher-order coefficients for all the subjects. Second, there is some substantial variation between subjects for the 2nd order coefficients, while this is not the case for the higher order terms. This is mainly due to the fact that the amount of the sphero-cylindrical error (determining the 2nd order terms) varies between the subjects, while there is little variation in the magnitude of 3rd and 4th order coefficients. Third, the scaling method improves the signal-to-noise ratio in all cases, but this is more profound in the higher-order terms, as also displayed in figure 3.

Figure 5

Signal-to-noise ratios for different radial order. Signal-to-noise ratios (mean / SD) for different radial orders of the wavefront aberration at 3 mm pupil as calculated by the direct and scaling methods. Data for all subjects are shown.

The improvement of S/N for each expansion coefficient is depicted in figure 6, which uses the standard graph that describes the radial and angular modes of the Zernike polynomial expansion. Determination of coefficients having values smaller than two times the noise level, (S/N < 2) are below the detection limit, regarded imprecise and are therefore displayed by red color. This graph corresponds to measurements of subject OL. We derived similar results for the other subjects.

Figure 6

Signal-to-noise ratio for individual Zernike coefficients. Signal-to-noise (S/N) ratio in the standard pyramidical layout of Zernike expansion coefficients for one subject (OL) as calculated by the two methods. Green colours show high S/N ratio, red colours show low S/N ratio.

Discussion

The aim of this study was to use a Shack-Hartmann clinical aberrometer (COAS, Wavefront Sciences Ltd) to evaluate the variability of low- and higher-order aberration measurement of the eye. Using the wave aberration polynomial determined for a full-size pupil we compared the Zernike expansion coefficients for a smaller, 3 mm pupil derived by two methods: first, by re-calculating the wave aberration coefficients to the reduced sampled points corresponding to that pupil ("direct" method), and, second, by using a matrix method to reconstruct a new set of coefficients appropriate for the reduced pupil ("scaling" method) (see figure 1).

Our results suggest that, for full-size pupil, the efficiency of the measurements varies across x and y pupil position: where the wavefront is larger, measurement variance is higher, especially near the margins of the pupil, where increased standard deviation results to higher wave-<span class="Disease">aberration error. Some of this increased variance may have been due to poorer image quality in parts of the Shack-Hartmann images[36]. Also, it may be partially attributed to the fact that saccadic movements, during the time required for data collection, lead to alignment errors that continuously change the set of sensor elements contributing to wavefront sensing. Although such dispn>lacements in respn>ect to the optical axis of the instrument cannot justify significant fluctuations of the wavefront aberration[37], we cannot exclude the possibility that the algorithm employed in COAS software may generate the increased noise in periphery during wavefront expansion, since pupil translation magnitude (~100 μm) is comparable to lenslet array spacing (as magnified by the conjugating optics) in the particular instrument.

The increased standard deviation of wavefront aberration at the periphery has implications in calculating wavefront-guided ablation patterns. An error of 0.45 μm in the measured wavefront aberration at the periphery of the treatment zone may lead to a substantial error in the calculated shot pattern depending on the laser beam delivery (scanning) method as well as beam parameters and compensation for <span class="Disease">corneal curvature[38].

As pupil becomes smaller, the magnitude of wavefront aberrations decreases. At 3 mm pupil, the "direct" method (employed by COAS) induces considerable variance in the measurement of higher-order aberration coefficients attributed to inherent <span class="Disease">fitting error. This results from the small number of sensor elements involved in the wavefront inclination measurement. Moreover, there is a clear shift in the magnitude of each coefficient to higher values, which may lead to inaccurate determination of higher-order terms. In contrast, the "scaling" method produces coefficients of higher variability, and this is not surprising since it allows the use of the information from a larger set of sensor elements, reducing instrument noise.

On the other hand, second-order terms, and especially C20, are measured with higher variability (see figure 5), and this is an interesting feature as the coefficients C20, <span class="Gene">C2-2, C2+2 can be used in the calculation of the conventional sphero-cylindrical correction[33]. This observation, in conjuction with the high accuracy of the refraction estimated objectively from wavefront aberration data when small (~3 mm) pupils are used[13,39,40], is of considerable importance, as this means that the COAS clinical aberrometer may be used as a reliable and accurate autorefractor. However, care must be taken when large pupils are tested, as the objective estimation of refractive error may lead to ambiguous results, due to the influence of higher-order aberrations on the determination of correction[40,41].

Another finding is that the variability of C20 (corresponding to defocus) improves only slightly when the "scaling" method is used. This probably occurs because the defocus term is mostly contained in the Shack-Hartmann spots at the centre of pupil. In contrast, spherical aberration, for example, depends on the 4th power of pupil radius, which means that most of the information is outside the central 3 mm and what is measured with a small pupil is mostly noise. Another reason may be the fact that the variance in C20 is not related to the inherent noise of the instrument itself, but to the dynamic characteristics of the <span class="Species">human visual system, such as accommodation micro-fluctuations[42,43], tear film changes[44], and eye movements leading to alignment errors[32,45].

Although instruments based on the Shack – Hartmann sensor have been extensively validated for experimental work[11,13,31,32,40], our results indicate that special care should be taken when measurement of aberration is used in clinical applications, such as refractive surgery, either decision making or outcome evaluation. The diversity of the measured values of various coefficients suggest that a number of measurements should be taken and averaged for each subject in order to calculate coefficients of higher efficiency[32]. This is of major importance in customized laser surgery, where aberration data are used to correct higher-order aberrations for the potential enhancement of visual performance [46-48].

Conclusions

Wavefront aberration of the eye, as derived from Shack-Hartmann images is determined with a certain degree of accuracy that varies considerably with pupil position. Zernike expansion coefficients are determined with less accuracy when re-calculated at "cropped" pupils with the use of the algorithm employed in COAS. This study shows that these errors attributed to the reduced number of sensor elements could be, at least partially, overcome using an appropriate algorithm that calculates the aberration coefficients for smaller pupils based on full-size pupil set of data. Moreover, micro – fluctuations observed in the C20 corresponding to defocus, are probably inherent characteristics of the eye and therefore show no improvement when the algorithm is applied.

Consequently, it must be emphasised that wavefront aberration data used in clinical care should not be extracted from a single measurement, which represents only a static snapshot of a dynamically changing aberration pattern.

Appendix A

Correcting Zernike coefficients for chromatic aberration

All Zernike coefficients, except the coefficient of defocus (C20), were adjusted to appropriate values for 550 nm, using a chromatic correction factor, K (equation 1):

K = (n1-1)/(n2-1)     (1)

Where n1 and n2 are the refractive index values for 550 nm and 840 nm, respondingly. These were calculated using equation 2[49], where the wavelength λ is written in nm.

nλ = 1.320535 - (4.685/ (λ - 214.102))     (2)

The C20 was corrected for wavelength, using equations (3) and (4). This is necessary because infrared light used for measurements is not reflected from the photoreceptors plane, where the subjective focal plane lies, but passes through to deeper layers and is reflected from the choroid[50]. Equation (3) calculates the spherical equivalent power corrected at 550 nm (S550), assuming that the 840-nm light is reflected from a plane 0.125 mm posterior to the retina of a 60D model eye (with a focal length of 16.667 mm). Equation (4) is used to derive C20 (in OSA format) from the spherical equivalent (S) for a specific pupil diameter (d) at each wavelength.

Competing interests (medicine)

None declared.

Authors' contributions

HG participated in data analysis, manuscript preparation and study design. SP carried out measurements, performed data analysis and participated in manuscript and figure preparation. AP performed wavefront measurements, implemented the MATLAB code, participated in figure preparation and created the attached computer program. All authors read and approved the final manuscript.

Pre-publication history

The pre-publication history for this paper can be accessed here:

Additional File 1

This zip file containes two MATLAB (v 5.2) files used for Zernike coefficient scaling to different pupil sizes. To run these files, MATLAB needs to be installed. This was tested for versions 5.2 for Macintosh and 6.1 for windows. Instructions for use as well as author credits are included in the files.

Click here for file

Additional File 2

This zip file contains three files that form a Windows utility (an executable) used for Zernike coefficient scaling. Instructions for use, as well as author credits and a disclaimer, are included in the files.

Click here for file