 Research article
 Open access
 Published:
Nonlocal total variation based on symmetric KullbackLeibler divergence for the ultrasound image despeckling
BMC Medical Imaging volume 17, Article number: 57 (2017)
Abstract
Background
Ultrasound imaging is safer than other imaging modalities, because it is noninvasive and nonradiative. Speckle noise degrades the quality of ultrasound images and has negative effects on visual perception and diagnostic operations.
Methods
In this paper, a nonlocal total variation (NLTV) method for ultrasonic speckle reduction is proposed. A spatiogram similarity measurement is introduced for the similarity calculation between image patches. It is based on symmetric KullbackLeibler (KL) divergence and signaldependent speckle model for logcompressed ultrasound images. Each patch is regarded as a spatiogram, and the spatial distribution of each bin of the spatiogram is regarded as a weighted Gamma distribution. The similarity between the corresponding bins of the two spatiograms is computed by the symmetric KL divergence. The SplitBregman fast algorithm is then used to solve the adapted NLTV object function. KolmogorovSmirnov (KS) test is performed on synthetic noisy images and real ultrasound images.
Results
We validate our method on synthetic noisy images and clinical ultrasound images. Three measures are adopted for the quantitative evaluation of the despeckling performance: the signaltonoise ratio (SNR), structural similarity index (SSIM), and natural image quality evaluator (NIQE). For synthetic noisy images, when the noise level increases, the proposed algorithm achieves slightly higher SNRS than that of the other two algorithms, and the SSIMS yielded by the proposed algorithm is obviously higher than that of the other two algorithms. For liver, IVUS and 3DUS images, the NIQE values are 8.25, 6.42 and 9.01, all of which are higher than that of the other two algorithms.
Conclusions
The results of the experiments over synthetic and real ultrasound images demonstrate that the proposed method outperforms current stateoftheart despeckling methods with respect to speckle reduction and tissue texture preservation.
Background
Ultrasound imaging is one of the four modern medical imaging techniques that deliver crosssectional images of a patient’s anatomy and physiology in real time [1]. The use of ultrasound in the diagnosis and assessment of imaging organs and soft tissue structures, as well as human blood, is well established [2]. This technique is progressively achieving an important role in the assessment and characterization of cardiac imaging because it is noninvasive and has a strong ability to identify biological soft tissues [3]. However, ultrasound imaging generates images with low contrast resolution, because of speckle noise.
The speckle pattern, which contains typical light and dark spots in an image, is generated by the interference effect of the ultrasonic echoes and scattering of randomly distributed structure scatters. In fact, speckle is technically not a noise in the typical engineering sense because its texture often carries useful information on the image. Therefore, distinguishing between speckle generated from tissue and that from the received RF signal is essential [3]. The former refers to the image texture and the latter, the speckle noise. Speckle noise is a form of multiplicative noise subjected to Gamma distribution [4,5,6]. It is the primary factor that limits the contrast resolution in diagnostic ultrasound imaging, thereby limiting the detectability of small lowcontrast lesions and rendering the ultrasound images generally difficult for nonspecialist to interpret [3]. Therefore, reducing speckle noise is important for the interpretation and computeraided analysis of ultrasound images. The aim of despeckling is to eliminate the speckle noise and maintain the image boundaries and image texture.
Decreasing the frequency of the transducer can increase the signalnoiseratio (SNR) in ultrasound images. However, doing so decreases the resolution of the image, and consequently hinders the capturing of internal details of human organs. Meanwhile, highresolution ultrasound images usually indicate low SNR, which can be increased by despeckling post processing techniques. One of the classical despeckling methods is the Wiener filter, which is a linear filter based on local statistics [7]. It uses the weighted average calculation of subregion statistics to estimate the statistical measures over different pixel windows [3]. However, the linear filters have an inherent drawback of oversmoothing anatomical details and transitional boundaries [8]. Thus, in the last two decades, researchers have focused on nonlinear filtering methods to remove speckle noise in ultrasound images.
Similar to ultrasound images, the synthetic aperture radar (SAR) images are usually degraded by multiplicative speckle noise. Some local speckle statisticsbased adaptive filters proposed for SAR images are also used in ultrasound images [9,10,11,12,13,14,15]. These adaptive filters use local speckle statistics to improve smoothness in homogeneous regions where the speckle is fully developed and reduce smoothness appreciably in the other regions showing the useful details of an anatomical structure [3]. Adaptive filters include the Lee filter [9], Kuan filter [10], Frost filter [11], adaptive weighted median filter [12], speckle reducing bilateral filter [13], maximum likelihood filer [14] and MRGMAP filter [15]. These adaptive filters use different weighting coefficients calculated from the subregion statistics over different pixel windows to preserve important image structures while removing speckle noise. However, although these adaptive filters can retain the boundaries, they cannot eliminate speckle noise effectively and are sensitive to the size of the pixel window. When the window size is large, the image edges are blurred.
Anisotropic diffusion, proposed by [16], is an efficient nonlinear technique that can enhance contrast and reduce noise simultaneously [17,18,19]. It smooths homogeneous image regions and preserves image edges without requiring any information from the image power spectrum [3]. For example, the specklereducing anisotropic diffusion filtering has been proposed to enhance image contrast [3]. Wavelet filtering is another method that exploits the decomposition of an image into the wavelet basis and zeros out the wavelet coefficients to despeckle the image. Specklereduction filtering in the wavelet domain is based on the concept of the Daubechies Symlet wavelet and softthresholding denoising [3]. It was proposed first by Donoho [20] and investigated further by Zhong and Cherkassky [21] and Gupta et al. [22]. However, although diffusion and wavelet filters can suppress speckle noise considerably, they produce excessive smoothing details, particularly in the image boundaries and texture.
Nonlocal denoising methods that use the similarities among small patches have been reckoned to be an effective way to preserve details while decreasing noise [6, 23, 24]. A simple and classic method is the wellknown nonlocal means (NLM) filter [23], which measures the similarities among the patches centered at a given pixel in the image and prevents the smoothing of boundaries by assigning only high weights to pixels with similar local properties. However, it cannot be applied to speckle filtering directly since the speckle model of the ultrasound images is not subjected to the Gaussian distribution. The NLM filter has been adapted to Gammadistributed speckle pattern by introducing a Bayesian estimator to the weighting function and Rayleighdistributed speckle pattern through the introduction of maximumlikelihood estimation [25, 26]. These adapted methods demonstrate better results than the original NLM filter. The NLM filter generates better structural preservation than the diffusion filters. However, it also generates excessive smoothing textures. Based on the NLM filter, the nonlocal total variation (NLTV) minimization scheme analyzes the patterns around the pixels and performs pattern matching in the global image [27, 28]. The nonlocal total variation norm processes textures and repetitive structures effectively. However, although the NLTV filter performs well in Gaussian noise reduction and sharp boundaries preservation, it cannot be applied to logcompressed ultrasound images directly, because the speckle is not subjected to the Gaussian distribution. Thus, a Bayesian NLTV speckle filter (BNLTV) has been developed for ultrasound images and thus is capable of improving speckle suppression and edge enhancement, outperforming the adaptive filters, diffusion filters, wavelet filters, nonlocal denoising methods, and original NLTV filter [29]. However, the BNLTV filter oversmooth image texture and introduces some artificial traces, while it eliminates the speckle noise. It cannot determine whether or not the speckle received from RF signal.
Another nonlocal denoising method is the blocking matching 3D (BM3D) filter proposed by Dabov et al. [30]. This method combines nonlocal and transformdomain approaches and presents an effective denoising framework by grouping similar patches into a 3D array. It subsequently filters the 3D array by using sparse representation in the transform domain. Finally, it aggregates multiple estimates at each location. Based on the BM3D filter, the SARBM3D filter [31] was proposed for despeckling SAR images. The SARBM3D filter is generally used to despeckle breast ultrasound images and exhibited the impressive despeckling ability [32].
In this paper, an adapted NLTV speckle filter is proposed to address the problem on speckle noise. This filter can eliminate the speckle noise while maintaining the edges and image texture. For the improvement of the despeckling performance, spatiogram similarity measurement based on symmetric KL divergence is introduced to the similarity calculation between the image patches. The proposed method is validated by considering both real ultrasound images and synthetic noisy images. Before the experiment, a KolmogorovSmirnov (KS) test has been conducted on the experimental images to ensure that the images are subjected to Gamma distribution [5]. For the evaluation of the despeckling performance without the ground truth image, the natural image quality evaluator (NIQE) [33] is introduced. NIQE can predict the quality of distorted images (noise, ringing, blur, or blocking), even without any prior knowledge of the reference images or their distortions. The remainder of this paper is organized as follows: Section 2 provides the proposed NLTV filtering procedure. Section 3 reports the experimental results. Finally, Section 4 presents the conclusions.
Methods
NLTV algorithm
Based on nonlocal denoising methods, NLTV denoising is generally designed for the zero mean Gaussian noise. Dong et al. [28] and Wen et al. [29] proposed a nonlocal total variation noise removal model for multiplicative noise. In these models, the noisy image Y from a noisefree image X can be modeled as follows [28, 29]:
where η generally denotes a zero mean Gaussian noise with a variance σ ^{2}. γ is a factor that relies on ultrasound devices and additional image reconstruction processing, and when γ = 0.5, the formation process for the logcompressed ultrasound image can be approximated [12]. Denoising methods aim to restore X or find a good estimate \( \widehat{X} \) of X from Y.
The NLTV denoising method searches similar patches in a certain searching window centered at a given pixel and obtains the weight coefficients according to the similarities among the patches. The restored pixel is reconstructed by minimizing nonlocal TV model. This operation is repeated for each pixel in the sliding searching window fashion in the whole noisy image. The details of the NLTV denoising process are described in the subsequent text.
Given a p × p patch and a iw × iw searching window centered at a certain pixel, i, in the noisy image, m similar patches (including the given patch) are found. The distance of two patches can be formulated as follows [28]:
where G _{ α } is a Gaussian kernel with standard deviation α, X(i + ∙) and X(j + ∙) are the patches centered at i and j, respectively, in a given searching window. The similarity coefficients can be calculated as follows [28]:
where h is a filtering parameter. The smaller the distance, the greater the w(i, j), which means the two patches is more similar.
Then the NL gradient is defined based on the weighting of similar patches, which is determined as follows [28]:
Based on the Bayesian rational [29] and Eq. (4), the restored gray value of i can be obtained by minimizing the nonlocal TV model represented as follows [28]:
where λ > 0 and is a Lagrange parameter balancing the influence between the nonlocal regularization term ∑_{ i }∇ _{ NL } X _{ i } and data fidelity term.
The above operation is performed in each pixel in the sliding search window fashion in the whole noisy image. In this operation, an efficient minimizing algorithm is necessary to improve the performance at a relatively fast rate. Many efficient algorithms for the computation of the nonlocal TV problem have been proposed [34, 35]. Of these algorithms, the split Bregman method can solve the nonlocal TV minimization efficiently. It simplifies the nonlocal TV problem to two subproblems calculated by the Gauss–Seidel iterative scheme and softthresholding formula [34].
Symmetric KL divergence NLTV speckle filter
The NLTV method was originally proposed for additive Gaussian noise removal [27]. It uses the nonrobust Euclidean distance to measure the similarities among the patches. To improve the performance of the NLTV despeckling method, we propose to replace the similarity measurement based on Euclidean distance with a spatiogram similarity measurement based on symmetric KL divergence. This similarity measurement is performed on the noisy patches to determine the similarities between each pair of patches. Each patch is regarded as a spatiogram, and the spatial distribution of each bin of a spatiogram is regarded as a weighted Gamma distribution, where the weight is the probability of the corresponding bin; subsequently, the similarity in the corresponding bins of the two spatiograms is computed by the symmetric KL divergence with two weighted Gamma distributions [36].
Given a spatiogram,
where n _{ b } is the probability value of the bin, μ _{ b } and Σ_{ b } are the mean and covariance matrix, respectively, of the pixel coordinates in the bin, B is the number of bins. A weighted Gamma distribution is represented as follows:
where
where α = μ ^{2}/Σ, β = μ/Σ. The KL distance between the two weighted Gamma distributions, f _{ b }(x) and \( {f}_b^{\hbox{'}}(x) \), can be calculated as follows:
where
where d is the space dimensionality (for spatiogram, d = 2). The symmetric KL distance between the two weighted Gamma distributions is represented as follows:
Then, the weight coefficient between the two patches can be calculated as follows:
SplitBregman implementation
After obtaining the weight coefficients by Eq. (12), we use the Split Bregman iteration [28] to minimize the adapted NLTVfunctional Eq. (5).
To prevent the singularity and complexity of the numerical difference, Goldstein and Osher [35] proposed an alternating minimization scheme. They introduced an auxiliary variable to approximate the image gradient. The NLTV functional can be modified as the following minimization problem by introducing the auxiliary variable d instead of ∇ _{ NL } X, such that [28]
The above constrained problem can be converted into an unconstrained problem by using the quadratic penalty method [28].
Then, the Bregman iteration is used to solve the above minimization [28].
Where β is a positive constant and b is the Bregman parameter.
Apparently, Eq. (15) is a multivariate function and can be minimized by alternatively solving the two minimization subproblems with respect to X and d. Fixed d, the minimization over X is [28]
The optimal solution of X satisfies the following equation using the EulerLagrange formula [28]:
To get a fast solution for Eq. (17), we use the GaussianSeidel iteration, and the solution is represented as [28]
Fixed X, the minimization over d is [28]
According to the softthresholding formula [28], the solution of Eq. (19) is provided by
Finally, the Bregman variable b is updated as follows [28]:
The algorithm of NLTV model is summarized in Algorithm 1. The procedure consists of two steps: the NL weights computation step and NLTV minimization step. In the NL weights computation step, the KL distance is calculated using Eq. (11). It is then used to calculate the NL weights using Eq. (12). After the NL weights are calculated, the NLTV minimization step is performed using the SplitBregman scheme. In this algorithm, n represents the number of GaussianSeidel iteration and is used to determine a good approximation in Eq. (18). Meanwhile, k is the number of the overall iteration.
Experiments
The synthetic noisy images and clinical ultrasound images from [29] are studied to assess and compare the performances of the despeckling methods quantitatively. In [29], a 2D synthetic image “Modified SheppLogan Phantom” available in MATLAB is considered and corrupted with different noise levels. The speckle simulation for the synthetic image is based on the noisy model, the Eq. (1). Three levels of noise are tested by setting standard deviations sigma = {0.8,1.0,1.3}.The clinical ultrasound images include liver, intravascular, and 3D ultrasound images. Many experiments performed on logcompressed ultrasound images showed that the speckle noise is better described by the Gamma distribution [4,5,6, 37]. A simple experiment is then conducted to determine whether or not the speckle noise is subjected to Gamma distribution. In this experiment, the KolmogorovSmirnov (KS) test is performed on the experimental images [5]. For each image, the comparison between Gaussian, Gamma, and Rayleigh distribution is done. The formula of the KS test is as follows [5]:
where \( \widehat{F_n} \) is the empirical cumulative distribution function (CDF) of the observed images, F _{ X } is the CDF of either Gaussian, Gamma, − Rayleigh, or FisherTippett distribution. The GlivenkoGantelli theorem states that, if the samples are subjected to distribution F _{ X }, then D _{ KS } converges to 0 [38]. In Table 1, the values in the Gamma column are smaller than the other values, and this means the real ultrasound images and synthetic noisy image are tended to subject to the Gamma distribution.
Three measures are adopted for the quantitative evaluation of the despeckling performance: the signaltonoise ratio (SNR), structural similarity index (SSIM), and natural image quality evaluator (NIQE). The SNR is established on the variance ratio of the effective signal and noise using the despeckled and noisefree image such that
where M is the total number of pixels in the noisy image, \( \widehat{X_i} \) and X _{ i } are the restored value and true value at pixel i, respectively. The SSIM is used to measure the structural similarity between the images [39, 40], and is defined as follows:
where μ _{ X } and \( {\mu}_{\widehat{X}} \) are the mean of the ground truth X and denoised image \( \widehat{X} \), respectively, \( {\sigma}_{X\widehat{X}} \) is the covariance of X and \( \widehat{X} \), and c _{1} and c _{2} are constant parameters. In this research, the SSIM is locally measured with an 8 × 8 Gaussian kernel, and the mean SSIM (MSSIM) is estimated as a global SSIM by all the local SSIMs. The NIQE [33] is based on the construction of a “quality aware” collection of statistical features, which are based on a simple and successful space domain natural scene statistic model, represented as
where v _{1}, v _{2} and Σ_{1}, Σ_{2} are the mean vectors and covariance matrices, respectively, of the natural Multivariate Gaussian (MVG) model and MVG model of a distorted image.
The proposed algorithm is compared with two adaptive algorithms, namely, the BNLTV filter [29], which is based on Bayesian framework, and SARBM3D filter [31] to provide relevant comparisons. The results of the [29] are provided by the author. The parameters of proposed filter and SARBM3D filter are adjusted to obtain the best SNR scores. The detailed parameters for the proposed filter and SARBM3D filter are listed in Table 2. All the algorithms are implemented with MATLAB R2012a, and the computer is equipped with 2.80 GHz CPU (10core, E5) and 64 GB RAM.
Filter parameters
The performance of the proposed despeckled method depends on the setting of the two patchrelated parameters (i.e., the patch size p and the number of similar patches in a search window m) and two nonlocal total variation schemerelated parameters. The parameters are the Lagrange parameter λ, which balances the action between the nonlocal regularization and data fidelity terms, and the quadratic penalty parameter β, which must be large enough to ensure d is sufficiently close to ∇ _{ NL } X. In our despeckling experiments with synthetic noisy images, p and m are selected as 5 and 10, respectively, which yield higher denoising SNR. For the real ultrasound images, the proposed algorithm is implemented with parameters of p and m with values of 7 and 10, respectively. The suitable values for λ and β are then ascertained by implementing the proposed algorithm on the synthetic noisy images with varying λ and β values. The SNR values are used as quantitative indicators to evaluate the sensitivity to the setting of these two parameters. For the proposed algorithm, the curves as a function of the λ and β parameters are shown in Fig. 1. These two parameters have an apparent effect on the speckle suppression and edge sharpness of the despeckling image. As shown in Fig. 1, the speckle cannot be reduced because of the small Lagrange parameter λ and quadratic penalty parameter β. Meanwhile, large λ and β values result in a blurred edge transition region. Only the proper parameters can balance the performance between speckle reduction and boundary preservation. Based on the observation from Fig. 1, λ and β are fixed at 30 and 25 in the following experiments while using the clinical ultrasound images, respectively, which yield an improved SNR for data with varying noise levels, although these two values are not strictly optimal for each image.
The searching window size is another parameter that affects the despeckling effect. The effect of searching window size is illustrated in Fig. 2, where the experimental results demonstrate that larger search windows do not necessarily produce favorable results. When the searching window size increases, the searching time increases dramatically, whereas the SNR score has no significant difference. The result implies that a 40 × 40 searching window size is suitable for despeckling and can reduce the computation time.
Results
Despeckling of synthetic noisy images
The SNR and MSSIM of the BNLTV, SARBM3D, and proposed methods on the Shepp–Logan phantom image corrupted with Gammadistributed noise are shown in Tables 3 and 4. With regard to SNR, the proposed method significantly outperforms the other two methods in images under all noise levels (standard deviation (sigma) from 0.8 to 1.3). When the noise level increases, the proposed algorithm achieves slightly higher SNRS than the SARBM3D algorithm, and the SSIMS yielded by the proposed algorithm is obviously higher than that of the SARBM3D algorithm. The proposed algorithm outperforms SARBM3D algorithm with SNR improvements ranging from 2.02 dB to 4.51 dB and SSIM improvements from 0.0036 to 0.0063. Compared with the BNLTV algorithm, SNR improvements of the proposed method are ranging from 0.54 dB to 4.93 dB and SSIM improvements from 0.0041 to 0.005.
Figure 3 provides a visual evaluation of the despeckling results from the Shepp–Logan phantom image corrupted with Gammadistributed noise corresponding to sigma of 1.3. The proposed algorithm shows better performance in removes speckle without considerable detail loss and edge blurring and thus exhibits better performance than the other despeckling algorithm. The BNLTV result has the edge preservation but with the texture information loss. The proposed method generates fewer intensity oscillations in the homogeneous region than the SARBM3D method.
The Shepp–Logan phantom image corrupted by the Rayleighdistributed speckle with sigma of 1.3 is used to evaluate the ability of the despeckling filters to manage speckles with different distributions. The results obtained from the SARBM3D, BNLTV and proposed method on the Rayleigh corrupted image are shown in Fig. 4. Table 5 presents the SNR and MSSIM scores in despeckling the Rayleigh corrupted image. Among the methods, the proposed method has the highest SNR and MSSIM scores for the Rayleighbased image.
Despeckling of field II kidney phantom noisy image
The performance of the proposed algorithm on the image is verified by performing Field II simulation of speckle noise. The Bmode image of a synthetic kidney is shown in Fig. 5a. The results despeckled by the three methods are given in the Fig. 5b–d. All the algorithms exhibit improved speckle noise reduction performance. Notably, the result of the BNLTV filter has the most texture information loss and artificial traces compared with the other two despeckling methods. Compared with the SARBM3D filter, the proposed method shows stronger speckle noise reduction ability.
As no ground truth image for the simulated kidney image is present, a blind or noreference image quality assessment metric, natural image quality evaluator (NIQE), is introduced. The quantitative results are listed in Table 6. The results show that the proposed method achieves the highest NIQE value, because it can maintain the texture of the image while removing the speckle noise.
Despeckling of clinical ultrasound image
The effectiveness of the proposed approach on the real ultrasound images is verified. The experiments are conducted on real 2D liver ultrasound image, intravascular ultrasound (IVUS) image, and 3D ultrasound image. Compared with the previous synthetic images, these clinical images contain more fruitful structure information, such as the edges, tiny features, and uniform areas.
The despeckled results of these compared methods on the real 2D liver ultrasound, IVUS, and 3D ultrasound images are presented in Figs. 6, 7, and 8, respectively. The proposed method reduces the speckle noise considerably without losing substantial texture information and smooths the details. As shown by the enlarged images in Figs. 9, 10, and 11, the BNLTV filter slightly smooths the image texture information slightly. Notably, the proposed filter generates the most visually pleasant results. These observations are consistent with the results obtained from the synthetic data.
The quantitative perceptual results of NIQE scoring are provided in Table 7. The results show that the proposed method has the highest NIQE values for 2D liver, IVUS and 3D ultrasound images.
Discussion
The excellent despeckling performance of the proposed filter can be attributed to the spatiogram similarity based on symmetric KL divergence. This similarity measurement is more adaptable to the complex speckle noise, because the KL divergence can measure the difference between two probability distributions and the spatiogram can consider the spatial information of the image patches. Therefore, it demonstrates better performance than the similarity computation based on the nonrobust Euclidean distance, because Euclidean distance uses only the intensity information of the pixels in the patch.
Finally, the limitation of the proposed filter is that its filtering parameters (m, p, λ, and β) are currently manually determined by experience. The optimal parameter setting may change with noise levels and image characteristics. Although the experimental results in this study show that the proposed filter can achieve stateoftheart performance with experiential parameters, the automatic determination of filtering parameters with theoretical foundations is warranted in a future study.
Conclusion
In this paper, an adapted NLTVbased speckle filter has been presented to despeckle the ultrasound images. The filter exploits the spatiogram similarity based on the symmetric KL divergence for the similarity calculation between image patches. As a result, the performance of NLTV is improved. The adapted filter is applied on synthetic and real ultrasound images and then compared with current stateoftheart methods. The results demonstrate that the proposed filter outperforms current stateoftheart filters with regard to ultrasound despeckling. We believe that our method can improve the quality of ultrasound images and has benefit effects on visual perception and diagnostic operations.
Abbreviations
 BM3D:

Blocking matching 3D filter
 BNLTV:

Bayesian NLTV speckle filter
 KL:

Kullbackleibler
 KS test:

Kolmogorovsmirnov test
 NIQE:

Natural image quality evaluator
 NLM:

Nonlocal means
 NLTV:

Nonlocal total variation
 SAR:

Synthetic aperture radar
 SNR:

Signalto noise ratio
 SSIM:

Structural similarity index
References
Gao HT, Huang QH, XM X, Li XL. Wireless and sensorless 3D ultrasound imaging. Neurocomputing. 2016;195:159–71.
Achim A, Bezerianos A, Tsakalides P. Novel Bayesian multiscale method for speckle removal in medical ultrasound images. IEEE Trans Med Imaging. 2001;20:772.
Loizou C, Pattichic C. Despeckle filtering algorithms and software for ultrasound Imagimg. Synth Lect Algorithms Softw Eng. 2015;7:1–180.
Zhong T, Tagare HD, Beaty JD. Evaluation of four probability distribution models for speckle in clinical cardiac ultrasound images. IEEE Trans Med Imaging. 2006;25:1483–91.
VegasSánchezFerrero G, et al. A gamma mixture model for IVUS imaging. Ultrasound Imaging. 2012:25–47.
Sudeep PV, et al. Speckle reduction in medical ultrasound images using an unbiased nonlocal means method. Biomed Signal Process Control. 2016;28:1–8.
Lou L. Study on image noise elimination algorithms based on MATLAB. Journal of Xi'an Shiyou University(Natural Science Edition). 2004; 19:76–80.
Balocco S, Gatta C, Ferre JM, Radeva P. Ultrasound Despeckle Methods. Springer US 2012:4971.
Lee JS. Digital image enhancement and noise filtering by use of local statistics. IEEE Trans Pattern Anal Mach Intell. 1980;2:165–8.
Kuan DT, Sawchuk AA, Strand TC, Chavel P. Adaptive noise smoothing filter for images with signaldependent noise. IEEE Trans Pattern Anal Mach Intell. 1985;7:165–77.
Frost VS, Stiles JA, Shanmugan KS, Holtzman JCA. Model for radar images and its application to adaptive digital filtering of multiplicative noise. IEEE Trans Pattern Anal Mach Intell. 1982;4:157–66.
Loupas T, Mcdicken WN, Allan PL. An adaptive weighted median filter for speckle suppression in medical ultrasonic images. IEEE Trans Circuits Syst. 1989;36:129–35.
Balocco S, Gatta C, Pujol O, Mauri J, Radeva PSRBF. Speckle reducing bilateral filtering. Ultrasound Med Biol. 2010;36:1353–63.
Aysal TC, Barner KE. Rayleighmaximumlikelihood filtering for speckle reduction of ultrasound images. IEEE Trans Med Imaging. 2007;26:712.
Baraldi A, Parmiggiani F. A refined gamma MAP SAR speckle filter with improved geometrical adaptivity. IEEE Trans Geosci Remote Sens. 1995;33(5):1245–57.
Perona P, Malik J. Scalespace and edge detection using anisotropic diffusion. IEEE Trans Pattern Anal Mach Intell. 1990;12:629–39.
YJ Y, Acton ST. Speckle reducing anisotropic diffusion. IEEE Trans Image Process. 2002;11:1260–70.
Krissian K, Westin CF, Kikinis R, Vosburgh GV. Oriented speckle reducing anisotropic diffusion. IEEE Trans Image Process. 2007;16:1412–24.
Shao, Dangguo, et al. Ultrasound speckle reduction based on fractional order differentiation. J Med Ultrasonics. 2017; 44 3: 227237.
Donoho DL. DENOISING BY SOFTTHRESHOLDING. IEEE Trans Inf Theory. 1995;41:613–27.
Zhong S, Cherkassky V. Image denoising using wavelet thresholding and model selection. Int Conf Image Process. 2000;3:262–5.
Gupta S, Chaouhan RC, Sexana SC. Waveletbased statistical approach for speckle reduction in medical ultrasound images. Med Biol Eng Comput. 2004;42:189–92.
Buades A, Coll B, Morel JMA. Review of image Denoising algorithms, with a new one. Multiscale Model Simul. 2005;4:490–530.
Gokul J, Nair MS, Rajan J. Guided SAR image despeckling with probabilistic non local weights. Comput Geosci. 2017; https://doi.org/10.1016/j.cageo.2017.07.004.
Coupe P, Hellier P, Kervrann C, Barillot C. Nonlocal meansbased speckle filtering for ultrasound images. IEEE Trans Image Process. 2009;18:2221–9.
Guo Y, Wang Y, Hou T. Speckle filtering of ultrasonic images using a modified non localbased algorithm. Biomed Signal Process Control. 2011;6:129–38.
Gilboa G, Osher S. Nonlocal operators with applications to image processing. Multiscale Model Simul. 2008;7:1005–28.
Dong FF, Zhang HL, Kong DX. Nonlocal total variation models for multiplicative noise removal using split Bregman iteration. Math Comput Model. 2012;55:939–54.
Wen TX, Gu J, Li L, Qin WJ, Wang L, Xie YQ. Nonlocal Totalvariationbased speckle filtering for ultrasound images. Ultrason Imaging. 2016;38:254–75.
Dabov K, Foi A, Katkovnik V, Egiazarian K. Image denoising by sparse 3D transformdomain collaborative filtering. IEEE Trans Image Process. 2007;16:2080–95.
Parrilli S, Poderico M, Angelino CV, Verdoliva LA, Nonlocal SAR. Image Denoising algorithm based on LLMMSE wavelet shrinkage. IEEE Trans Geosci Remote Sens. 2012;50:606–16.
Zhang J, Wang C, Cheng Y. Comparison of Despeckle filters for breast ultrasound images. Circuits Syst Signal Process. 2014;34:185–208.
Mittal A, Soundararajan R, Bovik AC. Making a completely blind image quality analyzer. IEEE Signal Process Lett. 2013;20:209–12.
Zhang XQ, Burger M, Bresson X, Osher S. Bregmanized nonlocal regularization for deconvolution and sparse reconstruction. SIAM J Imaging Sciences. 2010;3:253–76.
Goldstein T, Osher S. The split Bregman method for L1regularized problems. SIAM J Imaging Sciences. 2009;2:323–43.
Yao ZJ, Liu J, Zhou Y, Liu WY. Similarity measure method using symmetric KL divergence. J Huazhong Univ Sci Tech. 2011;39:1–4.
Eltoft E. Modeling the amplitude statistics of ultrasonic images. IEEE Trans Med Imaging. 2006;25:229–40.
Dudley RM. Uniform Central Limit Theorems. Cambridge University Press 1999; 57:509534.
Manjon JV, Coupe P, Buades A, Louis CD, Robles M. New methods for MRI denoising based on sparseness and selfsimilarity. Med Image Anal. 2012;16:18–27.
Wang Z, Bovik AC, Sheikh HR, Simoncelli EP. Image quality assessment from error visibility to structural similarity. IEEE Trans Image Process. 2004;13:600–12.
Acknowledgements
The author(s) would like to thank the reviewers for their fruitful comments.
Funding
This study was supported by the National Natural Science Foundation of China (Grants 661,771,233,61,271,155, 61,401,451, 61,372,007, 61,571,193) and Guangdong Provincial Key Laboratory of Medical Image Processing (no. 2014B030301042).
Availability of data and materials
The data used to draw the conclusion has been described in “Experiments” section. No further material will be provided.
Author information
Authors and Affiliations
Contributions
SL carried out the study, analyzed the results and drafted the article. FY revised the article and offered useful discussion on validation experiments. TW revised the article and offered the experimental results of the BNLTV filter for comparison. ZY revised the article and offered useful discussion on the proposed method. QH revised the article, offered useful suggestion on the algorithm architecture and the design of experiment. CY revised the article and sorted out the results of experiment. All authors read and approved the final version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
These studies are subject to all necessary approvals from local ethics research committee.
Consent for publication
Informed consent for publication was obtained from all the participants.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Liang, S., Yang, F., Wen, T. et al. Nonlocal total variation based on symmetric KullbackLeibler divergence for the ultrasound image despeckling. BMC Med Imaging 17, 57 (2017). https://doi.org/10.1186/s1288001702317
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1288001702317