OSA's Digital Library

Biomedical Optics Express

Biomedical Optics Express

  • Editor: Joseph A. Izatt
  • Vol. 4, Iss. 10 — Oct. 1, 2013
  • pp: 1890–1908
« Show journal navigation

Multi-MHz retinal OCT

Thomas Klein, Wolfgang Wieser, Lukas Reznicek, Aljoscha Neubauer, Anselm Kampik, and Robert Huber  »View Author Affiliations


Biomedical Optics Express, Vol. 4, Issue 10, pp. 1890-1908 (2013)
http://dx.doi.org/10.1364/BOE.4.001890


View Full Text Article

Acrobat PDF (2478 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We analyze the benefits and problems of in vivo optical coherence tomography (OCT) imaging of the human retina at A-scan rates in excess of 1 MHz, using a 1050 nm Fourier-domain mode-locked (FDML) laser. Different scanning strategies enabled by MHz OCT line rates are investigated, and a simple multi-volume data processing approach is presented. In-vivo OCT of the human ocular fundus is performed at different axial scan rates of up to 6.7 MHz. High quality non-mydriatic retinal imaging over an ultra-wide field is achieved by a combination of several key improvements compared to previous setups. For the FDML laser, long coherence lengths and 72 nm wavelength tuning range are achieved using a chirped fiber Bragg grating in a laser cavity at 419.1 kHz fundamental tuning rate. Very large data sets can be acquired with sustained data transfer from the data acquisition card to host computer memory, enabling high-quality averaging of many frames and of multiple aligned data sets. Three imaging modes are investigated: Alignment and averaging of 24 data sets at 1.68 MHz axial line rate, ultra-dense transverse sampling at 3.35 MHz line rate, and dual-beam imaging with two laser spots on the retina at an effective line rate of 6.7 MHz.

© 2013 OSA

1. Introduction

1.1. Importance of acquisition speed in retinal OCT

Optical coherence tomography is a noninvasive, three-dimensional imaging modality [1

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 254(5035), 1178–1181 (1991). [CrossRef] [PubMed]

]. OCT has found widespread use in many applications, especially in ophthalmology for cross-sectional investigation of the human retina [2

2. M. R. Hee, J. A. Izatt, E. A. Swanson, D. Huang, J. S. Schuman, C. P. Lin, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography of the human retina,” Arch. Ophthalmol. 113(3), 325–332 (1995). [CrossRef] [PubMed]

]. Over the last years, OCT imaging speed has been increased drastically by the introduction of Fourier-domain OCT (FD-OCT) [3

3. M. Wojtkowski, T. Bajraszewski, I. Gorczyńska, P. Targowski, A. Kowalczyk, W. Wasilewski, and C. Radzewicz, “Ophthalmic imaging by spectral optical coherence tomography,” Am. J. Ophthalmol. 138(3), 412–419 (2004). [CrossRef] [PubMed]

]. Nowadays, the fastest research systems operate at axial line rates of several MHz [4

4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

]. In contrast, commercial ophthalmic systems usually operate at tens of kHz axial line rate, about two orders of magnitude slower. Since sample motion restricts the acquisition time for in-vivo imaging, imaging speed effectively determines the total number of A-scans in each data set. Hence, either sampling density or field of view is restricted. For retinal imaging, snapshot-like acquisition times of well below 1 s are necessary for a high probability that the 3D data sets do not contain microsaccades [5

5. B. Povazay, B. Hofer, C. Torti, B. Hermann, A. R. Tumlinson, M. Esmaeelpour, C. A. Egan, A. C. Bird, and W. Drexler, “Impact of enhanced resolution, speed and penetration on three-dimensional retinal optical coherence tomography,” Opt. Express 17(5), 4134–4150 (2009). [CrossRef] [PubMed]

]. Therefore, commercial systems typically have a variety of acquisition modes in order to accomplish different goals. In “hi-res” modes, sampling density is increased and/or multiple frames are averaged in order to produce high-quality single B-frames. Frame averaging is especially important to reduce speckle noise, which obstructs fine anatomical detail [6

6. J. M. Schmitt, S. H. Xiang, and K. M. Yung, “Speckle in Optical Coherence Tomography,” J. Biomed. Opt. 4(1), 95–105 (1999). [CrossRef] [PubMed]

, 7

7. M. Wojtkowski, B. Kaluzny, and R. J. Zawadzki, “New Directions in Ophthalmic Optical Coherence Tomography,” Optom. Vis. Sci. 89(5), 524–542 (2012). [CrossRef] [PubMed]

]. In contrast, sampling density is typically reduced for acquisition of 3D data sets, which provide important diagnostic data, for instance thickness maps. It is usually not possible to have both the highest quality B-frames and high-resolution 3D capability in a single data set. This is why methods based on the field of compressive sensing are being investigated currently [8

8. L. Fang, S. Li, Q. Nie, J. A. Izatt, C. A. Toth, and S. Farsiu, “Sparsity based denoising of spectral domain optical coherence tomography images,” Biomed. Opt. Express 3(5), 927–942 (2012). [CrossRef] [PubMed]

]. Due to the small field of view, 3D data sets are usually restricted to either the macula or the optical nerve head (ONH) area, and the operator needs to decide on the scanned region prior to data acquisition. Stitching or “mosaicing” of multiple data sets can be applied to increase the effective field of view, but this techniques increase acquisition time considerably [9

9. R. J. Zawadzki, S. S. Choi, A. R. Fuller, J. W. Evans, B. Hamann, and J. S. Werner, “Cellular resolution volumetric in vivo retinal imaging with adaptive optics-optical coherence tomography,” Opt. Express 17(5), 4084–4094 (2009). [CrossRef] [PubMed]

11

11. S. Makita, F. Jaillon, M. Yamanari, M. Miura, and Y. Yasuno, “Comprehensive in vivo micro-vascular imaging of the human eye by dual-beam-scan Doppler optical coherence angiography,” Opt. Express 19(2), 1271–1283 (2011). [CrossRef] [PubMed]

]. Moreover, all the volumes need to be acquired before fusing them, which results in significant delay between capturing images and displaying large field of view.

For these reasons, many systems employ an eye tracker, which corrects for involuntary eye motion and thus allows for an extended measurement time [12

12. R. D. Ferguson, D. X. Hammer, L. A. Paunescu, S. Beaton, and J. S. Schuman, “Tracking optical coherence tomography,” Opt. Lett. 29(18), 2139–2141 (2004). [CrossRef] [PubMed]

]. Eye trackers make the system more complex, since they need to incorporate an additional imaging mechanism with a frame rate which is higher than the OCT volume acquisition rate. Typically a scanning laser ophthalmoscope (SLO) with a separate beam steering mechanism is used, which makes it difficult to ensure pixel-to-pixel correlation between OCT and SLO data over the entire imaging volume. This issue may be overcome if most of the optical path is used for both the SLO and the OCT channel [13

13. Y. K. Tao, S. Farsiu, and J. A. Izatt, “Interlaced spectrally encoded confocal scanning laser ophthalmoscopy and spectral domain optical coherence tomography,” Biomed. Opt. Express 1(2), 431–440 (2010). [CrossRef] [PubMed]

]. Moreover, since tracking accuracy is limited, it is usually not possible to generate high-resolution en-face images directly from the tracked OCT data set [14

14. G. Gregori and P. J. Rosenfeld, “Using OCT Fundus Images to Evaluate the Performance of the Spectralis OCT Eye Tracking System (Poster 1032/A434),” in ARVO Annual Meeting 2011, (Fort Lauderdale, Florida, 2011).

], although some research systems have shown higher tracking accuracy [15

15. K. V. Vienola, B. Braaf, C. K. Sheehy, Q. Yang, P. Tiruveedhula, D. W. Arathorn, J. F. de Boer, and A. Roorda, “Real-time eye motion compensation for OCT imaging with tracking SLO,” Biomed. Opt. Express 3(11), 2950–2963 (2012). [CrossRef] [PubMed]

]. Hence usually only the SLO image is used for a transverse view of the retina, not the reconstructed OCT en-face image. However, SLO images represent only the total integrated intensity from all retinal layers, whereas OCT can generate depth-resolved en-face views with very high depth resolution [16

16. A. G. Podoleanu and R. B. Rosen, “Combinations of techniques in imaging the retina with high resolution,” Prog. Retin. Eye Res. 27(4), 464–499 (2008). [CrossRef] [PubMed]

]. These depth-resolved en-face views are a unique feature of OCT [17

17. S. Jiao, R. Knighton, X. Huang, G. Gregori, and C. Puliafito, “Simultaneous acquisition of sectional and fundus ophthalmic images with spectral-domain optical coherence tomography,” Opt. Express 13(2), 444–452 (2005). [CrossRef] [PubMed]

], and facilitate the visualization of many pathologies [18

18. U. Schmidt-Erfurth, R. A. Leitgeb, S. Michels, B. Považay, S. Sacu, B. Hermann, C. Ahlers, H. Sattmann, C. Scholda, A. F. Fercher, and W. Drexler, “Three-Dimensional Ultrahigh-Resolution Optical Coherence Tomography of Macular Diseases,” Invest. Ophthalmol. Vis. Sci. 46(9), 3393–3402 (2005). [CrossRef] [PubMed]

20

20. I. Gorczynska, V. J. Srinivasan, L. N. Vuong, R. W. S. Chen, J. J. Liu, E. Reichel, M. Wojtkowski, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Projection OCT fundus imaging for visualising outer retinal pathology in non-exudative age-related macular degeneration,” Br. J. Ophthalmol. 93(5), 603–609 (2009). [CrossRef] [PubMed]

]. For instance, depth-resolved OCT allows the visualization and three-dimensional characterization of the choroid [21

21. T. Agawa, M. Miura, Y. Ikuno, S. Makita, T. Fabritius, T. Iwasaki, H. Goto, K. Nishida, and Y. Yasuno, “Choroidal thickness measurement in healthy Japanese subjects by three-dimensional high-penetration optical coherence tomography,” Graefes Arch. Clin. Exp. Ophthalmol. 249(10), 1485–1492 (2011). [CrossRef] [PubMed]

, 22

22. B. Povazay, B. Hermann, B. Hofer, V. Kajić, E. Simpson, T. Bridgford, and W. Drexler, “Wide-Field Optical Coherence Tomography of the Choroid In Vivo,” Invest. Ophthalmol. Vis. Sci. 50(4), 1856–1863 (2008). [CrossRef] [PubMed]

]. Due to the limited 3D data set size, these depth-resolved en-face views usually cover a smaller area than a typical SLO scan.

The problem of limited data set size can be overcome with MHz-OCT, i.e. with OCT systems capable of axial line rates in excess of 1 MHz. With MHz-OCT, a single data set can be acquired in short time that spans a very large field of view. These large data sets cover both the macula and ONH area, and enable ultrawide-field depth-resolved fundus views over an area comparable or even larger than covered by typical SLO scans [23

23. T. Klein, W. Wieser, C. M. Eigenwillig, B. R. Biedermann, and R. Huber, “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express 19(4), 3044–3062 (2011). [CrossRef] [PubMed]

]. Moreover, if the OCT imaging speed is increased up into the MHz and multi-MHz range, efficient speckle reduction by frame averaging can be used without loss of wide-field 3D and functional imaging capabilities, as will be shown in this paper. Other particularly promising applications of high-speed 3D OCT are the tomographic analysis of the posterior pole [24

24. A. A. Ellabban, A. Tsujikawa, A. Matsumoto, K. Yamashiro, A. Oishi, S. Ooto, I. Nakata, Y. Akagi-Kurashige, M. Miyake, H. S. Elnahas, T. M. Radwan, K. A. Zaky, and N. Yoshimura, “Three-Dimensional Tomographic Features of Dome-Shaped Macula by Swept-Source Optical Coherence Tomography,” Am. J. Ophthalmol. 155(2), 320–328, e2 (2013). [CrossRef] [PubMed]

], volumetric adaptive optics OCT [9

9. R. J. Zawadzki, S. S. Choi, A. R. Fuller, J. W. Evans, B. Hamann, and J. S. Werner, “Cellular resolution volumetric in vivo retinal imaging with adaptive optics-optical coherence tomography,” Opt. Express 17(5), 4084–4094 (2009). [CrossRef] [PubMed]

] and functional imaging [25

25. S. Makita, Y. Hong, M. Yamanari, T. Yatagai, and Y. Yasuno, “Optical coherence angiography,” Opt. Express 14(17), 7821–7840 (2006). [CrossRef] [PubMed]

29

29. B. White, M. Pierce, N. Nassif, B. Cense, B. Park, G. Tearney, B. Bouma, T. Chen, and J. de Boer, “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography,” Opt. Express 11(25), 3490–3497 (2003). [CrossRef] [PubMed]

]. For instance, the system configuration described in this paper has also been used to generate in-vivo single-shot wide-field microangiography images of the retina for the first time [30

30. C. Blatter, T. Klein, B. Grajciar, T. Schmoll, W. Wieser, R. Andre, R. Huber, and R. A. Leitgeb, “Ultrahigh-speed non-invasive widefield angiography,” J. Biomed. Opt. 17(7), 070505 (2012). [CrossRef] [PubMed]

].

1.2. Pathways to higher speed in OCT

There are two approaches to achieve effective axial scan rates of more than one MHz: Increased single-beam scan-rate or parallelization of illumination and/or data acquisition. Increased single-beam scan rate may seem straightforward, and relies on faster swept lasers or higher spectrometer readout rates. However, for a variety of technical reasons, only a few systems have demonstrated speeds in excess of 200 kHz: The fastest single-camera spectral domain OCT (SD-OCT) system has been demonstrated by Potsaid et. al., running at up to 312.5 kHz line rate with reduced image quality [31

31. B. Potsaid, I. Gorczynska, V. J. Srinivasan, Y. L. Chen, J. Jiang, A. Cable, and J. G. Fujimoto, “Ultrahigh speed Spectral / Fourier domain OCT ophthalmic imaging at 70,000 to 312,500 axial scans per second,” Opt. Express 16(19), 15149–15169 (2008). [CrossRef] [PubMed]

]. More complex system architectures with parallelized data acquisition allow for a further increase in SD-OCT speed: With two spectrometers, axial line rates of 500 kHz have been demonstrated, with a ~3 dB sensitivity penalty compared to a single spectrometer [32

32. L. An, P. Li, T. T. Shen, and R. Wang, “High speed spectral domain optical coherence tomography for retinal imaging at 500,000 A‑lines per second,” Biomed. Opt. Express 2(10), 2770–2783 (2011). [CrossRef] [PubMed]

]. The standard spectrometer design with a single camera can also be replaced with an optical demultiplexer and hundreds of single photodiodes. This constitutes an interesting approach, but currently suffers from very high cost and system complexity [33

33. D. Choi, H. Hiro-Oka, H. Furukawa, R. Yoshimura, M. Nakanishi, K. Shimizu, and K. Ohbayashi, “Fourier domain optical coherence tomography using optical demultiplexers imaging at 60,000,000 lines/s,” Opt. Lett. 33(12), 1318–1320 (2008). [CrossRef] [PubMed]

, 34

34. D. H. Choi, H. Hiro-Oka, K. Shimizu, and K. Ohbayashi, “Spectral domain optical coherence tomography of multi-MHz A-scan rates at 1310 nm range and real-time 4D-display up to 41 volumes/second,” Biomed. Opt. Express 3(12), 3067–3086 (2012). [CrossRef] [PubMed]

].

In swept source OCT (SS-OCT), also known as optical frequency-domain imaging (OFDI) [35

35. S. Yun, G. Tearney, J. de Boer, N. Iftimia, and B. Bouma, “High-speed optical frequency-domain imaging,” Opt. Express 11(22), 2953–2963 (2003). [CrossRef] [PubMed]

], axial line rate is given by the sweep speed of the laser light source. In short cavity lasers, sweep rate is physically limited by the buildup time for laser activity [36

36. R. Huber, M. Wojtkowski, K. Taira, J. G. Fujimoto, and K. Hsu, “Amplified, frequency swept lasers for frequency domain reflectometry and OCT imaging: design and scaling principles,” Opt. Express 13(9), 3513–3528 (2005). [CrossRef] [PubMed]

]. Due to this limit, it seems unlikely that speed can be increased over about 200 kHz to 400 kHz with short cavity lasers [37

37. B. Potsaid, B. Baumann, D. Huang, S. Barry, A. E. Cable, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Ultrahigh speed 1050nm swept source/Fourier domain OCT retinal and anterior segment imaging at 100,000 to 400,000 axial scans per second,” Opt. Express 18(19), 20029–20048 (2010). [CrossRef] [PubMed]

, 38

38. W.-Y. Oh, B. J. Vakoc, M. Shishkov, G. J. Tearney, and B. E. Bouma, “>400 kHz repetition rate wavelength-swept laser and application to high-speed optical frequency domain imaging,” Opt. Lett. 35(17), 2919–2921 (2010). [CrossRef] [PubMed]

]. Another concept, which might achieve higher OCT imaging rates in the future, is based on Vernier-tuned distributed Bragg reflector (VT-DBR) technology. A VT-DBR laser has only been used for OCT imaging at 1550nm with limited axial resolution [39

39. M. P. Minneman, J. Ensher, M. Crawford, and D. Derickson, “All-semiconductor high-speed akinetic swept-source for OCT,” in Optical Sensors and Biophotonics III, Proc. SPIE 8311 (SPIE, 2011), 831116–831116.

]. Currently, the physical limit to sweep speed can be overcome in four completely different ways: (1) Extremely short cavities in vertical-cavity surface-emitting lasers (VCSELs) [40

40. B. Potsaid, V. Jayaraman, J. G. Fujimoto, J. Jiang, P. J. S. Heim, and A. E. Cable, “MEMS tunable VCSEL light source for ultrahigh speed 60kHz - 1MHz axial scan rate and long range centimeter class OCT imaging,” in (SPIE, 2012), 82130M.

], (2) extremely long cavities in Fourier-domain mode-locked (FDML) lasers [41

41. R. Huber, M. Wojtkowski, and J. G. Fujimoto, “Fourier Domain Mode Locking (FDML): A new laser operating regime and applications for optical coherence tomography,” Opt. Express 14(8), 3225–3237 (2006). [CrossRef] [PubMed]

], (3) time-stretching of short pulses outside the cavity [42

42. S. Moon and D. Y. Kim, “Ultra-high-speed optical coherence tomography with a stretched pulse supercontinuum source,” Opt. Express 14(24), 11575–11584 (2006). [CrossRef] [PubMed]

], (4) or no cavity at all in wavelength swept amplified spontaneous emission (ASE) sources [43

43. C. M. Eigenwillig, B. R. Biedermann, W. Wieser, and R. Huber, “Wavelength swept amplified spontaneous emission source,” Opt. Express 17(21), 18794–18807 (2009). [CrossRef] [PubMed]

]:

(1) VCSELs can operate in a single longitudinal mode due to their short cavity length, which is on the size of about a few wavelengths. Hence, very long instantaneous coherence lengths can be achieved [44

44. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, J. Jiang, J. G. Fujimoto, and A. E. Cable, “High-precision, high-accuracy ultralong-range swept-source optical coherence tomography using vertical cavity surface emitting laser light source,” Opt. Lett. 38(5), 673–675 (2013). [CrossRef] [PubMed]

]. However, it should be pointed out that this ultra-long coherence length is not necessary for retinal imaging, since typically only ~2 mm depth range are necessary. Retinal imaging with a VCSEL has been performed at up to 580 kHz axial line rate [45

45. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, C. D. Lu, J. Jiang, A. E. Cable, J. S. Duker, and J. G. Fujimoto, “Retinal, anterior segment and full eye imaging using ultrahigh speed swept source OCT with vertical-cavity surface emitting lasers,” Biomed. Opt. Express 3(11), 2733–2751 (2012). [CrossRef] [PubMed]

].

(2) In Fourier domain mode-locked (FDML) lasers [41

41. R. Huber, M. Wojtkowski, and J. G. Fujimoto, “Fourier Domain Mode Locking (FDML): A new laser operating regime and applications for optical coherence tomography,” Opt. Express 14(8), 3225–3237 (2006). [CrossRef] [PubMed]

], the cavity round trip time is synchronized to the intracavity wavelength filter frequency. Therefore, the fundamental limit to laser tuning rate can be overcome, and single beam scan rates of up to 5.2 MHz were demonstrated at 1300 nm wavelength [4

4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

]. With chromatic dispersion management [46

46. B. R. Biedermann, W. Wieser, C. M. Eigenwillig, T. Klein, and R. Huber, “Dispersion, coherence and noise of Fourier domain mode locked lasers,” Opt. Express 17(12), 9947–9961 (2009). [CrossRef] [PubMed]

, 47

47. S. Todor, B. Biedermann, W. Wieser, R. Huber, and C. Jirauschek, “Instantaneous lineshape analysis of Fourier domain mode-locked lasers,” Opt. Express 19(9), 8802–8807 (2011). [CrossRef] [PubMed]

], the coherence length of FDML lasers was enhanced considerably to more than 20 mm [48

48. D. C. Adler, W. Wieser, F. Trepanier, J. M. Schmitt, and R. A. Huber, “Extended coherence length Fourier domain mode locked lasers at 1310 nm,” Opt. Express 19(21), 20930–20939 (2011). [CrossRef] [PubMed]

]. FDML lasers can also have high output power directly from the cavity, on the order of tens of milliwatts [4

4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

, 23

23. T. Klein, W. Wieser, C. M. Eigenwillig, B. R. Biedermann, and R. Huber, “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express 19(4), 3044–3062 (2011). [CrossRef] [PubMed]

, 49

49. S. Marschall, T. Klein, W. Wieser, B. R. Biedermann, K. Hsu, K. P. Hansen, B. Sumpf, K. H. Hasler, G. Erbert, O. B. Jensen, C. Pedersen, R. Huber, and P. E. Andersen, “Fourier domain mode-locked swept source at 1050 nm based on a tapered amplifier,” Opt. Express 18(15), 15820–15831 (2010). [CrossRef] [PubMed]

].

(3) Since a sweep corresponds to an extremely chirped pulse, very high sweep speeds can also be generated by temporally dispersing ultrashort pulses [42

42. S. Moon and D. Y. Kim, “Ultra-high-speed optical coherence tomography with a stretched pulse supercontinuum source,” Opt. Express 14(24), 11575–11584 (2006). [CrossRef] [PubMed]

]. To date, record sweep speeds of 90.9 MHz were demonstrated, albeit only at a sensitivity of 42 dB, unsuitable for imaging of biological tissue [50

50. K. Goda, A. Fard, O. Malik, G. Fu, A. Quach, and B. Jalali, “High-throughput optical coherence tomography at 800 nm,” Opt. Express 20(18), 19612–19617 (2012). [CrossRef] [PubMed]

].

(4) In wavelength swept ASE sources, light is repeatedly filtered and amplified by a series of filter and gain elements, which form no laser cavity [43

43. C. M. Eigenwillig, B. R. Biedermann, W. Wieser, and R. Huber, “Wavelength swept amplified spontaneous emission source,” Opt. Express 17(21), 18794–18807 (2009). [CrossRef] [PubMed]

]. Hence, the instantaneous linewidth is directly given by the number of consecutive filter events. Using two filter events, retinal imaging has been performed with a fast sweep rate of 340 kHz [51

51. C. M. Eigenwillig, T. Klein, W. Wieser, B. R. Biedermann, and R. Huber, “Wavelength swept amplified spontaneous emission source for high speed retinal optical coherence tomography at 1060 nm,” J Biophotonics 4(7-8), 552–558 (2011). [CrossRef] [PubMed]

].

Due to the intrinsic limitations of all spectrometer solutions and other swept laser sources that have been used for OCT imaging, currently swept-source OCT with either VCSEL or FDML laser sources seems to be the most promising way for (multi-) MHz speeds in retinal OCT.

1.3. Challenges of MHz speeds

As discussed in the first section, increased speed offers a multitude of advantages over slower systems for retinal imaging. However, line rates in the MHz and multi-MHz range also raise a couple of challenges. First, the high speed imposes strict requirements on many parts of the OCT system, not only on the light source. Second, increased speed goes in hand with decreased sensitivity, which may limit the reasonable speed to only a few MHz maximum. We will briefly discuss both of these challenges for SS-OCT, starting with sensitivity.

For in-vivo retinal imaging of humans, power incident on the patient is limited by international standards, such that the sample arm power cannot be arbitrarily increased. Since shot noise limited sensitivity scales inversely with the optical power incident on the sample, sensitivity decreases with increased speed. For retinal imaging, the incident power level is limited by maximum permissible exposure (MPE) values of standards such as ANSI Z136.1, ISO 15004-2 or IEC 60825-1. For instance, ISO 15004-2 limits power of continuous non-scanned illumination incident on the cornea in typical OCT configurations to 1.4 mW at 1050 nm, while the MPE is slightly higher for the other standards.

With 1.4 mW on the cornea, even shot noise limited sensitivity falls below a reasonable mark of ~90 dB when coupling and interferometer losses are taken into account (Fig. 1(a)
Fig. 1 Sensitivity vs. axial line rate for retinal imaging at 1050 nm (0.7A/W photodetector responsivity) for a) 1.4 mW and b) 3.5 mW sample arm power. Blue shaded regions: Shot noise limited sensitivity with 100% collection efficiency. Red shaded regions: sensitivity with a 3 dB penalty (for instance, 65% coupling efficiency and additional 1.13 dB penalty). Horizontal black lines indicate the limits for 95 dB (“clinical imaging”) and 90 dB (“research”) sensitivity.
). The calculation is based on the following formula for shot-noise limited sensitivity S [52

52. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 28(21), 2067–2069 (2003). [CrossRef] [PubMed]

]:
S=10log(ρPTe)IL,
(1)
where ρ=0.7A/Wis a typical detector responsivity at 1060 nm, Pis the optical power incident on the sample, Tis the duration of an A-scan (inverse A-scan rate), eis the elementary charge of the electron, and ILis loss on the way back from the sample to the detector. The formula assumes that apart fromIL, all power returning from the sample is incident on the photodetector. For clinical imaging, the requirements on sensitivity may even be higher than in a laboratory setting, since the collected signal intensity may be severely lowered by media opacities and optical aberrations of the patient eye. Assuming a minimum required sensitivity of 95 dB and a IL=3dBin the setup, the maximum axial line rate is below 1 MHz for 1.4 dB sample arm power. Hence, for fast systems it is very important to assure a high collection efficiency, and a suitable interferometer design for sensitivity close to the shot noise limit [53

53. A. M. Rollins and J. A. Izatt, “Optimal interferometer designs for optical coherence tomography,” Opt. Lett. 24(21), 1484–1486 (1999). [CrossRef] [PubMed]

]. Additionally, there are various pathways to acceptable signal strength:

First, higher optical power levels may be applicable when beam scanning and shorter exposure times are taken into account. However, this requires additional system safety mechanisms to assure that power is always below MPE, even if the beam scanning system fails or exposure time is extended. Figure 1(b) shows sensitivity vs. speed for an optical power of 3.5 mW, which has already been applied for in-vivo human imaging at 1050 nm [54

54. B. Považay, B. Hermann, A. Unterhuber, B. Hofer, H. Sattmann, F. Zeiler, J. E. Morgan, C. Falkner-Radler, C. Glittenberg, S. Blinder, and W. Drexler, “Three-dimensional optical coherence tomography at 1050 nm versus 800 nm in retinal pathologies: enhanced performance and choroidal penetration in cataract patients,” J. Biomed. Opt. 12(4), 041211 (2007). [CrossRef] [PubMed]

]. At this power level, multi-MHz retinal imaging with high sensitivity is feasible even in the presence of interferometer losses.

Second, effective sensitivity is increased by averaging. As already discussed, averaging is already commonly employed to reduce speckle noise. In fact, the perceived OCT image quality depends so critically on speckle reduction that, in the long term, it is likely that most retinal OCT systems will employ some type of averaging. Hence, the signal level is automatically increased, which readily alleviates the detrimental effects of higher speeds. Since SNR increases with the square root of the number of averaged frames N for incoherent averaging [55

55. M. Szkulmowski and M. Wojtkowski, “Averaging techniques for OCT imaging,” Opt. Express 21(8), 9757–9773 (2013). [CrossRef] [PubMed]

], even a moderate amount of 10x averaging will theoretically increase SNR by 5 dB. We will show later in this paper that high image quality can be achieved at multi-MHz speeds and low sample arm power, using different averaging schemes.

Third, optical power may be even further increased with multi-beam approaches, in which each beam illuminates a different sample region [56

56. N. Suehira, S. Ooto, M. Hangai, K. Matsumoto, N. Tomatsu, T. Yuasa, K. Yamada, and N. Yoshimura, “Three-beam spectral-domain optical coherence tomography for retinal imaging,” J. Biomed. Opt. 17(10), 106001 (2012). [CrossRef] [PubMed]

]. Hence, the effective line rate increases with the number of spatially separated beams, whereas the sensitivity for each beam remains constant. This multi-beam approach is especially promising for swept-source OCT, since the detection scheme can be more easily extended to a multi-channel setup than in spectral domain OCT. Multi-beam swept-source imaging has been demonstrated at 1300nm with 20 MHz line rate [4

4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

]. For retinal imaging, two beams have been used to achieve 400 kHz line rate in a SS-OCT setup at 1050 nm [37

37. B. Potsaid, B. Baumann, D. Huang, S. Barry, A. E. Cable, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Ultrahigh speed 1050nm swept source/Fourier domain OCT retinal and anterior segment imaging at 100,000 to 400,000 axial scans per second,” Opt. Express 18(19), 20029–20048 (2010). [CrossRef] [PubMed]

], and three beams at 840 nm for 210 kHz line rate [56

56. N. Suehira, S. Ooto, M. Hangai, K. Matsumoto, N. Tomatsu, T. Yuasa, K. Yamada, and N. Yoshimura, “Three-beam spectral-domain optical coherence tomography for retinal imaging,” J. Biomed. Opt. 17(10), 106001 (2012). [CrossRef] [PubMed]

].

Fourth, the amount of collected optical power also depends on the collection efficiency of the system. Due to optical aberrations, the collection efficiency reaches a maximum for rather small numerical apertures in retinal imaging. Adaptive optics OCT and joint-aperture OCT (JA-OCT) have been demonstrated to increase the collection efficiency of retinal OCT. JA-OCT additionally provides intrinsic speckle reduction while maintaining MHz line rates [57

57. T. Klein, R. André, W. Wieser, T. Pfeiffer, and R. Huber, “Joint aperture detection for speckle reduction and increased collection efficiency in ophthalmic MHz OCT,” Biomed. Opt. Express 4(4), 619–634 (2013). [CrossRef] [PubMed]

].

Besides these limitations considering OCT system sensitivity, the system requirements for (multi-) MHz OCT currently push the limits of the detection and the beam deflection systems. Fortunately, available electronic bandwidth, analog-to-digital (AD) sampling rate, and computational power have increased considerably over the last years, such that sustained acquisition of large data sets as demonstrated in this paper is possible. However, the mechanical specifications of the beam deflection system have not evolved accordingly. Standard galvanometer scanners restrict the frame rate to less than ~1 kHz, which is not enough to take full advantage of multi-MHz OCT systems. While resonant scanners provide high speed, they are not adjustable in terms of frame repetition rate, which is inconvenient for many applications. Both the sensitivity and beam steering issue may be overcome with parallel illumination. With line and full-field illumination, optical power on the sample can be increased considerably [5

5. B. Povazay, B. Hofer, C. Torti, B. Hermann, A. R. Tumlinson, M. Esmaeelpour, C. A. Egan, A. C. Bird, and W. Drexler, “Impact of enhanced resolution, speed and penetration on three-dimensional retinal optical coherence tomography,” Opt. Express 17(5), 4134–4150 (2009). [CrossRef] [PubMed]

]. Very high speeds can be achieved, but crosstalk between pixels may lead to a severe degradation of image quality compared to standard confocal OCT. Maintaining confocal beam delivery, multi-beam OCT also reduces requirements on the beam deflection system.

2. Dispersion-compensated 1050nm FDML laser

2.1. Laser Setup

2.2. Laser characterization

The laser spectrum and coherence properties were characterized directly after the cavity output, i.e. before the buffer stage (see Fig. 2), using an optical spectrum analyzer (OSA, Yokogawa AQ6370B), a motorized Mach-Zehnder interferometer (MZI, home built), a 1 GHz dual-balanced photodetector (BPD, Wieserlabs WL-BPD1GA), and a fast digital oscilloscope (Tektronix DPO 5104).

Figure 3
Fig. 3 Laser spectra and sensitivity roll-off at axial line rates of a) 1.68 MHz and b) 3.35 MHz. Note that for long imaging depths, sensitivity decay is dominated by the detection system’s 1 GHz electronic bandwidth.
(left) shows the laser output spectrum at both 1.68 MHz and 3.35 MHz, with a sweep range of 72 nm at both speeds. The sweep range of 72 nm corresponds to the 3 dB reflection bandwidth of the chirped fiber Bragg grating. Sensitivity decay with imaging depth was measured with the MZI, the BPD and the oscilloscope at 10 GS/s sampling rate (Fig. 3 right). The roll-off measurement is dominated by the electronic bandwidth of the detection system (1 GHz for both the BPD and the oscilloscope). Due to AC-coupling of the BPD, maximum sensitivity is reached a couple of hundred microns away from zero delay.

The signal remains almost constant up to 2 mm / 1 mm at 1.68 MHz / 3.35 MHz, respectively, and drops by about 6 dB (20*log scale) at around 3.8 mm (1.68 MHz) / 1.6 mm (3.35 MHz). Our 3.8 mm roll-off at 1.68 MHz sweep rate compares well with the 2.2 mm reported for a 1050 nm VCSEL based OCT at only 580 kHz [45

45. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, C. D. Lu, J. Jiang, A. E. Cable, J. S. Duker, and J. G. Fujimoto, “Retinal, anterior segment and full eye imaging using ultrahigh speed swept source OCT with vertical-cavity surface emitting lasers,” Biomed. Opt. Express 3(11), 2733–2751 (2012). [CrossRef] [PubMed]

].

3. OCT imaging setup, data acquisition and processing

The OCT imaging setup was designed to enable sustained multi-MHz single-beam and dual-beam retinal imaging over a wide field of view. These requirements can be divided into two groups: First, the optical system in the sample arm has to be designed for wide-field retinal imaging. Since alignment is critical in wide-field retinal imaging, a fixation target and visual alignment control should be available. Second, the detection system needs to offer enough bandwidth (analog electronic and digital sampling rate), while at the same time allowing detection close to shot noise. As will be shown later, it is also beneficial to provide a means for sustained acquisition of very large data sets, on the order of tens of gigabytes. The basic interferometer layout and optical design has been adapted from [23

23. T. Klein, W. Wieser, C. M. Eigenwillig, B. R. Biedermann, and R. Huber, “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express 19(4), 3044–3062 (2011). [CrossRef] [PubMed]

] (Fig. 4
Fig. 4 Dual-beam OCT imaging setup. BPD: Balanced photo-receiver; PC: polarization controller; DC: dispersion compensating glass blocks; Recal.: recalibration arm; Ref.: reference arm; Galvo: galvanometer mirrors (x and y), L1/L2: relay lenses; BS: beam-splitter; TG: target; DM: Dichroic mirror. All fiber couplers (blue) have a 50/50 coupling ratio.
), with improvements already used in [30

30. C. Blatter, T. Klein, B. Grajciar, T. Schmoll, W. Wieser, R. Andre, R. Huber, and R. A. Leitgeb, “Ultrahigh-speed non-invasive widefield angiography,” J. Biomed. Opt. 17(7), 070505 (2012). [CrossRef] [PubMed]

]. The fiber output was collimated with an aspheric lens of 11 mm focal length (Thorlabs C220TME-C). A combination of 8 spherical two-inch lenses, arranged in two groups (L1 and L2), was used to relay the scanner pivot point into the human pupil, where the beam had a diameter of approximately 1.1 mm (1/e2 intensity). A three inch dichroic mirror (Layertec GmbH, reflectivity > 99.9% from 1000 nm to 1090 nm) placed between the two relay lens elements coupled a fixation target and a camera for pupil alignment into the sample arm. For wide-field retinal imaging, especially in the non-mydriatic eye, it is important that the incoming OCT beam passes through the center of the pupil to reduce tilt of the retina in the OCT images. Moreover, we experienced that, for our setup, best image quality can be achieved when the pivot point of the imaging beam is located approximately at the nodal point of the eye, close to the posterior surface of the crystalline lens. Axial alignment to the nodal point also reduces retinal curvature in the OCT images [64

64. R. J. Zawadzki, A. R. Fuller, S. S. Choi, D. F. Wiley, B. Hamann, and J. S. Werner, “Correction of motion artifacts and scanning beam distortions in 3D ophthalmic optical coherence tomography imaging,” Proc. SPIE 64260, 64260 (2007).

].

It should be emphasized that this alignment strategy is different to the common approach, where the axial position of the pivot point is placed close to the pupil, in order to avoid clipping of the extreme rays [64

64. R. J. Zawadzki, A. R. Fuller, S. S. Choi, D. F. Wiley, B. Hamann, and J. S. Werner, “Correction of motion artifacts and scanning beam distortions in 3D ophthalmic optical coherence tomography imaging,” Proc. SPIE 64260, 64260 (2007).

, 65

65. A. Podoleanu, I. Charalambous, L. Plesea, A. Dogariu, and R. Rosen, “Correction of distortions in optical coherence tomography imaging of the eye,” Phys. Med. Biol. 49(7), 1277–1294 (2004). [CrossRef] [PubMed]

].

The OCT signal was detected with a 1GHz dual-balanced photodiode (Wieserlabs WL-BPD1GA), low-pass filtered (BLP-600 + , Minicircuits) and digitized by an 8 bit PCIe data acquisition card at 1.5 GS/s (A/D card, Signatec PX1500-4). At this sampling rate, the A/D card enabled sustained data transfer to host computer (PC) memory with 85% duty cycle. Hence, our data set size is no longer limited by the onboard memory of the digitizer, which is usually much more expensive than PC RAM. In our case, more than 44 GB of data can be acquired using 48 GB of PC memory. After data acquisition, the OCT data set was stored on the computer hard drive and processed using standard OCT signal processing, including numerical dispersion compensation. Due to the ultrahigh frame rate, the fast galvanometer scanner was driven with a sinusoidal signal. Both scan directions were used with an 85% duty cycle, and images were corrected for nonlinear scanning via linear interpolation in post-processing [4

4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

, 63

63. T. Klein, W. Wieser, R. Andre, T. Pfeiffer, C. M. Eigenwillig, and R. Huber, “Multi-MHz FDML OCT: snapshot retinal imaging at 6.7 million axial-scans per second,” in (SPIE, 2012), 82131E–82136.

]. In our processing chain it appeared advantageous to carry out interpolation on the logarithmized OCT intensity images to minimize interpolation artifacts caused by noise. However, the physically correct way to average sample reflectivity is to average OCT signal intensity, i.e. the squared OCT signal amplitude. Hence, we performed a multi-step procedure for averaging multiple B-frames: First, linear intensity (squared amplitude) values from each scan direction were averaged before logarithmization and interpolation. Subsequently, the interpolated logarithmic images from both scan directions were averaged. Cut-levels were used to restrict the dynamic range to 25 dB in the logarithmic images, with the lower cut-level automatically set 2 dB lower than the mean noise floor for each data set.

4. MHz and multi-MHz retinal imaging

4.1. Overview

In-vivo non-mydriatic retinal imaging was performed on a healthy volunteer in accordance with the tenets of the Declaration of Helsinki. To assure comparability of the different imaging speeds and protocols, the same eye of the same subject was imaged for all data sets presented in this paper.

Three imaging modes were investigated: Sustained multi-data set acquisition at 1.68 MHz, ultra-dense sampling at 3.35 MHz and multispot imaging at an effective axial line rate of 6.7 MHz. The imaging modes were chosen to investigate the benefits and challenges of (multi-) MHz OCT: At 1.68 MHz, speed is already high enough to acquire a wide-field view of the fundus in less than a second. Hence, the resulting data sets have a low probability of distortion. By acquisition of many data sets the remaining distortion can be compensated by simple alignment and averaging. We demonstrate a combination of multi-data set and B-frame averaging for speckle reduction and higher signal levels. At 3.35 MHz, the sensitivity of our system dropped below 90 dB, resulting in low signal strength. We will show that massive frame averaging from an ultra densely sampled data set yields high image quality and sufficient signal at this multi-MHz axial line rate. Finally, we demonstrate two-beam imaging at 6.7 MHz as a means to further increase speed. Table 1

Table 1. Imaging parameters

table-icon
View This Table
gives an overview of the imaging system configuration and performance parameters for each of the three modes. Note that the optical power incident on the cornea was lowered to 1.4 mW at 1.68 MHz.

4.2. 1.68 MHz imaging – sustained acquisition of 24 data sets

At 1.68 MHz laser sweep rate, ultra wide-field fundus coverage with dense isotropic sampling is possible with acquisition times of less than a second. Here, we acquired more than a million axial scans in 0.83 s, including galvanometer dead time. This short acquisition time already yields undistorted data sets with high probability. Moreover, even the slow galvanometer axis moves with an angular velocity of about 60 °/s, on the order of the typical saccade speed [68

68. S. Martinez-Conde, S. L. Macknik, X. G. Troncoso, and D. H. Hubel, “Microsaccades: a neurophysiological analysis,” Trends Neurosci. 32(9), 463–475 (2009). [CrossRef] [PubMed]

]. This means that distortions by saccades and other involuntary patient movements are spread out over large areas. Hence more B-frames are affected, and the relative displacement between successive B-frames is reduced, leading to a smoother transition between undistorted and distorted areas. Figure 5
Fig. 5 Retinal imaging at 1.68 MHz axial line rate. More than one million axial scans were acquired in 0.83 s. a) En-face fundus view, reconstructed from the OCT data set. b) Unaveraged B-frame, at the position indicated by the red bar in the en-face view. Media 1 shows an unaveraged flight-through of the entire data set at a playback speed 50 frames per second, i.e. much slower than the OCT acquisition frame rate of 1310 frames per second.
shows an en-face view of a data set consisting of 10882 A-scans, and a representative, unaveraged B-frame. The total imaging depth of 3.4 mm in air is sufficient for wide-field imaging, covering the macula, the optical nerve head (ONH) and parts of the periphery. Apart from one horizontal line in the middle of the B-frame, no artifacts are visible in the image. Sensitivity was measured to be 91 dB with a neutral density filter and a test mirror, with 1.4 mW sample arm power. Taking into account 3 dB loss in the interferometer and 65% coupling efficiency, this measured sensitivity value is shot-noise limited. We chose a power of 1.4 mW to allow for uninterrupted sequential acquisition of many data sets. As already discussed, this acquisition of tens of gigabytes of data was enabled by data sustained data transfer from the A/D card to host computer memory. With a total of 24 acquired data sets, new averaging protocols can be applied to reduce the impact of involuntary eye motion and to reduce speckle noise.

4.3 The simple “acquire-align-average” (AAA) processing approach

There have been various publications about averaging of several OCT volumes [15

15. K. V. Vienola, B. Braaf, C. K. Sheehy, Q. Yang, P. Tiruveedhula, D. W. Arathorn, J. F. de Boer, and A. Roorda, “Real-time eye motion compensation for OCT imaging with tracking SLO,” Biomed. Opt. Express 3(11), 2950–2963 (2012). [CrossRef] [PubMed]

, 67

67. B. Braaf, K. A. Vermeer, V. A. D. P. Sicam, E. van Zeeburg, J. C. van Meurs, and J. F. de Boer, “Phase-stabilized optical frequency domain imaging at 1-µm for the measurement of blood flow in the human choroid,” Opt. Express 19(21), 20886–20903 (2011). [CrossRef] [PubMed]

, 69

69. M. F. Kraus, B. Potsaid, M. A. Mayer, R. Bock, B. Baumann, J. J. Liu, J. Hornegger, and J. G. Fujimoto, “Motion correction in optical coherence tomography volumes on a per A-scan basis using orthogonal scan patterns,” Biomed. Opt. Express 3(6), 1182–1199 (2012). [CrossRef] [PubMed]

71

71. S. Ricco, M. Chen, H. Ishikawa, G. Wollstein, and J. Schuman, “Correcting Motion Artifacts in Retinal Spectral Domain Optical Coherence Tomography via Image Registration,” in Medical Image Computing and Computer-Assisted Intervention – MICCAI 2009, G.-Z. Yang, D. Hawkes, D. Rueckert, A. Noble, and C. Taylor, eds. (Springer Berlin Heidelberg, 2009), pp. 100–107.

]. Most of them are either complex, or slow, or require a tracking system during data acquisition. Here we demonstrate a simple and very fast approach, which still shows remarkably good results. Our approach has some similarities to automatic volume-stitching (mosaicing) methods. However, instead of acquiring many volumes from different locations of the retina, the same volume is repeatedly scanned to form a series of data sets. Afterwards, these data sets are automatically aligned to each other and averaged, an approach that we call “acquire-align-average (AAA)”, which produces single highly averaged images (Fig. 6
Fig. 6 a) Multi-data set averaging at 1.68 MHz with the acquire-align-average (AAA) scheme. After acquisition of N data sets, the data sets are aligned and averaged. Thus, distortion per frame is “averaged out” with increasing N. b) Result of AAA approach to imaging at 1.68 MHz. 4 frames from each of the 24 data sets were averaged (i.e. a total of 96 frames) to yield strong speckle reduction. c) Enlarged region from the region indicated by the white frame in b). All retinal layers including the ELM are clearly visible, indicating that alignment worked well. Note that image displayed in the right column of a) is the averaged en-face image after AAA processing of the 24 en-face images. Media 2 shows the en-face reconstructions of all acquired data sets before registration at a playback speed of 2 volumes per second. Slight “zipper” artifacts can be observed due to instabilities in the galvanometer scanners operating at very high speed. A blinking artifact is clearly visible in the first data set.
). Since image distortions are randomly spread over all data sets, distortions will be statistically removed from the average of all data sets. In AAA, transverse alignment based on the OCT en-face projections is performed first. Second, B-frames from the same location in each of these transversally aligned data sets are aligned to each other and averaged. This scheme yields highly averaged images for both en-face views and B-frames. It is conceptually similar to extracting motion information from an SLO channel in post-processing – since the OCT volume rate is fast enough, the OCT en-face image can be used instead of the SLO image. The AAA approach does obviously not produce a single corrected 3D data set, since axial shift will still be observed in the direction of the slow scan axis. However, in almost all cases (depth-resolved) en-face views and B-frames are the only views used anyway. Moreover, the additional 3D information is mainly used to generate thickness maps, which should not be affected by the improper 3D alignment along the slow axis. Thus, all information that is used in current clinical retinal OCT systems will benefit from the averaging-based signal improvement in AAA. In future, full 3D correction may be carried out, for instance by using a reference B-scan along the slow axis [64

64. R. J. Zawadzki, A. R. Fuller, S. S. Choi, D. F. Wiley, B. Hamann, and J. S. Werner, “Correction of motion artifacts and scanning beam distortions in 3D ophthalmic optical coherence tomography imaging,” Proc. SPIE 64260, 64260 (2007).

]. Other methods, which have been used for mosaicing, may be employed as well [10

10. Y. Li, G. Gregori, B. L. Lam, and P. J. Rosenfeld, “Automatic montage of SD-OCT data sets,” Opt. Express 19(27), 26239–26248 (2011). [CrossRef] [PubMed]

, 72

72. M. Niemeijer, M. K. Garvin, K. Lee, B. van Ginneken, M. D. Abràmoff, and M. Sonka, “Registration of 3D spectral OCT volumes using 3D SIFT feature point matching ” in Proc. of SPIE, (2009), p. 72591I.

74

74. B. Antony, M. D. Abràmoff, L. Tang, W. D. Ramdas, J. R. Vingerling, N. M. Jansonius, K. Lee, Y. H. Kwon, M. Sonka, and M. K. Garvin, “Automated 3-D method for the correction of axial artifacts in spectral-domain optical coherence tomography images,” Biomed. Opt. Express 2(8), 2403–2416 (2011). [CrossRef] [PubMed]

], but without the need for special scan patterns [69

69. M. F. Kraus, B. Potsaid, M. A. Mayer, R. Bock, B. Baumann, J. J. Liu, J. Hornegger, and J. G. Fujimoto, “Motion correction in optical coherence tomography volumes on a per A-scan basis using orthogonal scan patterns,” Biomed. Opt. Express 3(6), 1182–1199 (2012). [CrossRef] [PubMed]

, 70

70. H. C. Hendargo, R. Estrada, S. J. Chiu, C. Tomasi, S. Farsiu, and J. A. Izatt, “Automated non-rigid registration and mosaicing for robust imaging of distinct retinal capillary beds using speckle variance optical coherence tomography,” Biomed. Opt. Express 4(6), 803–821 (2013). [CrossRef] [PubMed]

].

After acquisition of all data sets, transverse alignment was performed using the reconstructed en-face images. Here, registration was performed with the StackReg plug-in in FiJi [75

75. P. Thévenaz, U. E. Ruttimann, and M. Unser, “A pyramid approach to subpixel registration based on intensity,” IEEE Trans. Image Process. 7(1), 27–41 (1998). [CrossRef] [PubMed]

, 76

76. J. Schindelin, I. Arganda-Carreras, E. Frise, V. Kaynig, M. Longair, T. Pietzsch, S. Preibisch, C. Rueden, S. Saalfeld, B. Schmid, J.-Y. Tinevez, D. J. White, V. Hartenstein, K. Eliceiri, P. Tomancak, and A. Cardona, “Fiji: an open-source platform for biological-image analysis,” Nat. Methods 9(7), 676–682 (2012). [CrossRef] [PubMed]

], using the rigid body model algorithm that allows translation and rotation only. After registration, the transverse displacement information from the en-face image should be applied to the entire 3D data sets. In our case, the data sets showed only little translation and rotation with respect to each other.

So surprisingly, MHz OCT enables fast averaging of frames from multi-volumes even without the information from transverse alignment. The AAA approach exhibits high resolution and low speckle noise without the need of a hardware eye tracker and without the need of very accurate alignment.

4.4. 3.35 MHz imaging – very densely sampled data sets

First, the galvanometer scanners did not work reliably at larger scan angles at these speeds anymore. Second, the imaging range of our system was reduced to 1.7 mm due to the finite sampling range and electronic bandwidth. Due to the natural curvature of the fundus, the retina moved out of the imaging range at the periphery, such that field of view was limited by A/D sampling rate (see Media 3).

To illustrate how low signal strength can be compensated by averaging, an extremely dense sampling pattern was applied, with 1900 A-scans per B-frame and 9500 B-frames per volume acquired in 6.3 s. While this acquisition time is too long for clinical application, it enables the controlled investigation of averaging very large amounts of adjacent B-frames. Figure 8
Fig. 8 a) 24x and b) 48x average of adjacent B-frames 3.35 MHz axial line rate. Note that, even though the speckle pattern is still partially developed in the 48x averaged image, some anatomical features are already distorted. For instance, the red circle in b) highlights an image detail, which is not visible in the 24x averaged frame. c,d): Enlarged view at the position indicated by the white box in a).
shows the average of 24 and 48 B-frames taken around the location of the B-frame in Fig. 7(b). Note that in this case, no alignment was necessary for averaging. Two points are noteworthy: First, signal intensity is increased to a level where all fundus layers up to the choroidal-scleral interface can be clearly identified. Moreover, fine structures such as retinal capillaries can be readily observed. Second, the B-frame sampling density is so high, that the speckle pattern is still partially visible despite the very large amount of up to 48 averaged frames. This may not be surprising, since the distance between neighboring B-frames is only about 1.3 µm (using a scaling factor of 288 µm/degree, and 9500 B-frames over 42°).

Hence, the 24 averaged frames span a distance of only slightly more than the spot size on the retina. Since the speckle pattern should remain fairly stable within a spot size, the speckle pattern is still visible even with 24x averaging. While speckle reduction is more pronounced with 48x averaging, the larger area covered already causes blurs in steep anatomical structures such as in the optical nerve head.

Furthermore, new anatomical features such as vitreoretinal adhesion are suddenly visible in the 48x averaged frame (see red circle in Fig. 8(b)). As already noted in [57

57. T. Klein, R. André, W. Wieser, T. Pfeiffer, and R. Huber, “Joint aperture detection for speckle reduction and increased collection efficiency in ophthalmic MHz OCT,” Biomed. Opt. Express 4(4), 619–634 (2013). [CrossRef] [PubMed]

], the optimum interframe distance for speckle reduction (see for instance [77

77. M. Szkulmowski, I. Gorczynska, D. Szlag, M. Sylwestrzak, A. Kowalczyk, and M. Wojtkowski, “Efficient reduction of speckle noise in Optical Coherence Tomography,” Opt. Express 20(2), 1337–1359 (2012). [CrossRef] [PubMed]

]) may be too large to keep sufficient transverse resolution. So in summary, the low image quality of 3.35 MHz OCT can be highly improved by averaging. A value of 24 averaged frames yield good image quality especially considering the fact, that 24x averaging of a 3.35 MHz OCT system still results in an effective 140 kHz A-scan rate system.

4.5. 6.7 MHz dual-spot imaging: FDML without dispersion compensation

As a proof-of principle experiment, we performed dual-spot imaging at 2x 3.35 MHz axial line rate with an FDML laser without dispersion compensation. The lower laser performance compared to the dispersion compensated laser highlights the importance of a well designed laser for swept-source multi-MHz imaging. The reduced axial resolution and shorter sensitivity roll-off decay are clearly visible in the resulting OCT images (Fig. 9
Fig. 9 Two-beam imaging at an effective axial line rate of 6.7 MHz, using only 0.85 mW optical power for each beam. a) Reconstructed en-face view: White dotted outline shows the area scanned by beam 1. There is around 10% overlap between the areas scanned by the two beams. b,c) B-frames from beams 1 and 2, as indicated by the red bars in the en-face view. The frames were 2x decimated in the transverse direction. d) 6x B-frame average at the same position as the frame shown in c), with enlarged region e).
). This is also due to the fact that we only used 0.85 mW sample arm optical power for each of the two beams, i.e. the same combined power as for single-beam imaging. Nevertheless, a moderate amount of 6x averaging increases the image quality considerably.

Effective axial line rate with 6x averaging is 1.1 MHz, still higher than the raw speed demonstrated with any non-FMDL light source. The en-face image also shows ciliary shadowing, a typical problem encountered in wide-field retinal imaging [78

78. A. S. Neubauer, M. Kernt, C. Haritoglou, S. G. Priglinger, A. Kampik, and M. W. Ulbig, “Nonmydriatic screening for diabetic retinopathy by ultra-widefield scanning laser ophthalmoscopy (Optomap),” Graefes Arch. Clin. Exp. Ophthalmol. 246(2), 229–235 (2008). [CrossRef] [PubMed]

].

5. Summary and outlook

We presented in-vivo OCT of the posterior pole at axial line rates of multiple MHz for the first time. At these ultrahigh speeds, image quality is affected by several challenges, such as limited sensitivity and imaging range. We have demonstrated three ways to overcome some of these issues: First, already moderate MHz speeds are sufficient for successive acquisition of many data sets. With the resulting large amount of acquired data, high image quality can be achieved. Here, 24 data sets have been acquired in less than a second, with low probability of distortion in each single data set. Hence, averaging reduces distortions and image artifacts are reduced, which are spread out statistically over all data sets. The presented alignment and averaging scheme corresponds to a software based eye-tracking solution. Hence system complexity is considerably reduced compared to commonly used hardware (SLO) based tracking schemes, since no further hardware is required. Second, we have shown that even a further increase to multi-MHz speeds is feasible even with relatively low optical power levels. Sufficient image quality can be obtained simply by averaging a correspondingly larger amount of frames. It might be argued that higher speeds make no sense if the necessary number of averaged increases accordingly, lowering effective line rate again. However, frame averaging will usually be performed anyway for a different reason, namely speckle reduction. We believe that perceived image quality of a fast system with 10x frame averaging will be higher than the image quality of a 10x slower system without averaging, assuming that the initial sensitivity is not too low. Suitable sensitivity values larger than ~90 dB should be easily achievable in ophthalmic multi-MHz OCT, by using higher MPE values for scanned systems. At the same time, even with averaging, the effective line rates in multi-MHz OCT exceed the line rates of commercial systems. For instance, 24x averaging at 3.35 MHz, as demonstrated in this paper, still yields an effective line rate of 140 kHz. Hence, very fast OCT systems are able to perform frame averaging for speckle reduction without reduction of 3D capabilities. Third, multi-beam imaging may be a solution for further increase in speed, especially if the total sample arm power is increased. With faster resonant beam scanners, the demonstrated speed of 6.7 MHz is sufficient for video-rate 3D OCT with reasonably dense sampling. For instance, a volume of 477 x 477 A-scans could be acquired with a volume rate of 25 Hz, using only 85% scanner duty cycle. Real-time video rate ophthalmic 3D OCT may be useful in applications such as OCT for surgical guidance, which until now is restricted to 2D live display or small 3D sampling densities only [79

79. Y. K. Tao, J. P. Ehlers, C. A. Toth, and J. A. Izatt, “Intraoperative spectral domain optical coherence tomography for vitreoretinal surgery,” Opt. Lett. 35(20), 3315–3317 (2010). [CrossRef] [PubMed]

, 80

80. J. U. Kang, Y. Huang, K. Zhang, Z. Ibrahim, J. Cha, W. P. A. Lee, G. Brandacher, and P. L. Gehlbach, “Real-time three-dimensional Fourier-domain optical coherence tomography video image guided microsurgeries,” J. Biomed. Opt. 17(8), 081403 (2012). [CrossRef] [PubMed]

].

Acknowledgments

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 254(5035), 1178–1181 (1991). [CrossRef] [PubMed]

2.

M. R. Hee, J. A. Izatt, E. A. Swanson, D. Huang, J. S. Schuman, C. P. Lin, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography of the human retina,” Arch. Ophthalmol. 113(3), 325–332 (1995). [CrossRef] [PubMed]

3.

M. Wojtkowski, T. Bajraszewski, I. Gorczyńska, P. Targowski, A. Kowalczyk, W. Wasilewski, and C. Radzewicz, “Ophthalmic imaging by spectral optical coherence tomography,” Am. J. Ophthalmol. 138(3), 412–419 (2004). [CrossRef] [PubMed]

4.

W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express 18(14), 14685–14704 (2010). [CrossRef] [PubMed]

5.

B. Povazay, B. Hofer, C. Torti, B. Hermann, A. R. Tumlinson, M. Esmaeelpour, C. A. Egan, A. C. Bird, and W. Drexler, “Impact of enhanced resolution, speed and penetration on three-dimensional retinal optical coherence tomography,” Opt. Express 17(5), 4134–4150 (2009). [CrossRef] [PubMed]

6.

J. M. Schmitt, S. H. Xiang, and K. M. Yung, “Speckle in Optical Coherence Tomography,” J. Biomed. Opt. 4(1), 95–105 (1999). [CrossRef] [PubMed]

7.

M. Wojtkowski, B. Kaluzny, and R. J. Zawadzki, “New Directions in Ophthalmic Optical Coherence Tomography,” Optom. Vis. Sci. 89(5), 524–542 (2012). [CrossRef] [PubMed]

8.

L. Fang, S. Li, Q. Nie, J. A. Izatt, C. A. Toth, and S. Farsiu, “Sparsity based denoising of spectral domain optical coherence tomography images,” Biomed. Opt. Express 3(5), 927–942 (2012). [CrossRef] [PubMed]

9.

R. J. Zawadzki, S. S. Choi, A. R. Fuller, J. W. Evans, B. Hamann, and J. S. Werner, “Cellular resolution volumetric in vivo retinal imaging with adaptive optics-optical coherence tomography,” Opt. Express 17(5), 4084–4094 (2009). [CrossRef] [PubMed]

10.

Y. Li, G. Gregori, B. L. Lam, and P. J. Rosenfeld, “Automatic montage of SD-OCT data sets,” Opt. Express 19(27), 26239–26248 (2011). [CrossRef] [PubMed]

11.

S. Makita, F. Jaillon, M. Yamanari, M. Miura, and Y. Yasuno, “Comprehensive in vivo micro-vascular imaging of the human eye by dual-beam-scan Doppler optical coherence angiography,” Opt. Express 19(2), 1271–1283 (2011). [CrossRef] [PubMed]

12.

R. D. Ferguson, D. X. Hammer, L. A. Paunescu, S. Beaton, and J. S. Schuman, “Tracking optical coherence tomography,” Opt. Lett. 29(18), 2139–2141 (2004). [CrossRef] [PubMed]

13.

Y. K. Tao, S. Farsiu, and J. A. Izatt, “Interlaced spectrally encoded confocal scanning laser ophthalmoscopy and spectral domain optical coherence tomography,” Biomed. Opt. Express 1(2), 431–440 (2010). [CrossRef] [PubMed]

14.

G. Gregori and P. J. Rosenfeld, “Using OCT Fundus Images to Evaluate the Performance of the Spectralis OCT Eye Tracking System (Poster 1032/A434),” in ARVO Annual Meeting 2011, (Fort Lauderdale, Florida, 2011).

15.

K. V. Vienola, B. Braaf, C. K. Sheehy, Q. Yang, P. Tiruveedhula, D. W. Arathorn, J. F. de Boer, and A. Roorda, “Real-time eye motion compensation for OCT imaging with tracking SLO,” Biomed. Opt. Express 3(11), 2950–2963 (2012). [CrossRef] [PubMed]

16.

A. G. Podoleanu and R. B. Rosen, “Combinations of techniques in imaging the retina with high resolution,” Prog. Retin. Eye Res. 27(4), 464–499 (2008). [CrossRef] [PubMed]

17.

S. Jiao, R. Knighton, X. Huang, G. Gregori, and C. Puliafito, “Simultaneous acquisition of sectional and fundus ophthalmic images with spectral-domain optical coherence tomography,” Opt. Express 13(2), 444–452 (2005). [CrossRef] [PubMed]

18.

U. Schmidt-Erfurth, R. A. Leitgeb, S. Michels, B. Považay, S. Sacu, B. Hermann, C. Ahlers, H. Sattmann, C. Scholda, A. F. Fercher, and W. Drexler, “Three-Dimensional Ultrahigh-Resolution Optical Coherence Tomography of Macular Diseases,” Invest. Ophthalmol. Vis. Sci. 46(9), 3393–3402 (2005). [CrossRef] [PubMed]

19.

A. G. Podoleanu, G. M. Dobre, R. G. Cucu, R. Rosen, P. Garcia, J. Nieto, D. Will, R. Gentile, T. Muldoon, J. Walsh, L. A. Yannuzzi, Y. Fisher, D. Orlock, R. Weitz, J. A. Rogers, S. Dunne, and A. Boxer, “Combined multiplanar optical coherence tomography and confocal scanning ophthalmoscopy,” J. Biomed. Opt. 9(1), 86–93 (2004). [CrossRef] [PubMed]

20.

I. Gorczynska, V. J. Srinivasan, L. N. Vuong, R. W. S. Chen, J. J. Liu, E. Reichel, M. Wojtkowski, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Projection OCT fundus imaging for visualising outer retinal pathology in non-exudative age-related macular degeneration,” Br. J. Ophthalmol. 93(5), 603–609 (2009). [CrossRef] [PubMed]

21.

T. Agawa, M. Miura, Y. Ikuno, S. Makita, T. Fabritius, T. Iwasaki, H. Goto, K. Nishida, and Y. Yasuno, “Choroidal thickness measurement in healthy Japanese subjects by three-dimensional high-penetration optical coherence tomography,” Graefes Arch. Clin. Exp. Ophthalmol. 249(10), 1485–1492 (2011). [CrossRef] [PubMed]

22.

B. Povazay, B. Hermann, B. Hofer, V. Kajić, E. Simpson, T. Bridgford, and W. Drexler, “Wide-Field Optical Coherence Tomography of the Choroid In Vivo,” Invest. Ophthalmol. Vis. Sci. 50(4), 1856–1863 (2008). [CrossRef] [PubMed]

23.

T. Klein, W. Wieser, C. M. Eigenwillig, B. R. Biedermann, and R. Huber, “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express 19(4), 3044–3062 (2011). [CrossRef] [PubMed]

24.

A. A. Ellabban, A. Tsujikawa, A. Matsumoto, K. Yamashiro, A. Oishi, S. Ooto, I. Nakata, Y. Akagi-Kurashige, M. Miyake, H. S. Elnahas, T. M. Radwan, K. A. Zaky, and N. Yoshimura, “Three-Dimensional Tomographic Features of Dome-Shaped Macula by Swept-Source Optical Coherence Tomography,” Am. J. Ophthalmol. 155(2), 320–328, e2 (2013). [CrossRef] [PubMed]

25.

S. Makita, Y. Hong, M. Yamanari, T. Yatagai, and Y. Yasuno, “Optical coherence angiography,” Opt. Express 14(17), 7821–7840 (2006). [CrossRef] [PubMed]

26.

T. Schmoll, C. Kolbitsch, and R. A. Leitgeb, “In vivo functional retinal optical coherence tomography,” J. Biomed. Opt. 15(4), 041513 (2010). [CrossRef] [PubMed]

27.

L. An, T. T. Shen, and R. K. Wang, “Using ultrahigh sensitive optical microangiography to achieve comprehensive depth resolved microvasculature mapping for human retina,” J. Biomed. Opt. 16(10), 106013 (2011). [CrossRef] [PubMed]

28.

T. Torzicky, S. Marschall, M. Pircher, B. Baumann, M. Bonesi, S. Zotter, E. Götzinger, W. Trasischker, T. Klein, W. Wieser, B. Biedermann, R. Huber, P. Andersen, and C. K. Hitzenberger, “Retinal polarization-sensitive optical coherence tomography at 1060 nm with 350 kHz A-scan rate using an Fourier domain mode locked laser,” J. Biomed. Opt. 18(2), 026008 (2013). [CrossRef] [PubMed]

29.

B. White, M. Pierce, N. Nassif, B. Cense, B. Park, G. Tearney, B. Bouma, T. Chen, and J. de Boer, “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography,” Opt. Express 11(25), 3490–3497 (2003). [CrossRef] [PubMed]

30.

C. Blatter, T. Klein, B. Grajciar, T. Schmoll, W. Wieser, R. Andre, R. Huber, and R. A. Leitgeb, “Ultrahigh-speed non-invasive widefield angiography,” J. Biomed. Opt. 17(7), 070505 (2012). [CrossRef] [PubMed]

31.

B. Potsaid, I. Gorczynska, V. J. Srinivasan, Y. L. Chen, J. Jiang, A. Cable, and J. G. Fujimoto, “Ultrahigh speed Spectral / Fourier domain OCT ophthalmic imaging at 70,000 to 312,500 axial scans per second,” Opt. Express 16(19), 15149–15169 (2008). [CrossRef] [PubMed]

32.

L. An, P. Li, T. T. Shen, and R. Wang, “High speed spectral domain optical coherence tomography for retinal imaging at 500,000 A‑lines per second,” Biomed. Opt. Express 2(10), 2770–2783 (2011). [CrossRef] [PubMed]

33.

D. Choi, H. Hiro-Oka, H. Furukawa, R. Yoshimura, M. Nakanishi, K. Shimizu, and K. Ohbayashi, “Fourier domain optical coherence tomography using optical demultiplexers imaging at 60,000,000 lines/s,” Opt. Lett. 33(12), 1318–1320 (2008). [CrossRef] [PubMed]

34.

D. H. Choi, H. Hiro-Oka, K. Shimizu, and K. Ohbayashi, “Spectral domain optical coherence tomography of multi-MHz A-scan rates at 1310 nm range and real-time 4D-display up to 41 volumes/second,” Biomed. Opt. Express 3(12), 3067–3086 (2012). [CrossRef] [PubMed]

35.

S. Yun, G. Tearney, J. de Boer, N. Iftimia, and B. Bouma, “High-speed optical frequency-domain imaging,” Opt. Express 11(22), 2953–2963 (2003). [CrossRef] [PubMed]

36.

R. Huber, M. Wojtkowski, K. Taira, J. G. Fujimoto, and K. Hsu, “Amplified, frequency swept lasers for frequency domain reflectometry and OCT imaging: design and scaling principles,” Opt. Express 13(9), 3513–3528 (2005). [CrossRef] [PubMed]

37.

B. Potsaid, B. Baumann, D. Huang, S. Barry, A. E. Cable, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Ultrahigh speed 1050nm swept source/Fourier domain OCT retinal and anterior segment imaging at 100,000 to 400,000 axial scans per second,” Opt. Express 18(19), 20029–20048 (2010). [CrossRef] [PubMed]

38.

W.-Y. Oh, B. J. Vakoc, M. Shishkov, G. J. Tearney, and B. E. Bouma, “>400 kHz repetition rate wavelength-swept laser and application to high-speed optical frequency domain imaging,” Opt. Lett. 35(17), 2919–2921 (2010). [CrossRef] [PubMed]

39.

M. P. Minneman, J. Ensher, M. Crawford, and D. Derickson, “All-semiconductor high-speed akinetic swept-source for OCT,” in Optical Sensors and Biophotonics III, Proc. SPIE 8311 (SPIE, 2011), 831116–831116.

40.

B. Potsaid, V. Jayaraman, J. G. Fujimoto, J. Jiang, P. J. S. Heim, and A. E. Cable, “MEMS tunable VCSEL light source for ultrahigh speed 60kHz - 1MHz axial scan rate and long range centimeter class OCT imaging,” in (SPIE, 2012), 82130M.

41.

R. Huber, M. Wojtkowski, and J. G. Fujimoto, “Fourier Domain Mode Locking (FDML): A new laser operating regime and applications for optical coherence tomography,” Opt. Express 14(8), 3225–3237 (2006). [CrossRef] [PubMed]

42.

S. Moon and D. Y. Kim, “Ultra-high-speed optical coherence tomography with a stretched pulse supercontinuum source,” Opt. Express 14(24), 11575–11584 (2006). [CrossRef] [PubMed]

43.

C. M. Eigenwillig, B. R. Biedermann, W. Wieser, and R. Huber, “Wavelength swept amplified spontaneous emission source,” Opt. Express 17(21), 18794–18807 (2009). [CrossRef] [PubMed]

44.

I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, J. Jiang, J. G. Fujimoto, and A. E. Cable, “High-precision, high-accuracy ultralong-range swept-source optical coherence tomography using vertical cavity surface emitting laser light source,” Opt. Lett. 38(5), 673–675 (2013). [CrossRef] [PubMed]

45.

I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, C. D. Lu, J. Jiang, A. E. Cable, J. S. Duker, and J. G. Fujimoto, “Retinal, anterior segment and full eye imaging using ultrahigh speed swept source OCT with vertical-cavity surface emitting lasers,” Biomed. Opt. Express 3(11), 2733–2751 (2012). [CrossRef] [PubMed]

46.

B. R. Biedermann, W. Wieser, C. M. Eigenwillig, T. Klein, and R. Huber, “Dispersion, coherence and noise of Fourier domain mode locked lasers,” Opt. Express 17(12), 9947–9961 (2009). [CrossRef] [PubMed]

47.

S. Todor, B. Biedermann, W. Wieser, R. Huber, and C. Jirauschek, “Instantaneous lineshape analysis of Fourier domain mode-locked lasers,” Opt. Express 19(9), 8802–8807 (2011). [CrossRef] [PubMed]

48.

D. C. Adler, W. Wieser, F. Trepanier, J. M. Schmitt, and R. A. Huber, “Extended coherence length Fourier domain mode locked lasers at 1310 nm,” Opt. Express 19(21), 20930–20939 (2011). [CrossRef] [PubMed]

49.

S. Marschall, T. Klein, W. Wieser, B. R. Biedermann, K. Hsu, K. P. Hansen, B. Sumpf, K. H. Hasler, G. Erbert, O. B. Jensen, C. Pedersen, R. Huber, and P. E. Andersen, “Fourier domain mode-locked swept source at 1050 nm based on a tapered amplifier,” Opt. Express 18(15), 15820–15831 (2010). [CrossRef] [PubMed]

50.

K. Goda, A. Fard, O. Malik, G. Fu, A. Quach, and B. Jalali, “High-throughput optical coherence tomography at 800 nm,” Opt. Express 20(18), 19612–19617 (2012). [CrossRef] [PubMed]

51.

C. M. Eigenwillig, T. Klein, W. Wieser, B. R. Biedermann, and R. Huber, “Wavelength swept amplified spontaneous emission source for high speed retinal optical coherence tomography at 1060 nm,” J Biophotonics 4(7-8), 552–558 (2011). [CrossRef] [PubMed]

52.

J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 28(21), 2067–2069 (2003). [CrossRef] [PubMed]

53.

A. M. Rollins and J. A. Izatt, “Optimal interferometer designs for optical coherence tomography,” Opt. Lett. 24(21), 1484–1486 (1999). [CrossRef] [PubMed]

54.

B. Považay, B. Hermann, A. Unterhuber, B. Hofer, H. Sattmann, F. Zeiler, J. E. Morgan, C. Falkner-Radler, C. Glittenberg, S. Blinder, and W. Drexler, “Three-dimensional optical coherence tomography at 1050 nm versus 800 nm in retinal pathologies: enhanced performance and choroidal penetration in cataract patients,” J. Biomed. Opt. 12(4), 041211 (2007). [CrossRef] [PubMed]

55.

M. Szkulmowski and M. Wojtkowski, “Averaging techniques for OCT imaging,” Opt. Express 21(8), 9757–9773 (2013). [CrossRef] [PubMed]

56.

N. Suehira, S. Ooto, M. Hangai, K. Matsumoto, N. Tomatsu, T. Yuasa, K. Yamada, and N. Yoshimura, “Three-beam spectral-domain optical coherence tomography for retinal imaging,” J. Biomed. Opt. 17(10), 106001 (2012). [CrossRef] [PubMed]

57.

T. Klein, R. André, W. Wieser, T. Pfeiffer, and R. Huber, “Joint aperture detection for speckle reduction and increased collection efficiency in ophthalmic MHz OCT,” Biomed. Opt. Express 4(4), 619–634 (2013). [CrossRef] [PubMed]

58.

R. Huber, D. C. Adler, V. J. Srinivasan, and J. G. Fujimoto, “Fourier domain mode locking at 1050 nm for ultra-high-speed optical coherence tomography of the human retina at 236,000 axial scans per second,” Opt. Lett. 32(14), 2049–2051 (2007). [CrossRef] [PubMed]

59.

V. J. Srinivasan, D. C. Adler, Y. L. Chen, I. Gorczynska, R. Huber, J. S. Duker, J. S. Schuman, and J. G. Fujimoto, “Ultrahigh-Speed Optical Coherence Tomography for Three-Dimensional and En Face Imaging of the Retina and Optic Nerve Head,” Invest. Ophthalmol. Vis. Sci. 49(11), 5103–5110 (2008). [CrossRef] [PubMed]

60.

W. Wieser, G. Palte, C. M. Eigenwillig, B. R. Biedermann, T. Pfeiffer, and R. Huber, “Chromatic polarization effects of swept waveforms in FDML lasers and fiber spools,” Opt. Express 20(9), 9819–9832 (2012). [CrossRef] [PubMed]

61.

W. Wieser, T. Klein, D. C. Adler, F. Trépanier, C. M. Eigenwillig, S. Karpf, J. M. Schmitt, and R. Huber, “Extended coherence length megahertz FDML and its application for anterior segment imaging,” Biomed. Opt. Express 3(10), 2647–2657 (2012). [CrossRef] [PubMed]

62.

R. Huber, D. C. Adler, and J. G. Fujimoto, “Buffered Fourier domain mode locking: unidirectional swept laser sources for optical coherence tomography imaging at 370,000 lines/s,” Opt. Lett. 31(20), 2975–2977 (2006). [CrossRef] [PubMed]

63.

T. Klein, W. Wieser, R. Andre, T. Pfeiffer, C. M. Eigenwillig, and R. Huber, “Multi-MHz FDML OCT: snapshot retinal imaging at 6.7 million axial-scans per second,” in (SPIE, 2012), 82131E–82136.

64.

R. J. Zawadzki, A. R. Fuller, S. S. Choi, D. F. Wiley, B. Hamann, and J. S. Werner, “Correction of motion artifacts and scanning beam distortions in 3D ophthalmic optical coherence tomography imaging,” Proc. SPIE 64260, 64260 (2007).

65.

A. Podoleanu, I. Charalambous, L. Plesea, A. Dogariu, and R. Rosen, “Correction of distortions in optical coherence tomography imaging of the eye,” Phys. Med. Biol. 49(7), 1277–1294 (2004). [CrossRef] [PubMed]

66.

Y. Chen, D. M. de Bruin, C. Kerbage, and J. F. de Boer, “Spectrally balanced detection for optical frequency domain imaging,” Opt. Express 15(25), 16390–16399 (2007). [CrossRef] [PubMed]

67.

B. Braaf, K. A. Vermeer, V. A. D. P. Sicam, E. van Zeeburg, J. C. van Meurs, and J. F. de Boer, “Phase-stabilized optical frequency domain imaging at 1-µm for the measurement of blood flow in the human choroid,” Opt. Express 19(21), 20886–20903 (2011). [CrossRef] [PubMed]

68.

S. Martinez-Conde, S. L. Macknik, X. G. Troncoso, and D. H. Hubel, “Microsaccades: a neurophysiological analysis,” Trends Neurosci. 32(9), 463–475 (2009). [CrossRef] [PubMed]

69.

M. F. Kraus, B. Potsaid, M. A. Mayer, R. Bock, B. Baumann, J. J. Liu, J. Hornegger, and J. G. Fujimoto, “Motion correction in optical coherence tomography volumes on a per A-scan basis using orthogonal scan patterns,” Biomed. Opt. Express 3(6), 1182–1199 (2012). [CrossRef] [PubMed]

70.

H. C. Hendargo, R. Estrada, S. J. Chiu, C. Tomasi, S. Farsiu, and J. A. Izatt, “Automated non-rigid registration and mosaicing for robust imaging of distinct retinal capillary beds using speckle variance optical coherence tomography,” Biomed. Opt. Express 4(6), 803–821 (2013). [CrossRef] [PubMed]

71.

S. Ricco, M. Chen, H. Ishikawa, G. Wollstein, and J. Schuman, “Correcting Motion Artifacts in Retinal Spectral Domain Optical Coherence Tomography via Image Registration,” in Medical Image Computing and Computer-Assisted Intervention – MICCAI 2009, G.-Z. Yang, D. Hawkes, D. Rueckert, A. Noble, and C. Taylor, eds. (Springer Berlin Heidelberg, 2009), pp. 100–107.

72.

M. Niemeijer, M. K. Garvin, K. Lee, B. van Ginneken, M. D. Abràmoff, and M. Sonka, “Registration of 3D spectral OCT volumes using 3D SIFT feature point matching ” in Proc. of SPIE, (2009), p. 72591I.

73.

S. Lee, M. Young, M. V. Sarunic, and M. F. Beg, “End-to-end pipeline for spectral domain optical coherence tomography and morphometric analysis of human optic nerve head,” J. Med. Biol. Eng. 31, 9 (2010).

74.

B. Antony, M. D. Abràmoff, L. Tang, W. D. Ramdas, J. R. Vingerling, N. M. Jansonius, K. Lee, Y. H. Kwon, M. Sonka, and M. K. Garvin, “Automated 3-D method for the correction of axial artifacts in spectral-domain optical coherence tomography images,” Biomed. Opt. Express 2(8), 2403–2416 (2011). [CrossRef] [PubMed]

75.

P. Thévenaz, U. E. Ruttimann, and M. Unser, “A pyramid approach to subpixel registration based on intensity,” IEEE Trans. Image Process. 7(1), 27–41 (1998). [CrossRef] [PubMed]

76.

J. Schindelin, I. Arganda-Carreras, E. Frise, V. Kaynig, M. Longair, T. Pietzsch, S. Preibisch, C. Rueden, S. Saalfeld, B. Schmid, J.-Y. Tinevez, D. J. White, V. Hartenstein, K. Eliceiri, P. Tomancak, and A. Cardona, “Fiji: an open-source platform for biological-image analysis,” Nat. Methods 9(7), 676–682 (2012). [CrossRef] [PubMed]

77.

M. Szkulmowski, I. Gorczynska, D. Szlag, M. Sylwestrzak, A. Kowalczyk, and M. Wojtkowski, “Efficient reduction of speckle noise in Optical Coherence Tomography,” Opt. Express 20(2), 1337–1359 (2012). [CrossRef] [PubMed]

78.

A. S. Neubauer, M. Kernt, C. Haritoglou, S. G. Priglinger, A. Kampik, and M. W. Ulbig, “Nonmydriatic screening for diabetic retinopathy by ultra-widefield scanning laser ophthalmoscopy (Optomap),” Graefes Arch. Clin. Exp. Ophthalmol. 246(2), 229–235 (2008). [CrossRef] [PubMed]

79.

Y. K. Tao, J. P. Ehlers, C. A. Toth, and J. A. Izatt, “Intraoperative spectral domain optical coherence tomography for vitreoretinal surgery,” Opt. Lett. 35(20), 3315–3317 (2010). [CrossRef] [PubMed]

80.

J. U. Kang, Y. Huang, K. Zhang, Z. Ibrahim, J. Cha, W. P. A. Lee, G. Brandacher, and P. L. Gehlbach, “Real-time three-dimensional Fourier-domain optical coherence tomography video image guided microsurgeries,” J. Biomed. Opt. 17(8), 081403 (2012). [CrossRef] [PubMed]

OCIS Codes
(120.3890) Instrumentation, measurement, and metrology : Medical optics instrumentation
(140.3510) Lasers and laser optics : Lasers, fiber
(170.3880) Medical optics and biotechnology : Medical and biological imaging
(170.4460) Medical optics and biotechnology : Ophthalmic optics and devices
(170.4500) Medical optics and biotechnology : Optical coherence tomography

ToC Category:
Optical Coherence Tomography

History
Original Manuscript: May 29, 2013
Revised Manuscript: July 17, 2013
Manuscript Accepted: July 18, 2013
Published: August 30, 2013

