Next Article in Journal
Stabilization of All Bell States in a Lossy Coupled-Cavity Array
Next Article in Special Issue
The Solutions to the Uncertainty Problem of Urban Fractal Dimension Calculation
Previous Article in Journal
Research on Fault Extraction Method of Variational Mode Decomposition Based on Immunized Fruit Fly Optimization Algorithm
Previous Article in Special Issue
Explicit Lump Solitary Wave of Certain Interesting (3+1)-Dimensional Waves in Physics via Some Recent Traveling Wave Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Complexity Estimation Based Filtering Method in Wavelet Domain for Magnetic Resonance Imaging Denoising

1
Escuela Superior de Ingeniería Mecánica y Eléctrica, Instituto Politécnico Nacional Av. IPN s/n, Edificio Z, acceso 3, 3er piso; SEPI-Electrónica, Col. Lindavista, 07738 Ciudad de México, Mexico
2
Instituto Nacional de Neurología y Neurocirugía, Av. Insurgentes Sur 3877, Col. La Farma, 14269 Ciudad de México, Mexico
*
Author to whom correspondence should be addressed.
Submission received: 1 March 2019 / Revised: 22 March 2019 / Accepted: 4 April 2019 / Published: 16 April 2019
(This article belongs to the Collection Wavelets, Fractals and Information Theory)

Abstract

:
In this paper, we propose the local complexity estimation based filtering method in wavelet domain for MRI (magnetic resonance imaging) denoising. A threshold selection methodology is proposed in which the edge and detail preservation properties for each pixel are determined by the local complexity of the input image. In the proposed filtering method, the current wavelet kernel is compared with a threshold to identify the signal- or noise-dominant pixels in a scale providing a good visual quality avoiding blurred and over smoothened processed images. We present a comparative performance analysis with different wavelets to find the optimal wavelet for MRI denoising. Numerical experiments and visual results in simulated MR images degraded with Rician noise demonstrate that the proposed algorithm consistently outperforms other denoising methods by balancing the tradeoff between noise suppression and fine detail preservation. The proposed algorithm can enhance the contrast between regions allowing the delineation of the regions of interest between different textures or tissues in the processed images. The proposed approach produces a satisfactory result in the case of real MRI denoising by balancing the detail preservation and noise removal, by enhancing the contrast between the regions of the image. Additionally, the proposed algorithm is compared with other approaches in the case of Additive White Gaussian Noise (AWGN) using standard images to demonstrate that the proposed approach does not need to be adapted specifically to Rician or AWGN noise; it is an advantage of the proposed approach in comparison with other methods. Finally, the proposed scheme is simple, efficient and feasible for MRI denoising.

1. Introduction

Magnetic resonance imaging (MRI) is a powerful medical imaging modality used to produce detailed images of soft tissues and anatomical body structures that can be visualized non-invasively at the millimeter scale [1,2]. MRI processing provides detailed quantitative brain analysis for accurate disease diagnosis [3,4] (i.e., brain tumor diagnosis [5], Alzheimer’s disease (AD), Parkinson’s disease, multiple sclerosis [6], dementia, schizophrenia, brain disorder identification and whole brain analysis of traumatic injury), detection, treatment planning and classification of abnormalities (i.e., extracting tissues like white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF)) [3].
In clinical evaluation and neuroscience research, MRI images are often corrupted by several artifact sources, such as intensity inhomogeneity, abnormal tissues with heterogeneous signal intensities, non-ideal hardware characteristics and the poor choice of scanning parameters [2,3,7]. In order to improve the quality of noisy MRI images to facilitate clinical diagnosis, the MRI pre-processing operations are introduced to improve the qualities of other MRI applications such as segmentation [8], detection [9] and classification [2,10].
Image denoising is a standard pre-processing task for MRI to precisely delineate regions of interest between different brain tissues, to enhance the contrast between regions and to reduce noise, while preserving, as much as possible, the image features as well as structural details [2,10].
Many denoising methods for MRI have been proposed in the literature, these methods can be divided into three major classes [11,12]: (1) filtering techniques include linear filters (i.e., spatial and temporal methods) and non-linear filters (i.e., anisotropic diffusion filtering (ADF) -based methods) [10], 4th order partial differential equation (PDE) –based methods, non-local means (NLM) –based methods [13] and combination of domain and range filters (i.e., bilateral and trilateral filters); (2) transform domain methods, this class consider the curvelet and the contourlet transforms [14,15] and the wavelet transform based methods (i.e., wavelet thresholding, wavelet domain filter, wavelet packet analysis, adaptive multiscale product thresholding, multiwavelet and undecimated wavelet) [7,12,16]; (3) Statistical methods such as maximum likelihood estimation approach [17], Bayesian approach [18], linear minimum mean square error estimation approach, phase error estimation approach, nonparametric neighborhood statistics/estimation approach and singularity function analysis [11,18,19]. Additionally, there exist some hybrid methodologies that belong to both NLM-based methods and Statistical approaches [20,21].
The spectrum of applications in medicine and biology of the wavelet transform has been extremely large, it includes the analysis of the electrocardiogram (ECG) and imaging modalities such as positron emission tomography (PET) and MRI [22]. The main difficulty in dealing with biomedical objects is the variability of the signals and the necessity to operate on a case by case basis [22]. On the other hand, the wavelet decomposition is determined by one mother wavelet function and its dilation and shift versions [23]. There are a lot of wavelet families published in the literature, but researchers commonly have difficulty selecting an optimal wavelet for a specific image processing application [23]. The choice of the optimal wavelet function depends on different criteria in several applications and in some of the distinctive properties (i.e., region of support and the number of vanishing moments) of the wavelet function [23,24].
In this paper, we propose the local complexity estimation based filtering method in wavelet domain for MRI denoising. A threshold selection methodology is proposed in which the edge and detail preservation properties for each pixel are determined by the local complexity of the input image. Statistics of standard deviation select the pixels whose values can be changed since low-energy wavelet coefficients correspond to the smooth regions and high-energy wavelet coefficients are in agreement with the signal features of sharp variation (i.e., edges and textures). In the proposed filtering method, the current wavelet kernel is compared with a threshold to identify the signal- or noise-dominant pixels in a scale providing a good visual quality avoiding blurred and over smoothened processed images. We present a comparative performance analysis with different wavelets to find the optimal wavelet for MRI denoising. The purpose of this research is to eliminate the noise in the MR image as much as possible without losing the details corresponding to image features as the structural details, which will be of highly useful in the quantitative brain analysis for accurate disease diagnosis. Numerical experiments and visual results in simulated MR images degraded with different percentages of Rician noise demonstrate that the proposed algorithm consistently outperforms other denoising methods by balancing the tradeoff between noise suppression and fine detail preservation. The proposed algorithm can enhance the contrast between regions allowing to delineate the regions of interest between different textures or tissues in the processed images. The proposed approach shows a satisfactory result in the case of real MRI denoising by balancing the detail preservation and noise removal, with enhancing the contrast between the regions of the image; otherwise, the comparative methods produce smooth results or limited denoising effectiveness. Additionally, the proposed algorithm is compared with other approaches using standard images degraded with different standard deviations of Additive White Gaussian Noise (AWGN) to demonstrate that the proposed approach does not need to be adapted specifically to Rician or AWGN noise, it is an advantage of the proposed approach in the denoising task of both AWGN and Rician noises against other methods. Finally, the proposed scheme is simple, efficient and feasible for the MRI denoising, the obtained results suggest that the application of the proposed method can benefit many quantitative techniques (i.e., segmentation, tractography or relaxometry) that can take advantage from the denoising and enhanced data produced for the application of the proposed method.
The paper is organized as follows. Section 2 designs the proposed filtering algorithm to MRI denoising. Section 3 presents the performance results in image filtering. Finally, Section 4 concludes the paper.

