Translator Disclaimer
1 July 2008 High-speed processing architecture for spectral-domain optical coherence microscopy
Author Affiliations +
Abstract
Optical coherence microscopy (OCM) is an interferometric technique that combines principles of confocal microscopy and optical coherence tomography (OCT) to obtain high-resolution en face images. Axial and lateral resolutions of several microns can be achieved using OCM depending on the numerical aperture (NA) of the objective and sample properties. We address the computational complexity that is inherent in spectral-domain OCM systems that limits its real-time capability as a microscope. An architecture that will improve the efficiency of the computation involved is presented. Currently, spectral-domain OCM images are obtained by individually taking the Fourier transform of each axial scan in cross-sectional frames and computationally slicing them to generate en face images. The real-time architecture presented here relies on the fact that only one Fourier domain point of a given axial scan needs to be computed rather than computing all the Fourier domain points, which can frequently require a significant amount of time to compute. This new realization has been shown to reduce the processing time to obtain the en face OCM images by a factor of 30.

Optical coherence tomography (OCT) is a biomedical imaging technique that is gaining wide interest as it performs noninvasive high-resolution cross-sectional imaging of internal microstructures in materials and biological systems by measuring backscattered light. Although OCT has achieved superior resolution in comparison to many other biomedical imaging techniques such as computed tomography, magnetic resonance imaging, and ultrasound, the resolution of OCT is not always sufficient enough to view living cells ( 10μm in size). The idea of combining principles from confocal microscopy with OCT has improved the resolution of OCT to view living cells. This technique is called optical coherence microscopy (OCM).1, 2, 3, 4, 5 and takes advantage of very high NA (> 0.9) focused illumination with spatially filtered as well as coherence-gated detection.

A simplified OCM system involves an interferometric setup as shown in Fig. 1 . A source signal is split into two arms, namely sample and reference arms, and interfered as backscattered signals return from the two arms. High-speed instrumentation that allows en face OCM imaging has been of great interest since real-time imaging of live specimen is often desirable. A combined acquisition of optical coherence and two-photon-excited (TPE) fluorescence microscopy images was reported to have an acquisition speed of 2.8s .2 A piezoelectric translator was used in Ref. 2 at the reference arm to displace the reference mirror, and XY raster scanning was performed using mirrors mounted on galvanometers in the sample arm. Time-domain (TD) OCM systems have been reported to have signal-to-noise ratios over 100dB , often utilizing high-power delivery into the sample. This is not often desirable, as high powers (> 10mW) can damage the sample, particularly when using the higher NA objectives for OCM. Also, high-speed imaging can compromise the sensitivity of TD-OCM systems. High sensitivity, high-speed acquisition, and the availability of phase and spectroscopic information have directed research investigations toward spectral-domain OCM.6, 7 Spectral-domain OCM systems have been reported to have sensitivity over 100dB with only a few milliwatts at the sample.8 In spectral-domain OCM, the interfered signal is detected using a spectral-domain detection system that is comprised of a collimator, grating, and a linear CCD array. The signal, which contains depth-resolved amplitude information, is processed and displayed using a personal computer. Moreover, no movement in the reference arm is needed and galvanometer-controlled mirrors are often used to rapidly scan the beam over the sample.

Fig. 1

Representative setup for optical coherence microscopy.

044013_1_036804jbo1.jpg

In spectral-domain OCM, as the sample is scanned, conventionally raw cross-sectional frames are first obtained. Then, a one-dimensional Fourier transform (FT) on each axial line is performed on each raw cross-sectional frame and the processed cross-sectional images are obtained. In order to obtain en face OCM images, these processed cross-sectional images are computationally sliced at the location in the images that correspond to the focal plane of the objective. Thus, the number of Fourier domain points that need to be computed to obtain the en face OCM image is very large (over a million points for a typical image). Computing a large number of Fourier domain points is therefore time-consuming and greatly impedes the real-time capability of an SD-OCM microscope. An illustration of the conventional architecture is shown in Fig. 2 .

Fig. 2

Conventional architecture to obtain en face OCM images.

044013_1_036804jbo2.jpg

A new real-time architecture that improves the efficiency of the computation involved in obtaining the en face OCM image is needed. To address the aforementioned computational complexity in the conventional architecture, we investigated if all the Fourier domain points corresponding to an axial line need to be computed in a raw cross-sectional frame to obtain the en face OCM image. Since en face OCM imaging is at an en face plane (i.e., z=p ), we demonstrate that only the one Fourier domain point at that plane [i.e., G(p) ] is needed. This plane usually corresponds to the focal plane of the objective lens. The new real-time architecture is based on this idea and is illustrated schematically in Fig. 3 .

Fig. 3

Real-time architecture to obtain en face OCM images.

044013_1_036804jbo3.jpg

The discrete Fourier transform of a discrete signal s(n) of width size 2048 is given by

Eq. 1

G(k)=n=0N1s(n)e2πiNkn,whereN=2048.
Conventionally, Fourier domain points from G(0) through G(2047) are computed for each axial line in a raw cross-sectional frame. In the real-time architecture, the Fourier domain point that corresponds to the focal plane is computed as shown in Eq. 2:

Eq. 2

{G(0)=G(1)=G(2)=}Fourierdomainpointsoutsidethefocalplanedonotneedtobecomputedsincethesedonotcontributetotheenfaceplanefeatures.G(p)=s(0)+s(1)ej2πp2048+s(2)ej4πp2048++s(2047)ej4094πp2048.{G(1024)=G(2047)=}Fourierdomainpointsoutsidethefocalplanedonotneedtobecomputedsincethesedonotcontributetotheenfaceplanefeatures.
The exponential coefficients in Eq. 2 can be precomputed and stored rather than being computed them for each individual axial line. The primary disadvantage of using this equation and approach is that the spectral information contained within the fully acquired signal is not obtained, as all the Fourier domain points are not computed.

Due to dispersion of the light during its pass through the various optical components, the OCM signal received at the detector is skewed. That is, the frequency associated with a certain depth falls onto a pixel that is associated with a different depth on the CCD array. Therefore, the expression for G(p) should be modified using a resampling array to digitally compensate for the dispersion.9 The resampling array is given by

Eq. 3

in=n+β2(nNωctr)2+β3(nNωctr)3,
where N=axial array size (e.g., N=2048 ), β2=second -order dispersion, β3=third -order dispersion, ωctr=center wavelength. Therefore, Eq. 2 is modified to have the final expression

Eq. 4

G(p)=s(0)ej2πp*i02048+s(1)ej2πp*i12048+s(2)ej2πp*i22048++s(2047)ej2πp*i20472048.

From Eq. 4, it can be seen that the processor has to perform several floating-point multiplications and additions to compute the Fourier domain point G(p) . Multiplying a single pair of floating-point values significantly contributes to the complexity of this computation when real-time en face OCM images need to be displayed. For this reason, the “SSE” capability of the processors is utilized here. SSE stands for “streaming SIMD (single-instruction, multiple-data) extensions.” SSE Instructions are a set of instructions that are supported by processors such as Intel® Pentium® III, IV, etc. Using SSE, four pairs of floating-point (32-bit) values can be multiplied in a single instruction, as illustrated in Fig. 4 . This feature of the processor is added in the new real-time architecture to compute the Fourier domain point G(p) . This further reduces the complexity of the computation involved.

Fig. 4

Differences between conventional and SSE multiplications.

044013_1_036804jbo4.jpg

This real-time architecture was implemented on our integrated microscope that enables simultaneous acquisition of OCM and multiphoton fluorescence images.4 An en face OCM image of TiO2 particles ( 5-μm diameter) in silicone was acquired and is shown in Fig. 5 . Computation times were measured using a 3.66-GHz Pentium® IV processor. The conventional architecture required 0.225ms to process one pixel of a 512×512-pixel en face OCM image, corresponding to a total processing time of 59s . Using the expression in Eq. 4, the real-time architecture required 0.0078ms to process one pixel of a 512×512-pixel en face OCM image, corresponding to a total processing time of 2s . Thus, there was an increase in the image processing speed by a factor of 30 .

Fig. 5

SD-OCM image (512×512pixels) of TiO2 particles ( 5-μm diameter) suspended in silicone acquired using the real-time architecture. The SD-OCM image acquired using the conventional architecture is identical to the image shown.

044013_1_036804jbo5.jpg

Since the real-time architecture uses the same reconstruction equation as the conventional architecture, except that the real-time architecture computes only one point at the focal plane, there was no reduction in the signal-to-noise ratio, contrast, or resolution. An image of the TiO2 particles captured using the conventional architecture showed no observable differences from the image shown in Fig. 5, which uses the real-time architecture. The observation made here is that not all the Fourier domain points computed in the conventional architecture contribute to the en face OCM image features. In this real-time architecture, only those Fourier domain points that actually contribute to the features of the en face OCM image are computed.

The principles of this architecture may be partially extended to spectroscopic spectral-domain OCM (SD-OCM). Spectroscopic SD-OCM analyzes the spectra as well as the intensity of the backscattered light in a depth-resolved manner. In spectroscopic SD-OCM, raw cross-sectional frames are demodulated to baseband using the fast Fourier transform (FFT) to obtain the depth-dependent signal. After windowing the confocal gate position using a Gaussian window, an inverse FFT is performed. This provides the spectral information corresponding to the confocal gate. This spectral information can be used to study absorption and scattering profiles of particles within the focal volume.10 Certain aspects of the real-time architecture presented here can be used in obtaining a real time spectroscopic OCM image. Instead of performing two FFTs as mentioned above, an infinite impulse response (IIR) filter can be applied on the raw cross-sectional frames to obtain the real-time spectroscopic SD-OCM. The application of the IIR filter on the raw cross-sectional frames can be expedited using the SSE instructions.

