You are here: Home > Research > Optical Computational Imaging

Adaptive Spectral Apodization

Adaptive Spectal Adopization In optical coherence tomography (OCT), as in other imaging modalities, one desires to achieve the highest possible resolution given instrument limitations. The bandwidth of the OCT source chiefly determines the useful resolution. However, if the spectrum is not smooth, then the point response in the interferogram will have large sidelobes that cause a degradation of effective resolution and introduce artifacts. As a result, a smooth, Gaussian-like spectrum is often employed to minimize sidelobes. However, there are many sources that do not produce smooth spectra but still produce a wide bandwidth, such as Ti-sapphire oscillators, nonlinear supercontinuum generation optical fibers including microstructured and tapered fibers, and ultrahigh-numerical-aperture fiber sources. We have created a method that produces an estimate of object scattering density (or reflectivity profile) that minimizes sidelobes while not overly increasing the contribution of noise. This method can make broad spectrum sources much more useful for high-resolution OCT.

Spectrum of source (solid line) and sidelobe-corrected instrument response function (dashed line). The OCT system used a UHNA3 fiber as a source, and spectra were obtained from the reflections off of a microscope slide.
Point response of original spectrum (top) and sidelobe-corrected instrument response function (bottom) obtained with a microscope slide as a test sample.
Marks DL, Carney PS, Boppart SA. Adaptive spectral apodization for side-lobe suppression in optical coherence tomography. J. Biomed. Opt. 9(6), 1281-1287, 2004. PDF

Fast FPGA/DSP We have implemented a real-time, multi-dimensional, all digital, OCT acquisition and processing system for use with OCT and doppler imaging. The system uses a high speed A/D converter to acquire data from OCT optical systems that employ a high-speed delay line. Using a floating-point DSP processor and a field programmable gate array, the acquisition and processing system demodulates, filters, and computes the magnitude and phase shifts of the input signal. The data is read by a host machine and displayed on screen at real-time rates. This system offers flexible acquisition and processing parameters for a wide range of optical microscopy techniques.

One implementation of FPGA for Doppler Imaging.

Dispersion Compensation Algorithms Practical clinical optical coherence tomography systems require automatic tools for identifying and correcting flaws in OCT images. One type of flaw is the loss of image detail owing to the dispersion of the medium, which in most cases is unknown. We have developed an autofocus algorithm for estimating the delay line and material dispersion from OCT reflectance data. This autofocus algorithm can be used in conjunction with a high-speed, digital-signal-processor-based OCT acquisition system for rapid image correction.

On the left is a cross section schematic of a PDMS microfluidic structure. On the right are autofocus digitally corrected reflections off of interfaces of the microfluidic structure. Plots (a), (c), (e), and (g) correspond to the uncorrected reflectance functions, whereas (b), (d), (f), and (h) are the corrected point spread functions.

D.L. Marks, A.L. Oldenburg, J.J. Roynolds, S.A. Boppart, Autofocus algorithm for dispersion correction in optical coherence tomography. Applied Optics. 42:3038-3046, 2003. PDF
Marks DL, Oldenburg AL, Reynolds JJ, Boppart SA. Digital algorithm for dispersion correction in optical coherence tomography for homogeneous and stratified media. Applied Optics 42:204-217, 2003. PDF

Projected Index Computed Tomography (PICT) Projected index computed tomography (PICT) is an imaging technique that provides a computed reconstruction of the index of refraction of a sample. PICT makes use of data from standard optical coherence tomography images taken from several view angles to determine a mapping of the refractive indices of the sample. A rectilinear propagation model is assumed, so the data are understood to be related to the line integral of the refractive index in the beam paths. These data thus provide a set of angular projections of the sample. The spatial distribution of the index of the object may then be reconstructed by use of standard filtered backprojection techniques. The resultant PICT images are free of the spatial distortion that is inherent in standard optical cross-sectional images and correspond well to the manufactured dimensions of specific samples. PICT has the potential to produce high-resolution, highly sensitive, and spatially accurate images of a variety of biological and material samples.

Glass capillary tube sample (left) is imaged by PICT method (right), which shows a spatial map refractive index. Linear artifacts in the image stem from the sharp discontinuities in projected data.

Zysk AM, Reynolds JJ, Marks DL, Carney PS, Boppart SA. Projected index computed tomography. Optics Letters. 28:701-703, 2003. PDF

Image Processing It is possible to achieve Doppler imaging with OCT when imaging blood vessels or any other movement. When blood flows towards or away from the transducer, the frequency of the backscattered light varies depending on the velocity of the fluid the light is incident upon. This shift can then be used to calculate the velocity and direction of the blood. This technology is especially useful in the study of microfluidics. Below is an example of Doppler imaging in microfluidics. The left image is a cross-section of a tubular microfluidic channel after Doppler OCT processing. On the right is a graph showing the flow velocity measured vs. the ideal flow profile.

Schaefer AW, Reynolds JJ, Marks DL, Boppart SA. Real-time digital signal processing-based optical coherence tomography and Doppler optical coherence tomography. IEEE Trans. Biomed. Engr. 51:186-190, 2004 PDF

Speckle Reduction in OCT Images In OCT, ultrasound, synthetic-aperture radar, and other coherent ranging methods, speckle can cause spurious details that detract from the utility of the image.  Speckle is a problem inherent to imaging densely scattering objects with limited bandwidth.  We have developed methods to minimize speckle and improve image quality while maintaining features that are consistent with the known data.

OCT images of structures in a Xenopus laevis (African frog) tadpole before and after application of the I-divergence regularization method to remove speckle.

Marks DL, Ralston TS, Boppart SA. Speckle reduction by I-divergence regularization in optical coherence tomography. JOSA A, 22:2366-2371, 2005. PDF

Interferometric Synthetic Aperture Microscopy (ISAM) State-of-the-art methods in high-resolution three-dimensional optical microscopy requires that the focus be scanned through the entire region of interest. However, an analysis of the physics of the light-sample interaction reveals that the Fourier space coverage is independent of depth. We have shown that be solving the inverse scattering problem for coherence microscopy, computed reconstruction yields volumes with a resolution in all planes that is equivalent to the resolution achieved only at the focal plane for conventional high-resolution microscopy. The entire illuminated volume has spatially-invariant resolution, thus eliminating the compromise between resolution and depth-of-field. This novel computational image-formation technique is called ISAM, and has the potential to broadly impact real-time 3-D microscopy and analysis in the fields of cell and tumor biology, as well as in clinical diagnosis where in vivo imaging is preferable to biopsy. We have analyzed and demonstrated this technique for rectangular Cartesian-coordinate imaging, for rotational catheter-based imaging, and for full-field imaging.

Original object model of point scatterers (left) and simulated OCT image (right) showing the central in-focus region with blurring above and below the focal region.

Unfiltered reconstruction of data shown above (left) and Tikhonov regularized solution of data shown above (right).  Note that spatially-invariant resolution is apparent, both inside and outside of the confocal parameter of the focused Gaussian beam.

Ralston TS, Marks DL, Carney PS, Boppart SA.  Inverse scattering for optical coherence tomography.  JOSA A, 23(5):1027-1037, 2006. PDF

Simulated (top) and reconstructed (bottom) radial OCT catheter images of randomly distributed point scatterers.

Marks DL, Ralston TS, Carney PS, Boppart SA.  Inverse scattering for rotationally scanned optical coherence tomography.  J. Opt. Soc. Am. A, 23(10):2433-2439, 2006. PDF

Simulated volume of point scatterers imaged with OCT to demonstrate spatially-invariant resolution. Projection of (a) unprocessed and (b) processed data on 2-D plane, and full 3-D rendering of (c) unprocessed and (d) processed volume of point scatterers. All axes are labeled in units of wavelength. Arrows in (c) indicate wavefront interference which is resolved in (d).

Ralston TS, Marks DL, Boppart SA, Carney PS. Inverse scattering for high-resolution interferometric microscopy. Optics Letters, In press, 2006.

Questions or comments? Please contact us
Web site designed and built by Simon Schlachter and Freddy Nguyen. Special thanks to Ron Stack.
Copyright © 2005 Stephen A. Boppart, Biophotonics Imaging Laboratory