Research - International Journal of Medical Research & Health Sciences ( 2021) Volume 10, Issue 2
Analysis of Filtering Approaches to Brain MRI in Spatial DomainVedant Shukla1*, Prasad Khandekar2 and Arti Khaparde1
2Department of Engineering and Technology, Dr. Vishwanath Karad MIT World Peace University, Pune, India
Vedant Shukla, School of Electronics and Communication Engineering, Dr. Vishwanath Karad MIT World Peace University, Pune, India, Email: email@example.com
Received: 30-Dec-2020 Accepted Date: Feb 11, 2021 ; Published: 18-Feb-2021
Non-linear filters are used to filter out the artifacts and noise present in MR data. The balance between signal preservation and noise reduction makes restoration of MR data a complex task. Application of non-linear filters such as median and Non-Local Means Filter (NLM) filter converts the right-skewed Rician Distribution into un-skewed Gaussian distribution. NLM filter gives better results than Bilateral and Median filter. As the distribution is un-skewed after the application of non-linear filters, standard linear filters such as Gaussian and Wiener filters were applied and results were drawn. A linear combination of NLM and Gaussian filter gives satisfactory results. The experimentation was performed on 40 clinical images and the NLM filter was found to have the best results. The image quality indices used for comparison are Peak Signal-to-Noise Ratio (PSNR), Root Mean Squared Error (RMSE), Structural Similarity Index (SSIM), and entropy. The experimentation was performed on MATLAB 2020a.
Rician distribution, Bilateral filter, NLM filter, PSNR, RMSE, SSIM, Entropy
MR signal is prone to signal-dependent Rician noise. Non-linear filters prove to be the most optimum model to filter such images. Since decades of research, several filters are applied and tested on MR Data. In 1986, Pitas I, et al. applied non-linear filters to multiple data, to understand the properties of noise present in the data . In 2008, Fernendez SA, et al. developed Linear Minimum Mean Squared Error (LMMSE) filter for the restoration of MR brain images . In 2015, Mohan J, et al. published a survey of MR data de-noising . A fourth-order PDE is one non-linear filter that can be applied to MR data as suggested by Yu-Li, et al. .
In this article, multiple filtering approaches were discussed, such as linear and non-linear filtering, transform domain as well as statistical approach. A combination of linear and non-linear filters can also be implemented; one such approach is discussed by Meenakshy, et al. . A bilateral filter is one non-linear filter that smoothens the images and preserves the edges as well. Chaudhary K, et al. from IISc Bangalore modified the standard bilateral filter to make it more accurate and fast . Though trilateral filters are also developed, which an extended Bilateral filters by Chang H, et al. in 2018 . In 2020, Chen KL, et al. proposed a novel filter that was a combination of Fuzzy c-means and Adaptive Non-Local Means .
The synthetic data on which the filtering models were performed are referred from the BrainWeb database . National Library of Medicine, USA provides MR images for research and analysis, which were used in this research article . Clinical Brain MR Images were obtained from BRATS 2012 database .
MRI Data Noise Distribution
MR Data is mainly distorted due to thermal noise from the scanned object. Thermal noise is in turn caused due to electronic components, coil, etc. The raw data obtained is complex representing Fourier Transform, and inverse Fourier Transform needs to be computed. The magnitude image is thus non-linear and the PDF of such images changes. The Probability Distribution Function (PDF) of a Rician distributed image is expressed as:
Here, I0 (∙) is the 0th order Bessel Function of the first kind and u (∙) is the Heaviside step function.
Application of Standard Non-linear Filters
Non-linear filters are applied to Rician distributed data. Some standard non-linear filters are the Bilateral filter, Median filter, and the Non-Local Means filter. Although the NLM filter should show better-filtering results, the distribution of the filtered image should be of importance. In Figure 1, the histogram of the Rician corrupted image and the median filtered image is shown. It can be observed that although the distribution in Rician corrupted image is Rician, the distribution of the median filtered image is Gaussian i.e. unskewed .
A similar change is observed in the distribution of in NLM filtered image as shown in Figure 1. Although the change in distribution occurs when input noise is high, hence it can be stated that the non-linear filters convert Rician intensity distributed data into Gaussian.
Application of Linear Filters to Gaussian Distributed Data
From the intensity plot of non-linear filtered image, it is evident that the distribution of image is Gaussian. For such distribution, linear filters such as Gaussian filter and Wiener filter can be applied. Linear filters are low-pass filters, i.e. they help in image smoothing. In frequency domain filtering, the homomorphic filter is one well-known filter for converting noise from signal-dependent to signal independent [12-17].
The filters were applied on MR data including synthetic T1 weighted brain MRI and pelvis subset MRI as well as clinical brain MRI. The block diagram of the experimentation is demonstrated in Figure 2.
Application of Non-linear Filters to Rician Distributed MR Images
The very general non-linear filters applied to MR data are Bilateral, Median, and Non-Local Means filters. In Table 1, we compare the filters using PSNR, RMSE, SSIM, and entropy.
|%N||s||PSNR 1||PSNR 2||PSNR 3||PSNR 4||SSIM 1||SSIM 2||SSIM 3||SSIM 4||RMSE 1||RMSE 2||RMSE 3||RMSE 4||Entropy 1||Entropy 2||Entropy 3||Entropy 4|
Table 1 :Image quality indices comparison.
In Table 1, PSNR of Rician corrupted and the original image is PSNR 1, bilateral filtered and the original image is PSNR 2, median filtered and the original image is PSNR 3 and that of NLM filtered and the original image is PSNR 4. Similar terminologies are used for RMSE, SSIM, and entropy. It can be noted that the NLM filter gives better results than the other non-linear filters.
In Figure 3, it can be observed that the Rician distribution upon application of non-linear filters tends to unskewed normal distribution. Since the distribution is now normal, we apply linear filters such as Gaussian and Wiener filters.
Linear Filters to Gaussian Distributed Data to Synthetic Brain MR Images
Linear filters such as Gaussian and Wiener filters are applied to non-linear filtered MR images. In Table 2, PSNR 5, PSNR 6, PSNR 7, PSNR 8 represent the Peak Signal to Noise Ratio (PSNR) of Gaussian filter applied to Median filtered image, Gaussian filter applied to NLM filtered image, Wiener filter applied to Median filtered image, Wiener filter applied to NLM filtered image concerning the original image. Other parameters such as SSIM, RMSE, and Entropy follow the same terminology. The corresponding filtered image can be seen in Figure 3e, Figure 3f, Figure 3g, and Figure 3h respectively. In Table 2, the combination of an application of Gaussian filter on NLM filter gives the best results. This can be observed in RMSE, SSIM, PSNR, and Entropy.
|%N||s||PSNR 5||PSNR 6||PSNR 7||PSNR 8||SSIM 5||SSIM 6||SSIM 7||SSIM 8||RMSE 5||RMSE 6||RMSE 7||RMSE 8||Entropy 5||Entropy 6||Entropy 7||Entropy 8|
Table 2: Image Quality Indices Comparison
It should be noted that SSIM 6 in Table 2 shows better results than SSIM 4 in Table 1. Hence, it can be concluded that a linear combination of Linear Gaussian Filter and non-linear NLM filter gives the best results.
Application of Linear Combination of Gaussian and NLM Filter to Clinical MR Images
In the last subsection, it was concluded that a linear combination of non-linear NLM Filter and linear Gaussian filter gives best better results for synthetic MR images. This combination of filters is now applied to clinical images. As a demonstration, all the non-linear filters and combination is applied to the preliminary clinical image (Figure 4).
The study was performed on 40 clinical images. The Structural Similarity Index (SSIM) of randomly selected 10 images is tabulated in Table 3. The terminologies are the same as in Tables 1 and 2. It can be noted that similar observations can be seen on clinical images, as on synthetic images. Experimentation on clinical images is challenging because of the absence of ground truth. The non-linear NLM filter works best as compared to the Bilateral and Median filter, while a linear combination of NLM and Gaussian filter gives satisfactory results. This can be interpreted using SSIM values in Table 3 and visually analysed using the corresponding line plot as in Figure 5.
|#||SSIM 2||SSIM 3||SSIM 4||SSIM 5||SSIM 6||SSIM 7||SSIM 8|
Table 3 :Noise level vs. SSIM of filtered brain images
Due to signal dependency of Rician noise, non-linear filtering is preferred. In this research work, standard non-linear filters such as Bilateral filter, Median filter, and NLM filters are applied to synthetic as well as clinical brain MR data. It is observed that the NLM filter works best when compared to other non-linear filters. It was also observed that upon application of non-linear filters, the distribution of the MR data tends to be Gaussian. Hence further experiments were performed by applying a linear combination of non-linear and linear filters and it was observed that NLM and Gaussian filter in series gives SSIM values closest to unity. But in the case of clinical MR images, the NLM filter gives the best results while the NLM-Gaussian filter gives acceptable results. This was validated by computing the standard image quality parameters.
Conflict of Interest
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
- Pitas, Ioannis, and A. Venetsanopoulos. "Nonlinear mean filters in image processing."IEEE Transactions on Acoustics, Speech, and Signal Processing,Vol. 34, No. 3, 1986, pp. 573-84.
- Aja-Fernández, Santiago, Carlos Alberola-López, and Carl-Fredrik Westin. "Noise and signal estimation in magnitude MRI and Rician distributed images: A LMMSE approach."IEEE Transactions on Image Processing,Vol. 17, No. 8, 2008, pp. 1383-98.
- Mohan, J., V. Krishnaveni, and Yanhui Guo. "A survey on the magnetic resonance image denoising methods."Biomedical Signal Processing and Control,Vol. 9, 2014, pp. 56-69.
- You, Yu-Li, et al. "Behavioral analysis of anisotropic diffusion in image processing."IEEE Transactions on Image Processing,Vol. 5, No. 11, 1996, pp. 1539-53.
- Meenakshy, K., and C. Shwetha. "Combined NLM-Weiner filter based image denoising."2015 International Conference on Power, Instrumentation, Control and Computing (PICC). IEEE, 2015.
- Chaudhury, Kunal N., and Swapnil D. Dabhade. "Fast and provably accurate bilateral filtering."IEEE Transactions on Image Processing,Vol. 25, No. 6, 2016, pp. 2519-28.
- Chang, Herng-Hua, Cheng-Yuan Li, and Audrey Haihong Gallogly. "Brain MR image restoration using an automatic trilateral filter with GPU-based acceleration."IEEE Transactions on Biomedical Engineering,Vol. 65, No. 2, 2017, pp. 400-13.
- Chen, Kaixin, et al. "An enhanced adaptive non-local means algorithm for Rician noise reduction in magnetic resonance brain images."BMC Medical Imaging,Vol. 20, No. 1, 2020, pp. 1-9.
- BrainWeb: Simulated brain database for synthetic MR images. http://www.bic.mni.mcgill.ca/brainweb/
- National Library of Medicine, USA. The Visible Human Data. https://www.nlm.nih.gov/research/visible/visible_human.html
- BRATS. "Multimodal brain tumor segmentation challenge." 2012. http://www.imm.dtu.dk/projects/BRATS2012/
- Aarya, Isshaa, Danchi Jiang, and Timothy Gale. "Signal dependent Rician noise denoising using nonlinear filter."Lecture Notes on Software Engineering,Vol. 1, No. 4, 2013, pp. 344-49.
- Saleh, Sami Abdulla Mohsen, and Haidi Ibrahim. "Mathematical equations for homomorphic filtering in frequency domain: A literature survey."Proceedings of the International Conference on Information and Knowledge Management, Vol. 45, 2012, pp. 74-77.
- Rajan, Jeny, K. Kannan, and M. R. Kaimal. "An improved hybrid model for molecular image denoising."Journal of Mathematical Imaging and Vision,Vol. 31, No. 1, 2008, pp. 73-79.
- Liu, X. Y., C-H. Lai, and Kyriacos A. Pericleous. "A fourth-order partial differential equation denoising model with an adaptive relaxation method."International Journal of Computer Mathematics,Vol. 92, No. 3, 2015, pp. 608-22.
- Lai, Rui, et al. "Local and nonlocal steering kernel weighted total variation model for image denoising."Symmetry,Vol. 11, No. 3, 2019, p. 329.
- Khaparde, Arti, and Vaidehi Deshmukh. "Optimized multi-focus image fusion using genetic algorithm."Advances in Science Technology and Engineering Systems Journal,Vol. 2, No. 1, 2017, pp. 51-56.