Pulsed illumination spectral-domain optical coherence tomography for human retinal imaging

Similar documents
High-speed imaging of human retina in vivo with swept-source optical coherence tomography

Ultrahigh speed volumetric ophthalmic OCT imaging at 850nm and 1050nm

High-speed spectral-domain optical coherence tomography at 1.3 µm wavelength

Polarization-sensitive spectral-domain optical coherence tomography using a single line scan camera

Full-range k -domain linearization in spectral-domain optical coherence tomography

Ultrahigh-resolution high-speed retinal imaging using spectral-domain optical coherence tomography

Improved lateral resolution in optical coherence tomography by digital focusing using twodimensional numerical diffraction method

Pulsed-source and swept-source spectral-domain optical coherence tomography with reduced motion artifacts

Pulsed-source spectral-domain optical coherence tomography with reduced motion artifacts

Optical frequency domain imaging with a rapidly swept laser in the nm range

Volumetric microvascular imaging of human retina using optical coherence tomography with a novel motion contrast technique

Visualization of human retinal micro-capillaries with phase contrast high-speed optical coherence tomography

Phase-resolved optical frequency domain imaging

Ultrahigh Speed Spectral / Fourier Domain Ophthalmic OCT Imaging

High-speed optical frequency-domain imaging

Simultaneous acquisition of the real and imaginary components in Fourier domain optical coherence tomography using harmonic detection

MEMS tunable VCSEL light source for ultrahigh speed 60kHz - 1MHz axial scan rate and long range centimeter class OCT imaging

Single camera spectral domain polarizationsensitive optical coherence tomography using offset B-scan modulation

Removing the depth-degeneracy in optical frequency domain imaging with frequency shifting

60 MHz A-line rate ultra-high speed Fourier-domain optical coherence tomography

some aspects of Optical Coherence Tomography

In vivo three-dimensional microelectromechanical endoscopic swept source optical coherence tomography

Temporal coherence characteristics of a superluminescent diode system with an optical feedback mechanism

Medical imaging has long played a critical role in diagnosing

Optical coherence tomography

Spectral domain optical coherence tomography with balanced detection using single line-scan camera and optical delay line

Characteristics of point-focus Simultaneous Spatial and temporal Focusing (SSTF) as a two-photon excited fluorescence microscopy

Photoacoustic ophthalmoscopy for in vivo retinal imaging

Phase-stabilized optical frequency domain imaging at 1-µm for the measurement of blood flow in the human choroid

Axsun OCT Swept Laser and System

Improved spectral optical coherence tomography using optical frequency comb

doi: /OE

Isolator-Free 840-nm Broadband SLEDs for High-Resolution OCT

Common Path Side Viewing Monolithic Ball Lens Probe for Optical Coherence Tomography

Optimization for Axial Resolution, Depth Range, and Sensitivity of Spectral Domain Optical Coherence Tomography at 1.3 µm

Frequency comb swept lasers

Heterodyne swept-source optical coherence tomography for complete complex conjugate ambiguity removal

Frequency comb swept lasers for optical coherence tomography

Low-noise broadband light generation from optical fibers for use in high-resolution optical coherence tomography

Single-shot two-dimensional full-range optical coherence tomography achieved by dispersion control

Optical Coherence Tomography Systems and signal processing in SD-OCT

Citation for published version (APA): Nguyen, D. V. (2013). Integrated-optics-based optical coherence tomography

University of Lübeck, Medical Laser Center Lübeck GmbH Optical Coherence Tomography

Gabor fusion technique in a Talbot bands optical coherence tomography system

Frequency comb swept lasers

PHASE DETECTION WITH SUB-NANOMETER SEN- SITIVITY USING POLARIZATION QUADRATURE EN- CODING METHOD IN OPTICAL COHERENCE TOMOG- RAPHY

SUPPLEMENTARY INFORMATION

μoct imaging using depth of focus extension by self-imaging wavefront division in a commonpath fiber optic probe

Optical coherence angiography

Adaptive ranging for optical coherence tomography

University of Oulu, Finland.

Retinal nerve fiber layer retardation measurements using a polarization-sensitive fundus camera

Contents. Acknowledgments. iii. 1 Structure and Function 1. 2 Optics of the Human Eye 3. 3 Visual Disorders and Major Eye Diseases 5

doi: /OE

Lecture 25 Optical Coherence Tomography

Monte Carlo simulation of an optical coherence tomography signal in homogeneous turbid media

Transmission- and side-detection configurations in ultrasound-modulated optical tomography of thick biological tissues

GRENOUILLE.

Fourier Domain (Spectral) OCT OCT: HISTORY. Could OCT be a Game Maker OCT in Optometric Practice: A THE TECHNOLOGY BEHIND OCT

FIRST REPORTED in the field of fiber optics [1], [2],

UC Davis UC Davis Previously Published Works

OCT mini-symposium. Presenters. Donald Miller, Indiana Univ. Joseph Izatt, Duke Univ. Thomas Milner, Univ. of Texas at Austin Jay Wei, Zeiss Meditec

Glaucoma Advanced, LAbel-free High resolution Automated OCT Diagnostics GALAHAD

Radial Polarization Converter With LC Driver USER MANUAL

GPU accelerated real-time multi-functional spectral-domain optical coherence tomography system at 1300nm

Numerical analysis of gradient index lens based optical coherence tomography imaging probes

OCT Spectrometer Design Understanding roll-off to achieve the clearest images

Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second

1272. Phase-controlled vibrational laser percussion drilling

Supplementary Materials

Talbot bands in the theory and practice of optical coherence tomography

A broadband achromatic metalens for focusing and imaging in the visible

Optical design of a dynamic focus catheter for high-resolution endoscopic optical coherence tomography

Analysis of posterior retinal layers in spectral optical coherence tomography images of the normal retina and retinal pathologies

Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, 77. Table of Contents 1

Micromachining with tailored Nanosecond Pulses

All fiber optics circular-state swept source polarization-sensitive optical coherence tomography

Use of Computer Generated Holograms for Testing Aspheric Optics

Ultra High Speed Space Division Multiplexing OCT

Characterization of a fibre optic swept laser source at 1!m for optical coherence tomography imaging systems

High-Coherence Wavelength Swept Light Source

101 W of average green beam from diode-side-pumped Nd:YAG/LBO-based system in a relay imaged cavity

Performance Comparison of Spectrometers Featuring On-Axis and Off-Axis Grating Rotation

Bias errors in PIV: the pixel locking effect revisited.

Improvement of terahertz imaging with a dynamic subtraction technique

Multi-functional angiographic OFDI using frequency-multiplexed dual-beam illumination

Supplementary Figure 1. GO thin film thickness characterization. The thickness of the prepared GO thin

EE119 Introduction to Optical Engineering Spring 2003 Final Exam. Name:

Improving efficiency of CO 2

Spatially Resolved Backscatter Ceilometer

Extended coherence length megahertz FDML and its application for anterior segment imaging

7 CHAPTER 7: REFRACTIVE INDEX MEASUREMENTS WITH COMMON PATH PHASE SENSITIVE FDOCT SETUP

Physics 431 Final Exam Examples (3:00-5:00 pm 12/16/2009) TIME ALLOTTED: 120 MINUTES Name: Signature:

Development of a new multi-wavelength confocal surface profilometer for in-situ automatic optical inspection (AOI)

CHAPTER 5 FINE-TUNING OF AN ECDL WITH AN INTRACAVITY LIQUID CRYSTAL ELEMENT

Extremely simple device for measuring 1.5-µm ultrashort laser pulses

Transferring wavefront measurements to ablation profiles. Michael Mrochen PhD Swiss Federal Institut of Technology, Zurich IROC Zurich

Digital Wavefront Sensors Measure Aberrations in Eyes

Fast Lateral Line-scan Model with Doublet-cylinder-lens in Intelligent OCT

Off-axis full-field swept-source optical coherence tomography using holographic refocusing

Transcription:

Pulsed illumination spectral-domain optical coherence tomography for human retinal imaging Jang-Woo You 1, 1) Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology 373-1, Kusong-dong, Yusong-gu, Daejeon, 305-701, Republic of Korea Teresa C. Chen Massachusetts Eye and Ear Infirmary and Harvard Medical School, 43 Charles Street, Boston, Massachusetts 0114 Mircea Mujat, B. Hyle Park, and Johannes F. de Boer ) Wellman Center for Photomedicine, Massachusetts General Hospital, Harvard Medical School deboer@helix.mgh.harvard.edu Abstract: We present pulsed illumination spectral-domain optical coherence tomography (SD-OCT) for in vivo human retinal imaging. We analyze the signal-to-noise (SNR) for continuous wave (CW) and pulsed illumination SD-OCT. The lateral beam scan motion is responsible for a SNR drop due to lateral scanning induced interference fringe washout. Pulsed illumination can reduce the SNR drop by shorter sample illumination time during the integration time of a camera. First, we demonstrate the SNR benefit of pulsed illumination over CW as function of lateral scan speed for a paper sample. Finally, we show better SNR in retinal images of a normal subject with pulsed illumination SD-OCT over CW at high lateral scanning speed. 006 Optical Society of America OCIS codes: (170.3880) Medical and Biological Imaging; (170.4470) Ophthalmology; (170.4500) Optical Coherence Tomography; References and links 1. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, Optical Coherence Tomography, Science 54, 1178-1181, (1991). G. J. Tearney, M. E. Brezinski, B. E. Bouma, S. A. Boppart, C. Pitris, J. F. Southern, and J. G. Fujimoto, In vivo endoscopic optical biopsy with optical coherence tomography, Science 76, 037-039, (1997) 3. H. Ozdemir, S. A. Karacorlu, and M. Karacorlu, Early Optical Coherence Tomography Changes After Photodynamic Therapy in Patients With Age-Related Macular Degeneration, Am. J. Ophthalmol. 141, 574-576, (006) 4. A. Salinas-Alaman, A. Garcia-Layana, M. J. Maldonado, C. Sainz-Gomez, and A. Alvarez-Vidal, Using optical coherence tomography to monitor photodynamic therapy in age related macular degeneration, Am. J. Ophthalmol. 140, 3-8, (005) 5. W. Goebel and R. Franke, Retinal thickness in diabetic retinopathy - Comparison of optical coherence tomography, the retinal thickness analyzer, and fundus photography, Retina, 6, 49-57, (006) 6. A. Polito, M. Del Borrello, G. Polini, F. Furlan, M. Isola, and F. Bandello, Diurnal variation in clinically significant diabetic macular edema measured by the stratus OCT, Retina, 6, 14-0, (006) 7. A. J. Witkin, T. H. Ko, J. G. Fujimoto, J. S. Schuman, C. R. Baumal, A. H. Rogers, E. Reichel, and J. S. Duker, Redefining lamellar holes and the vitreomacular interface: An ultrahigh-resolution optical coherence tomography study, Ophthalmol. 113, 388-397, (006) 8. M. L. Subramanian, S. N. Truong, A. H. Rogers, J. S. Duker, E. Reichel, and C. R. Baumal, Vitrectomy for stage 1 macular holes identified by optical coherence tomography, Ophthalmic Surg. Lasers 37, 4-46, (006) (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6739

9. C. L. Shields, M. A. Materin, and J. A. Shields, Review of optical coherence tomography for intraocular tumors, Curr. Opin. Ophthalmol. 16, 141-154, (005) 10. B. Cense, T. C. Chen, B. H. Park, M. C. Pierce, and J. F. de Boer, Thickness and birefringence of healthy retinal nerve fiber layer tissue measured with polarization sensitive optical coherence tomography, Invest. Ophthalmol. Visual Sci. 45, 606-61, (004) 11. G. Wollstein, J. S. Schuman, L. L. Price, A. Aydin, P. C. Stark, E. Hertzmark, E. Lai, H. Ishikawa, C. Mattox, J. G. Fujimoto, and L. A. Paunescu, Optical coherence tomography longitudinal evaluation of retinal nerve fiber layer thickness in glaucoma, Arch. Ophthalmol. 13, 464-470, (005) 1. F. A. Medeiros, L. M. Zangwill, C. Bowd, R. M. Vessani, R. Susanna, and R. N. Weinreb, Evaluation of retinal nerve fiber layer, optic nerve head, and macular thickness measurements for glaucoma detection using optical coherence tomography, Am. J. Ophthalmol. 139, 44-55, (005) 13. J. M. J. Roodhooft, Leading causes of blindness worldwide, Bull. Soc. belge. Ophtalmol. 83, 19-5, (00) 14. A. F. Fercher, C. K. Hitzenberger, G. Kamp, and S. Y. Elzaiat, Measurement of Intraocular Distances by Backscattering Spectral Interferometry, Opt. Commun. 117, 43-48, (1995) 15. S. R. Chinn, E. A. Swanson, and J. G. Fujimoto, Optical coherence tomography using a frequency-tunable optical source, Opt. Lett., 340-34, (1997) 16. S. H. Yun, G. J. Tearney, J. F. de Boer, N. Iftimia, and B. E. Bouma, High-speed optical frequencydomain imaging, Opt. Express 11, 953-963, (003). 17. R. Leitgeb, C. K. Hitzenberger, and A. F. Fercher, Performance of fourier domain vs. time domain optical coherence tomography, Opt. Express 11, 889-894, (003). 18. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, Improved signal-tonoise ratio in spectral-domain compared with time-domain optical coherence tomography, Opt. Lett. 8, 067-069, (003) 19. M. A. Choma, M. V. Sarunic, C. H. Yang, and J. A. Izatt, Sensitivity advantage of swept source and Fourier domain optical coherence tomography, Opt. Express 11, 183-189, (003). 0. N. Nassif, B. Cense, B. H. Park, S. H. Yun, T. C. Chen, B. E. Bouma, G. J. Tearney, and J. F. de Boer, In vivo human retinal imaging by ultrahigh-speed spectral domain optical coherence tomography, Opt. Lett. 9, 480-48, (004) 1. W. Y. Oh, S. H. Yun, B. J. Vakoc, G. J. Tearney, and B. E. Bouma, Ultrahigh-speed optical frequency domain imaging and application to laser ablation monitoring, Appl. Phys. Lett. 88, (006). B. Cense, N. Nassif, T. C. Chen, M. C. Pierce, S. H. Yun, B. H. Park, B. E. Bouma, G. J. Tearney, and J. F. de Boer, Ultrahigh-Resolution High-Speed Retinal Imaging Using Spectral-Domain Optical Coherence Tomography, Opt. Express 1, 435-447, (004). 3. M. Wojtkowski, V. J. Srinivasan, T. H. Ko, J. G. Fujimoto, A. Kowalczyk, and J. S. Duker, Ultrahighresolution, high-speed, Fourier domain optical coherence tomography and methods for dispersion compensation, Opt. Express 1, 404-4, (004). 4. T. C. Chen, B. Cense, M. C. Pierce, N. Nassif, B. H. Park, S. H. Yun, B. R. White, B. E. Bouma, G. J. Tearney, and J. F. de Boer, Spectral domain optical coherence tomography - Ultra-high speed, ultra-high resolution ophthalmic imaging, Arch. Ophthalmol. 13, 1715-170, (005) 5. S. H. Yun, G. J. Tearney, J. F. de Boer, and B. E. Bouma, Motion artifacts in optical coherence tomography with frequency-domain ranging, Opt. Express 1, 977-998, (004). 6. S. H. Yun, G. J. Tearney, J. F. de Boer, and B. E. Bouma, Pulsed-source and swept-source spectraldomain optical coherence tomography with reduced motion artifacts, Opt. Express 1, 5614-564, (004). 7. G. Moneron, A. C. Boccara, and A. Dubois, Stroboscopic ultrahigh-resolution full-field optical coherence tomography, Opt. Lett. 30, 1351-1353, (005) 8. N. Nassif, B. Cense, B. H. Park, M. C. Pierce, S. H. Yun, B. E. Bouma, G. J. Tearney, T. C. Chen, and J. F. de Boer, In vivo high-resolution video-rate spectral-domain optical coherence tomography of the human retina and optic nerve, Opt. Express, 1 367-376, (004). 9. American National Standards Institute, American National Standard for Safe Use of Lasers Z136.1. 000: Orlando. 30. M. Mujat, R. C. Chan, B. Cense, B. H. Park, C. Joo, T. Akkin, T. C. Chen, and J. F. de Boer, Retinal nerve fiber layer thickness map determined from optical coherence tomography images, Opt. Express 13, 9480-9491, (005). 31. S. L. Jiao, R. Knighton, X. R. Huang, G. Gregori, and C. A. Puliafito, Simultaneous acquisition of sectional and fundus ophthalmic images with spectral-domain optical coherence tomography, Opt. Express 13, 444-45, (005), 3. S. H. Yun, G. J. Tearney, B. E. Bouma, B. H. Park, and J. F. de Boer, High-speed spectral-domain optical coherence tomography at 1.3 mu m wavelength, Opt. Express 11, 3598-3604, (003). 33. I. Escudero-Sanz and R. Navarro, Off-axis aberrations of a wide-angle schematic eye model, J. Opt. Soc. Am. A-Opt. Image Sci. Vis. 16, 1881-1891, (1999) 34. D. A. Goss and R. W. West, Introduction to the Optics of the Eye, (Butterworth-Heinemann, 00), Chap. 3. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6740

1. Introduction OCT has emerged as a promising technology [1, ] for in-vivo imaging in a wide range of medical applications. Current ophthalmic applications of OCT have included diseases such as age-related macular degeneration, [3, 4] diabetic retinopathy, [5, 6] macular holes, [7, 8] intraocular tumors, [9] and glaucoma. [10-1] In particular, glaucoma is the third leading cause of worldwide blindness, after cataract and trachoma. [13] The evolution from timedomain [1] to spectral-domain [14] and optical frequency domain imaging (OFDI), [15, 16] has made simultaneous improvement in sensitivity [17-19], high speed, [0, 1] and axial resolution feasible in many applications, such as retinal imaging. [-4] As OCT utilizes lateral point-scanning, motion of the sample or scanning beam during the measurement causes SNR reduction and image degradation in SD-OCT and OFDI. [5] Yun et al. theoretically investigated axial and lateral motion artifacts in CW SD-OCT and sweptsource OFDI, and experimentally demonstrated reduced axial and lateral motion artifacts using a pulsed source and a swept source in endoscopic imaging of biological tissue. [5, 6] Stroboscopic illumination in full field OCT was demonstrated, resulting in reduced motion artifacts for in-vivo measurement. [7] In ophthalmic applications of SD-OCT, SNR reduction caused by high speed lateral scanning of the beam over the retina may be dominant over axial patient motion. Using pulsed illumination can reduce lateral motion artifacts. We analyzed the SNR benefit of pulsed-illumination over CW SD-OCT, demonstrating that pulsed illumination provides a better SNR for in-vivo high speed human retinal imaging. The structure of the paper is as follows: the SNR benefit of pulsed-source over CW SD- OCT theoretically analyzed for different lateral scanning speeds. Then, a paper sample was measured to show the validity of the theory. Finally, the SNR for in-vivo human retinal imaging with pulsed-illumination and CW SD-OCT was analyzed and compared with the theory.. Benefits of pulsed illumination SD-OCT.1. Theory Yun et al. [5] derived the equations for the SNR decrease due to lateral motion in CW SD- OCT. For a CW source, the SNR decrease is given by, 5log 10 1 0.5, Δx SNR decrease + (1) w o where x is the scanning distance during the camera integration time and w o denotes fullwidth-half-maximum (FWHM) of the beam profile. The normalized displacement is defined as x/w o. For pulsed illumination, x is replaced by (τpulse/τcamera) x, where τpulse is the pulse width in time and τcamera is the integration time of the camera for a single A-line. As can be seen in Fig. 1, pulsed-illumination SD-OCT has a significant SNR benefit over CW SD-OCT when the displacement x is larger than the FWHM w o of the beam profile under the same average power. For example, at a normalized displacement of 4, a 0% duty cycle pulsed illumination has a 4. db better SNR than CW illumination. Figure 1 also shows that for CW illumination, a normalized displacement of 1 gives 1 db decrease in SNR. At a source center wavelength of 841 nm, an axial motion of 50 nm per A- line would also result in a 1 db SNR reduction. [5] At an A-line rate of 9.3 khz, this corresponds to an axial velocity of about 1.5 mm/s, which for retinal imaging is, in general, much higher than encountered in practice. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6741

1 0-1 SNR decrease (db) - -3-4 -5-6 -7-8 CW 10 % 0 % 50 % -9 0 4 6 8 Normalized displacement Fig. 1. Theoretical SNR decrease due to lateral motion for a CW source and for pulsed illumination with different duty cycles... Experimental setup Our measurement system has been previously described by Nassif et.al. [8] For the measurements presented here, the light source was a high power SLD (SLD-371-HP3, Superlum, Rusia) with a spectral bandwidth of 48.8 nm FWHM at a center wavelength of λo=843.1 nm. The current driver to the source could be operated in CW, and pulsed mode triggered by a signal from the computer. The maximum output power of the SLD was 3 mw and the measured axial resolution was 7.9 µm in air. The maximum average output power of the SLD in pulsed mode was determined by the pulse width and the pulse repetition rate. The pulse repetition rate was the same as the spectrometer triggering frequency, i.e. 9.3 khz. The pulse width was set to 7 µs (~0% duty cycle), giving 400 µw in the sample arm after the slit lamp. For CW operation, a neutral density filter was inserted in the source arm to adjust the average power to the same value as in pulsed operation, such that the SNR decrease for pulsed and CW illumination could be compared directly. The collimated beam exiting the slit lamp was focused by a 30 mm focal length biconvex lens on a paper sample placed in the focal plane of the lens as illustrated in Fig.. The FWHM beam diameter w o was calculated as 8. µm, and the lateral scanning range was 16 mm as measured on the sample. To control the normalized displacement, the lateral scan range remained unchanged but the number of A-lines per one complete lateral scan was decreased, resulting in increased lateral scan velocity. The scan range, the number of A-lines and the spot size were used to calculate the normalized displacement x/w o. The paper sample was measured in two different configurations: parallel to the focal plane, and at a 4.1 tilt as illustrated at the top of Fig. 3. The number of A-lines per image was reduced from 1000 to 300 in 100 A-line decrements. The SNR was measured in each image with CW and pulsed illumination. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 674

lens Paper sample Fig.. Slit lamp integrated with a biconvex lens and a paper sample.3. Results and discussion For pulsed source illumination, the pulse width was set to 7 µs (0% duty cycle), which is equivalent to changing x to 0. x in Eq. (1). The average power incident on the paper sample was kept constant for CW and pulsed illumination. The SNR was calculated from the average of the maximum signal in each A-line near the center of each image. 80 % and 0 % of the A-lines near the center in each image were considered to calculate SNR for horizontal and slanted flat paper, respectively. The SNR from experiment and theory for pulsed and CW illumination are presented in Fig. 3. 80 % 0 % 1000 A-lines 500 A-lines 0 0-1 -1 SNR decrease (db) - -3-4 -5 pulsed experiment pulsed theory cw experiment cw theory SNR decrease (db) - -3-4 -5 pulsed experiment pulsed theory cw experiment cw theory -6-6 -7-7 3 4 5 6 7 3 4 5 6 7 Normalized displacement Normalized displacement Fig. 3. SNR decrease for the flat paper sample parallel to the focal plane (left) and with 4.1 tilt (right) for CW and pulsed illumination compared with theory. 1000 and 500 A-lines are equivalent to.3 and 4.3 normalized displacements, respectively. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6743

The best fit for w o was 7 µm, in good agreement with the calculated value of 8. µm. The data was normalized to give 0 db SNR at the smallest normalized displacement of the pulsed illumination. Figure 3 shows a significantly reduced SNR decrease for pulsed illumination. The equal SNR decrease between horizontal and slanted surfaces shows that a stationary slanted surface does not generate a further SNR decrease, despite that at a 4.1 slant, a scan length of 16 mm and 1000 A-lines per scan, the path length to the sample surface increases by 1.1 µm per A-line. For retinal imaging, this implies that the retinal curvature in an image is not expected to contribute to a SNR decrease. The experimental data agrees very well with the theory and demonstrates that the SNR for pulsed illumination is higher than for CW. 3. Application to human retinal imaging 3.1. Experimental setup We applied, for the first time to our knowledge, pulsed illumination to human retinal imaging. We analyzed the SNR by changing the normalized displacement, and we compared SNR with CW illumination. The pulse width was 8 µs and the pulse repetition rate was 9.3 KHz, synchronized with the integration time of the high-speed line scan camera. The voltage input for horizontal scan was the same as in the paper sample experiment, ±.5 V, however the actual scan range was 8.6 mm on the retina, which is about half the scan range on the paper sample. This is due to the higher refractive power of the eye as compared to that of the 30 mm bi-convex lens, which also will result in a smaller spot size on the retina. The normalized displacement was adjusted based on the actual scan range on the retina. According to the ANSI standard for safe use of lasers, the maximum permissible exposure expressed as the average power of a pulse train increases with the pulse repetition rate up to a frequency of 55 khz, where the limit to continuous wave exposure is reached. The MPE sp for a single pulse between 1 ns to 18 µs is given by, MPE sp 7 = 5.0 C 10 [ J / cm ] () a with C A =1.66 for 810 nm. The MPE of a pulse in the pulse train is given by the single pulse MPE, multiplied by a correction factor determined by the number of pulses n during a time interval of at most 10 seconds for small source ocular exposure. The maximum average power into the eye is then given by this MPE of a pulse in the pulse train, multiplied by the repetition frequency f and the limiting aperture of 0.385 cm. rep MPE = 0.5 ( n ) = MPE sp MPE 10 sp f 0.5 rep f 0.385cm 0.75 rep 0.385 cm [ μw ] Figure 4 shows the maximum permissible average power in the pulse train as a function of repetition frequency, and the maximum permissible power of CW illumination (which can not be exceeded). Although the maximum average power increases up to a frequency of 55 khz, the energy in a single pulse, which ultimately determines the SNR or sensitivity of a single A- line, decreases. The optimal sensitivity would be achieved for a single pulse every 10 sec, which is incompatible with high-speed imaging. The maximum permissible exposure (MPE) for human retinal imaging at a wavelength of 810 nm is 43.5 µw for pulsed illumination with 9.3 khz repetition rate and 8 µs pulse length, and 640 µw for CW illumination. [9] (3) (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6744

800 700 CW illumination 600 500 MPE (μw) 400 300 00 Pulsed illumination 100 0 0 10 0 30 40 50 60 70 pulse repetition frequency (khz) Fig. 4. Maximum permissible exposure expressed as average power into the eye, according to the ANSI standards for CW and pulsed illumination. Red line: CW power limitation. Black line: Pulsed source power limitation calculated from Eq. (3) as a function of pulse repetition frequency. The smaller MPE of CW and pulsed illumination at specific frequency determines the power limitation to the retina. B C A (a) (b) D (c) (d) Fig. 5. Integrated reflectance images of in-vivo human retina centered on the optic nerve head for two kinds of illumination (CW and pulsed) at two different lateral scanning speeds. (a) CW, 1000 A-lines (b) pulsed, 1000 A-lines (c) CW, 500 A-lines (d) pulsed, 500 A-lines. Four regions of interest are indicated (zones A to D), where the SNR was analyzed in each image. Some discontinuities in (a) and (b) were caused by eye motion during the two times longer measurement time as compared to (c) and (d). The lateral scan angle of 30 in the experiment results in 8.6 mm lateral scan range. Vertical scan range is 5. mm on the retina. The sample arm power after the slit lamp of CW and pulsed illumination was set to 600 µw and 385 µw respectively. The 385 µw average power for pulsed illumination will cause a 1.9 db smaller SNR compared to the 600 µw for CW. This has been accounted for in the SNR analysis comparing pulsed and CW illumination. As shown in Fig. 5, integrated reflectance images of the retina [30, 31] of a healthy volunteer at two different lateral scanning speeds characterized by 1000 and 500 A-lines per image for CW and pulsed llumination were measured. We can see some discontinuities in the 1000 A-line images because of the two times longer measurement time relative to the 500 A-line images. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6745

3.. Results and discussion The SNR in each measurement was calculated and analyzed as follows. First, we selected 4 regions of interest based on the integrated reflectance image as shown in Fig. 5 to compare the SNR in the same region of the retina for CW and pulsed illumination. Zone A is comprised of 5 frames and the other zones contain 9 frames. Even if we selected the same regions of the retina based on Fig. 5, the depth location of the retina in each frame was slightly different for the pulsed and CW case. The sensitivity of SD-OCT decreases with depth [17, 3], and the normalized sensitivity decay was measured in our previous work. [8] Each original A-line was divided by the normalized sensitivity of the system as a weighting factor to correct for the depth dependent sensitivity. Since the depth dependent sensitivity has the smallest roll-off close to the zero-delay position of the reference arm, the regions of interest were selected to the left side of the optical nerve head (ONH) because they are closest to the zero-delay position of the reference arm, as shown in Fig. 6. (a) (b) (c) (d) Fig. 6. ONH images from zone C in Fig 4(a) ~ 4(d) at the same location, gray scale coded over the same dynamic range. The ONH is in axially different locations. SNR was compared after correcting for the depth dependent sensitivity decrease. The horizontal dimension is 8.6 mm and vertical dimension is 1.7 mm. The highest signal in each A-line was detected and averaged over the region of interest (ROI) in each image and the SNR of each image was averaged again in each zone. Table 1 shows the SNR of each zone for different experimental conditions. Table 1. SNR calculated in each zone of the in-vivo human retina data shown in Fig. 6 after correcting for the depth dependent sensitivity decrease. SNR (db) CW, 1000 A-line pulsed, 1000 A-line CW, 500 A-line pulsed, 500 A-line zone A 9.0 8.7 5.1 8.1 zone B 9.8 9.4 6.8 8.9 zone C 31.3 30.6 8.9 30.3 zone D 31.1 30.5 8.0 30. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6746

The data was fitted to the theory in the following way: since the reflectivity of a particular location on the retina varies, SNR values in the 4 zones could not be compared directly. However, SNR values for pulsed and CW illumination within the same zone could be compared directly. The same offset was subtracted from the SNR values within a zone. The offset varied by zone, and was determined by a fit to the theory. The beam diameter w o in each zone was also determined by a fit to the theory, giving 3.8, 4.6, 5.7, and 4.7 µm for zones A to D, respectively. The smallest RMS spot size on the retina layer can be simulated with Zemax as ~ 3 µm with Navarro s model [33] and ~5 µm with Gullstrand Schematic Eye No.1 [34]. The normalized displacement is slightly different in each zone due to the variation of w o. The SNR values in each zone and the theory are shown in Fig. 7. 1 0 1000 Alines/image 500 Alines/image cw pulsed Relative SNR (db) -1 - -3-4 A B C D -5 0 1 3 4 5 6 7 Normalized displacement Fig. 7. Relative SNR fitted with w o and offset as parameters. The theoretical fit for CW illumination (continuous line) starts 1.9 db above the theoretical fit for pulsed illumination (dotted line) because of the higher average CW power. Despite the lower average power for pulsed illumination, the SNR is better for normalized displacements larger than. The 1.9 db better relative SNR for CW illumination over pulsed illumination at zero normalized displacement in Fig. 7 reflects the higher average power permitted under the ANSI standards. We can see that the SNR was almost the same at 1000 A-lines per image with CW and pulsed illumination but it was 1.4 to 3.0 db higher with pulsed illumination for larger normalized displacement. The variation in beam diameter is attributed to a difference in curvature between the retina and the focal plane of the imaging system, where the best overlap between retina and focal plane was realized in zone A, good overlap was realized in zone B and D, and the worst overlap was realized in zone C. One consideration that should be taken into account is that shortening the pulse length while maintaining average power and pulse repetition rate will decrease the lateral motion artifact but also increases the ratio of RIN over shot noise. For a shot noise limited system, this ratio needs to be smaller than 1. Using Eq. () in ref [0] and replacing Pref with J p τ i, with J p the energy of the pulse and τ i the pulse duration, gives, σ noise ηe J p ηej p = σ r + d + + τ [ e iτ coh ] E ν Eντ (4) i (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6747

The reference arm pulse energy J p to saturate the spectrometer to 90% of the full well capacity is independent of pulse length. The pulse energy for which the shot noise still Evτ i dominates over RIN noise is given by J p <, which shows that reducing the pulse ητ coh length also reduces this pulse energy. In our case of an 8 µs pulse length, the pulse energy is still more than a factor of two below the value where RIN noise is equal to shot noise. [0] 4. Conclusions We examined the SNR decrease due to lateral motion for CW and pulsed illumination and demonstrated good correlation between experimental data and theory for two different samples: a) paper, and b) in-vivo human retina. With equal incident power, depending on the normalized displacement of the beam across a paper sample, a ~3 to ~6 db better SNR was found for pulsed over CW illumination. In ophthalmic applications, lateral motion from high speed beam scanning can be regarded as a dominant factor contributing to SNR decrease. As mentioned in section 3.1, according to the ANSI standards, the MPE for pulsed illumination in our setup is 43.5 µw, as compared to 640 µw for CW. Therefore, based on the MPE, CW should have a 1.8 db better SNR over pulsed illumination. However, due to the SNR decrease resulting from fast lateral scanning, the SNR was almost the same between CW and pulsed illumination with 1000 A-lines per image over an 8. mm scan range. The measured SNR for pulsed illumination was 1.4 to 3.0 db higher than for the CW source at a two times larger normalized displacement and agreed well with the theory. In conclusion, with CW illumination under our experimental conditions we have demonstrated that doubling the scan speed or normalized displacement resulted in up to.5 db of SNR reduction. With pulsed illumination the lateral scan speed could be doubled at a cost of up to 0.5 db SNR reduction. Pulsed illumination permits higher scan speeds and reduced scan times for large area retinal scans at a much smaller SNR penalty than CW illumination. Acknowledgments This research was supported by research grants from the National Institute of Health (R01 EY 14975 and R01 RR 19768) and the Department of Defense (F4 960-01-1-0014). JWY would like to thank Dr. R. Navarro for his advice on modeling the human eye for numerical simulation with Zemax. (C) 006 OSA 4 July 006 / Vol. 14, No. 15 / OPTICS EXPRESS 6748