Graphics Processing Unit-Based Ultrahigh Speed Real-Time Fourier Domain Optical Coherence Tomography

This paper describes a series of novel graphics processing unit (GPU)-based image reconstruction and visualization methods that we developed for realizing ultrahigh speed, real-time Fourier domain optical coherence tomography (FD-OCT). Several GPU-based algorithms including high-speed linear/cubic interpolation, gridding-based nonuniform fast Fourier transform (NUFFT), numerical dispersion compensation, and multi-GPU implementation were developed to improve the point-spread function, SNR roll-off, and stability of the system. Full-range complex-conjugate-free FD-OCT was also implemented on the GPU architecture to double the imaging range and to improve the SNR. The maximum processing speed of <formula formulatype="inline"><tex Notation="TeX">$>$</tex></formula>3 Gigavoxels/s ( <formula formulatype="inline"><tex Notation="TeX">$>$</tex></formula>6 Mega-A-scans/s of 1024-pixel FD-OCT) was achieved using the latest GPU module from NVIDIA. The GPU-based volume rendering enabled real-time four-dimensional (4-D), i.e., 3-D + time, FD-OCT imaging, and a 5 volume/s 4-D FD-OCT system was demonstrated via in vivo tissue imaging. These GPU technologies were highly effective in circumventing the well-known imaging reconstruction and visualization bottlenecks existing among current ultrahigh speed FD-OCT systems and could open the door to realize interventional OCT imaging for applications in guiding microsurgery.

[1]  T. Yatagai,et al.  Simultaneous B-M-mode scanning method for real-time full-range Fourier domain optical coherence tomography. , 2006, Applied optics.

[2]  J. Duker,et al.  Ultrahigh-resolution, high-speed, Fourier domain optical coherence tomography and methods for dispersion compensation. , 2004, Optics express.

[3]  Shoude Chang,et al.  Swept Source Optical Coherence Tomography with Nonuniform Frequency Domain Sampling , 2008 .

[4]  Jing Xu,et al.  Performance and scalability of Fourier domain optical coherence tomography acceleration using graphics processing units. , 2011, Applied optics.

[5]  Kate Sugden,et al.  Processing and rendering of Fourier domain optical coherence tomography images at a line rate over 524 kHz using a graphics processing unit. , 2011, Journal of biomedical optics.

[6]  J. Fujimoto,et al.  Ultrahigh speed spectral / Fourier domain OCT ophthalmic imaging at 70,000 to 312,500 axial scans per second. , 2008, Optics express.

[7]  Y. Yasuno,et al.  Full-range, high-speed, high-resolution 1 microm spectral-domain optical coherence tomography using BM-scan for volumetric imaging of the human posterior eye. , 2008, Optics express.

[8]  K. Cleary,et al.  Image Guided Interventions. , 2020, Biomedizinische Technik. Biomedical engineering.

[9]  Zhilin Hu,et al.  Fourier domain optical coherence tomography with a linear-in-wavenumber spectrometer. , 2007, Optics letters.

[10]  Kang Zhang,et al.  Real-time 4D signal processing and visualization using graphics processing unit on a regular nonlinear-k Fourier-domain OCT system , 2010, Optics express.

[11]  J. Duker,et al.  Ultrahigh speed 1050nm swept source/Fourier domain OCT retinal and anterior segment imaging at 100,000 to 400,000 axial scans per second. , 2010, Optics express.

[12]  R. Huber,et al.  K-space linear Fourier domain mode locked laser and applications for optical coherence tomography. , 2008, Optics express.

[13]  Gangjun Liu,et al.  Real-time polarization-sensitive optical coherence tomography data processing with parallel computing. , 2009, Applied optics.

[14]  K. Cleary,et al.  Image-guided interventions : technology and applications , 2008 .

[15]  Uretz J. Oliphant,et al.  Optical Coherence Tomography: The Intraoperative Assessment of Lymph Nodes in Breast Cancer , 2010, IEEE Engineering in Medicine and Biology Magazine.

[16]  J. Fujimoto,et al.  Intraoperative assessment of microsurgery with three-dimensional optical coherence tomography. , 1998, Radiology.

[17]  Yuuki Watanabe,et al.  Real-time processing for full-range Fourier-domain optical-coherence tomography with zero-filling interpolation using multiple graphic processing units. , 2010, Applied optics.

[18]  Yuuki Watanabe,et al.  Real-time display on Fourier domain optical coherence tomography system using a graphics processing unit. , 2009, Journal of biomedical optics.

[19]  Kang Zhang,et al.  Real-time intraoperative 4D full-range FD-OCT based on the dual graphics processing units architecture for microsurgery guidance , 2011, Biomedical optics express.

[20]  Leslie Greengard,et al.  Accelerating the Nonuniform Fast Fourier Transform , 2004, SIAM Rev..

[21]  Kang Zhang,et al.  Real-time numerical dispersion compensation using graphics processing unit for Fourier-domain optical coherence tomography , 2011 .

[22]  Ruikang K. Wang,et al.  Use of a scanner to modulate spatial interferograms for in vivo full-range Fourier-domain optical coherence tomography. , 2007, Optics letters.

[23]  Sébastien Vergnole,et al.  Experimental validation of an optimized signal processing method to handle non-linearity in swept-source optical coherence tomography. , 2010, Optics express.

[24]  Kang Zhang,et al.  Graphics processing unit accelerated non-uniform fast Fourier transform for ultrahigh-speed, real-time Fourier-domain OCT , 2010, Optics express.

[25]  Daniel X Hammer,et al.  Real-time processing for Fourier domain optical coherence tomography using a field programmable gate array. , 2008, The Review of scientific instruments.

[26]  Adrian Bradu,et al.  Real-time resampling in Fourier domain optical coherence tomography using a graphics processing unit. , 2010, Journal of biomedical optics.

[27]  Wolfgang Wieser,et al.  Multi-megahertz OCT: High quality 3D imaging at 20 million A-scans and 4.5 GVoxels per second. , 2010, Optics express.

[28]  Marc Levoy,et al.  Display of surfaces from volume data , 1988, IEEE Computer Graphics and Applications.

[29]  J. Fujimoto,et al.  Buffered Fourier domain mode locking: Unidirectional swept laser sources for optical coherence tomography imaging at 370,000 lines/s. , 2006, Optics letters.

[30]  Adrien E. Desjardins,et al.  Real-Time FPGA Processing for High-Speed Optical Frequency Domain Imaging , 2009, IEEE Transactions on Medical Imaging.

[31]  Kai Wang,et al.  Development of a non-uniform discrete Fourier transform based high speed spectral domain optical coherence tomography system. , 2009, Optics express.

[32]  R. Leitgeb,et al.  Complex ambiguity-free Fourier domain optical coherence tomography through transverse scanning. , 2007, Optics letters.

[33]  R. Huber,et al.  Megahertz OCT for ultrawide-field retinal imaging with a 1050 nm Fourier domain mode-locked laser. , 2011, Optics express.

[34]  Michael Pircher,et al.  Full range complex spectral domain optical coherence tomography without additional phase shifters. , 2007, Optics express.

[35]  B. Vakoc,et al.  >400 kHz repetition rate wavelength-swept laser and application to high-speed optical frequency domain imaging. , 2010, Optics letters.

[36]  J. Fujimoto,et al.  Three-dimensional endomicroscopy using optical coherence tomography , 2007 .

[37]  J. Schuman,et al.  Optical coherence tomography. , 2000, Science.

[38]  Sébastien Vergnole,et al.  Artifact removal in Fourier-domain optical coherence tomography with a piezoelectric fiber stretcher. , 2008, Optics letters.