Citation
Thomas Klein, Wolfgang Wieser, Lukas Reznicek, Aljoscha Neubauer, Anselm Kampik, and Robert Huber, "Multi-MHz retinal OCT," Biomed. Opt. Express 4, 1890-1908 (2013)
http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-4-10-1890


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  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,” Science254(5035), 1178–1181 (1991). [CrossRef] [PubMed]
  2. M. R. Hee, J. A. Izatt, E. A. Swanson, D. Huang, J. S. Schuman, C. P. Lin, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography of the human retina,” Arch. Ophthalmol.113(3), 325–332 (1995). [CrossRef] [PubMed]
  3. M. Wojtkowski, T. Bajraszewski, I. Gorczyńska, P. Targowski, A. Kowalczyk, W. Wasilewski, and C. Radzewicz, “Ophthalmic imaging by spectral optical coherence tomography,” Am. J. Ophthalmol.138(3), 412–419 (2004). [CrossRef] [PubMed]
  4. W. Wieser, B. R. Biedermann, T. Klein, C. M. Eigenwillig, and R. Huber, “Multi-Megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second,” Opt. Express18(14), 14685–14704 (2010). [CrossRef] [PubMed]
  5. B. Povazay, B. Hofer, C. Torti, B. Hermann, A. R. Tumlinson, M. Esmaeelpour, C. A. Egan, A. C. Bird, and W. Drexler, “Impact of enhanced resolution, speed and penetration on three-dimensional retinal optical coherence tomography,” Opt. Express17(5), 4134–4150 (2009). [CrossRef] [PubMed]
  6. J. M. Schmitt, S. H. Xiang, and K. M. Yung, “Speckle in Optical Coherence Tomography,” J. Biomed. Opt.4(1), 95–105 (1999). [CrossRef] [PubMed]
  7. M. Wojtkowski, B. Kaluzny, and R. J. Zawadzki, “New Directions in Ophthalmic Optical Coherence Tomography,” Optom. Vis. Sci.89(5), 524–542 (2012). [CrossRef] [PubMed]
  8. L. Fang, S. Li, Q. Nie, J. A. Izatt, C. A. Toth, and S. Farsiu, “Sparsity based denoising of spectral domain optical coherence tomography images,” Biomed. Opt. Express3(5), 927–942 (2012). [CrossRef] [PubMed]
  9. R. J. Zawadzki, S. S. Choi, A. R. Fuller, J. W. Evans, B. Hamann, and J. S. Werner, “Cellular resolution volumetric in vivo retinal imaging with adaptive optics-optical coherence tomography,” Opt. Express17(5), 4084–4094 (2009). [CrossRef] [PubMed]
  10. Y. Li, G. Gregori, B. L. Lam, and P. J. Rosenfeld, “Automatic montage of SD-OCT data sets,” Opt. Express19(27), 26239–26248 (2011). [CrossRef] [PubMed]
  11. S. Makita, F. Jaillon, M. Yamanari, M. Miura, and Y. Yasuno, “Comprehensive in vivo micro-vascular imaging of the human eye by dual-beam-scan Doppler optical coherence angiography,” Opt. Express19(2), 1271–1283 (2011). [CrossRef] [PubMed]
  12. R. D. Ferguson, D. X. Hammer, L. A. Paunescu, S. Beaton, and J. S. Schuman, “Tracking optical coherence tomography,” Opt. Lett.29(18), 2139–2141 (2004). [CrossRef] [PubMed]
  13. Y. K. Tao, S. Farsiu, and J. A. Izatt, “Interlaced spectrally encoded confocal scanning laser ophthalmoscopy and spectral domain optical coherence tomography,” Biomed. Opt. Express1(2), 431–440 (2010). [CrossRef] [PubMed]
  14. G. Gregori and P. J. Rosenfeld, “Using OCT Fundus Images to Evaluate the Performance of the Spectralis OCT Eye Tracking System (Poster 1032/A434),” in ARVO Annual Meeting 2011, (Fort Lauderdale, Florida, 2011).
  15. K. V. Vienola, B. Braaf, C. K. Sheehy, Q. Yang, P. Tiruveedhula, D. W. Arathorn, J. F. de Boer, and A. Roorda, “Real-time eye motion compensation for OCT imaging with tracking SLO,” Biomed. Opt. Express3(11), 2950–2963 (2012). [CrossRef] [PubMed]
  16. A. G. Podoleanu and R. B. Rosen, “Combinations of techniques in imaging the retina with high resolution,” Prog. Retin. Eye Res.27(4), 464–499 (2008). [CrossRef] [PubMed]
  17. S. Jiao, R. Knighton, X. Huang, G. Gregori, and C. Puliafito, “Simultaneous acquisition of sectional and fundus ophthalmic images with spectral-domain optical coherence tomography,” Opt. Express13(2), 444–452 (2005). [CrossRef] [PubMed]
  18. U. Schmidt-Erfurth, R. A. Leitgeb, S. Michels, B. Považay, S. Sacu, B. Hermann, C. Ahlers, H. Sattmann, C. Scholda, A. F. Fercher, and W. Drexler, “Three-Dimensional Ultrahigh-Resolution Optical Coherence Tomography of Macular Diseases,” Invest. Ophthalmol. Vis. Sci.46(9), 3393–3402 (2005). [CrossRef] [PubMed]
  19. A. G. Podoleanu, G. M. Dobre, R. G. Cucu, R. Rosen, P. Garcia, J. Nieto, D. Will, R. Gentile, T. Muldoon, J. Walsh, L. A. Yannuzzi, Y. Fisher, D. Orlock, R. Weitz, J. A. Rogers, S. Dunne, and A. Boxer, “Combined multiplanar optical coherence tomography and confocal scanning ophthalmoscopy,” J. Biomed. Opt.9(1), 86–93 (2004). [CrossRef] [PubMed]
  20. I. Gorczynska, V. J. Srinivasan, L. N. Vuong, R. W. S. Chen, J. J. Liu, E. Reichel, M. Wojtkowski, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Projection OCT fundus imaging for visualising outer retinal pathology in non-exudative age-related macular degeneration,” Br. J. Ophthalmol.93(5), 603–609 (2009). [CrossRef] [PubMed]
  21. T. Agawa, M. Miura, Y. Ikuno, S. Makita, T. Fabritius, T. Iwasaki, H. Goto, K. Nishida, and Y. Yasuno, “Choroidal thickness measurement in healthy Japanese subjects by three-dimensional high-penetration optical coherence tomography,” Graefes Arch. Clin. Exp. Ophthalmol.249(10), 1485–1492 (2011). [CrossRef] [PubMed]
  22. B. Povazay, B. Hermann, B. Hofer, V. Kajić, E. Simpson, T. Bridgford, and W. Drexler, “Wide-Field Optical Coherence Tomography of the Choroid In Vivo,” Invest. Ophthalmol. Vis. Sci.50(4), 1856–1863 (2008). [CrossRef] [PubMed]
  23. T. Klein, W. Wieser, C. M. Eigenwillig, B. R. Biedermann, and R. Huber, “Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser,” Opt. Express19(4), 3044–3062 (2011). [CrossRef] [PubMed]
  24. A. A. Ellabban, A. Tsujikawa, A. Matsumoto, K. Yamashiro, A. Oishi, S. Ooto, I. Nakata, Y. Akagi-Kurashige, M. Miyake, H. S. Elnahas, T. M. Radwan, K. A. Zaky, and N. Yoshimura, “Three-Dimensional Tomographic Features of Dome-Shaped Macula by Swept-Source Optical Coherence Tomography,” Am. J. Ophthalmol.155(2), 320–328, e2 (2013). [CrossRef] [PubMed]
  25. S. Makita, Y. Hong, M. Yamanari, T. Yatagai, and Y. Yasuno, “Optical coherence angiography,” Opt. Express14(17), 7821–7840 (2006). [CrossRef] [PubMed]
  26. T. Schmoll, C. Kolbitsch, and R. A. Leitgeb, “In vivo functional retinal optical coherence tomography,” J. Biomed. Opt.15(4), 041513 (2010). [CrossRef] [PubMed]
  27. L. An, T. T. Shen, and R. K. Wang, “Using ultrahigh sensitive optical microangiography to achieve comprehensive depth resolved microvasculature mapping for human retina,” J. Biomed. Opt.16(10), 106013 (2011). [CrossRef] [PubMed]
  28. T. Torzicky, S. Marschall, M. Pircher, B. Baumann, M. Bonesi, S. Zotter, E. Götzinger, W. Trasischker, T. Klein, W. Wieser, B. Biedermann, R. Huber, P. Andersen, and C. K. Hitzenberger, “Retinal polarization-sensitive optical coherence tomography at 1060 nm with 350 kHz A-scan rate using an Fourier domain mode locked laser,” J. Biomed. Opt.18(2), 026008 (2013). [CrossRef] [PubMed]
  29. B. White, M. Pierce, N. Nassif, B. Cense, B. Park, G. Tearney, B. Bouma, T. Chen, and J. de Boer, “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography,” Opt. Express11(25), 3490–3497 (2003). [CrossRef] [PubMed]
  30. C. Blatter, T. Klein, B. Grajciar, T. Schmoll, W. Wieser, R. Andre, R. Huber, and R. A. Leitgeb, “Ultrahigh-speed non-invasive widefield angiography,” J. Biomed. Opt.17(7), 070505 (2012). [CrossRef] [PubMed]
  31. B. Potsaid, I. Gorczynska, V. J. Srinivasan, Y. L. Chen, J. Jiang, A. Cable, and J. G. Fujimoto, “Ultrahigh speed Spectral / Fourier domain OCT ophthalmic imaging at 70,000 to 312,500 axial scans per second,” Opt. Express16(19), 15149–15169 (2008). [CrossRef] [PubMed]
  32. L. An, P. Li, T. T. Shen, and R. Wang, “High speed spectral domain optical coherence tomography for retinal imaging at 500,000 A‑lines per second,” Biomed. Opt. Express2(10), 2770–2783 (2011). [CrossRef] [PubMed]
  33. D. Choi, H. Hiro-Oka, H. Furukawa, R. Yoshimura, M. Nakanishi, K. Shimizu, and K. Ohbayashi, “Fourier domain optical coherence tomography using optical demultiplexers imaging at 60,000,000 lines/s,” Opt. Lett.33(12), 1318–1320 (2008). [CrossRef] [PubMed]
  34. D. H. Choi, H. Hiro-Oka, K. Shimizu, and K. Ohbayashi, “Spectral domain optical coherence tomography of multi-MHz A-scan rates at 1310 nm range and real-time 4D-display up to 41 volumes/second,” Biomed. Opt. Express3(12), 3067–3086 (2012). [CrossRef] [PubMed]
  35. S. Yun, G. Tearney, J. de Boer, N. Iftimia, and B. Bouma, “High-speed optical frequency-domain imaging,” Opt. Express11(22), 2953–2963 (2003). [CrossRef] [PubMed]
  36. R. Huber, M. Wojtkowski, K. Taira, J. G. Fujimoto, and K. Hsu, “Amplified, frequency swept lasers for frequency domain reflectometry and OCT imaging: design and scaling principles,” Opt. Express13(9), 3513–3528 (2005). [CrossRef] [PubMed]
  37. B. Potsaid, B. Baumann, D. Huang, S. Barry, A. E. Cable, J. S. Schuman, J. S. Duker, and J. G. Fujimoto, “Ultrahigh speed 1050nm swept source/Fourier domain OCT retinal and anterior segment imaging at 100,000 to 400,000 axial scans per second,” Opt. Express18(19), 20029–20048 (2010). [CrossRef] [PubMed]
  38. W.-Y. Oh, B. J. Vakoc, M. Shishkov, G. J. Tearney, and B. E. Bouma, “>400 kHz repetition rate wavelength-swept laser and application to high-speed optical frequency domain imaging,” Opt. Lett.35(17), 2919–2921 (2010). [CrossRef] [PubMed]
  39. M. P. Minneman, J. Ensher, M. Crawford, and D. Derickson, “All-semiconductor high-speed akinetic swept-source for OCT,” in Optical Sensors and Biophotonics III, Proc. SPIE 8311 (SPIE, 2011), 831116–831116.
  40. B. Potsaid, V. Jayaraman, J. G. Fujimoto, J. Jiang, P. J. S. Heim, and A. E. Cable, “MEMS tunable VCSEL light source for ultrahigh speed 60kHz - 1MHz axial scan rate and long range centimeter class OCT imaging,” in (SPIE, 2012), 82130M.
  41. R. Huber, M. Wojtkowski, and J. G. Fujimoto, “Fourier Domain Mode Locking (FDML): A new laser operating regime and applications for optical coherence tomography,” Opt. Express14(8), 3225–3237 (2006). [CrossRef] [PubMed]
  42. S. Moon and D. Y. Kim, “Ultra-high-speed optical coherence tomography with a stretched pulse supercontinuum source,” Opt. Express14(24), 11575–11584 (2006). [CrossRef] [PubMed]
  43. C. M. Eigenwillig, B. R. Biedermann, W. Wieser, and R. Huber, “Wavelength swept amplified spontaneous emission source,” Opt. Express17(21), 18794–18807 (2009). [CrossRef] [PubMed]
  44. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, J. Jiang, J. G. Fujimoto, and A. E. Cable, “High-precision, high-accuracy ultralong-range swept-source optical coherence tomography using vertical cavity surface emitting laser light source,” Opt. Lett.38(5), 673–675 (2013). [CrossRef] [PubMed]
  45. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, C. D. Lu, J. Jiang, A. E. Cable, J. S. Duker, and J. G. Fujimoto, “Retinal, anterior segment and full eye imaging using ultrahigh speed swept source OCT with vertical-cavity surface emitting lasers,” Biomed. Opt. Express3(11), 2733–2751 (2012). [CrossRef] [PubMed]
  46. B. R. Biedermann, W. Wieser, C. M. Eigenwillig, T. Klein, and R. Huber, “Dispersion, coherence and noise of Fourier domain mode locked lasers,” Opt. Express17(12), 9947–9961 (2009). [CrossRef] [PubMed]
  47. S. Todor, B. Biedermann, W. Wieser, R. Huber, and C. Jirauschek, “Instantaneous lineshape analysis of Fourier domain mode-locked lasers,” Opt. Express19(9), 8802–8807 (2011). [CrossRef] [PubMed]
  48. D. C. Adler, W. Wieser, F. Trepanier, J. M. Schmitt, and R. A. Huber, “Extended coherence length Fourier domain mode locked lasers at 1310 nm,” Opt. Express19(21), 20930–20939 (2011). [CrossRef] [PubMed]
  49. S. Marschall, T. Klein, W. Wieser, B. R. Biedermann, K. Hsu, K. P. Hansen, B. Sumpf, K. H. Hasler, G. Erbert, O. B. Jensen, C. Pedersen, R. Huber, and P. E. Andersen, “Fourier domain mode-locked swept source at 1050 nm based on a tapered amplifier,” Opt. Express18(15), 15820–15831 (2010). [CrossRef] [PubMed]
  50. K. Goda, A. Fard, O. Malik, G. Fu, A. Quach, and B. Jalali, “High-throughput optical coherence tomography at 800 nm,” Opt. Express20(18), 19612–19617 (2012). [CrossRef] [PubMed]
  51. C. M. Eigenwillig, T. Klein, W. Wieser, B. R. Biedermann, and R. Huber, “Wavelength swept amplified spontaneous emission source for high speed retinal optical coherence tomography at 1060 nm,” J Biophotonics4(7-8), 552–558 (2011). [CrossRef] [PubMed]
  52. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett.28(21), 2067–2069 (2003). [CrossRef] [PubMed]
  53. A. M. Rollins and J. A. Izatt, “Optimal interferometer designs for optical coherence tomography,” Opt. Lett.24(21), 1484–1486 (1999). [CrossRef] [PubMed]
  54. B. Považay, B. Hermann, A. Unterhuber, B. Hofer, H. Sattmann, F. Zeiler, J. E. Morgan, C. Falkner-Radler, C. Glittenberg, S. Blinder, and W. Drexler, “Three-dimensional optical coherence tomography at 1050 nm versus 800 nm in retinal pathologies: enhanced performance and choroidal penetration in cataract patients,” J. Biomed. Opt.12(4), 041211 (2007). [CrossRef] [PubMed]
  55. M. Szkulmowski and M. Wojtkowski, “Averaging techniques for OCT imaging,” Opt. Express21(8), 9757–9773 (2013). [CrossRef] [PubMed]
  56. N. Suehira, S. Ooto, M. Hangai, K. Matsumoto, N. Tomatsu, T. Yuasa, K. Yamada, and N. Yoshimura, “Three-beam spectral-domain optical coherence tomography for retinal imaging,” J. Biomed. Opt.17(10), 106001 (2012). [CrossRef] [PubMed]
  57. T. Klein, R. André, W. Wieser, T. Pfeiffer, and R. Huber, “Joint aperture detection for speckle reduction and increased collection efficiency in ophthalmic MHz OCT,” Biomed. Opt. Express4(4), 619–634 (2013). [CrossRef] [PubMed]
  58. R. Huber, D. C. Adler, V. J. Srinivasan, and J. G. Fujimoto, “Fourier domain mode locking at 1050 nm for ultra-high-speed optical coherence tomography of the human retina at 236,000 axial scans per second,” Opt. Lett.32(14), 2049–2051 (2007). [CrossRef] [PubMed]
  59. V. J. Srinivasan, D. C. Adler, Y. L. Chen, I. Gorczynska, R. Huber, J. S. Duker, J. S. Schuman, and J. G. Fujimoto, “Ultrahigh-Speed Optical Coherence Tomography for Three-Dimensional and En Face Imaging of the Retina and Optic Nerve Head,” Invest. Ophthalmol. Vis. Sci.49(11), 5103–5110 (2008). [CrossRef] [PubMed]
  60. W. Wieser, G. Palte, C. M. Eigenwillig, B. R. Biedermann, T. Pfeiffer, and R. Huber, “Chromatic polarization effects of swept waveforms in FDML lasers and fiber spools,” Opt. Express20(9), 9819–9832 (2012). [CrossRef] [PubMed]
  61. W. Wieser, T. Klein, D. C. Adler, F. Trépanier, C. M. Eigenwillig, S. Karpf, J. M. Schmitt, and R. Huber, “Extended coherence length megahertz FDML and its application for anterior segment imaging,” Biomed. Opt. Express3(10), 2647–2657 (2012). [CrossRef] [PubMed]
  62. R. Huber, D. C. Adler, and J. G. Fujimoto, “Buffered Fourier domain mode locking: unidirectional swept laser sources for optical coherence tomography imaging at 370,000 lines/s,” Opt. Lett.31(20), 2975–2977 (2006). [CrossRef] [PubMed]
  63. T. Klein, W. Wieser, R. Andre, T. Pfeiffer, C. M. Eigenwillig, and R. Huber, “Multi-MHz FDML OCT: snapshot retinal imaging at 6.7 million axial-scans per second,” in (SPIE, 2012), 82131E–82136.
  64. R. J. Zawadzki, A. R. Fuller, S. S. Choi, D. F. Wiley, B. Hamann, and J. S. Werner, “Correction of motion artifacts and scanning beam distortions in 3D ophthalmic optical coherence tomography imaging,” Proc. SPIE64260, 64260 (2007).
  65. A. Podoleanu, I. Charalambous, L. Plesea, A. Dogariu, and R. Rosen, “Correction of distortions in optical coherence tomography imaging of the eye,” Phys. Med. Biol.49(7), 1277–1294 (2004). [CrossRef] [PubMed]
  66. Y. Chen, D. M. de Bruin, C. Kerbage, and J. F. de Boer, “Spectrally balanced detection for optical frequency domain imaging,” Opt. Express15(25), 16390–16399 (2007). [CrossRef] [PubMed]
  67. B. Braaf, K. A. Vermeer, V. A. D. P. Sicam, E. van Zeeburg, J. C. van Meurs, and J. F. de Boer, “Phase-stabilized optical frequency domain imaging at 1-µm for the measurement of blood flow in the human choroid,” Opt. Express19(21), 20886–20903 (2011). [CrossRef] [PubMed]
  68. S. Martinez-Conde, S. L. Macknik, X. G. Troncoso, and D. H. Hubel, “Microsaccades: a neurophysiological analysis,” Trends Neurosci.32(9), 463–475 (2009). [CrossRef] [PubMed]
  69. M. F. Kraus, B. Potsaid, M. A. Mayer, R. Bock, B. Baumann, J. J. Liu, J. Hornegger, and J. G. Fujimoto, “Motion correction in optical coherence tomography volumes on a per A-scan basis using orthogonal scan patterns,” Biomed. Opt. Express3(6), 1182–1199 (2012). [CrossRef] [PubMed]
  70. H. C. Hendargo, R. Estrada, S. J. Chiu, C. Tomasi, S. Farsiu, and J. A. Izatt, “Automated non-rigid registration and mosaicing for robust imaging of distinct retinal capillary beds using speckle variance optical coherence tomography,” Biomed. Opt. Express4(6), 803–821 (2013). [CrossRef] [PubMed]
  71. S. Ricco, M. Chen, H. Ishikawa, G. Wollstein, and J. Schuman, “Correcting Motion Artifacts in Retinal Spectral Domain Optical Coherence Tomography via Image Registration,” in Medical Image Computing and Computer-Assisted Intervention – MICCAI 2009, G.-Z. Yang, D. Hawkes, D. Rueckert, A. Noble, and C. Taylor, eds. (Springer Berlin Heidelberg, 2009), pp. 100–107.
  72. M. Niemeijer, M. K. Garvin, K. Lee, B. van Ginneken, M. D. Abràmoff, and M. Sonka, “Registration of 3D spectral OCT volumes using 3D SIFT feature point matching ” in Proc. of SPIE, (2009), p. 72591I.
  73. S. Lee, M. Young, M. V. Sarunic, and M. F. Beg, “End-to-end pipeline for spectral domain optical coherence tomography and morphometric analysis of human optic nerve head,” J. Med. Biol. Eng.31, 9 (2010).
  74. B. Antony, M. D. Abràmoff, L. Tang, W. D. Ramdas, J. R. Vingerling, N. M. Jansonius, K. Lee, Y. H. Kwon, M. Sonka, and M. K. Garvin, “Automated 3-D method for the correction of axial artifacts in spectral-domain optical coherence tomography images,” Biomed. Opt. Express2(8), 2403–2416 (2011). [CrossRef] [PubMed]
  75. P. Thévenaz, U. E. Ruttimann, and M. Unser, “A pyramid approach to subpixel registration based on intensity,” IEEE Trans. Image Process.7(1), 27–41 (1998). [CrossRef] [PubMed]
  76. J. Schindelin, I. Arganda-Carreras, E. Frise, V. Kaynig, M. Longair, T. Pietzsch, S. Preibisch, C. Rueden, S. Saalfeld, B. Schmid, J.-Y. Tinevez, D. J. White, V. Hartenstein, K. Eliceiri, P. Tomancak, and A. Cardona, “Fiji: an open-source platform for biological-image analysis,” Nat. Methods9(7), 676–682 (2012). [CrossRef] [PubMed]
  77. M. Szkulmowski, I. Gorczynska, D. Szlag, M. Sylwestrzak, A. Kowalczyk, and M. Wojtkowski, “Efficient reduction of speckle noise in Optical Coherence Tomography,” Opt. Express20(2), 1337–1359 (2012). [CrossRef] [PubMed]
  78. A. S. Neubauer, M. Kernt, C. Haritoglou, S. G. Priglinger, A. Kampik, and M. W. Ulbig, “Nonmydriatic screening for diabetic retinopathy by ultra-widefield scanning laser ophthalmoscopy (Optomap),” Graefes Arch. Clin. Exp. Ophthalmol.246(2), 229–235 (2008). [CrossRef] [PubMed]
  79. Y. K. Tao, J. P. Ehlers, C. A. Toth, and J. A. Izatt, “Intraoperative spectral domain optical coherence tomography for vitreoretinal surgery,” Opt. Lett.35(20), 3315–3317 (2010). [CrossRef] [PubMed]
  80. J. U. Kang, Y. Huang, K. Zhang, Z. Ibrahim, J. Cha, W. P. A. Lee, G. Brandacher, and P. L. Gehlbach, “Real-time three-dimensional Fourier-domain optical coherence tomography video image guided microsurgeries,” J. Biomed. Opt.17(8), 081403 (2012). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Supplementary Material


» Media 1: MOV (47421 KB)     
» Media 2: MOV (57376 KB)     
» Media 3: MOV (5873 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited