0
JBO Letters

Speckle reduction in optical coherence tomography using angular compounding by B-scan Doppler-shift encoding

[+] Author Affiliations
Hui Wang, Andrew M. Rollins

Case Western Reserve University, Department of Biomedical Engineering, 10900 Euclid Avenue, Cleveland, Ohio 44120

J. Biomed. Opt. 14(3), 030512 (June 23, 2009). doi:10.1117/1.3155523
History: Received February 04, 2009; Revised April 16, 2009; Accepted May 07, 2009; Published June 23, 2009
Text Size: A A A

Open Access Open Access

* Address all correspondence to Hui Wang; Tel: 216–368–0609; E-mail: Hxw26@case.edu

We propose a novel method for speckle reduction for optical coherence tomography based on angular compounding by B-scan Doppler-shift encoding (AngularCBD). By decentering the probe beam from the pivot of a scanning mirror, the illumination angle represented by different components of the beam can be encoded in Doppler shift. Compounding multiple images reconstructed from different Doppler-shift bands, we can suppress speckle without sacrificing image acquisition speed. Speckle reduction with AngularCBD is demonstrated by imaging a phantom and tissue sample in vitro and in vivo.

Figures in this Article

In optical coherence tomography (OCT) imaging, speckle arises from multiple scattering in the sample and degrades the visibility of fine structures in images. However, speckle is also the carrier of the information of tissue structure. Many efforts have been made to reduce speckle in OCT. In general, noise can be suppressed by averaging if successive acquisitions contain the same signal but uncorrelated noise. Speckle can be suppressed in this way by using frequency compounding, polarization compounding, or angular compounding.15 Speckle reduction by angular compounding has been realized by multiple methods that involve acquiring multiple images with slightly different beam incidence angles.2,67 This has required sacrificing imaging speed or use of bulky optical components.4,69 Angular compounding by path-length encoding (ACPE) can suppress speckle without reducing imaging speed.2 However, ACPE requires a long imaging range in order to separate the images from different angles. Thus, it is difficult to achieve this with spectral-domain OCT (SDOCT). Also, the necessary delay element introduces extra dispersion, which could be a problem using ultrahigh-resolution OCT. Here, we introduce a novel and simple method to suppress speckle, angular compounding by B-scan Doppler-shift encoding (AngularCBD).

A schematic of AngularCBD is shown in Fig. 1, which represents a typical sample arm used in a bench-top OCT. The collimated beam in the sample arm is aligned offset from the pivot of the scanning mirror, so that a Doppler shift is introduced during the course of a B-scan. This principle has been previously proposed for generating a carrier signal for enface time-domain OCT imaging10 and has more recently been used to resolve the complex-conjugate ambiguity in SDOCT.1113 Because the sample beam is incident on a range of positions on the scanning mirror before the objective lens, different incidence angles are encoded by different Doppler shifts. It has been shown that the Doppler shift is linearly proportional to the offset from the pivot of the scanning mirror.11

Grahic Jump LocationF1 :

(a) Schematic of the scanner used for AngularCBD and angle encoded Doppler shift spectrum by physically blocking part of the beam. (b) Measured angle encoded Doppler shift spectrums at 1300 nm, which are illustrated in (a). The spectrum was measured with a uniform phantom described in the context.

To illustrate the concept, an experiment was conducted by physically blocking halves of the beam alternatively. The corresponding Doppler shifts are shown in Fig. 1. A phantom made by mixing 2% intralipid and agarose gel was imaged with SDOCT at 1300nm. The axial resolution of the system was 14μm in air. The beam in the sample arm, with a beam diameter of 6mm, was focused by an achromatic lens with a focal length of 35mm to a spot of 10μm diam. With a power at the sample of 10mW, the sensitivity was measured to be 110dB with 6dB decay at 2mm. The center of the beam was offset from the pivot by 4.2mm, which introduce a Doppler shift of 8kHz during a B-scan. All images were acquired at 40k A-scans/s. Each frame consists of 1000 A-scans across 2.5mm on the sample, resulting in a transverse imaging pitch of 2.5μm. Half-beams 1 and 2 were blocked alternately as illustrated in Fig. 1. The Dopplershift acquired over the course of a B-scan was calculated by Fourier transform of the data at one pixel of the photodiode array. With half-beam 1 blocked, only the signal from path 2-2 was acquired and the observed Doppler shift was concentrated in the higher frequency range, as shown in Figs. 1. Similarly, the observed Doppler shift was in the lower frequency range when half-beam 2 was blocked. The gap between the 1-1 and 2-2 represents the Doppler shift from beams 1-2 and 2-1, which was lost by blocking the half-beam.

The angle-encoded Doppler shift can be used to reconstruct multiple OCT images representing different illumination/collection angles. We used a moving window to bandpass filter the B-scan Doppler shift spectrum. Separated images were achieved by reconstructing OCT images from the filtered Doppler shift bands, which were then summed to form the final image. Because the separated images are angle dependent, they contain different speckle patterns. Therefore, speckle reduction can be expected by image compounding. The same phantom employed in Fig. 1 was also used to quantify the speckle contrast after speckle reduction. A moving Tukey window was used to filter the B-scan Doppler shift spectrum as shown in Fig. 1 without beam blocking. Figure 2 shows the five filtered B-scan Doppler shift subbands with the moving window. Figure 2 plots images reconstructed from subbands 1, 2, 3, and 4. To quantify the correlation of the speckle patterns between different subband images, the cross-correlation coefficients were calculated following a previously described method,3 where the exponential decay with depth was first corrected and then an area of about 200×200μm2 was selected for calculation. The cross-correlation between the neighbor images in Fig. 2 is <0.1. This confirms that each image predominantly represents independent speckle patterns resulting from different illumination/collection angles. Therefore, speckle reduction can be expected after compounding the five images. For comparison, the original image and the image after speckle reduction are shown in Fig. 3. The speckle contrast C is defined as the ratio of the standard deviation to the mean of the pixel values in the selected region. The speckle contrast of Fig. 3 is 0.57, which is very close to the theoretical value of 0.52 for fully developed speckle.3 After image compounding, the speckle contrast is 0.31 measured from the same region. Thus, speckle contrast has been reduced by a factor of 1.84. A reduction factor of 2.2 is expected for five independent angle-related images.3 The difference can be attributed to nonuniform light power distribution between the five filtered subbands and residual overlap between the subbands. Significant reduction of variation caused by speckle along single A-scan is shown in Figs. 3 as well as shown in the images of Figs. 3.

Grahic Jump LocationF2 :

(a) Five B-scan Doppler shift subbands filtered with a moving Tukey window with a width of 150 points. (b) The corresponding images of subbands 1 2 3 and 4 (50×50μm2). The image pixel values have been normalized for display and are shown in logarithmic scale.

Grahic Jump LocationF3 :

(a) A segment of the original image without speckle reduction, (b) the image after AngularCBD speckle reduction, and (c) a image after a 7-A-scan rolling average. (d–f) Single A-scans of the entire depth of phantom images. Images are shown in logarithmic scale. Scale bar: 50μm.

Although significant speckle reduction is achieved, implementation of AngularCBD requires a compromise between speckle reduction and lateral resolutionBecause only part of the beam is used for each subimage, the effective imaging N.A is reduced by the bandpass filter. To measure the transverse resolution reduction factor, a USAF target with multiple parallel bars covered by a phantom with a thickness of 600μm was imaged, and the slope of the edges of the bars across a B-scan was measured. The lateral resolution was found to be reduced by a factor of 1.84. A simple blurring function can also be expected to reduce speckle contract. Therefore, in order to demonstrate that the AngularCBD method is not simply equivalent to a lateral blur, a rolling average of seven neighboring A-scans was computed on the original image, and the result is shown in Fig. 3. As a result, the lateral resolution was reduced by a factor of 1.72, closely matching the resolution reduction resulting from the AngularCBD processing. The speckle contrast was reduced by the blur to 0.49, which represents a factor of 1.16 reduction, much less than the reduction factor achieved by AngularCBD.

We demonstrated the technique by imaging biological tissue samples in vivo and in vitro. Speckle reduction is apparent and fine structures are visualized with better contrast after AngularCBD, as shown in Fig. 4 compared to the original image [Fig. 4] by imaging a living quail embryo in vivo. Figures 4 show images of a porcine esophagus in vitro before and after AngularCBD. The layered structures are better defined after AngularCBD.

Grahic Jump LocationF4 :

In vivo image of quail embryo (a) before and (b) after AngularCBD. In vitro porcine esophagus (c) before and (d) after AngularCBD. Contrast of fine structures and layers are clearly enhanced after speckle reduction. Images are shown in logarithmic scale. Scale bar: 1mm.

Further speckle reduction is possible by choosing a smaller filtering window to increase the number of angle-related subband images. However, the lateral resolution will be decreased further. Just increasing the number of subband images with a fixed width window does not improve the speckle reduction further due to increased subband overlap. The AngularCBD method is very flexible because all processing is carried out after image acquisition. The width and the number of the windows can be chosen based on the sample and desired speckle reduction factor. As other B-scan Doppler methods, transverse sampling must be sufficient to sample the full speckle pattern.

In practical implementation, the method is straightforward to optimize. In order for the entirety of the B-scan Doppler spectrum to be shifted to one side of dc, the entirety of the beam should be offset to one side of the mirror pivot. It will be good to choose a larger beam to compensate the loss of the lateral resolution after ABCD.

If there is flow or motion in the tissue being imaged, an additional Doppler shift will be imposed on the sample beam, which will upshift or downshift the B-scan Doppler frequency at that location. Thus, the flow signal will be translated to different locations in the B-scan Doppler spectrum. This effect has been used for Doppler OCT flow imaging,12 but should not affect the AngularCBD method as long as the range of the B-scan Doppler spectrum analyzed encompasses the shifted signal components. Therefore, AngularCBD can be used simultaneously with deconjugation of the mirror images in SDOCT and flow tracking.1114

In conclusion, we have introduced a novel method, AngularCBD, to reduce speckle in OCT images. We demonstrated speckle reduction by a factor of 1.84 by incoherently adding five Doppler-shift angle-encoded images. AngularCBD is compatible with most typical OCT scanners making use of tilting mirrors. Important features of AngularCBD are that it does not interfere with the frame rate nor introduce any optical nor mechanical components nor extra dispersion. The only needed modification to a conventional scanner is to decenter the beam from the pivot of the scanning mirror, which is typically readily realizable.

Acknowledgments

The project described was supported by the National Institutes of Health (NIH) (grant Nos. R24CA110943 and R01 CA114276–01).

Schmitt  J. M., , Xiang  S. H., , and Yung  K. M., “ Speckle in optical coherence tomography. ,” J. Biomed. Opt..  1083-3668 4, (1 ), 95–105  ((1999)).
Iftimia  N., , Bouma  B. E., , and Tearney  G. J., “ Speckle reduction in optical coherence tomography by “path length encoded” angular compounding. ,” J. Biomed. Opt..  1083-3668 8, (2 ), 260–263  ((2003)).
Pircher  M., , Gotzinger  E., , Leitgeb  R., , Fercher  A. F., , and Hitzenberger  C. K., “ Speckle reduction in optical coherence tomography by frequency compounding. ,” J. Biomed. Opt..  1083-3668 8, (3 ), 565–569  ((2003)).
Kim  J., , Miller  D. T., , Kim  E., , Oh  S., , Oh  J., , and Milner  T. E., “ Optical coherence tomography speckle reduction by a partially spatially coherent source. ,” J. Biomed. Opt..  1083-3668 10, (6 ), 064034–064039  ((2005)).
Desjardins  A. E., , Vakoc  B. J., , Tearney  G. J., , and Bouma  B. E., “ Speckle reduction in OCT using massively-parallel detection and frequency-domain ranging. ,” Opt. Express.  1094-4087 14, (11 ), 4736–4745  ((2006)).
Schmitt  J. M., “ Array detection for speckle reduction in optical coherence microscopy. ,” Phys. Med. Biol..  0031-9155 42, (7 ), 1427−1439  ((1997)).
Desjardins  A. E., , Vakoc  B. J., , Oh  W. Y., , Motaghiannezam  S. M., , Tearney  G. J., , and Bouma  B. E., “ Angle-resolved optical coherence tomography with sequential angular selectivity for speckle reduction. ,” Opt. Express.  1094-4087 15, (10 ), 6200–6209  ((2007)).
Bashkansky  M., and Reintjes  J., “ Statistics and reduction of speckle in optical coherence tomography. ,” Opt. Lett..  0146-9592 25, (8 ), 545–547  ((2000)).
Popescu  D. P., , Hewko  M. D., , and Sowa  M. G., “ Speckle noise attenuation in optical coherence tomography by compounding images acquired at different positions of the sample. ,” Opt. Commun..  0030-4018 269, (1 ), 247–251  ((2007)).
Podoleanu  A. G., , Dobre  G. M., , and Jackson  D. A., “ En face coherence imaging using galvanometer scanner modulation. ,” Opt. Lett..  0146-9592 23, (3 ), 147–149  ((1998)).
An  L., and Wang  R. K., “ Use of a scanner to modulate spatial interferograms for in vivo full-range Fourier-domain optical coherence tomography. ,” Opt. Lett..  0146-9592 32, (23 ), 3423–3425  ((2007)).
Wang  R. K., , Jacques  S. L., , Ma  Z., , Hurst  S., , Hanson  S. R., , and Gruber  A., “ Three dimensional optical angiography. ,” Opt. Express.  1094-4087 15, (7 ), 4083–4097  ((2007)).
Leitgeb  R. A., , Michaely  R., , Lasser  T., , and Sekhar  S. C., “ Complex ambiguity-free Fourier domain optical coherence tomography through transverse scanning. ,” Opt. Lett..  0146-9592 32, (23 ), 3453–3455  ((2007)).
Baumann  B., , Pircher  M., , Guzinger  E., , and Hitzenberger  C. K., “ Full range complex spectral domain optical coherence tomography without additional phase shifters. ,” Opt. Express.  1094-4087 15, (20 ), 13375–13387  ((2007)).
© 2009 Society of Photo-Optical Instrumentation Engineers

Citation

Hui Wang and Andrew M. Rollins
"Speckle reduction in optical coherence tomography using angular compounding by B-scan Doppler-shift encoding", J. Biomed. Opt. 14(3), 030512 (June 23, 2009). ; http://dx.doi.org/10.1117/1.3155523


Access This Article

Figures

Grahic Jump LocationF1 :

(a) Schematic of the scanner used for AngularCBD and angle encoded Doppler shift spectrum by physically blocking part of the beam. (b) Measured angle encoded Doppler shift spectrums at 1300 nm, which are illustrated in (a). The spectrum was measured with a uniform phantom described in the context.

Grahic Jump LocationF2 :

(a) Five B-scan Doppler shift subbands filtered with a moving Tukey window with a width of 150 points. (b) The corresponding images of subbands 1 2 3 and 4 (50×50μm2). The image pixel values have been normalized for display and are shown in logarithmic scale.

Grahic Jump LocationF3 :

(a) A segment of the original image without speckle reduction, (b) the image after AngularCBD speckle reduction, and (c) a image after a 7-A-scan rolling average. (d–f) Single A-scans of the entire depth of phantom images. Images are shown in logarithmic scale. Scale bar: 50μm.

Grahic Jump LocationF4 :

In vivo image of quail embryo (a) before and (b) after AngularCBD. In vitro porcine esophagus (c) before and (d) after AngularCBD. Contrast of fine structures and layers are clearly enhanced after speckle reduction. Images are shown in logarithmic scale. Scale bar: 1mm.

Tables

References

Schmitt  J. M., , Xiang  S. H., , and Yung  K. M., “ Speckle in optical coherence tomography. ,” J. Biomed. Opt..  1083-3668 4, (1 ), 95–105  ((1999)).
Iftimia  N., , Bouma  B. E., , and Tearney  G. J., “ Speckle reduction in optical coherence tomography by “path length encoded” angular compounding. ,” J. Biomed. Opt..  1083-3668 8, (2 ), 260–263  ((2003)).
Pircher  M., , Gotzinger  E., , Leitgeb  R., , Fercher  A. F., , and Hitzenberger  C. K., “ Speckle reduction in optical coherence tomography by frequency compounding. ,” J. Biomed. Opt..  1083-3668 8, (3 ), 565–569  ((2003)).
Kim  J., , Miller  D. T., , Kim  E., , Oh  S., , Oh  J., , and Milner  T. E., “ Optical coherence tomography speckle reduction by a partially spatially coherent source. ,” J. Biomed. Opt..  1083-3668 10, (6 ), 064034–064039  ((2005)).
Desjardins  A. E., , Vakoc  B. J., , Tearney  G. J., , and Bouma  B. E., “ Speckle reduction in OCT using massively-parallel detection and frequency-domain ranging. ,” Opt. Express.  1094-4087 14, (11 ), 4736–4745  ((2006)).
Schmitt  J. M., “ Array detection for speckle reduction in optical coherence microscopy. ,” Phys. Med. Biol..  0031-9155 42, (7 ), 1427−1439  ((1997)).
Desjardins  A. E., , Vakoc  B. J., , Oh  W. Y., , Motaghiannezam  S. M., , Tearney  G. J., , and Bouma  B. E., “ Angle-resolved optical coherence tomography with sequential angular selectivity for speckle reduction. ,” Opt. Express.  1094-4087 15, (10 ), 6200–6209  ((2007)).
Bashkansky  M., and Reintjes  J., “ Statistics and reduction of speckle in optical coherence tomography. ,” Opt. Lett..  0146-9592 25, (8 ), 545–547  ((2000)).
Popescu  D. P., , Hewko  M. D., , and Sowa  M. G., “ Speckle noise attenuation in optical coherence tomography by compounding images acquired at different positions of the sample. ,” Opt. Commun..  0030-4018 269, (1 ), 247–251  ((2007)).
Podoleanu  A. G., , Dobre  G. M., , and Jackson  D. A., “ En face coherence imaging using galvanometer scanner modulation. ,” Opt. Lett..  0146-9592 23, (3 ), 147–149  ((1998)).
An  L., and Wang  R. K., “ Use of a scanner to modulate spatial interferograms for in vivo full-range Fourier-domain optical coherence tomography. ,” Opt. Lett..  0146-9592 32, (23 ), 3423–3425  ((2007)).
Wang  R. K., , Jacques  S. L., , Ma  Z., , Hurst  S., , Hanson  S. R., , and Gruber  A., “ Three dimensional optical angiography. ,” Opt. Express.  1094-4087 15, (7 ), 4083–4097  ((2007)).
Leitgeb  R. A., , Michaely  R., , Lasser  T., , and Sekhar  S. C., “ Complex ambiguity-free Fourier domain optical coherence tomography through transverse scanning. ,” Opt. Lett..  0146-9592 32, (23 ), 3453–3455  ((2007)).
Baumann  B., , Pircher  M., , Guzinger  E., , and Hitzenberger  C. K., “ Full range complex spectral domain optical coherence tomography without additional phase shifters. ,” Opt. Express.  1094-4087 15, (20 ), 13375–13387  ((2007)).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement


Buy this article ($18 for members, $25 for non-members).
Sign In