2. Proposed Method

Discrete wavelet transform (DWT) is an implementation of the wavelet transform using a discrete set of wavelet scales and translations [7]. DWT decomposes an image in different (approximation and detail) sub-bands at different frequencies (scales) with the help of high pass and low pass filters [25]. Figure 1a presents the DWT scheme using high pass filters to extract the high frequency information (i.e., edges and fine details of the image) and low pass filters to obtain the low frequency information (i.e., the low pass representation or the approximation of the image), these filters are first applied in one dimension and then in another one [25]; and Figure 1b depicts the decomposition of a noisy image using the DWT in four wavelet sub-bands labeled as the low-low (LLs) sub-band correspond to the approximation sub-band and the low-high (LHs), high-low (HLs) and high-high (HHs) sub-bands correspond to horizontal, vertical and diagonal details of the image, respectively, s = { 1 , 2 , , S } is the scale and S represents the coarsest scale [12,25].
In the DWT implementation, a standard decimated filterbank algorithm is used (see Figure 1a) [22], a high pass filter g[n] and a low pass filter h[n] are applied to a noisy signal y[n] in the following way [26]
y h i g h [ k ] = n y [ n ]   g [ 2 k n ] y l o w [ k ] = n y [ n ]   h [ 2 k n ]
where y h i g h [ k ] and y l o w [ k ] are the outputs of the high pass and low pass filters, respectively.
In the wavelet thresholding methods, the detail coefficients are processed with soft or hard thresholding to estimate the signal components [7]. The DWT denoising procedure depends upon the usage of wavelet function and thresholding [12]. The wavelet functions are used for estimating the noiseless coefficients from noisy wavelet coefficients in wavelet domain. Various threshold selection methodologies have been proposed to minimize the contribution of noise such as VisuShrink, SureShrink, BayesShrink and NeighShrink [7,12].
We propose a threshold selection methodology in which the edge and detail preservation properties for each pixel are determined by the local complexity of the input image. In the proposed method the current wavelet kernel is compared with a threshold to identify the signal- or noise-dominant pixels in a scale providing a good visual quality avoiding blurred and over smoothened processed images. The steps of the proposed algorithm are given as follows.
Step 1. Apply the DWT. Let obtain the decomposition of the noisy image using the DWT and choose the sub-band HH1 to realize the next steps.
Step 2. Compute the standard deviation. The standard deviation of wavelet coefficients shows the corresponding energy of wavelet coefficients (i.e., low-energy wavelet coefficients appertain to the smooth regions and high-energy wavelet coefficients appertain to the edges and textures). Let compute the standard deviation σ p in the sub-band HH1 where p is the current kernel. The standard deviation σ p is computed using a 3 × 3 kernel according with Figure 2
σ p = m p = 1 n p ( y m p y ¯ ) 2 / n p
where y m p is the mp-th element of the current kernel p = { 1 , 2 , , N } , N is the total number of kernels in the sub-band HH1, y ¯ = m p = 1 n p y m p / n p is the mean value of the current kernel and n p = 9 is the number of elements in the kernel.
Step 3. Compute the threshold. The pixels are classified using a threshold based on the local values of the standard deviations of all kernels in the sub-band HH1. The median value of the standard deviations has been chosen for this purpose. The median is used as robust estimation of the energy of the wavelet kernel coefficients given by its local standard deviation [27,28]. The threshold T selects the pixels whose values are considered as noisy
T = MED { σ 1 , σ 2 , σ 3 , , σ N }
where MED is the median.
Step 4. Apply condition to the current kernel. The proposed condition provides good noise removal, while the edges and the fine details are preserved. The proposed condition to provide denoising is given as follows,
w c o   = { 0 , α p > T w c , o t h e r w i s e
where w c o is the output of proposed procedure, w c is the original wavelet coefficient, α p = | σ p c MED { σ p 1 , σ p 2 , , σ p n } | is the proposed noise estimation parameter, it can be used as an impulsive noise detector when the impulsive noise levels are high [28], it verifies the difference between the value of the median of coefficients and the central coefficient in terms of standard deviation values; σ p c is the standard deviation located in the center of current kernel p , this is, each kernel provides such estimation; and { σ p 1 , σ p 2 , , σ p n } are the standard deviations contained in the current kernel p .
We note that the high-energy wavelet coefficients in the sub-band HH1 involve noise and the edges and textures. The proposed condition (4) distinguishes when a wavelet coefficient (pixel) is noisy or is a detail (edge or texture) in the following way: If the value of the proposed noise estimation parameter α p is bigger than the threshold T , then the current kernel p in the sub-band HH1 is classified as noisy and in such positions the values of the wavelet coefficients are setting in zero (see Figure 2). Otherwise, the wavelet coefficients of this kernel are classified as details and these are unaltered.
Step 5. Compute the Inverse Discrete Wavelet Transform (IDWT). We obtain the restored image applying the IDWT according to the wavelet decomposition.

3. Simulation Results

The proposed local complexity estimation based filtering method in wavelet domain is compared with some reference approaches commonly used in the literature in terms of objective performances given by PSNR (Peak Signal-to-Noise Ratio) and SSIM (Structural Similarity Index) [29] and subjective visual denoising results. The methods used to compare our approach were computed and used in accord with their references. Also, the parameters required by each comparative algorithm are set equal to the values assumed in such references. The reason for choosing these methods to compare them with the proposed one is that their performances have been compared with various known methods and their advantages have been demonstrated.
Four tests have been proposed to determine the performance of the proposed approach. First, a comparative performance analysis in a MRI database is done using the DWT with different wavelets to find the optimal wavelet for MRI denoising; Second, the proposed algorithm is compared with other approaches using standard images degraded with different standard deviations of Additive White Gaussian Noise (AWGN); Third, comparative results in simulated MR images degraded with Rician noise are obtained to evaluate our proposal; and Fourth, a real case of MRI denoising is shown to demonstrate the capabilities of noise filtering of the proposed approach against other methods.
We note that the use of AWGN with different standard deviations is proposed to demonstrate the robustness of the proposed approach in the denoising of standard images in comparison with other methods published recently. In the case of simulated MRI, the Rician noise is built from white Gaussian noise in the complex domain. The proposed approach does not need to be adapted specifically to Rician noise, it is an advantage of the proposed approach in the denoising task of AWGN and Rician noises against other methods. For this reason, we implement these test to determine the performance of the proposed approach.

3.1. Comparative Performance of Different Wavelets

In order to analyze different wavelets for MRI denoising, we utilize a database provided by the National Institute of Neurology and Neurosurgery of Mexico [30]. The real dataset has been recorded by using a Philips Achieva MRI 1.5T scanner with the following parameters: Echo and Repetition Times equal to 102 and 5000 ms, respectively, the Field Of View is 276 × 270 mm and image size of 512 × 512 pixels. This dataset has 900 MRI of three patients (300 MRI for each patient) in a DICOM (Digital Imaging and Communications in Medicine) format. We evaluate the wavelets Haar, Daubechis 2 (DB2), Daubechis 4 (DB4), Symlets 2 (SYM2), Symlets 4 (SYM4), Coiflets 1 (COIF1) y Coiflets 2 (COIF2) in the DWT to realize the MRI denoising. During the wavelet decomposition process, the detail coefficients can be processed with soft or hard thresholding to estimate the signal components for effective denoising [7,12]. Our aim is to find the optimal wavelet according to the best PSNR and SSIM values for a hard threshold of T = 0 , and with this, all high frequency information (the noise, edges and fine details of the image) of the horizontal (LH1), vertical (HL1), and/or diagonal (HH1) details is eliminated. After numerous simulations, we decide to apply this procedure in the wavelet coefficients of HH1 sub-band, the value of T = 0 was chosen only to find the optimal wavelet. With this wavelet, in Section 3.2 we apply the proposed threshold selection methodology to preserve the edges and fine details of the image.
Table 1 presents the average performance results in MRI denoising on the MRI database in terms of PSNR (Peak Signal-to-Noise Ratio) and SSIM (Structural Similarity Index). From Table 1, one can see that the best average PSNR and SSIM performances are given for DB4 and Haar wavelets, respectively. The differences between the results obtained using the two objective quality measures are given because the PSNR is sensitive to the energy of errors instead of real information loss in spite of it is still employed “universal” regardless of its questionable performance in several image applications and SSIM is designed to model any image distortion as a combination of the loss of correlation, luminance distortion and contrast distortion factors, it is applicable to different image processing applications because it does not depend on the images being tested, the viewing conditions or the individual observers [29]. For these reasons, there are no coincidences between both quality measures.
Figure 3 depicts the visual results applying the DWT with different wavelets in a MRI image in terms of PSNR and SSIM. The visual results reveal that the best performances of noise suppression and image distortion are given for DB4 and Haar wavelets, respectively. These results are in concordance with the PSNR and SSIM performances of Table 1.

3.2. Comparative Performance in Standard Images

To evaluate the proposed algorithm in the task of AWGN denoising, we apply the test presented in Reference [12] considered the same data and conditions. For this purpose, we use ten standard images (Lena, Jetplane, Mandrill, House, Boat, Lake, Peppers, Barbara, Pirate and Texture) of size 512 × 512 pixels degraded with the standard deviation σ = { 10 , 20 , 30 , 40 , 50 } of AWGN with zero mean. These images present natural noise, artifacts (noise, intensity, color inhomogeneity in the regions, regions with similar textures, shadows, object reflections, etc.) and diverse content such as fine structures (parallel edges), homogenous areas, texture details and structural information [12]. Comparative performance analysis is carried out for a) wavelet-based approaches such as, VisuShrink with hard threshold [31], BayesShrink [32] and NeighSureShrink [33] and b) NLM (nonlocal means) -based approaches, such as, the standard NLM [34], NLM-DCT (NLM-Discrete Cosine Transform) [35] and NLM-DCT-WEIGHTED (NLM-DCT-Weighted) [12].
Table 2 shows PSNR and SSIM performances for the proposed method with the use of different wavelets in the standard images degraded with a different standard deviation of AWGN. From Table 2, we observe that the best PSNR performance is for the proposed method with DB4 wavelet and in the case of SSIM performance is in favor of the DB4 wavelet in the most of cases ( σ = { 10 , 20 , 30 } ) followed by DB2 and Haar wavelets for σ = { 40 , 50 } .
We note that Table 1 and Table 2 show PSNR and SSIM performances but the first one presents the performances on MRI images using the DWT with a hard threshold of T = 0 to find the best wavelet and the second one shows the performance results for the proposed denoising method on standard images. From these results, we can conclude that with the best (Haar and DB4) wavelets can denoise images of different kind degraded with Rician noise (MRI images) and AWGN (standard images). Section 3.3 will confirm the findings of Table 1 but using the proposed denoising method instead of the hard threshold.
Figure 4 presents the PSNR and SSIM performance analysis for the proposed method and other ones used as comparative in the ten standard images degraded with σ = { 10 , 20 , 30 , 40 , 50 } of AWGN. We show experimental results in the images Lena and House, we observe in Figure 4a,b that in the case of image Lena the best PSNR performance is for the proposed method with DB4 wavelet and for the SSIM performance the proposed method outperforms other methods in the case of σ = { 20 , 30 , 40 } of AWGN; and Figure 4c,d shows that the best PSNR and SSIM performances are for the proposed method with DB4 wavelet in the image House for all standard deviations of AWGN. Then, we provide the average PSNR and SSIM performances for each standard deviation of AWGN using the ten standard images, these results are given in Figure 4e,f where the proposed method with DB4 wavelet provides the best results in terms of PSNR and SSIM performances for each noise level followed by the proposal with Haar wavelet. Figure 4g,h presents the average, minimum and maximum PSNR and SSIM values computed for each denoising method using the ten images. Finally, the results reveal that the proposed method outperforms other denoising methods used as comparative, in the case of average PSNR is in favor of proposed method from 1.06 to 2.4 dB in comparison with the best comparison method (NLM-DCT-Weighted) for the five levels of AWGN and the average SSIM changes from 0.007 to 0.042 in favor of proposed method in comparison with NLM-DCT-Weighted.
Figure 5 depicts the visual results obtained with different denoising algorithms in the images Lena, Mandrill, Lake, Pirate and Texture, degraded with noise level σ = 20 according to Figure 4. The denoised images obtained with the proposed method (DB4 wavelet) have better visual qualities in terms of denoising and fine detail preservation in comparison with other algorithms used as comparative. Moreover, the proposed methodology has the best capability for preserving edges and fine structural details and it enhances the contrast between regions of different texture. It is due to the localization property of wavelets and the proposed condition used to classify the pixels as noisy or details.

3.3. Comparative Performance in Simulated MRI

In this subsection, we realize two tests using simulated MR images from the BrainWeb database [36] and we compare our proposal with different state-of-art denoising methods using different percentages levels of Rician noise.
Test 1: We implement the test realized in Reference [37] with the same data and under the same conditions. In this case, we compare the proposed method with the standard NLM [34], UNLM (Unbiased NLM) [38] and UNLMDCT (UNLM Discrete Cosine Transform) [37] denoising algorithms using three images of 217 × 181 pixels simulated from the BrainWeb database [36] and degraded with 3%, 6%, 9%, 12%, 15% and 18% of Rician noise: (a) T1-weighted MR image, (b) T2-weighted MR image and (c) proton density-weighted (PD-weighted) MR image. Figure 6 shows the PSNR performance of various denoising methods in the simulated MR images, these results reveal that the proposed method with DB4 wavelet provides better PSNR performance for all percentages levels of Rician noise in comparison with other methods used as comparative, this is, the PSNR changes in favor of the proposed method from 0.16 to 1.97 dB in comparison with the best comparison method (UNLMDCT) for the six levels of Rician noise in the three tested images. Figure 7 depicts the visual results in the case of the T1-weighted and PD-weighted MR images degraded with 6% of Rician noise. This Figure shows that in the case of the T1-weighted MR image, the denoised image with the proposed method provides better noise removal and fine detail preservation and allowing the enhancement between different regions corresponding to different tissues in comparison with other algorithms used as comparative. For the PD-weighted MR image, the visual results reveal that the best performance is provided by the proposed method.
Test 2: This test is realized according to Reference [39] with the same data and conditions. For this purpose, the proposed method is compared with the ADF (Anisotropic Diffusion Filter) [40], WIENER Filter [41], TV (Total Variation minimization) [42], standard NLM [34] and NLNS (Nonlocal Neutrosophic Set) [39] denoising algorithms using three images of 217 × 181 pixels simulated from the BrainWeb database [36]: (a) T1-weighted MR image degraded by 7% of Rician noise, (b) T2-weighted MR image degraded by 9% of Rician noise and (c) T1-weighted MR image with multiple sclerosis (MS) lesion degraded by 15% of Rician noise. Figure 8 presents the PSNR and SSIM performance for the three MR images, the PSNR results indicate that the best performance is a favor of the proposed method, this is, the PSNR changes in favor of proposed method from 1.89 to 2.39 dB in comparison with the best comparison method (NLNS) but the SSIM performance of proposal disappoint in comparison with the NLNS from 0.0435 to 0.1192. The SSIM behavior differs from the PSNR because the PSNR is an objective criterion measurement, whereas the SSIM better captures human perception. Figure 9 depicts the visual results for the three MR images according to the results presented in Figure 8. From Figure 9, one can see that the use of the proposed methodology appears to have better visual qualities in comparison with other algorithms used as comparative.

