Graphics processing unit accelerated non-uniform fast Fourier transform for ultrahigh-speed, real-time Fourier-domain OCT

We implemented fast Gaussian gridding (FGG)-based non-uniform fast Fourier transform (NUFFT) on the graphics processing unit (GPU) architecture for ultrahigh-speed, real-time Fourier-domain optical coherence tomography (FD-OCT). The Vandermonde matrix-based non-uniform discrete Fourier transform (NUDFT) as well as the linear/cubic interpolation with fast Fourier transform (InFFT) methods are also implemented on GPU to compare their performance in terms of image quality and processing speed. The GPU accelerated InFFT/NUDFT/NUFFT methods are applied to process both the standard half-range FD-OCT and complex full-range FD-OCT (C-FD-OCT). GPU-NUFFT provides an accurate approximation to GPU-NUDFT in terms of image quality, but offers >10 times higher processing speed. Compared with the GPU-InFFT methods, GPU-NUFFT has improved sensitivity roll-off, higher local signal-to-noise ratio and immunity to side-lobe artifacts caused by the interpolation error. Using a high speed CMOS line-scan camera, we demonstrated the real-time processing and display of GPU-NUFFT-based C-FD-OCT at a camera-limited rate of 122 k line/s (1024 pixel/A-scan).

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

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

[3]  M. Wojtkowski,et al.  Ultra high-speed swept source OCT imaging of the anterior segment of human eye at 200 kHz with adjustable imaging range. , 2009, Optics express.

[4]  Michael W. Jenkins,et al.  4D embryonic cardiography using gated optical coherence tomography. , 2006, Optics express.

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

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

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

[8]  Eva Lankenau,et al.  Optical coherence tomography with online visualization of more than seven rendered volumes per second. , 2010, Journal of biomedical optics.

[9]  Cha-Min Tang,et al.  Optical coherence tomography guided neurosurgical procedures in small rodents , 2009, Journal of Neuroscience Methods.

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

[11]  Peter Koch,et al.  Using nonequispaced fast Fourier transformation to process optical coherence tomography signals , 2009, European Conference on Biomedical Optics.

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

[13]  C. Dorrer,et al.  Spectral resolution and sampling issues in Fourier-transform spectral interferometry , 2000 .

[14]  David L Wilson,et al.  High temporal resolution OCT using image-based retrospective gating. , 2009, Optics express.

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

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

[17]  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.

[18]  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.

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

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

[21]  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.

[22]  David L Wilson,et al.  Denoising and 4D visualization of OCT images. , 2008, Optics express.

[23]  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.

[24]  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.

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

[26]  Iwona Gorczynska,et al.  Anterior segment imaging with Spectral OCT system using a high-speed CMOS camera. , 2009, Optics express.

[27]  Kang Zhang,et al.  Real-time complex optical coherence tomography using graphics processing unit for surgical intervention , 2010, 2010 IEEE Photinic Society's 23rd Annual Meeting.

[28]  Kang Zhang,et al.  A Surface Topology and Motion Compensation System for Microsurgery Guidance and Intervention Based on Common-Path Optical Coherence Tomography , 2009, IEEE Transactions on Biomedical Engineering.

[29]  A. Low,et al.  Technology Insight: optical coherence tomography—current status and future development , 2006, Nature Clinical Practice Cardiovascular Medicine.

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

[31]  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.

[32]  Lin An,et al.  Ultra-high speed full range complex spectral domain optical coherence tomography for volumetric imaging at 140,000 A scans per second , 2010, BiOS.

[33]  Jin U. Kang,et al.  All-fiber common-path optical coherence tomography: sensitivity optimization and system analysis , 2005, IEEE Journal of Selected Topics in Quantum Electronics.

[34]  James G. Fujimoto,et al.  Surgical Guidance and Intervention , 2001 .

[35]  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.

[36]  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.

[37]  Chul Gyu Song,et al.  Endoscopic Functional Fourier Domain Common-Path Optical Coherence Tomography for Microsurgery , 2010, IEEE Journal of Selected Topics in Quantum Electronics.

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