A real-time architecture to obtain high-speed en face SD-OCM images has been presented. The real-time architecture takes into account the material dispersion introduced by the sample and compensates for this using a resampling array. Moreover, the real-time architecture takes advantage of the advanced instructions supported by the processor to improve the processing speed. These processing speed enhancements can improve real-time imaging rates when the advantages of SD-OCM are preferred over TD-OCM implementations.

Acknowledgments

This work was supported in part by the National Institutes of Health (Roadmap Initiative, 1 R21 EB005321; NIBIB, 1 R01 EB005221, S.A.B.) and the National Science Foundation (BES03-47747, BES06-19257, S.A.B.). Additional information can be found at www.biophotonics.uiuc.edu.

References

1. 

J. A. Izatt, M. R. Hee, G. M. Owen, E. A. Swanson, and J. G. Fujimoto, “Optical coherence microscopy in scattering media,” Opt. Lett., 19 (8), 590 –592 (1994). 0146-9592 Google Scholar

2. 

E. Beaurepaire, L. Moreaux, F. Amblard, and J. Mertz, “Combined scanning optical coherence and two-photon-excited fluorescence microscopy,” Opt. Lett., 24 (14), 969 –971 (1999). 0146-9592 Google Scholar

3. 

J. P. Dunkers, M. T. Cocerone, and N. R. Washburn, “Collinear optical coherence and confocal fluorescence microscopies for tissue engineering,” Opt. Express, 11 (23), 3074 –3079 (2003). 1094-4087 Google Scholar

4. 

C. Vinegoni, T. S. Ralston, W. Tan, W. Luo, D. L. Marks, and S. A. Boppart, “Integrated structural and functional optical imaging combining spectral-domain optical coherence and multiphoton microscopy,” Appl. Phys. Lett., 88 (5), 053901 –1-3 (2006). https://doi.org/10.1063/1.2171477 0003-6951 Google Scholar

5. 

S. Tang, C.–H. Sun, T. B. Krasieva, Z. Chen, and B. J. Tromberg, “Imaging subcellular scattering contrast by using combined optical coherence and multiphoton microscopy,” Opt. Lett., 32 (5), 503 –505 (2007). https://doi.org/10.1364/OL.32.000503 0146-9592 Google Scholar

6. 

R. Leitgeb, C. K. Hitzenberger, and A. F. Fercher, “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express, 11 (8), 889 –891 (2003). 1094-4087 Google Scholar

7. 

R. A. Leitgeb, W. Drexler, A. Unterhuber, B. Hermann, T. Bajraszewski, T. Le, A. Stingl, and A. F. Fercher, “Ultrahigh resolution Fourier domain optical coherence tomography,” Opt. Lett., 12 (10), 2156 –2158 (2004). 0146-9592 Google Scholar

8. 

R. A. Leitgeb, M. Villiger, A. H. Bachmann, L. Steinmann, and T. Lasser, “Extended focus depth for Fourier domain optical coherence microscopy,” Opt. Lett., 31 (16), 2450 –2452 (2006). https://doi.org/10.1364/OL.31.002450 0146-9592 Google Scholar

9. 

D. L. Marks, A. L. Oldenburg, J. J. Reynolds, and S. A. Boppart, “Autofocus algorithm for dispersion correction in optical coherence tomography,” Appl. Opt., 42 (16), 3038 –3046 (2003). https://doi.org/10.1364/AO.42.003038 0003-6935 Google Scholar

10. 

C. Xu, C. Vinegoni, T. S. Ralston, W. Luo, W. Tan, S. A. Boppart, “Spectroscopic spectral-domain optical coherence microscopy,” Opt. Lett., 31 (8), 1079 –1081 (2006). https://doi.org/10.1364/OL.31.001079 0146-9592 Google Scholar
©(2008) Society of Photo-Optical Instrumentation Engineers (SPIE)
Robin G. Chelliyil, Tyler S. Ralston, Daniel L. Marks, and Stephen A. Boppart M.D. "High-speed processing architecture for spectral-domain optical coherence microscopy," Journal of Biomedical Optics 13(4), 044013 (1 July 2008). https://doi.org/10.1117/1.2960018
Published: 1 July 2008
JOURNAL ARTICLE
4 PAGES


SHARE
Advertisement
Advertisement
Back to Top