Journal of Healthcare Engineering
Volume 2018 (2018), Article ID 7329548, 7 pages
https://doi.org/10.1155/2018/7329548

Research Article

Efficient OCT Image Enhancement Based on Collaborative Shock Filtering

Guohua Liu,¹ Ziyu Wang,² Guoying Mu,³ and Peijin Li⁴

¹Department of Ophthalmology, Qilu Children’s Hospital of Shandong University, Jinan 250022, China
²Department of Radiology, Yidu Central Hospital of Weifang, Qingzhou 262500, China
³Department of Ophthalmology, Shandong Provincial Hospital, Jinan 250021, China
⁴School of Computer Science and Technology, Shandong University, Jinan 250101, China

Correspondence should be addressed to Guoying Mu; mgyeyes@163.com

Received 5 May 2017; Accepted 14 August 2017; Published 1 February 2018

Academic Editor: Chengpo Mu

Copyright © 2018 Guohua Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Efficient enhancement of noisy optical coherence tomography (OCT) images is a key task for interpreting them correctly. In this paper, to better enhance details and layered structures of a human retina image, we propose a collaborative shock filtering for OCT image denoising and enhancement. Noisy OCT image is first denoised by a collaborative filtering method with new similarity measure, and then the denoised image is sharpened by a shock-type filtering for edge and detail enhancement. For dim OCT images, in order to improve image contrast for the detection of tiny lesions, a gamma transformation is first used to enhance the images within proper gray levels. The proposed method integrating image smoothing and sharpening simultaneously obtains better visual results in experiments.

1. Introduction

Optical coherence tomography (OCT) is an important emerging noninvasive technique in medical imaging, which offers high-speed tomographic imaging of human tissues [1, 2]. Layered structures of the human retina provided by OCT imaging are commonly detected and measured for diagnosis and prognostic evaluation [3]. Thus, medical OCT imaging has a variety of clinical applications and is widely used in the ophthalmology for vision-related diseases. Since OCT imaging takes advantage of low-coherence interferometry, it inevitably suffers from speckle noise, visible as dark and bright spots on the imaged objects [4], which severely degrades OCT image contrast and quality and makes it particularly challenging to identify fine structures and features. Therefore, it is necessary to remove the speckle noise and to enhance OCT images in order to interpret OCT images correctly.

There are many approaches to reduce the speckle noise of OCT images for enhancing their image quality [5, 6]. For example, the geometry-oriented calculus of variations and partial differential equations have been widely used in image processing [7, 8]: the total variation regularization [9] and the anisotropic diffusion [10] for image denoising and enhancement, the shock filters [11, 12] for image sharpening, and so forth. However, in the process of image denoising and enhancement, these classic geometrical regularization methods, based on operators in differential geometry such as gradient, divergence, and directional derivative, often tend to modify the image towards a piecewise constant function and blur fine features of the image, particularly the image’s details and layered structures.

Other image enhancement methods include wavelet decomposition [13], local statistical [14], and self-similarity and sparsity-based methods [15–17]. These algorithms can reduce noise to a certain extent, but they also lead to the loss of image details and consequently the blur of important image features. Moreover, some algorithms are fairly time-consuming for them to enhance a noisy OCT image.

In 2005, Buades et al. proposed a nonlocal means (NLM) filter based on the self-similarity in the whole image, which can effectively preserve image details and textures [18, 19]. The nonlocal means filter produces impressive results in denoising textured patterns. A more complicated neighborhood filter based on the self-similarity of image blocks was introduced by Dabov et al. [20], which is actually also a weighted average algorithm, just after a transform-domain collaborative filtering. This algorithm achieves better denoising performance in terms of peak signal-to-noise ratio [21].

In this paper, to better preserve or enhance image details and layered structures, we propose an effective enhancement algorithm for OCT image denoising and enhancement based on collaborative shock filtering. The basic idea is that noisy OCT image is first denoised with a collaborative filtering of speckle noise following a gamma distribution, and then the denoised image is sharpened by a shock-type filtering for edge and detail enhancement. For dim OCT image, in order to improve the contrast of the retinal image for the detection of tiny lesions, a gamma transformation [22] is first used to enhance the image within proper gray levels.

Image smoothing and sharpening are two opposite operations in image processing. Generally speaking, image smoothing is to eliminate unnecessary and false discontinuous features (such as noise), while image sharpening is to produce or enhance some discontinuous features (such as edges and details) in proper positions of the image. In many cases, image denoising methods often inevitably blur image edges and details even if these methods are designed elaborately. For example, see following denoising experiments on noisy images. Thus, both image smoothing and sharpening are needed simultaneously even if only in an image denoising task. One will see that the methods integrating the above two operations simultaneously obtain better visual results in experiments [23].

The rest of this paper is organized as follows. In Section 2, we briefly review some classic works: shock filters and self-similarity filtering. Then, a collaborative shock-filtering algorithm integrating these classic methods is proposed in Section 3. Experiments on test images are described in Section 4. Finally, conclusions are drawn out in Section 5.

2. Shock Filters and Self-Similarity Filtering

In this section, we review some classic methods in image denoising and enhancement.

2.1. Shock Filter

In [24], some special ideas and techniques developed in numerical solutions of nonlinear hyperbolic equations were applied to feature-oriented image enhancement for the first time. Osher and Rudin introduced a novel image sharpening technique called the shock filter (SF) [11], which is based on a modification of the nonlinear Burgers’ equation and simulates the shock wave calculation in the computational fluid mechanics [25]. Different from the nonlinear parabolic equation of diffusion-type process, they proposed a hyperbolic one:

where u is the observed noisy image, sign is a sign function, and is the second directional derivative of the image along a local normal direction to the isophote line. It detects an image edge using the zero crossing of , where a shock is formed at the speed of .

Considering the image noise in the estimation of edges, Alvarez and Mazorra added a smoothing kernel and coupled the anisotropic diffusion with the shock filter (ADSF) [12] for noise elimination and edge sharpening:

where is a Gaussian kernel with standard deviation σ, is the second directional derivative of the image along a local tangent direction, and c is a constant to balance the anisotropic diffusion and the shock filter.

2.2. Self-Similarity Filtering

Another effective image denoising technique is neighborhood filters [8], which consider self-similarity between two pixels or blocks of the image both in spatial location and in gray level. Buades et al. proposed the following nonlocal means (NLM) filter [18, 19]:

where is a normalization factor, and h is a filtering parameter related to noise level, and are 2D spatial coordinates. The nonlocal means filter gives better results in denoising textured patterns.

A more complicated neighborhood filter called the Block-matching and 3D filtering (BM3D) algorithm was introduced by Dabov et al. [20], which has three main steps: grouping of similar image blocks, 3D collaborative filtering of these blocks in the spectrum domain, and aggregation of all local estimates. This algorithm includes basic and final estimations using the following weighted average:

where R and S denote the reference block and the similar block, respectively, w is the weight of the corresponding block, is the local block-wise estimate, and is a characteristic function. As they claim, this algorithm is currently one of best denoising methods in terms of peak signal-to-noise ratio [20, 21].

Finally, as one will see, the neighborhood filters are mainly used for image denoising, while the shock filters can be used for both image denoising and enhancement simultaneously, because of their inherent local backward diffusion [25].

3. Collaborative Shock Filtering

The main information of an image is encoded in its edges, details, and structures. These components need to be processed in special ways to enhance the image as well as possible.

In order to better preserve or enhance image details and layered structures, we propose an efficient collaborative shock-filtering (CSF) algorithm for image denoising and enhancement, which fuses different advantages of previously mentioned classic methods: powerful capabilities of the self-similarity filtering in image edge denoising and the shock filters in edge sharpening.

To be specific, as illustrated in Figure 1, for a given degraded OCT image, it is first denoised with a collaborative filtering method; then, the denoised image is sharpened by a shock-type filtering for edge and detail enhancement [11]. For dim OCT images, in order to improve the contrast of the retinal image for the detection of tiny lesions, a gamma transformation [22] is first used to enhance the image within proper gray levels.

Figure 1: A flowchart of the proposed collaborative shock-filtering algorithm. A noisy OCT image is enhanced through three steps in sequence: gamma manipulation (gamma), collaborative filtering, and shock filter, respectively.

Different from the BM3D algorithm for Gaussian noise removal [20], the block similarity measure is elaborated in the proposed collaborative filtering of speckle noise following a gamma distribution [26]. Based on a generalized likelihood ratio, the similarity criteria is defined as follows [27]:

Here, is the Hadamard product of two image matrixes (blocks) and is the function taking the mean of all elements of matrices. The larger value of S implies that image block is more similar to .

4. Experimental Results

In this section, we perform a lot of enhancement experiments to verify our method on various degraded OCT images. For the above methods, numerical schemes presented in original works are used.

All methods are implemented using the MATLAB programming on OCT images with grayscale ranging from 0 to 1. In the comparison of different methods, parameters in each method are selected such that optimal visual results are obtained. Through three sequential processings of the gamma manipulation (), the collaborative filtering () and the shock filter (), where n and denote the iteration times and the temporal step size, tiny details and layered structures are shown clearly due to image contrast improvement and noise removal. Moreover, our method produces fewer overshoot artifacts while avoiding noise magnification.

In Figure 2, a severely noisy OCT image of size 183 × 894 is enhanced without image gray adjustment (gamma transformation). Its better effects with clearer layered structures and effective noise removal can be observed more clearly in zoomed parts of results in Figure 3, compared with other methods such as ADSF, NLM, and BM3D. The ADSF method does not efficiently remove speckle noise and produces blurry edges and layered structures with some annoying artifacts in flat regions. The NLM method produces very blurred layered structures with broken nonsmooth linear edges. The BM3D method produces better enhancement results with a few slightly blurry edges and layered structures. Blurry and nonsmooth-layered structures are bad for one to measure the thickness of layered structures to predict and evaluate related disorders in ophthalmology.

(a)

(b)

(c)

(d)

(e)

Figure 2: Enhancement of noisy OCT image: (a–e) original image, ADSF, NLM, BM3D, and the proposed collaborative shock filtering, respectively.

(a)

(b)

(c)

(d)

(e)

Figure 3: Zoomed parts in Figure 2: (a–e) original image, ADSF, NLM, BM3D, and the proposed collaborative shock filtering, respectively.

In order to observe enhancement effects by these methods more clearly, local profiles (500th column) of different results are shown in Figure 4. One can see that, compared with related methods, the proposed method removes speckle noise more effectively preserving layered structures without producing annoying artifacts, providing a chance to detect and evaluate retinal lesions faithfully by image enhancement. In Figure 5, local profiles of a flat region (line 55–90, 500th column) and an edge region (line 105–135, 500th column) of different results are shown. Superior denoising and edge sharpening by the proposed CSF method can be clearly observed: sharper edges with a smaller edge width are helpful for an accurate measure of layered structures.

Figure 4: Enhancement of noisy OCT image: comparison of profiles (500th column) of original and enhanced images by ADSF, NLM, and the proposed collaborative shock filtering (CSF), respectively.

(a)

(b)

Figure 5: Enhancement of noisy OCT image: comparison of profiles of different regions of original and enhanced images by ADSF, NLM, and the proposed collaborative shock filtering (CSF): (a) flat region (line 55–90, 500th column) and (b) edge region (line 105–135, 500th column), respectively.

In Figure 6, a very dim and noisy OCT image of size 194 × 477 is enhanced with image gray adjustment. Clearer layered structures can be easily observed with effective noise removal.

(a)

(b)

(c)

(d)

Figure 6: Enhancement of dim noisy OCT image: (a–d) original image, gamma transformation, BM3D, and the proposed collaborative shock filtering, respectively.

In Figures 7–9, three much noisy OCT images of sizes (245 × 586, 187 × 341, and 250 × 453) are enhanced without image gray adjustment. Clearer layered structures can be easily observed with effective noise removal.

(a)

(b)

(c)

Figure 7: Enhancement of noisy OCT image: (a–c): original image, BM3D, and proposed collaborative shock filtering, respectively.

(a)

(b)

(c)

Figure 8: Enhancement of dim noisy OCT image (a–c) original image, BM3D, and the proposed collaborative shock filtering, respectively.

(a)

(b)

(c)

Figure 9: Enhancement of dim noisy OCT image: (a–c): original image, BM3D, and the proposed collaborative shock filtering, respectively.

(a)

(b)

(c)

Figure 10: Enhancement of blurred noisy OCT image: (a–c) original image, BM3D, the and proposed collaborative shock filtering, respectively.

(a)

(b)

Figure 11: Enhancement of noisy OCT images: (a) results by BM3D and (b) results by the proposed collaborative shock filtering, respectively.

Figure 12: Enhancement of noisy OCT image (first test image in Figure 2): comparison of profiles (line 105–135, 500th column) of enhanced images by BM3D and the proposed collaborative shock filtering (CSF), respectively.

In Figure 10, a blurred noisy OCT image of size 594 × 1263 is enhanced without image gray adjustment. In this case, more sharpening with the shock filter () is used to ensure clearer layered structures and effective noise removal.

Finally, we compare the proposed CSF method with the BM3D method for edge sharpening in OCT image enhancement. In Figure 11, typical results of previous experiments by two methods are shown to carefully observe the enhancement of layered structures in OCT images. Obviously, the CSF method produces sharper edges and clearer layered structures around the retinal fovea. In Figure 12, moreover, it is clearer for one to observe profiles of an edge region (line 105–135, 500th column) in enhancing the first test image in Figure 2: the proposed CSF method produces sharper edges with smaller edge width, which is consistent with the visual observation.

5. Conclusions and Future Work

To better enhance image details and layered structures of noisy optical coherence tomography (OCT) image, we propose a collaborative shock filtering for OCT image denoising and enhancement. Noisy OCT image is first denoised with a collaborative filtering method based on a new similarity measure, followed by a sharpening step by a shock-type filtering for edge and detail enhancement. Finally, in order to improve the contrast of dim OCT image, a gamma transformation is used to enhance the image within proper gray levels. Simultaneously, integrating the image smoothing and sharpening of the proposed method obtains better visual results in image enhancement experiments.

In future work, we will develop robust measurement tools for quantitative analysis of vision-related diseases based on this work.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The research has been supported in part by the National Natural Science Foundation of China (61272239 and 61671276) and the Science and Technology Development Project of Shandong Province of China (2014GGX101024).

References

D. Huang, E. A. Swanson, C. P. Lin et al., “Optical coherence tomography,” Science, vol. 254, no. 5035, pp. 1178–1181, 1991.

J. M. Schmitt, “Optical coherence tomography (OCT): a review,” IEEE Journal of Selected Topics in Quantum Electronics, vol. 5, no. 4, pp. 1205–1215, 1999.

J. G. Fujimoto, W. Drexler, J. S. Schuman, and C. K. Hitzenberger, “Optical coherence tomography (OCT) in ophthalmology: introduction,” Optics Express, vol. 17, no. 5, pp. 3978-3979, 2009.

J. M. Schmitt, S. Xiang, and K. M. Yung, “Speckle in optical coherence tomography,” Journal of Biomedical Optics, vol. 4, no. 1, pp. 95–105, 1999.

H. M. Salinas and D. C. Fernández, “Comparison of PDE-based nonlinear diffusion approaches for image enhancement and denoising in optical coherence tomography,” IEEE Transactions on Medical Imaging, vol. 26, no. 6, pp. 761–771, 2007.

P. Sudeep, S. I. Niwas, P. Palanisamy et al., “Enhancement and bias removal of optical coherence tomography images: an iterative approach with adaptive bilateral filtering,” Computers in Biology and Medicine, vol. 71, pp. 97–107, 2016.

T. Chan and J. Shen, Image Processing and Analysis: Variational, PDE, Wavelet, and Stochastic Methods, the Society for Industrial and Applied Mathematics, Philadelphia, PA, USA, 2005.

G. Aubert and P. Kornprobst, Mathematical Problems in Image Processing: Partial Differential Equations and the Calculus of Variations, Springer-Verlag, New York, NY, USA, 2nd edition, 2006.

L. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physica D: Nonlinear Phenomena, vol. 60, no. 1-4, pp. 259–268, 1992.

P. Perona and J. Malik, “Scale-space and edge detection using anisotropic diffusion,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 12, no. 7, pp. 629–639, 1990.

S. Osher and L. Rudin, “Feature-oriented image enhancement using shock filters,” SIAM Journal on Numerical Analysis, vol. 27, no. 4, pp. 919–940, 1990.

L. Alvarez and L. Mazorra, “Signal and image restoration using shock filters and anisotropic diffusion,” SIAM Journal on Numerical Analysis, vol. 31, no. 2, pp. 590–605, 1994.

Z. Jian, Z. Yu, L. Yu, B. Rao, Z. Chen, and B. J. Tromberg, “Speckle attenuation in optical coherence tomography by curvelet shrinkage,” Optics Letters, vol. 34, no. 10, pp. 1516–1518, 2009.

A. Wong, A. Mishra, K. Bizheva, and D. A. Clausi, “General bayesian estimation for speckle noise reduction in optical coherence tomography retinal imagery,” Optics Express, vol. 18, no. 8, pp. 8338–8352, 2010.

L. Fang, S. Li, Q. Nie, J. A. Izatt, C. A. Toth, and S. Farsiu, “Sparsity based denoising of spectral domain optical coherence tomography images,” Biomedical Optics Express, vol. 3, no. 5, pp. 927–942, 2012.

H. Yu, J. Gao, and A. Li, “Probability-based non-local means filter for speckle noise suppression in optical coherence tomography images,” Optics Letters, vol. 41, no. 5, pp. 994–997, 2016.

D. Thapa, K. Raahemifar, and V. Lakshminarayanan, “Reduction of speckle noise from optical coherence tomography images using multi-frame weighted nuclear norm minimization method,” Journal of Modern Optics, vol. 62, no. 21, pp. 1856–1864, 2015.

A. Buades, B. Coll, and J. Morel, “A review of image denoising algorithms, with a new one,” Multiscale Modeling & Simulation, vol. 4, no. 2, pp. 490–530, 2005.

A. Buades, B. Coll, and J.-M. Morel, “Nonlocal image and movie denoising,” International Journal of Computer Vision, vol. 76, no. 2, pp. 123–139, 2008.

K. Dabov, A. Foi, V. Katkovnik, and K. Egiazarian, “Image denoising by sparse 3D transform-domain collaborative filtering,” IEEE Transactions on Image Processing, vol. 16, no. 8, pp. 2080–2095, 2007.

V. Katkovnik, A. Foi, K. Egiazarian, and J. Astola, “From local kernel to nonlocal multiple-model image denoising,” International Journal of Computer Vision, vol. 86, no. 1, pp. 1–32, 2010.

R. Gonzalez and R. Woods, Digital Image Processing, Prentice-Hall, New York, NY, USA, 2nd edition, 2004.

L. Wang, G. Liu, S. Fu, L. Xu, K. Zhao, and C. Zhang, “Retinal image enhancement using robust inverse diffusion equation and self-similarity filtering,” PLoS One, vol. 11, no. 7, article e0158480, 2016.

L. Rudin, Shock Filters, Tech. Rep., Defence Advanced Research Projects Agency (DARPA), Rockwell International Science Center, United States, Thousand Oaks, 1984.

S. Fu, C. Zhang, and X. Tai, “Image denoising and deblurring: non-convex regularization, inverse diffusion and shock filter,” Science China Information Sciences, vol. 54, no. 6, pp. 1184–1198, 2011.

D. A. Jesus and D. R. Iskander, “Assessment of corneal properties based on statistical modeling of OCT speckle,” Biomedical Optics Express, vol. 8, no. 1, pp. 162–176, 2017.

C.-A. Deledalle, L. Denis, and F. Tupin, “How to compare noisy patches? Patch similarity beyond Gaussian noise,” International Journal of Computer Vision, vol. 99, no. 1, pp. 86–102, 2012.