3.4. Comparative Performance in Real MRI

Here, a real case of MRI denoising is presented using the dataset provided in Reference [18]. In this work, Baselice et al. reported comparative results in the real MR image shown in Figure 10a. Denoising visual image results are depicted in Figure 10b–f for the proposed method (DB4 wavelet) and the LMMSE (Linear Minimum Mean Squared Error) [43], BM3D (Block-Matching and 3D) [44], MAP (Maximum A Posteriori estimator) [18] and ADF (Anisotropic Diffusion Filter) [45] denoising algorithms, respectively. From this Figure, the denoising image provided by the proposed approach shows a satisfactory result by balancing the detail preservation and noise removal, by enhancing the contrast between the regions of the image. Otherwise, comparative methods produce smooth results or limited denoising effectiveness. Finally, the obtained results of the proposed approach suggest that it can use as pre-processing stage in MRI applications such as segmentation, detection, and/or classification that can take advantage from the denoising and enhanced data produced for the application of the proposed method.

4. Conclusions

We propose the local complexity estimation based filtering method in wavelet domain for MRI denoising. In the proposed methodology, the edge and detail preservation properties for each pixel are determined by the local complexity of the input image to identify the signal- or noise-dominant pixels in a scale providing a good visual quality avoiding blurred and over smoothened processed images. Numerical experiments and visual results in simulated MR images degraded with different percentages of Rician noise have demonstrated that the proposed denoising algorithm provides better image denoising while preserving image features as well as structural details in comparison with other denoising methods proposed in the literature in most cases. This is due to the proposed condition used to classify the pixels as either noisy or details. In the case of real MRI denoising, the proposed approach produces a satisfactory result by balancing detail preservation and noise removal with enhancing the contrast between the regions of the image; otherwise, the comparative methods produce smooth results or limited denoising effectiveness. Additionally, performance results in standard images degraded with different standard deviations of AWGN indicate that the proposed approach does not need to be adapted specifically to Rician or AWGN noise; it is an advantage of the proposed approach in the denoising task of both AWGN and Rician noises, compared with other methods. The main advantages of the proposed scheme for the MRI denoising and other kinds of images are: a) it is simple because in each iteration to decide if the current pixel is noisy or is a detail only compute one standard deviation and two median values, for this reason, we assume that the time complexity of the proposed approach is much less than other methods such as the NLM-based methods; b) it is efficient because the objective results in terms of PSNR and SSIM criteria and subjective results produced by the visual denoised images reveal that the proposed method provides better results in comparison with other methods; and c) it is feasible because the obtained results suggest that the application of the proposed method can benefit many quantitative techniques (i.e., segmentation, tractography or relaxometry) that gain an advantage from the denoising and enhanced data produced for the application of the proposed method.

Author Contributions

Methodology, I.Y.O.-F., F.J.G.-F. and A.A.-R.; formal analysis, I.Y.O.-F., F.J.G.-F. and A.A.-R.; data curation, I.Y.O.-F., F.J.G.-F. and A.A.-R.; writing—original draft preparation, I.Y.O.-F., F.J.G.-F. and A.A.-R.; writing—review and editing, I.Y.O.-F., F.J.G.-F. and A.A.-R.

Funding

This work is supported by Instituto Politécnico Nacional de México (IPN) and Consejo Nacional de Ciencia y Tecnología (Conacyt), project 240820.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Makropoulos, A.; Counsell, S.J.; Rueckert, D. A review on automatic fetal and neonatal brain MRI segmentation. NeuroImage 2018, 170, 231–248. [Google Scholar] [CrossRef] [PubMed]
  2. El Hassani, A.; Majda, A. Efficient image denoising method based on mathematical morphology reconstruction and the Non-Local Means filter for the MRI of the head. In Proceedings of the 4th IEEE International Colloquium on Information Science and Technology (CiSt), Tangier, Morocco, 24–26 October 2016; pp. 422–427. [Google Scholar]
  3. Dora, L.; Agrawal, S.; Panda, R.; Abraham, A. State-of-the-Art Methods for Brain Tissue Segmentation: A Review. IEEE Rev. Biomed. Eng. 2017, 10, 235–249. [Google Scholar] [CrossRef]
  4. Gordillo, N.; Montseny, E.; Sobrevilla, P. State of the art survey on MRI brain tumor segmentation. Magn. Reson. Imaging 2013, 31, 1426–1438. [Google Scholar] [CrossRef]
  5. El-Dahshan, E.A.; Mohsen, H.M.; Revett, K.; Salem, A.M. Computer-aided diagnosis of human brain tumor through MRI: A survey and a new algorithm. Expert Syst. Appl. 2014, 41, 5526–5545. [Google Scholar] [CrossRef]
  6. Malka, D.; Vegerhof, A.; Cohen, E.; Rayhshtat, M.; Libenson, A.; Shalev, M.A.; Zalevsky, Z. Improved Diagnostic Process of Multiple Sclerosis Using Automated Detection and Selection Process in Magnetic Resonance Imaging. Appl. Sci. 2017, 7, 831. [Google Scholar] [CrossRef]
  7. Anand, C.S.; Sahambi, J.S. Wavelet domain non-linear filtering for MRI denoising. Magn. Reson. Imaging 2010, 28, 842–861. [Google Scholar] [CrossRef]
  8. Xu, R.; Wunsch, D.C. Clustering Algorithms in Biomedical Research: A Review. IEEE Rev. Biomed. Eng. 2010, 3, 120–154. [Google Scholar] [CrossRef]
  9. Shabairou, N.; Cohen, E.; Wagner, O.; Malka, D.; Zalevsky, Z. Color image identification and reconstruction using artificial neural networks on multimode fiber images: Towards an all-optical design. Opt. Lett. 2018, 43, 5603–5606. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, J.; Li, M.; Wang, J.; Wu, F.; Liu, T.; Pan, Y. A survey of MRI-based brain tumor segmentation methods. Tsinghua Sci. Technol. 2014, 19, 578–595. [Google Scholar]
  11. Mohan, J.; Krishnaveni, V.; Guo, Y. A survey on the magnetic resonance image denoising methods. Biomed. Signal Process. Control 2014, 9, 56–69. [Google Scholar] [CrossRef]
  12. Singh, K.; Ranade, S.K.; Singh, C. Comparative performance analysis of various wavelet and nonlocal means based approaches for image denoising. Optik 2017, 131, 423–437. [Google Scholar] [CrossRef]
  13. Foi, A.; Boracchi, G. Foveated Nonlocal Self-Similarity. Int. J. Comput. Vis. 2016, 120, 78–110. [Google Scholar] [CrossRef]
  14. Hyder, S.A.; Sukanesh, R. An efficient algorithm for denoising MR and CT images using digital curvelet transform. Adv. Exp. Med. Biol. 2011, 696, 471–480. [Google Scholar]
  15. Anila, S.; Sivaraju, S.S.; Devarajan, N. A new contourlet based multiresolution approximation for MRI image noise removal. Natl. Acad. Sci. Lett. 2017, 40, 39–41. [Google Scholar] [CrossRef]
  16. Nowak, R.D. Wavelet-Based Rician Noise Removal for Magnetic Resonance Imaging. IEEE Trans. Image Process. 1999, 8, 1408–1419. [Google Scholar] [CrossRef]
  17. He, L.; Greenshields, I.R. A Nonlocal Maximum Likelihood Estimation Method for Rician Noise Reduction in MR Images. IEEE Trans. Med. Imaging 2009, 28, 165–172. [Google Scholar] [PubMed]
  18. Baselice, F.; Ferraioli, G.; Pascazio, V.; Sorriso, A. Bayesian MRI Denoising in complex domain. Magn. Reson. Imaging 2017, 38, 112–122. [Google Scholar] [CrossRef]
  19. Jiang, D.; Dou, W.; Vosters, L.; Xu, X.; Sun, Y.; Tan, T. Denoising of 3D magnetic resonance images with multichannel residual learning of convolutional neural network. Jpn. J. Radiol. 2018, 39, 566–574. [Google Scholar] [CrossRef]
  20. Rajan, J.; den Dekker, A.J.; Sijbers, J. A new non-local maximum likelihood estimation method for Rician noise reduction in magnetic resonance images using the Kolmogorov–Smirnov. Signal Process. 2014, 103, 16–23. [Google Scholar] [CrossRef]
  21. Baselice, F.; Ferraioli, G.; Pascazio, V.; Sorriso, A. Denoising of MR images using Kolmogorov-Smirnov distance in a Non Local Framework. Magn. Reson. Imaging 2019, 57, 176–193. [Google Scholar] [CrossRef] [PubMed]
  22. Unser, M.; Aldroubi, A. A review of wavelets in biomedical applications. Proc. IEEE 1996, 84, 626–638. [Google Scholar] [CrossRef]
  23. Lo, S.C.; Li, H.; Freedman, M.T. Optimization of Wavelet Decomposition for Image Compression and Feature Preservation. IEEE Trans. Med. Imaging 2003, 22, 1141–1151. [Google Scholar] [CrossRef]
  24. Maheswaran, R.; Khosa, R. Comparative study of different wavelets for hydrologic forecasting. Comput. Geosci. 2012, 46, 284–295. [Google Scholar] [CrossRef]
  25. Arora, M.; Bashani, S.; Gupta, K.K.; Mohammed, A.M. Wavelet denoising: Comparative analysis and optimization using machine learning. In Proceedings of the 2014 9th International Conference on Industrial and Information Systems (ICIIS), Gwalior, India, 15–17 December 2014; pp. 1–6. [Google Scholar]
  26. Mitra, S.K. Digital Signal Processing: A Computer Based Approach, 2nd ed.; McGraw-Hill International: Columbus, OH, USA, 2003. [Google Scholar]
  27. Carvajal-Gamez, B.E.; Gallegos-Funes, F.J.; Rosales-Silva, A.J. Color local complexity estimation based steganographic (CLCES) method. Expert Syst. Appl. 2013, 40, 1132–1142. [Google Scholar] [CrossRef]
  28. Gallegos-Funes, F.J.; Rosales-Silva, A.J.; Toledo-Lopez, A. Multichannel image processing by using the Rank M-type L-filter. J. Vis. Commun. Image Represent. 2012, 23, 323–330. [Google Scholar] [CrossRef]
  29. Wang, Z.; Bovik, A.C. A Universal Image Quality Index. IEEE Signal Process. Lett. 2002, 9, 81–84. [Google Scholar] [CrossRef]
  30. Arellano-Reynoso, A. Instituto Nacional de Neurología y Neurocirugía. Available online: http://www.innn.salud.gob.mx/ (accessed on 25 February 2019).
  31. Donoho, D.L.; Johnstone, I.M. Ideal spatial adaptation via wavelet shrinkage. Biometrika 1994, 81, 425–455. [Google Scholar] [CrossRef]
  32. Chang, S.G.; Yu, B.; Vetterli, M. Adaptive wavelet thresholding for image denoising and compression. IEEE Trans. Image Process. 2000, 9, 1532–1546. [Google Scholar] [CrossRef]
  33. Zhou, D.; Cheng, W. Image denoising with an optimal threshold and neighboring window. Pattern Recognit. Lett. 2008, 29, 1694–1697. [Google Scholar]
  34. Buades, A.; Coll, B.; Morel, J. A nonlocal algorithm for image denoising. In Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’05), San Diego, CA, USA, 20–25 June 2005; pp. 60–65. [Google Scholar]
  35. Hu, J.; Pu, Y.; Zhang, Y.; Liu, Y.; Zhou, J. A novel nonlocal means denoising method using the DCT. In Proceedings of the International Conference on Image Processing, Computer Vision, and Pattern Recognition IPCV’11, Las Vegas, NV, USA, 18–21 July 2011; pp. 1–5. [Google Scholar]
  36. BrianWeb. Available online: http://www.bic.mni.mcgill.ca/brainweb/ (accessed on 25 February 2019).
  37. Hu, J.; Pu, Y.; Wu, X.; Zhang, Y.; Zhou, J. Improved DCT-based non local means filter for MR images denoising. Comput. Math. Methods Med. 2012, 2012, 1–14. [Google Scholar] [CrossRef]
  38. Manjón, J.V.; Carbonell-Caballero, J.; Lull, J.J.; García-Martí, G.; Martí-Bonmatí, L.; Robles, M. MRI denoising using nonlocal means. Med Image Anal. 2008, 12, 514–523. [Google Scholar] [CrossRef]
  39. Mohan, J.; Krishnaveni, V.; Guo, Y. MRI denoising using nonlocal neutrosophic set approach of Wiener filtering. Biomed. Signal Process. Control 2013, 8, 779–791. [Google Scholar] [CrossRef]
  40. Perona, P.; Malik, J. Scale-space and edge detection using anisotropic diffusion. IEEE Trans. Pattern Anal. Mach. Intell. 1990, 12, 629–639. [Google Scholar] [CrossRef]
  41. Gonzalez, R.C.; Woods, R.E. Digital Image Processing, 4th ed.; Pearson: Upper Saddle River, NJ, USA, 2017. [Google Scholar]
  42. Rudin, L.I.; Osher, S.; Fatemi, E. Nonlinear total variation based noise removal algorithms. Phys. D 1992, 60, 259–268. [Google Scholar] [CrossRef]
  43. Aja-Fernandez, S.; Niethammer, M.; Kubicki, M.; Shenton, M.; Westin, C.-F. Restoration of DWI data using a Rician LMMSE estimator. IEEE Trans. Med. Imaging 2008, 27, 1389–1403. [Google Scholar] [CrossRef]
  44. Dabov, K.; Foi, A.; Katkovnik, V.; Egiazarian, K. Image denoising by sparse 3-D transform-domain collaborative filtering. IEEE Trans. Image Process. 2007, 16, 2080–2095. [Google Scholar] [CrossRef]
  45. Gerig, G.; Kubler, O.; Kikinis, R.; Jolesz, F. Nonlinear anisotropic filtering of MRI data. IEEE Trans. Med Imaging 1992, 11, 221–232. [Google Scholar] [CrossRef]
