OSA's Digital Library

Optics Letters

Optics Letters

| RAPID, SHORT PUBLICATIONS ON THE LATEST IN OPTICAL DISCOVERIES

  • Editor: Alan E. Willner
  • Vol. 37, Iss. 14 — Jul. 15, 2012
  • pp: 2871–2873
« Show journal navigation

Speckle reduction of endovascular optical coherence tomography using a generalized divergence measure

Kyle H. Y. Cheng, Edmund Y. Lam, Beau A. Standish, and Victor X. D. Yang  »View Author Affiliations


Optics Letters, Vol. 37, Issue 14, pp. 2871-2873 (2012)
http://dx.doi.org/10.1364/OL.37.002871


View Full Text Article

Acrobat PDF (488 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Endovascular optical coherence tomography (EV-OCT) is an emerging intravascular imaging technique for observing blood vessel walls. Fluctuating speckle noise, especially during rapid pull-back, can severely degrade the visibility of morphological structures. Moreover, the speckle pattern varies in different parts of the image due to beam divergence and is further complicated by interpolation through the coordinate transformation necessary for displaying the rotary scanning images, challenging the use of frequency domain analysis. In this study, a computationally efficient method using a generalized divergence regularization procedure is presented to suppress speckle noise in EV-OCT images. Results show substantial smoothing of the grainy appearance and enhanced visualization of deeper structures as demonstrated in porcine carotid arteries.

© 2012 Optical Society of America

Endovascular optical coherence tomography (EV-OCT) is an emerging intravascular imaging technique for major blood vessel wall visualization with 15μm resolution [1

1. A. F. Low, G. J. Tearney, B. E. Bouman, and I. Jang, Nat. Clin. Pract. Cardiovasc. Med. 3, 154 (2006). [CrossRef]

]. Numerous studies have demonstrated that EV-OCT can distinguish among the media, external elastic lamina, and the adventitia clearly in normal arteries [1

1. A. F. Low, G. J. Tearney, B. E. Bouman, and I. Jang, Nat. Clin. Pract. Cardiovasc. Med. 3, 154 (2006). [CrossRef]

]. In addition, EV-OCT has also been shown to delineate the fibrous cap, lipid deposits, and neovascularization in atherosclerotic plaques [1

1. A. F. Low, G. J. Tearney, B. E. Bouman, and I. Jang, Nat. Clin. Pract. Cardiovasc. Med. 3, 154 (2006). [CrossRef]

]. However, due to water absorption and tissue scattering, the penetration depth of EV-OCT is limited to 2mm. The presence of speckle noise further complicates visualization as its modulated appearance can further distort the faint morphological signals of deeply situated structures (e.g., adventitia). In between B-mode frames of an EV-OCT pullback image sequence, the speckle pattern fluctuates dramatically due to motion artifacts from various sources [2

2. A. Mariampillai, B. A. Standish, E. H. Moriyama, M. Khurana, N. R. Munce, M. K. K. Leung, J. Jiang, A. Cable, B. C. Wilson, I. A. Vitkin, and V. X. D. Yang, Opt. Lett. 33, 1530 (2008). [CrossRef]

,3

3. A. Mariampillai, M. K. K. Leung, M. Jarvi, B. A. Standish, K. Lee, B. C. Wilson, A. Vitkin, and V. X. D. Yang, Opt. Lett. 35, 1257 (2010). [CrossRef]

].

Common hardware based approaches include frequency [4

4. M. Pircher, E. Gotzinger, R. Leitgeb, A. F. Fercher, and C. K. Hitzenberger, J. Biomed. Opt. 8, 565 (2003). [CrossRef]

] and angular compounding [5

5. A. E. Desjardins, B. J. Vakoc, W. Y. Oh, S. M. R. Motaghiannezam, G. J. Tearney, and B. E. Bouma, Opt. Express 15, 6200 (2007). [CrossRef]

]. They are robust ways of speckle suppression as speckle properties vary across wavelengths or different illumination angles [6

6. J. W. Goodman, Speckle Phenomena in Optics (Roberts & Company, 2007).

]. However, since they require additional components, they may not be applicable to existing commercial EV-OCT systems.

Various studies have been devoted to speckle reduction using digital image processing techniques. Many of these are based on integral transforms and the manipulation of the associated coefficients [7

7. D. C. Adler, T. H. Ko, and J. G. Fujimoto, Opt. Lett. 29, 2878 (2004). [CrossRef]

9

9. Z. Jian, Z. Yu, L. Yu, B. Rao, Z. Chen, and B. J. Tromberg, Opt. Lett. 34, 1516 (2009). [CrossRef]

]. Others rely on various algebraic iterations [10

10. J. Rogowska and M. E. Brezinski, IEEE Trans. Med. Imag. 19, 1261 (2000). [CrossRef]

,11

11. P. Puvanathasan and K. Bizheva, Opt. Express 17, 733 (2009). [CrossRef]

]. Optical beam divergence, nonuniform spatial sampling of the rotary scanning mechanism, and the subsequent interpolation all contribute to the speckle pattern shape and orientation variation in different parts of the image. These represent challenges to coefficient based speckle reduction methods. Since a typical EV-OCT pullback image sequence consists of several hundred images, algorithms with complex algebraic operations may not be efficient to tackle the large amounts of data.

In this Letter, we propose a simple, fast converging iterative algorithm modified from [12

12. D. L. Marks, T. S. Ralston, and S. A. Boppart, J. Opt. Soc. Am. A 22, 2366 (2005). [CrossRef]

]. It is demonstrated that by regularizing the image during each iteration using a generalized divergence measure, β-divergence, speckle suppression and edge preservation can be achieved simultaneously. Moreover, one can freely adjust the trade-off between speckle suppression and edge preservation by adjusting a single parameter.

β-divergence was used for blind source separation [13

13. M. Mihoko and S. Eguchi, Neural Comput. 14, 1859 (2002). [CrossRef]

] and is defined as
dβ(x,u)=1β(β1)[xβ+(β1)uββxuβ1],
(1)
where x equals, in our case, the despeckled image, and u denotes a reference image. In this equation, β is a scalar parameter to be chosen. The β-divergence has the important property, where only when x=u will the divergence measure vanish. It does not subject the image to error by generating undefined pixels since there are no reciprocal or logarithmic terms of x or u for β>1. The measure was also shown to be robust against outliers [13

13. M. Mihoko and S. Eguchi, Neural Comput. 14, 1859 (2002). [CrossRef]

].

In this study, the β-divergence algorithm was developed to first enforce the least-square consistency of the despeckled image with the measured data, and then focus on the matching of details between the image and a reference image u using the β-divergence regularization measure. To achieve this, the following inverse problem was solved:
x^=argxminyPx2+λi1β(β1)[xi2β+(β1)ui2ββxi2ui2(β1)],
(2)
where y is the measured image, P is a downscaled point spread function, and the subscript i indicates the ith pixel of the image. By taking the derivative of Eq. (2) with respect to x, we arrive at the iterative scheme:
x^i(n+1)=x^inε{2(P*PxP*y)i+λβ(β1)[2β|xi|2β12β|xi||ui|2(β1)]},
(3)
where P* is the adjoint of P, ε is the iteration constant, and λ is the regularization constant. The image was normalized at each iteration to limit its bandwidth.

OCT images were obtained using a commercial Lightlab C7-XR Fourier Domain OCT system (Lightlab Imaging). To demonstrate the speckle suppression effect of the algorithm, an initial test image of human finger skin was obtained by pressing two fingertips against the imaging catheter tip. In vivo endovascular OCT imaging of the porcine carotid arterial wall image was obtained subsequently with eccentric placement of the imaging catheter to demonstrate the versatility of our algorithm to suppress speckle at different parts of the image. The details of animal imaging protocol have been previously described elsewhere [14

14. K. H. Y. Cheng, C. Sun, B. Vuong, K. K. C. Lee, A. Mariampillai, T. R. Marotta, J. Spears, W. J. Montanera, P. R. Herman, K. Rasmus, B. A. Standish, and V. X. D. Yang, Biomed. Opt. Express 3, 388 (2012). [CrossRef]

]. All animal procedures were approved by St. Michael’s Hospital (Toronto, Ontatio) Animal Care Committee.

In the algorithm, ε was chosen to be 0.75 and λ was set to 0.6 determined interactively by the OCT operator. The reference image u was generated by passing the unprocessed image to a 12×12 median filter and then filtered by a 5×5 averaging filter. Median filters have been used before for speckle suppression [15

15. A. Ozcan, A. Bilenca, A. E. Desjardins, B. E. Bouma, and G. J. Tearney, J. Opt. Soc. Am. A 24, 1901 (2007). [CrossRef]

], while averaging filters can effectively remove additive noise, since speckle in the EV-OCT image after logarithmic compression becomes additive. The algorithm was set to run for 40 iterations as the result started to converge to 0.002 mean square difference from u.

Figure 1(a) shows the original finger skin image, and Fig. 1(b) demonstrates the despeckled finger skin image processed with our proposed algorithm. Most of the features were preserved in the despeckled image due to the robustness of β-divergence against outliers, while the grainy appearance was mostly smoothed. Figures 2 and 3 show the original and despeckled (β=3) carotid arterial wall image, respectively. After speckle suppression, the grainy appearance of the tissue was clearly removed.

Fig. 1. (a) Original finger skin image; (b) despeckled finger skin image. (*) indicates the stratus corneum, (#) indicates the dermis, and the red arrows point to the sweat glands.
Fig. 2. Original porcine carotid artery EV-OCT image. The red ROI indicates the signal region and the yellow ROI indicates the noise region used in the metrics calculation. The green ROI indicates the zoomed region in Fig. 5. The three blue ROIS are used for ENL calculations.
Fig. 3. Despeckled porcine carotid artery EV-OCT image.

To quantitatively evaluate the performance of the algorithm on real endovascular images, a number of metrics were calculated for a region of interest (ROI) encompassing high signal regions (media) and low signal regions (adventitia) versus a noise background as depicted in Fig. 2. The signal-to-noise ratio (SNR) was defined as SNR=20log(xlin/σlin), where xlin is the maximum intensity in the ROI in linear scale. σlin is the standard deviation of the noise region in linear scale. The contrast-to-noise ratio (CNR) was defined as CNR=10log[(μxμb)/(σx2σb2)0.5], where μ and σ are the mean and standard deviation and subscripts x and b denote the ROI and noise region, respectively. The equivalent number of looks (ENL) was defined as ENL=μx2/σx2 and is averaged between three ROIs depicted in Fig. 2. An edge preservation parameter is defined previously [8

8. P. Puvanathasan and K. Bizheva, Opt. Express 15, 15747 (2007). [CrossRef]

]. Briefly, the larger the parameter, the more edges are preserved. Figures 4(a), 4(b), and 4(c) show the enlarged views of the ROIs without processing, processed by [12

12. D. L. Marks, T. S. Ralston, and S. A. Boppart, J. Opt. Soc. Am. A 22, 2366 (2005). [CrossRef]

] and processed by our proposed algorithm, respectively. As noted by the arrows in Fig. 4(c), visually these structures are very difficult to be discerned in Fig. 4(a). Figure 5 shows a region far from the catheter without and with speckle suppression. The result demonstrates that our algorithm can better delineate features, despite the fact that they are located further from the imaging catheter where a diverging optical beam provides less SNR. This also demonstrates our technique is robust against variations in speckle orientation.

Fig. 4. (a) Enlarged view of the red ROI in the original porcine arterial wall image; (b) enlarged view of the red signal ROI in the despeckled image processed by I-divergence algorithm; (c) enlarged view of the red signal ROI in the despeckled image processed by our proposed algorithm. The arrows indicate structures that are much less visible in the unprocessed image.
Fig. 5. (a) Enlarged view of the green ROI in the original image; (b) enlarged view of the green signal ROI in the despeckled image processed by I-divergence algorithm; (c) enlarged view of the green signal ROI in the despeckled image processed by our proposed algorithm. The red arrow indicates the clearer visualization of the external elastic lamina.

Table 1 shows the results when β=1.5, 3, and 4. In all cases, the SNR, CNR, and ENL all showed improvement compared to the original image. The result obtained by a 12×12 median filter was also listed. It produces very good SNR, but edges were poorly preserved. To further show how β affects the speckle suppression effect, CNR and edge preservation were plotted against β in Fig. 6. When β is less than 3, the CNR is higher but the edge preservation drops, and vice versa when β is larger than 3. Thus, β acts as an adjustable parameter that controls the trade-off between speckle suppression and edge preservation. Such trade-off adjustment is needed when one wants to see different features. For example, when one wants to evaluate large area atherosclerotic plaques, β should be set lower than 3 to apply stronger speckle suppression so as to clearly delineate different parts of the plaque. However, when one wants to seek neovascularization, β should be set higher than 3 in order to preserve the vessel edges.

Table 1. Metrics Comparison Between the Original Image and Processed Images Using Different β Values

table-icon
View This Table
Fig. 6. Trade-off between CNR and edge preservation by adjustment of β.

We note that one of the shortcomings is the lack of “gold standard” reference image u, which is not unique to our technique. Furthermore, the algorithm was run on an Intel Corei5 computer implemented in MATLAB. Speckle suppression for a 969×969 image took 7 seconds. Future implementation in C++ and/or GPU parallel processing may be required to perform real-time EV-OCT speckle reduction.

In summary, we proposed and demonstrated an effective speckle suppression method for EV-OCT images while preserving edge information. Furthermore, to our knowledge, we demonstrated for the first time how speckle reduction can potentially benefit EV-OCT by enhancing visibility of endovascular morphological structures with a computational method.

The work is supported by the National Engineering Research Council of Canada, the Canada Research Chairs, and Ryerson University.

References

1.

A. F. Low, G. J. Tearney, B. E. Bouman, and I. Jang, Nat. Clin. Pract. Cardiovasc. Med. 3, 154 (2006). [CrossRef]

2.

A. Mariampillai, B. A. Standish, E. H. Moriyama, M. Khurana, N. R. Munce, M. K. K. Leung, J. Jiang, A. Cable, B. C. Wilson, I. A. Vitkin, and V. X. D. Yang, Opt. Lett. 33, 1530 (2008). [CrossRef]

3.

A. Mariampillai, M. K. K. Leung, M. Jarvi, B. A. Standish, K. Lee, B. C. Wilson, A. Vitkin, and V. X. D. Yang, Opt. Lett. 35, 1257 (2010). [CrossRef]

4.

M. Pircher, E. Gotzinger, R. Leitgeb, A. F. Fercher, and C. K. Hitzenberger, J. Biomed. Opt. 8, 565 (2003). [CrossRef]

5.

A. E. Desjardins, B. J. Vakoc, W. Y. Oh, S. M. R. Motaghiannezam, G. J. Tearney, and B. E. Bouma, Opt. Express 15, 6200 (2007). [CrossRef]

6.

J. W. Goodman, Speckle Phenomena in Optics (Roberts & Company, 2007).

7.

D. C. Adler, T. H. Ko, and J. G. Fujimoto, Opt. Lett. 29, 2878 (2004). [CrossRef]

8.

P. Puvanathasan and K. Bizheva, Opt. Express 15, 15747 (2007). [CrossRef]

9.

Z. Jian, Z. Yu, L. Yu, B. Rao, Z. Chen, and B. J. Tromberg, Opt. Lett. 34, 1516 (2009). [CrossRef]

10.

J. Rogowska and M. E. Brezinski, IEEE Trans. Med. Imag. 19, 1261 (2000). [CrossRef]

11.

P. Puvanathasan and K. Bizheva, Opt. Express 17, 733 (2009). [CrossRef]

12.

D. L. Marks, T. S. Ralston, and S. A. Boppart, J. Opt. Soc. Am. A 22, 2366 (2005). [CrossRef]

13.

M. Mihoko and S. Eguchi, Neural Comput. 14, 1859 (2002). [CrossRef]

14.

K. H. Y. Cheng, C. Sun, B. Vuong, K. K. C. Lee, A. Mariampillai, T. R. Marotta, J. Spears, W. J. Montanera, P. R. Herman, K. Rasmus, B. A. Standish, and V. X. D. Yang, Biomed. Opt. Express 3, 388 (2012). [CrossRef]

15.

A. Ozcan, A. Bilenca, A. E. Desjardins, B. E. Bouma, and G. J. Tearney, J. Opt. Soc. Am. A 24, 1901 (2007). [CrossRef]

OCIS Codes
(100.2000) Image processing : Digital image processing
(100.2980) Image processing : Image enhancement
(100.3190) Image processing : Inverse problems
(170.1610) Medical optics and biotechnology : Clinical applications
(170.4500) Medical optics and biotechnology : Optical coherence tomography

ToC Category:
Medical Optics and Biotechnology

History
Original Manuscript: March 28, 2012
Revised Manuscript: May 25, 2012
Manuscript Accepted: May 29, 2012
Published: July 11, 2012

Virtual Issues
Vol. 7, Iss. 9 Virtual Journal for Biomedical Optics

Citation
Kyle H. Y. Cheng, Edmund Y. Lam, Beau A. Standish, and Victor X. D. Yang, "Speckle reduction of endovascular optical coherence tomography using a generalized divergence measure," Opt. Lett. 37, 2871-2873 (2012)
http://www.opticsinfobase.org/ol/abstract.cfm?URI=ol-37-14-2871


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. A. F. Low, G. J. Tearney, B. E. Bouman, and I. Jang, Nat. Clin. Pract. Cardiovasc. Med. 3, 154 (2006). [CrossRef]
  2. A. Mariampillai, B. A. Standish, E. H. Moriyama, M. Khurana, N. R. Munce, M. K. K. Leung, J. Jiang, A. Cable, B. C. Wilson, I. A. Vitkin, and V. X. D. Yang, Opt. Lett. 33, 1530 (2008). [CrossRef]
  3. A. Mariampillai, M. K. K. Leung, M. Jarvi, B. A. Standish, K. Lee, B. C. Wilson, A. Vitkin, and V. X. D. Yang, Opt. Lett. 35, 1257 (2010). [CrossRef]
  4. M. Pircher, E. Gotzinger, R. Leitgeb, A. F. Fercher, and C. K. Hitzenberger, J. Biomed. Opt. 8, 565 (2003). [CrossRef]
  5. A. E. Desjardins, B. J. Vakoc, W. Y. Oh, S. M. R. Motaghiannezam, G. J. Tearney, and B. E. Bouma, Opt. Express 15, 6200 (2007). [CrossRef]
  6. J. W. Goodman, Speckle Phenomena in Optics (Roberts & Company, 2007).
  7. D. C. Adler, T. H. Ko, and J. G. Fujimoto, Opt. Lett. 29, 2878 (2004). [CrossRef]
  8. P. Puvanathasan and K. Bizheva, Opt. Express 15, 15747 (2007). [CrossRef]
  9. Z. Jian, Z. Yu, L. Yu, B. Rao, Z. Chen, and B. J. Tromberg, Opt. Lett. 34, 1516 (2009). [CrossRef]
  10. J. Rogowska and M. E. Brezinski, IEEE Trans. Med. Imag. 19, 1261 (2000). [CrossRef]
  11. P. Puvanathasan, and K. Bizheva, Opt. Express 17, 733 (2009). [CrossRef]
  12. D. L. Marks, T. S. Ralston, and S. A. Boppart, J. Opt. Soc. Am. A 22, 2366 (2005). [CrossRef]
  13. M. Mihoko and S. Eguchi, Neural Comput. 14, 1859 (2002). [CrossRef]
  14. K. H. Y. Cheng, C. Sun, B. Vuong, K. K. C. Lee, A. Mariampillai, T. R. Marotta, J. Spears, W. J. Montanera, P. R. Herman, K. Rasmus, B. A. Standish, and V. X. D. Yang, Biomed. Opt. Express 3, 388 (2012). [CrossRef]
  15. A. Ozcan, A. Bilenca, A. E. Desjardins, B. E. Bouma, and G. J. Tearney, J. Opt. Soc. Am. A 24, 1901 (2007). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited