**Curvelet-based ****Bayesian ****Estimator ****for ****Speckle ****Suppression ****in ****Ultrasound ****Imaging**

**Abstract. **Ultrasound images are inherently affected by speckle noise, and thus the reduction of this noise is a crucial pre-processing step for their successful interpretation. Bayesian estimation is a powerful signal estimation technique used for speckle noise removal in images. In the Bayesian-based despeckling schemes, the choice of suitable statistical models and the development of a shrinkage function for estimation of the noise-free signal are the major concerns. In this paper, a novel curvelet-based Bayesian estimation scheme for despeckling of ultrasound images is developed. The curvelet coefficients of the multiplicative degradation model of the noisy ultrasound image are additively decomposed into noise-free and signal-dependent noise components. The Cauchy and two-sided exponential distributions are assumed to be statistical models for the noise-free and signal-dependent noise components of the observed curvelet coefficients, respectively, and an efficient low-complexity realization of the Bayesian estimator is proposed. The experimental results demonstrate the validity of the proposed despeckling scheme in providing a significant suppression of the speckle noise and simultaneously preserving the image details.

**Keywords:**Ultrasound imaging, curvelet transform, speckle noise, Bayesian estimation, statistical modeling.

- Introduction

Ultrasound imaging is important for medical diagnosis and has the advantages of cost effectiveness, port-ability, acceptability and safety [1]. However, ultrasound images are of relatively poor quality due to its contamination by the speckle noise, which considerably degrades the image quality and leads to a negative impact on the diagnostic task. Thus, reducing speckle noise while preserving anatomic information is necessary to better delineate the regions of interest in the ultrasound images.

In the work of speckle suppression in ultrasound images, many spatial-based techniques that employ either single-scale or multi-scale filtering have been developed in the literature [2-4]. Earlydeveloped single-scale spatial filtering [2] are limited in their capability for significantly reducing the speckle noise. More promising spatial single-scale techniques such as those using bilateral filtering [4] and nonlocal filtering [3] have been recently proposed.

This work was supported in part by the Natural Sciences and Engineering Research Council (NSERC) of Canada and in part by the Regroupement Strategique en Microelectronique du Quebec (ReSMiQ).

These techniques depend on the size of the filter window, and hence, for a satisfactory speckle suppression, they require large computational time. Alternatively, multi-scale spatial techniques [5], based on partial differential equations, have been investigated in the literature. These techniques are iterative and can produce smooth images with preserved edges. However, important structural details are unfortunately degraded during the iteration process. As an appropriate alternative to spatial-based speckle suppression in ultrasound images, many other despeckling techniques based on different transform domains, such as the ones of wavelet, contourlet, and curvelet, have been recently proposed in the literature [6-8]. Wavelet transform has a good reputation as a tool for noise reduction but has the drawback of poor directionality, which makes its usage limited in many applications. Using contourlet transform provides an improved noise reduction performance due to its property of fi‚exible directional decomposability. However, curvelet transform offers a higher directional sensitivity than that of contourlet transform and is more efficient in representing the curve-like details in images.

For the development of despeckling techniques based on transform domains, thresholding [7] has been presented as a technique to build linear estimators of the noise-free signal coefficients. However, the main drawback of this thresholding technique is in the difficulty of determining a suitable threshold value. To circumvent this problem, non-linear estimators [6] have been statistically developed based on Bayesian estimation formalism. For the development of an efficient Bayesian-based despeckling scheme, the choice of a suitable probability distribution to model the transform domain coefficients is a major concern. Also, while investigating a suitable statistical model, the complexity of the Bayesian estimation process should be taken into consideration. Consequently, special attention should be paid to the realization complexity of the Bayesian estimator that results from employing the selected probabilistic model in one of the Bayesian frameworks.

In this paper, to achieve a satisfactory performance for despeckling of ultrasound images at a lower computational effort, a new curvelet-based Bayesian scheme is proposed. The multiplicative degradation model representing an observed ultrasound image is decomposed into an additive model consisting of noise-free and signal-dependent noise components. Two-sided exponential distribution is used as a prior statistical model for the curvelet coefficients of the signal-dependant noise. This model, along with the Cauchy distribution is used to develop a low-complexity Bayesian estimator. The performance of the proposed Bayesian despeckling scheme is evaluated on both syntheticallyspeckled and real ultrasound images, and the results are compared to that of some other existing despeckling schemes.

- Modeling of Curvelet Coefficients

The multiplicative degradation model of a speckle-corrupted ultrasound image *g*(*i**,**j*) in the spatial domain is given by

*g*(*i,**j*) = *v*(*i,**j*)*s*(*i,**j*)(1)

where *v*(*i**,**j*) and *s*(*i,**j*) denote the noise-free image and the speckle noise, respectively. This model of the noisy observation of *v*(*i**,**j*) can be additively decomposed as a noise-free signal component and a signal-dependant noise:

*g*(*i**,**j*) = *v*(*i**,**j*) + (*s*(*i,**j*) *âˆ’*1)*v*(*i**,**j*)

= *v*(*i,**j*) + *u*(*i,**j*)(2)

where (*s*(*i,**j*) *âˆ’*1)*v*(*i**,**j*) represents the signal-dependant noise. Taking the curvelet transform of (2) at level *l*, we have

*y*[*l**,d*](*i**,**j*) = *x*[*l**,d*](*i**,**j*) + *n*[*l**,d*](*i**,**j*)(3)

where *y*[*l**,d*](*i**,**j*), *x*[*l**,d*](*i**,**j*) and *n*[*l**,d*](*i**,**j*) denote, respectively, the (*i,**j*)th curvelet coefficient of the observed image, the corresponding noise free image and the corresponding additive signal-dependant noise at direction *d*= 1*,*2*,*3*,**Â·**Â·**Â·**,**D*. In order to simplify the notation, we will henceforth drop both the superscripts *l*and *d*and the index (*i,**j*).

In this work, in order to reduce the noise inherited in ultrasound images, we propose exploiting the statistical characteristics of the curvelet coefficients in (3) to derive an efficient Bayesian estimator. Thus, one needs to provide a prior probabilistic model for the curvelet coefficients of *x*and

*n*. It has been shown that the distribution of the curvelet coefficients of noise-free images can be suitably modeled by the Cauchy distribution [9]. The zero-mean Cauchy distribution is given by

*p**x*(*x*) = (*Î³/Ï€*)(*x*2 + *Î³*2)(4)

where *Î³*is the dispersion parameter. The noisy observation is used to estimate the Cauchy distribution parameter *Î³*by minimizing the function

2

Â *Ï†*Ë†*y**y**t*

(*t*) *âˆ’**Ï†*(*t*) *e**âˆ’*

*dt*(5)

where *Ï†*Ë†*y*(*t*) is the empirical characteristic function corresponding to the curvelet coefficients *y*of

22

the noisy observation, *Ï†**y*(*t*) = *Ï†**x*(*t*)*Ï†**E*(*t*), *Ï†**x*(*t*) = *e**âˆ’**Î³**|**t**|*, and *Ï†**E*(*t*) = *e**âˆ’*(*Ïƒ**ï¿½**/*2)*|**t**|*

deviation *Ïƒ**E*obtained as

with the standard

*Ïƒ**E*=

MAD(*y*(*i,**j*))

0*.*6745

(6)

In (6), MAD denotes the median absolute deviation operation. Now, in order to formulate the Bayesian estimator, a prior statistical assumption for the curvelet coefficients of *n*of the signal dependant noise should also be assumed. From experimental observation, it is noticed that the tail part of the empirical distribution of *n*decays at a low rate. Hence, in this paper, we propose to use a two-sided exponential (TSE) distribution given by

1

*p**n*(*n*) =*e**âˆ’|**n**|**/Î²*

2*Î²*

(7)

where *Î²*is a positive real constant referred to as the scale parameter. The method of log-cummulants (MoLC) is adopted to estimate the parameter *Î²*, and thus the estimated *Î²*Ëœ is obtained by using the following expression:

*Î²*Ëœ = exp

1*N*1Â Â Â Â Â Â Â Â *N*2

log(*y*(*i,**j*))+ *Î¾*

(8)

*N*1*N*2

*i*=1*j*=1

where *Î¾*is the Euler-Mascheroni constant and *N*1 and *N*2 define the size *N*1 *Ã-**N*2 of the curvelet subband considered.

- Bayesian Estimator

Due to the fact that each of the Cauchy and TSE distributions has only one parameter, one could expect the process of Bayesian estimation to be of lower complexity. The values of the Bayes estimates *x*Ë† of the noise-free curvelet coefficients *x*of a subband under the quadratic loss function, which minimizes the mean square error (MSE), are given by the shrinkage function:

*x*Ë†(*y*) =*p**x**|**y*(*x**|**y*)*xdx*

P *p**n*(*y**âˆ’**x*)*p**x*(*x*)*xdx*

=P *p*(*y**âˆ’**x*)*p*(*x*)

(9)

It is noted that a closed-form expression for

*x*Ë†(*y*) given by the above equation does not exist.

Thus, in order to obtain the Bayesian estimates for the noise-free curvelet coefficients, the two integrations associated with (9) are numerically performed for each curvelet coefficient. Since this procedure requires an excessive computational effort, the bayseian estimates are obtained by replacing the associated integrals in (9) with infinite series as suggested in [10]. Accordingly, the Bayesian shrinkage function can be expressed as

*e**âˆ’**y/Î²*[*f *(*y*)*Î¶*] + *e**y/Î²*[ *f*(*y*) + *Î¶*]

*x*Ë†(*y*) =(10)

*e**âˆ’**y/Î²*[*f*21(*y*) + *Î¶*2] + *e**y/Î²*[*âˆ’**f*22(*y*) + *Î¶*2]

where

*f*11(*y*) = *f*12

(*âˆ’**y*) = sin(*Î³**/Î²*) Im *E*( *âˆ’**y*+ *j**Î³*)*âˆ’*Si(*Î³**/Î²*) + *Ï€*

1*Î²*2

__âˆ’____y____+____j____Î³__

*âˆ’*cos(*Î³/Î²*) Â ReÂ *E*1(*Î²*

+ Ci(*Î³**/Î²*) ,(11)

*f*(*y*) = *âˆ’**f*

1*âˆ’**y*+ *j**Î³*

(*âˆ’**y*) = *âˆ’ *sin(*Î³**/Î²*) Re *E*()+ Ci(*Î³**/Î²*)

2122*Î³*1*Î²*

__1____âˆ’____y____+____j____Î³____Ï€__

*âˆ’**Î³*cos(*Î³/Î²*) Â ImÂ *E*1(*Î²*

*âˆ’*Si(*Î³**/Î²*) + 2

,(12)

*Î¶*1 = lim *f*12

*y**â†’**âˆž*

(*y*) = sin(*Î³**/Î²*) *âˆ’*Si(*Î³**/Î²*) + *Ï€ **âˆ’*cos(*Î³**/Î²*)Ci(*Î³**/Î²*), and(13)

*Î¶*= lim *f*

11

(*y*) =sin(*Î³/Î²*)Ci(*Î³/Î²*) +cos(*Î³/Î²*)

*Ï€*

Si(*Î³**/Î²*) +

(14)

222

*y**â†’**âˆž*

In the equations above, *j*= *âˆš**âˆ’*1, Im*{Â·}*and Re*{Â·}*are the imaginary and real parts, respectively, of a complex argument, and *E*1(*Â·*), Si(*Â·*) and Ci(*Â·*) are, respectively, the exponential, sine and cosine integral functions obtained as in [10].

- Experimental Results

Extensive experimentations are carried out in order to study the performance of the proposed despeckling scheme. The results are compared with those of other existing despeckling schemes that use improved-Lee filtering [2], adaptive-wavelet shrinkage [6], and contourlet thresholding [7]. Performance evaluation of the various despeckling schemes is conducted on synthetically-speckled and real ultrasound images. In the implementation of the proposed speckling scheme, the 5-level decomposition of the curvelet transform is applied. From the experimental observation, applying a higher level of decomposition of the curvelet transform does not lead to any improvement in the despeckling performance. Since the curvelet transform is a shift-variant transform, the cycle spinning [11] is performed on the observed noisy image to avoid any possible pseudo-Gibbs artifacts in the neighborhood of discontinuities. In the proposed despeckling scheme, only the detail curvelet coefficients are despeckled using the Bayesian shrinkage function in (10).

The peak signal-to-noise ratio (PSNR) is used as a quantitative measure to assess the despeckling performance of the various schemes when applied on synthetically-speckled images. Table I gives the PSNR values obtained when applying the various schemes on two synthetically-speckled images of size 512*Ã-*512, namely, *Lena*and *Boat*. It is obviously seen from this table that, in all cases, the proposed despeckling scheme provides higher values of PSNR compared to that provided by the other schemes. To have a better insight on the despeckling performance of the various schemes, the results in Table 1 are visualized in Figure 1. It is obvious from this figure that the superiority of the proposed scheme over the other schemes is more evident when a higher level of speckle noise is introduced to the test images. In order to study the performances of the various despeckling schemes on real ultrasound images, two images obtained from [12] and shown in Figure 2 are used. Since the noise-free images cannot be made available, one can only give a subjective evaluation of the performance of the various despeckling schemes. From Figure 2, it is clearly seen that the schemes in [2] and [6] provide despeckled images that suffer from the presence of visually noticeable speckle noise. On the other hand, the scheme in [7] severely over-smooth the noisy images thus providing despeckled images in which some of the texture details are lost. However, the proposed despeckling scheme results in images with not only a significant reduction in the speckle noise but also a good preservation of the textures of the original images.

Table 1: The PSNR values obtained when applying the various despeckling schemes on *Lena*and

*Boat*images contaminated by speckle noise at different levels.

**34**

**[2]**

**32****[6]**

**30****[7]**

**Proposed**

**28**

**26**

**24**

**22**

**20**

**18**

**0.1****0.2****0.3****0.4****0.5****0.7****1**

**Standard deviation of noise**

(a)

**32**

**[2]**

**30****[6]**

**28****[7]**

**Proposed**

**26**

**24**

**22**

**20**

**18**

**16**

**0.1****0.2****0.3****0.4****0.5****0.7****1**

**Standard deviation of noise**

(b)

Fig. 1: Quantitative comparison between the various despeckling schemes in terms of PSNR values:

(a) *Lena*image; (b) *Boat*image.

- Conclusion

In this paper, a new curvelet-based scheme for suppressing the speckle noise in ultrasound images has been developed in the framework of Bayesian estimation. The observed ultrasound image is first additively decomposed into noise-free and signal-dependant noise components. The Cauchy and twosided exponential distributions have been used as probabilistic models for the curvelet coefficients of the noise-free and signal-dependant noise components, respectively, of the ultrasound image. The proposed probabilistic models of the curvelet coefficients of an observed ultrasound image has been employed to formulate a Bayesian shrinkage function in order to obtain the estimates of the noise-free curvelet coefficients. A low-complexity realization of this shrinkage function has been employed. Experiments have been carried out on both synthetically-speckled and real ultrasound images in order to demonstrate the performance of the proposed despeckling scheme. In comparison with some other existing despeckling schemes, the results have shown that the proposed scheme provides higher PSNR values and gives well-despeckled images with better diagnostic details.

- (b)

(c)(d)(e)(f)

(g)(h)(i)(j)

Fig. 2: Qualitative comparison between the various despeckling schemes. (a)(b) Noisy ultrasound images. Despeckled images obtained by applying the schemes in (c)(g) [2] ,(d)(h) [6] ,(e)(i) [7] and (f)(j) the proposed scheme.

References

- Dhawan, A.P.: Medical image analysis. Volume 31. John Wiley & Sons (2011)
- Loupas, T., McDicken, W., Allan, P.: Â An adaptive weighted median filter for speckle suppression in medical ultrasonic images. IEEE transactions on Circuits and Systems
**36**(1) (1989) 129-135 - CoupÂ´e, P., Hellier, P., Kervrann, C., Barillot, C.: Nonlocal means-based speckle filtering for ultrasound images. IEEE transactions on image processing
**18**(10) (2009) 2221-2229 - Sridhar, B., Reddy, K., Prasad, A.: An unsupervisory qualitative image enhancement using adaptive morphological bilateral filter for medical images. International Journal of Computer Applications
**10**(2i) (2014) 1 - Abd-Elmoniem, K.Z., Youssef, A.B., Kadah, Y.M.: Real-time speckle reduction and coherence enhancement in ultrasound imaging via nonlinear anisotropic diffusion. IEEE Transactions on Biomedical Engineering
**49**(9) (2002) 997-1014 - Swamy, M., Bhuiyan, M., Ahmad, M.: Spatially adaptive thresholding in wavelet domain for despeckling of ultrasound images. IET Image Process
**3**(3) (2009) 147-162 - Hiremath, P., Akkasaligar, P.T., Badiger, S.: Speckle reducing contourlet transform for medical ultrasound images. Int J Compt Inf Engg
**4**(4) (2010) 284-291 - Jian, Z., Yu, Z., Yu, L., Rao, B., Chen, Z., Tromberg, B.J.: Speckle attenuation in optical coherence tomography by curvelet shrinkage. Optics letters
**34**(10) (2009) 1516-1518 - Deng, C., Wang, S., Sun, H., Cao, H.: Multiplicative spread spectrum watermarks detection performance analysis in curvelet domain. In: 2009 International Conference on E-Business and Information System Security. (2009)
- Damseh, R.R., Ahmad, M.O.: A low-complexity mmse bayesian estimator for suppression of speckle in sar images. In: Circuits and Systems (ISCAS), 2016 IEEE International Symposium on, IEEE (2016) 1002-1005
- Temizel, A., Vlachos, T., Visioprime, W.: Wavelet domain image resolution enhancement using cycle-spinning. Electronics Letters
**41**(3) (2005) 119-121