Figure 1. Discrete wavelet transform (DWT): (a) DWT scheme using high pass and low pass filters and (b) Decomposition of a noisy image using the DWT in four wavelet sub-bands.
Figure 1. Discrete wavelet transform (DWT): (a) DWT scheme using high pass and low pass filters and (b) Decomposition of a noisy image using the DWT in four wavelet sub-bands.
Entropy 21 00401 g001
Figure 2. Proposed scheme to compute the standard deviation σ p in each kernel p = { 1 , 2 , , N } of the wavelet coefficients from the noisy color image.
Figure 2. Proposed scheme to compute the standard deviation σ p in each kernel p = { 1 , 2 , , N } of the wavelet coefficients from the noisy color image.
Entropy 21 00401 g002
Figure 3. Visual results on magnetic resonance image (MRI) image applying the discrete wavelet transform (DWT) with different wavelets: (a) Original MRI, (b) Haar (PSNR = 34.812, SSIM = 0.941), (c) DB2 (PSNR = 36.571, SSIM = 0.922), (d) DB4 (PSNR = 38.053, SSIM = 0.913), (e) SYM2 (PSNR = 34.598, SSIM = 0.862), (f) SYM4 (PSNR = 33.826, SSIM = 0.853), (g) COIF1 (PSNR = 34.766, SSIM = 0.852), (h) COIF2 (PSNR = 37.054, SSIM = 0.848). The best results are given in bold format.
Figure 3. Visual results on magnetic resonance image (MRI) image applying the discrete wavelet transform (DWT) with different wavelets: (a) Original MRI, (b) Haar (PSNR = 34.812, SSIM = 0.941), (c) DB2 (PSNR = 36.571, SSIM = 0.922), (d) DB4 (PSNR = 38.053, SSIM = 0.913), (e) SYM2 (PSNR = 34.598, SSIM = 0.862), (f) SYM4 (PSNR = 33.826, SSIM = 0.853), (g) COIF1 (PSNR = 34.766, SSIM = 0.852), (h) COIF2 (PSNR = 37.054, SSIM = 0.848). The best results are given in bold format.
Entropy 21 00401 g003
Figure 4. PSNR and SSIM performance analysis of various denoising methods in ten standard images with σ = { 10 , 20 , 30 , 40 , 50 } of AWGN: (a) PSNR performance in the image Lena, (b) SSIM performance in the image Lena, (c) PSNR performance in the image House, (d) SSIM performance in the image House, (e) Average PSNR performance using ten images, (f) Average SSIM performance using ten images, (g) Average, minimum and maximum PSNR values for each denoising method using ten images, (h) Average, minimum and maximum SSIM values for each denoising method using ten images.
Figure 4. PSNR and SSIM performance analysis of various denoising methods in ten standard images with σ = { 10 , 20 , 30 , 40 , 50 } of AWGN: (a) PSNR performance in the image Lena, (b) SSIM performance in the image Lena, (c) PSNR performance in the image House, (d) SSIM performance in the image House, (e) Average PSNR performance using ten images, (f) Average SSIM performance using ten images, (g) Average, minimum and maximum PSNR values for each denoising method using ten images, (h) Average, minimum and maximum SSIM values for each denoising method using ten images.
Entropy 21 00401 g004aEntropy 21 00401 g004b
Figure 5. Visual results for different denoising methods in the images Lena, Mandrill, Lake, Pirate and Texture: (a) Noisy images with noise level σ = 20, (b) Denoised images obtained with NeighSureShrink, (c) Denoised images obtained with NLM-DCT, (d) Denoised images obtained with NLM-DCT-WEIGHTED and (e) Denoised images obtained with proposed method (DB4 wavelet).
Figure 5. Visual results for different denoising methods in the images Lena, Mandrill, Lake, Pirate and Texture: (a) Noisy images with noise level σ = 20, (b) Denoised images obtained with NeighSureShrink, (c) Denoised images obtained with NLM-DCT, (d) Denoised images obtained with NLM-DCT-WEIGHTED and (e) Denoised images obtained with proposed method (DB4 wavelet).
Entropy 21 00401 g005
Figure 6. PSNR performance analysis of various denoising methods in simulated MR images degraded with 3%, 6%, 9%, 12%, 15% and 18% of Rician noise: (a) Original T1-weighted MR image, (b) PSNR performance in the T1-weighted MR image, (c) Original T2-weighted MR image, (d) PSNR performance in the T2-weighted MR image, (e) Original PD-weighted MR image, (f) PSNR performance in the PD-weighted MR image.
Figure 6. PSNR performance analysis of various denoising methods in simulated MR images degraded with 3%, 6%, 9%, 12%, 15% and 18% of Rician noise: (a) Original T1-weighted MR image, (b) PSNR performance in the T1-weighted MR image, (c) Original T2-weighted MR image, (d) PSNR performance in the T2-weighted MR image, (e) Original PD-weighted MR image, (f) PSNR performance in the PD-weighted MR image.
Entropy 21 00401 g006
Figure 7. Visual results for different denoising methods in the simulated MR images: (a) Original T1-weighted and PD-weighted MR images, (b) Noisy MR images degraded with 6% of Rician noise, (c) Denoised MR images obtained with NLM, (d) Denoised MR images obtained with UNLM, (e) Denoised MR images obtained with UNLMDCT and (f) Denoised MR images obtained with proposed method (DB4 wavelet).
Figure 7. Visual results for different denoising methods in the simulated MR images: (a) Original T1-weighted and PD-weighted MR images, (b) Noisy MR images degraded with 6% of Rician noise, (c) Denoised MR images obtained with NLM, (d) Denoised MR images obtained with UNLM, (e) Denoised MR images obtained with UNLMDCT and (f) Denoised MR images obtained with proposed method (DB4 wavelet).
Entropy 21 00401 g007
Figure 8. Performance analysis of various denoising methods in simulated MR images degraded with Rician noise: (a) Original T1-weighted MR image, (b) PSNR and SSIM performances in the T1-weighted MR image degraded by 7% of Rician noise, (c) Original T2-weighted MR image, (d) PSNR and SSIM performances in the T2-weighted MR image degraded by 9% of Rician noise, (e) Original T1-weighted MR image with MS lesion, (f) PSNR and SSIM performances in the T1-weighted MR image with MS lesion degraded by 15% of Rician noise.
Figure 8. Performance analysis of various denoising methods in simulated MR images degraded with Rician noise: (a) Original T1-weighted MR image, (b) PSNR and SSIM performances in the T1-weighted MR image degraded by 7% of Rician noise, (c) Original T2-weighted MR image, (d) PSNR and SSIM performances in the T2-weighted MR image degraded by 9% of Rician noise, (e) Original T1-weighted MR image with MS lesion, (f) PSNR and SSIM performances in the T1-weighted MR image with MS lesion degraded by 15% of Rician noise.
Entropy 21 00401 g008aEntropy 21 00401 g008b
Figure 9. Visual results for different denoising methods in the simulated MR images degraded with different Rician noise: (a) Original MR images, (b) Noisy MR images, (c) Denoised MR images obtained with ADF, (d) Denoised MR images obtained with WIENER, (e) Denoised MR images obtained with TV, (f) Denoised MR images obtained with NLM, (g) Denoised MR images obtained with NLNS and h) Denoised MR images obtained with proposed method (DB4 wavelet).
Figure 9. Visual results for different denoising methods in the simulated MR images degraded with different Rician noise: (a) Original MR images, (b) Noisy MR images, (c) Denoised MR images obtained with ADF, (d) Denoised MR images obtained with WIENER, (e) Denoised MR images obtained with TV, (f) Denoised MR images obtained with NLM, (g) Denoised MR images obtained with NLNS and h) Denoised MR images obtained with proposed method (DB4 wavelet).
Entropy 21 00401 g009
Figure 10. Visual results for different denoising methods in a real MR image: (a) Original MR images, (b) Denoised MR image obtained with proposed method (DB4 wavelet), (c) Denoised MR image obtained with LMMSE, (d) Denoised MR image obtained with BM3D, (e) Denoised MR image obtained with MAP, (f) Denoised MR image obtained with ADF.
Figure 10. Visual results for different denoising methods in a real MR image: (a) Original MR images, (b) Denoised MR image obtained with proposed method (DB4 wavelet), (c) Denoised MR image obtained with LMMSE, (d) Denoised MR image obtained with BM3D, (e) Denoised MR image obtained with MAP, (f) Denoised MR image obtained with ADF.
Entropy 21 00401 g010
Table 1. Average peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) performances obtained from different wavelets. The best results are given in bold format.
Table 1. Average peak signal-to-noise ratio (PSNR) and structural similarity index (SSIM) performances obtained from different wavelets. The best results are given in bold format.
WaveletPSNRSSIM
Haar31.600.941
DB233.500.928
DB435.800.936
SYM234.500.889
SYM432.650.891
COIF131.550.878
COIF234.100.870
Table 2. PSNR and SSIM performances for the proposed denoising method in different standard images with σ = { 10 , 20 , 30 , 40 , 50 } of Additive White Gaussian Noise (AWGN).
Table 2. PSNR and SSIM performances for the proposed denoising method in different standard images with σ = { 10 , 20 , 30 , 40 , 50 } of Additive White Gaussian Noise (AWGN).
ImageNoise(σ)Proposed Denoising Method
HAARDB2DB4SYM2SYM4COIF2COIF4
PSNRSSIMPSNRSSIMPSNRSSIMPSNRSSIMPSNRSSIMPSNRSSIMPSNRSSIM
Lena1034.160.88134.220.88334.980.88334.140.87833.920.87434.150.88034.080.878
2032.910.85432.460.85433.020.85732.190.85032.080.84832.380.85132.360.852
3029.900.81630.000.81930.070.82229.910.81229.860.81229.890.81429.880.813
4028.710.77928.700.77928.780.78828.660.77228.390.76428.450.76828.420.766
5027.360.70327.400.69927.430.70127.340.68927.290.68627.310.68627.330.684
Jetplane1035.050.91935.140.92235.170.92235.010.90634.890.90234.920.91034.990.904
2031.860.87131.870.87331.920.87431.820.86531.800.86231.840.86431.840.863
3029.920.81929.920.82129.960.82729.900.81629.870.81529.910.81729.900.815
4028.240.78828.250.78228.250.78428.200.76220.180.76128.180.76328.180.763
5027.090.74527.080.74427.110.74427.030.72626.980.72427.000.72427.010.726
Mandrill1034.390.92234.440.92234.460.92434.360.91334.350.90534.370.91734.360.909
2031.670.83731.710.84431.750.84731.630.83231.610.82631.630.83431.620.829
3029.560.77929.570.78629.610.79129.540.76329.530.76329.540.77129.560.773
4027.890.70527.930.70127.940.70427.860.69227.850.68827.880.69127.870.687
5027.080.66127.110.66127.150.66027.050.65027.010.63927.030.65427.030.638
House1039.210.93939.230.94439.230.94539.160.93639.150.92139.170.93739.180.924
2036.240.91136.300.91436.330.91936.220.90436.190.90136.220.90336.190.902
3034.660.86734.670.86334.740.86634.550.86934.520.86534.590.86734.580.866
4033.130.85133.120.84233.170.84733.030.84032.980.82732.990.83233.030.828
5031.670.83631.700.83431.720.83631.640.81931.640.81431.650.81831.670.818
Boat1034.000.89134.170.88334.190.88933.940.86933.910.86433.930.87133.940.866
2031.830.82131.850.82531.880.82931.790.81231.710.80431.770.81631.770.806
3029.670.76629.690.76929.690.76729.490.75829.480.75529.530.75929.580.759
4028.000.72328.040.72128.070.72327.920.70927.900.70027.970.70727.960.704
5026.770.65126.770.64726.790.64826.700.63826.690.62926.720.63626.760.633
Lake1033.070.88833.070.89733.080.89732.960.87632.95087032.980.87532.990.867
2030.290.82330.330.82930.370.82630.260.81930.260.80930.260.81730.280.811
3028.780.78028.790.78328.810.78428.730.76628.720.76628.760.77428.750.769
4027.340.73927.390.73727.440.73927.300.72427.280.72327.310.72027.330.721
5026.760.70126.840.69626.880.69926.770.68926.720.68926.740.69126.740.692
Peppers1034.790.87934.840.88834.850.88234.770.86834.760.85634.770.86334.780.863
2033.450.82133.460.82633.490.82733.440.81733.390.81133.420.80833.420.809
3031.370.78131.360.78131.400.78331.330.77931.280.76931.340.77231.330.770
4030.090.74630.100.74230.100.74629.910.73329.870.73229.880.73426.870.732
5028.870.71628.900.71328.940.71528.850.69828.820.68128.830.69728.810.694
Barbara1033.390.88633.390.88633.410.88833.340.87733.300.87533.350.87633.310.879
2029.020.82929.070.82929.170.83129.000.81728.970.81228.980.81329.010.815
3027.000.76127.020.76027.020.76326.890.75026.760.74426.720.74726.740.749
4025.610.72625.610.72425.650.72525.540.71925.520.71225.520.71425.500.714
5024.980.67825.020.67225.020.67424.970.65624.920.65524.950.65824.940.661
Pirate1034.100.88734.120.89034.170.89034.060.87634.000.87634.050.87934.060.875
2031.940.82732.010.83332.010.83531.860.82031.830.81431.880.80831.890.813
3029.460.77929.510.78129.520.78029.440.75529.410.75229.440.75629.420.752
4028.640.72228.690.72028.720.72028.610.71128.590.71028.610.71228.600.708
5026.780.68626.790.67726.790.68226.730.65926.710.66226.730.66326.720.664
Texture1032.350.93632.480.94332.490.94432.360.92832.350.92432.330.92832.340.928
2028.070.88528.110.88928.140.89328.040.87728.040.86928.060.87528.040.873
3026.280.82326.330.82226.360.82426.270.80726.250.80026.270.80326.260.807
4024.110.73724.150.73124.150.73324.050.72224.020.72224.040.71924.050.721
5022.880.68122.970.67423.030.67622.880.66322.830.66022.870.66422.860.661
Average1034.490.90234.510.90534.600.90734.410.89234.350.88634.400.89334.380.889
2031.720.84731.710.85131.810.85331.620.84131.580.83531.640.83831.640.837
3029.670.79729.680.79829.710.80026.660.78629.560.78429.590.78829.600.783
4028.280.75129.090.74929.120.75028.100.73827.250.73328.080.73627.780.734
5027.020.70527.050.70127.080.70326.990.68826.960.68326.980.68926.960.687
Note: The best results are given in bold format.

Share and Cite

MDPI and ACS Style

Orea-Flores, I.Y.; Gallegos-Funes, F.J.; Arellano-Reynoso, A. Local Complexity Estimation Based Filtering Method in Wavelet Domain for Magnetic Resonance Imaging Denoising. Entropy 2019, 21, 401. https://0-doi-org.brum.beds.ac.uk/10.3390/e21040401

AMA Style

Orea-Flores IY, Gallegos-Funes FJ, Arellano-Reynoso A. Local Complexity Estimation Based Filtering Method in Wavelet Domain for Magnetic Resonance Imaging Denoising. Entropy. 2019; 21(4):401. https://0-doi-org.brum.beds.ac.uk/10.3390/e21040401

Chicago/Turabian Style

Orea-Flores, Izlian Y., Francisco J. Gallegos-Funes, and Alfonso Arellano-Reynoso. 2019. "Local Complexity Estimation Based Filtering Method in Wavelet Domain for Magnetic Resonance Imaging Denoising" Entropy 21, no. 4: 401. https://0-doi-org.brum.beds.ac.uk/10.3390/e21040401

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop