Adaptive Optics Systems (AOSs) have been widely used for retinal imaging to obtain high resolution images of photoreceptor mosaics and retinal capillaries[
Chinese Optics, Volume. 15, Issue 1, 79(2022)
Longitudinal chromatic aberration compensation method for dual-wavelength retinal imaging adaptive optics systems
Dual-wavelength retinal imaging adaptive optics systems are suitable for high contrast and resolution imaging of retinal capillaries. The compensation of the Longitudinal Chromatic Aberrations (LCAs) in dual-wavelength adaptive systems is researched. The LCA is measured, the measured wavefronts are analyzed, and the arbitrary wavefront LCA compensation method is given. An adaptive correction experiment is carried out and the experimental results indicate that the root mean square error of the wavefront is reduced to 0.16 λ (λ=589 nm) and the retinal capillary resolution is improved to 6 μm. This work may be used for the clinical applications of retinal imaging.
1 Introduction
Adaptive Optics Systems (AOSs) have been widely used for retinal imaging to obtain high resolution images of photoreceptor mosaics and retinal capillaries[
Newton first described the LCA of the human eye hundreds of years ago[
2 Longitudinal chromatic aberration measurement
2.1 Optical layout
To compensate for the LCA, the AOS should be first measured and analyzed. An optical setup was established in the laboratory, as shown in Fig. 1. Two lasers with wavelengths of 589 and 808 nm were selected as the dual-wavelength optical sources. The collimated beam is reflected to the eye by BS3 and focused on the retina, then reflected out from the eye. The spatial coherence of the laser was decreased by the diffuser to reduce the speckle artifacts. An annular aperture was placed at the conjugated plane of the pupil to eliminate the reflected light from the cornea. The light reflected from the retina passed through L7 and L8 and went into a Shack Hartmann wavefront sensor (SH-WFS) for wavefront detection. The SH-WFS, fabricated in-house, had a microlens array of 20×20, was 3 mm in diameter, and had an acquistition frequency of 515 Hz. A pupil camera was used to acquire the correct position of the pupil. A target illuminated by LED was used for fixation.
Figure 1.Illustration of the optical setup for LCA measurement: two lasers are used with the wavelength of 589 nm and 808 nm
The LCA is the difference between two wavefronts acquired at two different wavelengths. Ocular aberration does not change within 50 ms [
Figure 2.Time sequence for system control
2.2 Measured results
Five subjects, between the ages of 24 and 28, with healthy eyes with myopia with degrees from 100 to 500, were selected as subjects to measure the LCA. A wavelength of 589 nm was used to image their retinal capillaries; thus, all the aberrations were expressed with λ=589 nm. Fig.3 (Color online) shows the measured wavefronts of subject A. The LCA was obtained by subtracting the aberration of 808 nm from the aberration of 589 nm and its Root Mean Square (RMS) value was 1.3 λ, as shown in Fig. 3(c). Fig. 3(d) shows the LCA without the defocus and its RMS value is 0.31 λ, which corresponds to 6 times the diffraction limit. Therefore, although the defocus is the main component of LCA, the system must compensate for the residual LCA to acquire high-resolution images of the retina. The LCAs of different subjects are shown in Fig.4. The magnitudes of the LCA of other subjects are similar to that of subject A. The defocus can be overcome by moving the imaging camera; thus, we only consider the compensation of LCA without the defocus in the following.
Figure 3.Measured wavefronts of subject A. (a) Wavelength of 589 nm; (b) wavelength of 808 nm; (c) LCA; (d) LCA without defocus
Figure 4.Measured LCA for different subjects, with and without defocus at λ=589 nm
3 Arbitrary wavefront LCA compensation method
For the dual-wavelength retinal imaging of the AOS, the aberration was detected at the 808 nm wavelength and corrected with the wavefront corrector (WFC). A light source with a wavelength of 589 nm was used for retinal capillary imaging and the LCA was produced and should be corrected by the WFC simultaneously and the LCA should be detected and added to the aberration detected at 808 nm. The LCA of the eye varies slowly; thus, it should be measured and selected correctly for effective compensation. Consequently, the variation of LCA with time should be considered first. We measured the LCA of subject A at time intervals of 5 minutes, 1 hour, 10 hours, 15 hours, 24 hours, 30 hours, and 36 hours, as shown in Fig. 5 (Color online). The wavefronts are similar to a certain extent, but the LCA changes at different times. The RMS of LCA as a function of the time is shown in Fig.6. It indicates that the RMS of LCA changes minimally with time. The mean value and standard deviation of LCA are 0.32 λ and 0.032 λ, respectively.
Figure 5.Wavefronts of LCA at different times for subject A. (a) Start; (b) 5 minutes later; (c) 1 hour later; (d) 10 hours later; (e) 15 hours later; (f) 24 hours later; (g) 30 hours later; (h) 36 hours later
Figure 6.Measured LCA at different times for subject A
Because the wavefronts of LCA measured at different times have a certain similarity (shown in Fig.5) and the RMS value of LCA almost has no change (shown as Fig. 6), we select the wavefront of aberration measured at an arbitrary time to represent the timely detected LCA. Thus, before adaptive correction, the LCA should be measured first; then, while the adaptive correction is started, the measured LCA is combined into the real-time detected aberration and will be corrected by the WFC. Hence, after adaptive correction, the residual of LCA may be calculated with:
Which we call the arbitrary wavefront method. To simulate the correction effect of LCA, the first wavefront in Fig. 5 was selected as the arbitrary wavefront and the compensation results are shown in Fig. 7 (Color online). It is shown that, after compensation, the RMS value of the residual light is between 0.13 λ and 0.18 λ, the mean value of the LCA is 0.16 λ with a standard deviation of 0.017 λ.
By following the same process, other wavefronts in Fig. 5 were chosen successively as the arbitrary wavefront and the mean value of the RMS of residual are shown in Fig. 8. It indicates that after compensation, the mean RMS of LCA decreased from 0.3 λ to about 0.15 λ. The mean RMS of LCA at different times was averaged again and the computed result shows that, the averaged RMS of LCA is decreased to 0.158±0.011 λ. The LCA compensations of other subjects were computed and their mean values are shown in Fig. 9. It indicates that, after compensation, the mean RMS values of LCA are reduced to 0.166±0.017 λ for different subjects. Consequently, the arbitrary wavefront compensation method is a simple and effective method of compensating for LCA in a dual-wavelength AOS.
Figure 7.Calculated compensation of LCA for subject A at different times while the first wavefront of Fig. 5 is chosen as the arbitrary wavefront. The mean value of LCA is 0.16 λ with the standard deviation of 0.017. (a) Start; (b) 5 minutes later; (c) 1 hour later; (d) 10 hours later; (e) 15 hours later; (f) 24 hours later; (g) 30 hours later; (h) 36 hours later
Figure 8.Compensation of LCA for subject A at different times. The mean value of the arbitrary wavefront is 0.158±0.011 λ.
Figure 9.Compensation of LCA for different subjects. The mean values of the arbitrary wavefront is 0.166±0.017 λ.
4 LCA compensation experiment
To validate the arbitrary wavefront compensation method of LCA, an AOS was built, as shown in Fig. 10 (Color online), and mainly consists of the illumination, aberration correction and imaging, pupil monitoring, and target fixation. The illumination system included two lasers (CNI Inc.) with wavelengths of 589 and 808 nm, which were used for imaging the retinal capillaries and for wavefront detection. An annular aperture was placed at the conjugated position of the pupil to eliminate the reflected light from the cornea. The incidence optical power was lower than the maximum permissible exposure of the ANSI[
Figure 10.Optical layout for the retinal imaging AOS: L1-L13, Lens 1- Lens 13; PBS, polarizing beam splitter; BS, beam splitter; an 808 nm laser is used for wavefront detection, tracing and positioning the capillary; a 589 nm laser is used for high contrast imaging of the capillary; the collimated beam comes from the illumination system and is reflected into the eye, and then reflected again out from the eye by the retina; this reflected light is detected and corrected by the adaptive optics system and imaged with an imaging camera; the pupil position is observed by the pupil monitoring system and the eye is fixed with the target staring system
Figure 11.Experimental configuration of the AOS on an optical flat
First, the LCA was measured with the method described in Section 2. Then, an adaptive correction experiment was performed with the compensation of the LCA. To do the comparison, an adaptive correction was done without the compensation of LCA. A Field Of View (FOV) of 250 μm was selected for retinal imaging, and the resolution of the optical system was approximately evaluated with the diameter of the capillary[
Figure 12.Experiment results of adaptive correction and LCA compensation for subjects A and C. (a) Measured aberration at 808 nm for subject A; (b) wavefront of LCA for subject A; (c) image of retinal capillary without LCA compensation for subject A; (d) image of retinal capillary after LCA compensation for subject A; (e) measured aberration at 808 nm for subject C; (f) wavefront of LCA for subject C; (g) image of retinal capillary without LCA compensation for subject C; (h) image of retinal capillary after LCA compensation for subject C
5 Conclusion
An LCA compensation method is proposed for a dual-wavelength retinal imaging adaptive optics system. First, the LCAs of the human eye were measured in 5 subjects. An optical setup was established in the laboratory with dual wavelengths of 589 and 808 nm. The measured results showed that the RMS of LCA without defocus was 0.31 λ or so. Then, subject A was chosen as an example for the analysis. The LCA of subject A was measured at different times and results showed that its RMS changed minimally with time and the wavefronts of LCA measured at different times showed similarities. Hence, arbitrary wavefronts were selected as the aberrations to do the compensation. Calculated compensation results indicated that the LCA is reduced from 0.32 λ to 0.158 λ. Moreover, the comparison was performed for different subjects and the results showed that the mean RMS values of LCA are reduced to 0.166±0.017 λ for different subjects.
At last, an adaptive correction experiment was done for two subjects using the arbitrary wavefront LCA compensation method. The results showed that with LCA compensation, the RMS error of the wavefront decreased from 0.31 λ to approximately 0.16 λ. The resolutions of the retinal capillaries were improved from 13 μm and 14 μm to 6 μm and 7 μm for subjects A and C, respectively. Therefore, the LCA may be counteracted effectively with the proposed method. Furthermore, the proposed method may be applied to multi-wavelength retinal imaging AOS as well. This work is helpful for the LCA compensation of multi-wavelength retinal imaging AOS and acquiring high-resolution and high-contrast retinal images for clinical applications.
[1] LIANG J ZH, WILLIAMS D R, MILLER D T. Supernormal vision and high-resolution retinal imaging through adaptive optics[J]. Journal of the Optical Society of America A, 14, 2884-2892(1997).
[2] ROORDA A, WILLIAMS D R. The arrangement of the three cone classes in the living human eye[J]. Nature, 397, 520-522(1999).
[3] JONNAL R S, RHA J, ZHANG Y, et al.
[4] LI K Y, ROORDA A. Automated identification of cone photoreceptors in adaptive optics retinal images[J]. Journal of the Optical Society of America A, 24, 1358-1363(2007).
[5] DUBRA A, SULAI Y, NORRIS J L, et al. Noninvasive imaging of the human rod photoreceptor mosaic using a confocal adaptive optics scanning ophthalmoscope[J]. Biomedical Optics Express, 2, 1864-1876(2011).
[6] CHUI T Y P, VANNASDALE D A, BURNS S A. The use of forward scatter to improve retinal vascular imaging with an adaptive optics scanning laser ophthalmoscope[J]. Biomedical Optics Express, 3, 2537-2549(2012).
[7] ZAYIT-SOUDRY S, DUNCAN J L, SYED R, et al. Cone structure imaged with adaptive optics scanning laser ophthalmoscopy in eyes with nonneovascular age-related macular degeneration[J]. Investigative Ophthalmology & Visual Science, 54, 7498-7509(2013).
[8] QUERQUES G, KAMAMI-LEVY C, GEORGES A, et al. Appearance of regressing drusen on adaptive optics in age-related macular degeneration[J]. Ophthalmology, 121, 611-612(2014).
[9] PAQUES M, BROLLY A, BENESTY J, et al. Venous nicking without arteriovenous contact: the role of the arteriolar microenvironment in arteriovenous nickings[J]. JAMA Ophthalmology, 133, 947-950(2015).
[10] MARTIN J A, ROORDA A. Direct and noninvasive assessment of parafoveal capillary leukocyte velocity[J]. Ophthalmology, 112, 2219-2224(2005).
[11] [11] HAN G Q. High contrast imaging study of retinal capillaries in human eyes[D]. Beijing: Changchun Institute of Optics, Fine Mechanics Physics University of Chinse Academy of Sciences, 2018: 2830. (in Chinese)
[12] FABER D J, AALDERS M C G, MIK E G, et al. Oxygen saturation-dependent absorption and scattering of blood[J]. Physical Review Letters, 93, 028102(2004).
[13] REINHOLZ F, ASHMAN R A, EIKELBOOM R H. Simultaneous three wavelength imaging with a scanning laser ophthalmoscope[J]. Cytometry Part A, 37, 165-170(1999).
[14] GRIEVE K, TIRUVEEDHULA P, ZHANG Y H, et al. Multi-wavelength imaging with the adaptive optics scanning laser ophthalmoscope[J]. Optics Express, 14, 12230-12242(2006).
[15] [15] NEWTON S I. Opticks, , A Treatise of the Reflections, Refractions, in Flections & Colours of Light[M]. 4th ed. London: G. Bell & Sons, Ltd. , 1931.
[16] WALD G, GRIFFIN D R. The change in refractive power of the human eye in dim and bright light[J]. Journal of the Optical Society of America, 37, 321-336(1947).
[17] THIBOS L N, BRADLEY A, STILL D L, et al. Theory and measurement of ocular chromatic aberration[J]. Vision Research, 30, 33-49(1990).
[18] FERNÁNDEZ E J, ARTAL P. Ocular aberrations up to the infrared range: from 632.8 to 1070 nm[J]. Optics Express, 16, 21199-21208(2008).
[19] RUCKER F J, OSORIO D. The effects of longitudinal chromatic aberration and a shift in the peak of the middle-wavelength sensitive cone fundamental on cone contrast[J]. Vision Research, 48, 1929-1939(2008).
[20] NAKAJIMA M, HIRAOKA T, HIROHARA Y, et al. Verification of the lack of correlation between age and longitudinal chromatic aberrations of the human eye from the visible to the infrared[J]. Biomedical Optics Express, 6, 2676-2694(2015).
[21] VINAS M, DORRONSORO C, CORTES D, et al. Longitudinal chromatic aberration of the human eye in the visible and near infrared from wavefront sensing, double-pass and psychophysics[J]. Biomedical Optics Express, 6, 948-962(2015).
[22] VINAS M, DORRONSORO C, GARZÓN N, et al. In vivo subjective and objective longitudinal chromatic aberration after bilateral implantation of the same design of hydrophobic and hydrophilic intraocular lenses[J]. Journal of Cataract & Refractive Surgery, 41, 2115-2124(2015).
[23] CHONG S P, ZHANG T W, KHO A, et al. Ultrahigh resolution retinal imaging by visible light OCT with longitudinal achromatization[J]. Biomedical Optics Express, 9, 1477-1491(2018).
[24] ZAWADZKI R J, CENSE B, ZHANG Y, et al. Ultrahigh-resolution optical coherence tomography with monochromatic and chromatic aberration correction[J]. Optics Express, 16, 8126-8143(2008).
[25] DUBRA A, SULAI Y. Reflective afocal broadband adaptive optics scanning ophthalmoscope[J]. Biomedical Optics Express, 2, 1757-1768(2011).
[26] JIANG X Y, KUCHENBECKER J A, TOUCH P, et al. Measuring and compensating for ocular longitudinal chromatic aberration[J]. Optica, 6, 981-990(2019).
[27] ZHENG X L, LIU R X, XIA M L, et al. Temporal properties study of ocular wave aberrations with high frequency sampling[J]. Chinese Optics, 34, 0733001(2014).
[28] MORGAN J I W, HUNTER J J, MASELLA B, et al. Light-induced retinal changes observed with high-resolution autofluorescence imaging of the retinal pigment epithelium[J]. Investigative Ophthalmology & Visual Science, 49, 3715-3729(2008).
[29] DELORI F C, WEBB R H, SLINEY D H. Maximum permissible exposures for ocular safety (ANSI 2000), with emphasis on ophthalmic devices[J]. Journal of the Optical Society of America A, 24, 1250-1265(2007).
[30] [30] Laser Institute of America. ANSI Z136.12014 American national stard f safe use of lasers[S]. lo: American National Stards Institute, 2007: 2229.
[31] LI CH, XIA M L, MU Q Q, et al. High-precision open-loop adaptive optics system based on LC-SLM[J]. Optics Express, 17, 10774-10781(2009).
[32] YANG L B, HU L F, LI D Y, et al. Determining the imaging plane of a retinal capillary layer in adaptive optical imaging[J]. Chinese Physics B, 25, 094219(2016).
[33] BURNS S A, ELSNER A E, CHUI T Y, et al.
Get Citation
Copy Citation Text
Qin-yu ZHU, Guo-qing HAN, Jian-tao PENG, Qi-long RAO, Yi-li SHEN, Mei-rui CHEN, Hui-juan SUN, Hong-min MAO, Guo-ding XU, Zhao-liang CAO, Li XUAN. Longitudinal chromatic aberration compensation method for dual-wavelength retinal imaging adaptive optics systems[J]. Chinese Optics, 2022, 15(1): 79
Category: Original Article
Received: Sep. 6, 2021
Accepted: --
Published Online: Jul. 27, 2022
The Author Email: Zhao-liang CAO (caozl@usts.edu.cn)