Variability of wavefront aberration measurements in small pupil sizes using a clinical Shack-Hartmann aberrometer
© Ginis et al 2004
Received: 02 October 2003
Accepted: 11 February 2004
Published: 11 February 2004
Skip to main content
© Ginis et al 2004
Received: 02 October 2003
Accepted: 11 February 2004
Published: 11 February 2004
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.
A clinical aberrometer (COAS Wavefront Scienses, Ltd) based on the Shack-Hartmann principle was employed in this study. Fifty consecutive measurements were perfomed 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.
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.
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.
Since the application of excimer laser technology for the correction of eye's simple refractive errors (i.e. 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, 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 ocular aberrations (e.g. the Spatially Resolved Refractometer (SRR)[5, 6], the Tscherning aberrometer and the Retinal ray-tracing, 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 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 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.
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 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.
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.
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. 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).
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 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 implemented in a MATLAB file, as previously described by Campbell.
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.
The two lower graphs depict similar plots for two of the subjects tested. The variance for the C2 0 is much higher in this case, as would be expected due to the dynamic nature of the 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 dispersion 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 C4 0 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 dispersion.
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-aberration error. Some of this increased variance may have been due to poorer image quality in parts of the Shack-Hartmann images. 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 displacements in respect to the optical axis of the instrument cannot justify significant fluctuations of the wavefront aberration, 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 corneal curvature.
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 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 C2 0, are measured with higher variability (see figure 5), and this is an interesting feature as the coefficients C2 0, C2 -2, C2 +2 can be used in the calculation of the conventional sphero-cylindrical correction. 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 C2 0 (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 C2 0 is not related to the inherent noise of the instrument itself, but to the dynamic characteristics of the human visual system, such as accommodation micro-fluctuations[42, 43], tear film changes, 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. 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].
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 C2 0 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.
All Zernike coefficients, except the coefficient of defocus (C2 0), 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, where the wavelength λ is written in nm.
nλ = 1.320535 - (4.685/ (λ - 214.102)) (2)
The C2 0 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. 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 C2 0 (in OSA format) from the spherical equivalent (S) for a specific pupil diameter (d) at each wavelength.
This work was supported by the EU, through the Research Training Network "Adaptive optics for retinal imaging and improved vision" under contract number HPRN-CT-2002-00301. Authors would like to thank Oleksandr Lutsenko for his involvement in data acquisition and Professor Ioannis Pallikaris for his valuable comments.
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.