Wavelet domain non-linear filtering for MRI denoising

Denoising of 3D magnetic resonance images via edge-enhanced low-rank tensor decomposition

Magnetic Resonance images (MRI) denoising is to obtain high quality image from infectant version. Recently, low-rank tensor (LRT) methods have been developed and attained resounding success in MRI denoising. However, these pure LRT models are incapable of utilizing the comprehensive inherent information of clean MRI. To overcome these drawbacks, we design a novel edge-enhanced low-rank tensor approximation (EELRTA) framework for Rician noise removal. The tensor gradient L0 norm regularization with describing the local structure information is incorporated into the weighted core tensor rank model for improving texture edge preservation. The application of weights can further preserve the potentially useful information distributed on the different core tensor coefficients with different physical meanings. What's more, non-local self-similarity tactic is employed for low-rank sparsity-encourage and enhancing anti-noise capability of EELRTA model. The proposed EELRTA method is tackled by an efficient alternating direction method of multipliers (ADMM). The Experiment results on simulation and multiple sclerosis lesion (MSL) data illustrate that the proposed method can effectively remove noise while reasonably retaining pathological structure information.

Artificial Intelligence for Neuroimaging in Pediatric Cancer

BACKGROUND/OBJECTIVES

Artificial intelligence (AI) is transforming neuroimaging by enhancing diagnostic precision and treatment planning. However, its applications in pediatric cancer neuroimaging remain limited. This review assesses the current state, potential applications, and challenges of AI in pediatric neuroimaging for cancer, emphasizing the unique needs of the pediatric population.

METHODS

A comprehensive literature review was conducted, focusing on AI's impact on pediatric neuroimaging through accelerated image acquisition, reduced radiation, and improved tumor detection. Key methods include convolutional neural networks for tumor segmentation, radiomics for tumor characterization, and several tools for functional imaging. Challenges such as limited pediatric datasets, developmental variability, ethical concerns, and the need for explainable models were analyzed.

RESULTS

AI has shown significant potential to improve imaging quality, reduce scan times, and enhance diagnostic accuracy in pediatric neuroimaging, resulting in improved accuracy in tumor segmentation and outcome prediction for treatment. However, progress is hindered by the scarcity of pediatric datasets, issues with data sharing, and the ethical implications of applying AI in vulnerable populations.

CONCLUSIONS

To overcome current limitations, future research should focus on building robust pediatric datasets, fostering multi-institutional collaborations for data sharing, and developing interpretable AI models that align with clinical practice and ethical standards. These efforts are essential in harnessing the full potential of AI in pediatric neuroimaging and improving outcomes for children with cancer.

Combining variational mode decomposition with regularisation techniques to denoise MRI data

In this paper, we propose a novel method for removing noise from MRI data by exploiting regularisation techniques combined with variational mode decomposition. Variational mode decomposition is a new decomposition technique for sparse decomposition of a 1D or 2D signal into a set of modes. In turn, regularisation is a method that can translate the ill-posed problem (e.g., image denoising) into a well-posed problem. The proposed method aims to remove the noise from the image in two steps. In the first step, the MR imaging data are decomposed by the 2D variational mode decomposition algorithm. In the second step, for effective suppression of Rician noise from MRI data, we used the fused lasso signal approximator with all modes acquired from the medical scan. The performance of the proposed approach was compared with state-of-the-art reference methods based on different metrics, that is, the peak signal-to-noise ratio, the structural similarity index metrics, the high-frequency error norm, the quality index based on local variance, and the sharpness index. The experiments were performed on the basis of both simulated and real images. The presented results prove the high denoising performance of the proposed algorithm; particularly under heavy noise conditions.

Automatic autism spectrum disorder detection using artificial intelligence methods with MRI neuroimaging: A review

Autism spectrum disorder (ASD) is a brain condition characterized by diverse signs and symptoms that appear in early childhood. ASD is also associated with communication deficits and repetitive behavior in affected individuals. Various ASD detection methods have been developed, including neuroimaging modalities and psychological tests. Among these methods, magnetic resonance imaging (MRI) imaging modalities are of paramount importance to physicians. Clinicians rely on MRI modalities to diagnose ASD accurately. The MRI modalities are non-invasive methods that include functional (fMRI) and structural (sMRI) neuroimaging methods. However, diagnosing ASD with fMRI and sMRI for specialists is often laborious and time-consuming; therefore, several computer-aided design systems (CADS) based on artificial intelligence (AI) have been developed to assist specialist physicians. Conventional machine learning (ML) and deep learning (DL) are the most popular schemes of AI used for diagnosing ASD. This study aims to review the automated detection of ASD using AI. We review several CADS that have been developed using ML techniques for the automated diagnosis of ASD using MRI modalities. There has been very limited work on the use of DL techniques to develop automated diagnostic models for ASD. A summary of the studies developed using DL is provided in the Supplementary Appendix. Then, the challenges encountered during the automated diagnosis of ASD using MRI and AI techniques are described in detail. Additionally, a graphical comparison of studies using ML and DL to diagnose ASD automatically is discussed. We suggest future approaches to detecting ASDs using AI techniques and MRI neuroimaging.

How Machine Learning is Powering Neuroimaging to Improve Brain Health

This report presents an overview of how machine learning is rapidly advancing clinical translational imaging in ways that will aid in the early detection, prediction, and treatment of diseases that threaten brain health. Towards this goal, we aresharing the information presented at a symposium, "Neuroimaging Indicators of Brain Structure and Function - Closing the Gap Between Research and Clinical Application", co-hosted by the McCance Center for Brain Health at Mass General Hospital and the MIT HST Neuroimaging Training Program on February 12, 2021. The symposium focused on the potential for machine learning approaches, applied to increasingly large-scale neuroimaging datasets, to transform healthcare delivery and change the trajectory of brain health by addressing brain care earlier in the lifespan. While not exhaustive, this overview uniquely addresses many of the technical challenges from image formation, to analysis and visualization, to synthesis and incorporation into the clinical workflow. Some of the ethical challenges inherent to this work are also explored, as are some of the regulatory requirements for implementation. We seek to educate, motivate, and inspire graduate students, postdoctoral fellows, and early career investigators to contribute to a future where neuroimaging meaningfully contributes to the maintenance of brain health.

A new unsupervised pseudo-siamese network with two filling strategies for image denoising and quality enhancement

Digital image noise may be introduced during acquisition, transmission, or processing and affects readability and image processing effectiveness. The accuracy of established image processing techniques, such as segmentation, recognition, and edge detection, is adversely impacted by noise. There exists an extensive body of work which focuses on circumventing such issues through digital image enhancement and noise reduction, but this work is limited by a number of constraints including the application of non-adaptive parameters, potential loss of edge detail information, and (with supervised approaches) a requirement for clean, labeled, training data. This paper, developed on the principle of Noise2Void, presents a new unsupervised learning approach incorporating a pseudo-siamese network. Our method enables image denoising without the need for clean images or paired noise images, instead requiring only noise images. Two independent branches of the network utilize different filling strategies, namely zero filling and adjacent pixel filling. Then, the network employs a loss function to improve the similarity of the results in the two branches. We also modify the Efficient Channel Attention module to extract more diverse features and improve performance on the basis of global average pooling. Experimental results show that compared with traditional methods, the pseudo-siamese network has a greater improvement on the ADNI dataset in terms of quantitative and qualitative evaluation. Our method therefore has practical utility in cases where clean images are difficult to obtain.

Deep Adaptive Blending Network for 3D Magnetic Resonance Image Denoising

The visual quality of magnetic resonance images (MRIs) is crucial for clinical diagnosis and scientific research. The main source of quality degradation is the noise generated during MRI acquisition. Although denoising MRI by deep learning methods shows great superiority compared with traditional methods, the deep learning methods reported to date in the literature cannot simultaneously leverage long-range and hierarchical information, and cannot adequately utilize the similarity in 3D MRI. In this paper, we address the two issues by proposing a deep adaptive blending network (DABN) characterized by a large receptive field residual dense block and an adaptive blending method. We first propose the large receptive field residual dense block that can capture long-range information and fuse hierarchical features simultaneously. Then we propose the adaptive blending method that produces denoised pixels by adaptively filtering 3D MRI, which explicitly utilizes the similarity in 3D MRI. Residual is also considered as a compensating item after adaptive filtering. The blending adaptive filter and residual are predicted by a network consisting of several large receptive field residual dense blocks. Experimental results show that the proposed DABN outperforms state-of-the-art denoising methods in both clinical and simulated MRI data.

Feasibility of Three-Dimensional Balanced Steady-State Free Precession Cine Magnetic Resonance Imaging Combined with an Image Denoising Technique to Evaluate Cardiac Function in Children with Repaired Tetralogy of Fallot

Objective

To investigate the feasibility of cine three-dimensional (3D) balanced steady-state free precession (b-SSFP) imaging combined with a non-local means (NLM) algorithm for image denoising in evaluating cardiac function in children with repaired tetralogy of Fallot (rTOF).

Materials and Methods

Thirty-five patients with rTOF (mean age, 12 years; range, 7–18 years) were enrolled to undergo cardiac cine image acquisition, including two-dimensional (2D) b-SSFP, 3D b-SSFP, and 3D b-SSFP combined with NLM. End-diastolic volume (EDV), end-systolic volume (ESV), stroke volume (SV), and ejection fraction (EF) of the two ventricles were measured and indexed by body surface index. Acquisition time and image quality were recorded and compared among the three imaging sequences.

Results

3D b-SSFP with denoising vs. 2D b-SSFP had high correlation coefficients for EDV, ESV, SV, and EF of the left (0.959–0.991; p < 0.001) as well as right (0.755–0.965; p < 0.001) ventricular metrics. The image acquisition time ± standard deviation (SD) was 25.1 ± 2.4 seconds for 3D b-SSFP compared with 277.6 ± 0.7 seconds for 2D b-SSFP, indicating a significantly shorter time with the 3D than the 2D sequence (p < 0.001). Image quality score was better with 3D b-SSFP combined with denoising than with 3D b-SSFP (mean ± SD, 3.8 ± 0.6 vs. 3.5 ± 0.6; p = 0.005). Signal-to-noise ratios for blood and myocardium as well as contrast between blood and myocardium were higher for 3D b-SSFP combined with denoising than for 3D b-SSFP (p < 0.05 for all but septal myocardium).

Conclusion

The 3D b-SSFP sequence can significantly reduce acquisition time compared to the 2D b-SSFP sequence for cine imaging in the evaluation of ventricular function in children with rTOF, and its quality can be further improved by combining it with an NLM denoising method.

Joint multi-field T1 quantification for fast field-cycling MRI

Purpose

Recent developments in hardware design enable the use of fast field‐cycling (FFC) techniques in MRI to exploit the different relaxation rates at very low field strength, achieving novel contrast. The method opens new avenues for in vivo characterizations of pathologies but at the expense of longer acquisition times. To mitigate this, we propose a model‐based reconstruction method that fully exploits the high information redundancy offered by FFC methods.

Methods

The proposed model‐based approach uses joint spatial information from all fields by means of a Frobenius ‐ total generalized variation regularization. The algorithm was tested on brain stroke images, both simulated and acquired from FFC patients scans using an FFC spin echo sequences. The results are compared to three non‐linear least squares fits with progressively increasing complexity.

Results

The proposed method shows excellent abilities to remove noise while maintaining sharp image features with large signal‐to‐noise ratio gains at low‐field images, clearly outperforming the reference approach. Especially patient data show huge improvements in visual appearance over all fields.

Conclusion

The proposed reconstruction technique largely improves FFC image quality, further pushing this new technology toward clinical standards.

A Survey on State-of-the-art Denoising Techniques for Brain Magnetic Resonance Images

The accuracy of the magnetic resonance (MR) image diagnosis depends on the quality of the image, which degrades mainly due to noise and artifacts. The noise is introduced because of erroneous imaging environment or distortion in the transmission system. Therefore, denoising methods play an important role in enhancing the image quality. However, a tradeoff between denoising and preserving the structural details is required. Most of the existing surveys are conducted on a specific MR image modality or on limited denoising schemes. In this context, a comprehensive review on different MR image denoising techniques is inevitable. This survey suggests a new direction in categorizing the MR image denoising techniques. The categorization of the different image models used in medical image processing serves as the basis of our classification. This study includes recent improvements on deep learning-based denoising methods alongwith important traditional MR image denoising methods. The major challenges and their scope of improvement are also discussed. Further, many more evaluation indices are considered for a fair comparison. An elaborate discussion on selecting appropriate method and evaluation metric as per the kind of data is presented. This study may encourage researchers for further work in this domain.

FFA-DMRI: A Network Based on Feature Fusion and Attention Mechanism for Brain MRI Denoising

Magnetic Resonance Imaging (MRI) is an indispensable tool in the diagnosis of brain diseases due to painlessness and safety. Nevertheless, Rician noise is inevitably injected during the image acquisition process, which leads to poor observation and interferes with the treatment. Owing to the complexity of Rician noise, using the elimination method of Gaussian to remove it does not perform well. Therefore, the feature fusion and attention network (FFA-DMRI) is proposed to separate noise from observed MRI. Inspired by the attention-guided CNN network (ADNet) and Convolutional block attention module (CBAM), a spatial attention mechanism has been specially designed to obtain the area of interest in MRI. Furthermore, the feature fusion block concatenates local with global information, which makes full use of the multilevel structure and boosts the expressive ability of network. The comprehensive experiments on Alzheimer's disease neuroimaging initiative dataset (ADNI) have demonstrated high effectiveness of FFA-DMRI with maintaining the crucial brain details. Moreover, in terms of visual inspections, the denoising results are also consistent with human perception.

Free-breathing liver fat and R 2 ∗ quantification using motion-corrected averaging based on a nonlocal means algorithm

PURPOSE

To propose a motion-robust chemical shift-encoded (CSE) method with high signal-to-noise (SNR) for accurate quantification of liver proton density fat fraction (PDFF) and R 2 ∗ .

METHODS

A free-breathing multi-repetition 2D CSE acquisition with motion-corrected averaging using nonlocal means (NLM) was proposed. PDFF and R 2 ∗ quantified with 2D CSE-NLM were compared to two alternative 2D techniques: direct averaging and single acquisition (2D 1ave) in a digital phantom. Further, 2D NLM was compared in patients to 3D techniques (standard breath-hold, free-breathing and navigated), and the alternative 2D techniques. A reader study and quantitative analysis (Bland-Altman, correlation analysis, paired Student's t-test) were performed to evaluate the image quality and assess PDFF and R 2 ∗ measurements in regions of interest.

RESULTS

In simulations, 2D NLM resulted in lower standard deviations (STDs) of PDFF (2.7%) and R 2 ∗ (8.2  s - 1 ) compared to direct averaging (PDFF: 3.1%, R 2 ∗ : 13.6  s - 1 ) and 2D 1ave (PDFF: 8.7%, R 2 ∗ : 33.2  s - 1 ). In patients, 2D NLM resulted in fewer motion artifacts than 3D free-breathing and 3D navigated, less signal loss than 2D direct averaging, and higher SNR than 2D 1ave. Quantitatively, the STDs of PDFF and R 2 ∗ of 2D NLM were comparable to those of 2D direct averaging (p>0.05). 2D NLM reduced bias, particularly in R 2 ∗ (-5.73 to -0.36  s - 1 ) that arises in direct averaging (-3.96 to 11.22  s - 1 ) in the presence of motion.

CONCLUSIONS

2D CSE-NLM enables accurate mapping of PDFF and R 2 ∗ in the liver during free-breathing.

MRI denoising using progressively distribution-based neural network

Magnetic Resonance (MR) images often suffer from noise pollution during image acquisition and transmission, which limits the accuracy of quantitative measurements from the data. Noise in magnitude MR images is usually governed by Rician distribution, due to the existence of uncorrelated Gaussian noise with zero-mean and equal variance in both the real and imaginary parts of the complex K-space data. Different from the existing MRI denoising methods that utilizing the spatial neighbor information around the pixels or patches, this work turns to capture the pixel-level distribution information by means of supervised network learning. A progressive network learning strategy is proposed via fitting the distribution of pixel-level and feature-level intensities. The proposed network consists of two residual blocks, one is used for fitting pixel domain without batch normalization layer and another one is applied for matching feature domain with batch normalization layer. Experimental results under synthetic, complex-valued and clinical MR brain images demonstrate great potential of the proposed network with substantially improved quantitative measures and visual inspections.

Denoising and contrast-enhancement approach of magnetic resonance imaging glioblastoma brain tumors

We investigate a new preprocessing approach for MRI glioblastoma brain tumors. Based on combined denoising technique (bilateral filter) and contrast-enhancement technique (automatic contrast stretching based on image statistical information), the proposed approach offers competitive results while preserving the tumor region's edges and original image's brightness. In order to evaluate the proposed approach's performance, quantitative evaluation has been realized through the Multimodal Brain Tumor Segmentation (BraTS 2015) dataset. A comparative study between the proposed method and four state-of-the art preprocessing algorithm attests that the proposed approach could yield a competitive performance for magnetic resonance brain glioblastomas tumor preprocessing. In fact, the result of this step of image preprocessing is very crucial for the efficiency of the remaining brain image processing steps: i.e., segmentation, classification, and reconstruction.

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

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.

MRI denoising by NeighShrink based on chi-square unbiased risk estimation

NeighShrink is an efficient image denoising algorithm for the reduction of additive white Gaussian noise. However, it does not perform well in terms of Rician noise removal for MRI (Magnetic Resonance Imaging). Allowing for the characteristics of squared-magnitude MR (Magnetic Resonance) images, which follow a non-central chi-square distribution, the CURE (Chi-Square Unbiased Risk Estimation) is used to determine an optimal threshold for NeighShrink. Therefore, we propose the NeighShrinkCURE denoising algorithm. Bilateral filtering and cycle spinning are used to further improve denoising performance. Experimental results show that the proposed algorithm is simple and efficient, and provides good noise reduction while preserving edges and details well. Compared with some similar MRI denoising algorithms, the proposed algorithm has improvements in PSNR (Peak Signal-to-Noise Ratio) and SSIM (Structural Similarity). Although running time of the proposed algorithm has some increment compared with some current MRI denoising algorithms, the comprehensive performance of the proposed algorithm is superior to several current MRI denoising algorithms.

Enhancement of the Nuclear Spin Noise Signal Using Wavelet Transform

Spin noise spectroscopy has attracted considerable attention recently owing partly to intrinsic interest in the phenomenon and partly to its significant application potential. Here, we address the inherent problem of low sensitivity of nuclear spin noise and examine the utility of wavelet transform to mitigate this problem by distinguishing real peaks from the noise contaminated data. Suppression of the random circuit noise and the consequent enhancement of the correlated nuclear spin noise signal have been demonstrated with discrete wavelet transform. Spectra of both ¹ H and ¹³ C nuclear spins have been considered and significant signal enhancements in both the cases have been observed. A detailed analysis of several possible wavelet, thresholding and decomposition solutions have been made to obtain the optimum condition for signal enhancement. It is observed that the application of wavelet transform leaves the spin noise signal line shape essentially unchanged, which is an advantage for several applications involving spin noise spectra.

Weighted Schatten p-norm minimization for 3D magnetic resonance images denoising

Magnetic resonance (MR) imaging plays an important role in clinical diagnosis and scientific research. A clean MR image can better provide patient's information to doctors or researchers for further treatment. However, in real life, MR images are inevitably corrupted by annoying Rician noise in the process of imaging. Aiming at the Rician noise of 3D MR images, a framework is proposed to suppress noise by low-rank matrix approximation (LRMA) with weighted Schatten p-norm minimization regularization (WSNMD-3D). The proposed method not only considers the importance of different rank components, but can also approximate the true rank of the latent low-rank matrix. This approach first groups similar non-local cubic patches extracted from the noisy 3D MR image into a matrix whose columns are vectorized patches. The above matrix can be modeled as a low-rank matrix approximate model. Then weighted Schatten p-norm minimization (WSNM) is applied to the model, which shrinks different rank components with different treatments. Finally, the denoised 3D MR image is acquired by aggregating all denoised patches with weighted averaging. Experimental results on synthetic and real 3D MR data show that the proposed method obtains better results than state-of-the-art methods, both visually and quantitatively.

De-noising of 3D multiple-coil MR images using modified LMMSE estimator

De-noising is a crucial topic in Magnetic Resonance Imaging (MRI) which focuses on less loss of Magnetic Resonance (MR) image information and details preservation during the noise suppression. Nowadays multiple-coil MRI system is preferred to single one due to its acceleration in the imaging process. Due to the fact that the model of noise in single-coil and multiple-coil MRI systems are different, the de-noising methods that mostly are adapted to single-coil MRI systems, do not work appropriately with multiple-coil one. The model of noise in single-coil MRI systems is Rician while in multiple-coil one (if no subsampling occurs in k-space or GRAPPA reconstruction process is being done in the coils), it obeys noncentral Chi (nc-χ). In this paper, a new filtering method based on the Linear Minimum Mean Square Error (LMMSE) estimator is proposed for multiple-coil MR Images ruined by nc-χ noise. In the presented method, to have an optimum similarity selection of voxels, the Bayesian Mean Square Error (BMSE) criterion is used and proved for nc-χ noise model and also a nonlocal voxel selection methodology is proposed for nc-χ distribution. The results illustrate robust and accurate performance compared to the related state-of-the-art methods, either on ideal nc-χ images or GRAPPA reconstructed ones.

A Bayesian Multiresolution Approach for Noise Removal in Medical Magnetic Resonance Images

A Bayesian approach using wavelet coefficient modeling is proposed for de-noising additive white Gaussian noise in medical magnetic resonance imaging (MRI). In a parallel acquisition process, the magnetic resonance image is affected by white Gaussian noise, which is additive in nature. A normal inverse Gaussian probability distribution function is taken for modeling the wavelet coefficients. A Bayesian approach is implemented for filtering the noisy wavelet coefficients. The maximum likelihood estimator and median absolute deviation estimator are used to find the signal parameters, signal variances, and noise variances of the distribution. The minimum mean square error estimator is used for estimating the true wavelet coefficients. The proposed method is simulated on MRI. Performance and image quality parameters show that the proposed method has the capability to reduce the noise more effectively than other state-of-the-art methods. The proposed method provides 8.83%, 2.02%, 6.61%, and 30.74% improvement in peak signal-to-noise ratio, structure similarity index, Pratt’s figure of merit, and Bhattacharyya coefficient, respectively, over existing well-accepted methods. The effectiveness of the proposed method is evaluated by using the mean squared difference (MSD) parameter. MSD shows the degree of dissimilarity and is 0.000324 for the proposed method, which is less than that of the other existing methods and proves the effectiveness of the proposed method. Experimental results show that the proposed method is capable of achieving better signal-to-noise ratio performance than other tested de-noising methods.

3D CMRO2 mapping in human brain with direct 17O MRI: Comparison of conventional and proton-constrained reconstructions

Oxygen metabolism is altered in brain tumor regions and is quantified by the cerebral metabolic rate of oxygen consumption (CMRO₂). Direct dynamic ¹⁷O MRI with inhalation of isotopically enriched ¹⁷O₂ gas can be used to quantify CMRO₂; however, pixel-wise CMRO₂ quantification in human brain is challenging due to low natural abundance of ¹⁷O isotope and, thus, the low signal-to-noise ratio (SNR) of ¹⁷O MR images. To test the feasibility CMRO₂ mapping at a clinical 3 T MRI system, a new iterative reconstruction was proposed, which uses the edge information contained in a co-registered ¹H gradient image to construct a non-homogeneous anisotropic diffusion (AD) filter. AD-constrained reconstruction of ¹⁷O MR images was compared to conventional Kaiser-Bessel gridding without and with Hanning filtering, and to iterative reconstruction with a total variation (TV) constraint. For numerical brain phantom and in two in vivo data sets of one healthy volunteer, AD-constrained reconstruction provided ¹⁷O images with improved resolution of fine brain structures and resulted in higher SNR. CMRO₂ values of 0.78 - 1.55µmol/g_tissue/min (white brain matter) and 1.03 - 2.01µmol/g_tissue/min (gray brain matter) as well as the CMRO₂ maps are in a good agreement with the results of ¹⁵O-PET and ¹⁷O MRI at 7 T and at 9.4 T. In conclusion, the proposed AD-constrained reconstruction enabled calculation of 3D CMRO₂ maps at 3 T MRI system, which is an essential step towards clinical translation of ¹⁷O MRI for non-invasive CMRO₂ quantification in tumor patients.

Noise Estimation and Reduction in Magnetic Resonance Imaging Using a New Multispectral Nonlocal Maximum-likelihood Filter

Denoising of magnetic resonance (MR) images enhances diagnostic accuracy, the quality of image manipulations such as registration and segmentation, and parameter estimation. The first objective of this paper is to introduce a new, high-performance, nonlocal filter for noise reduction in MR image sets consisting of progressively-weighted, that is, multispectral, images. This filter is a multispectral extension of the nonlocal maximum likelihood filter (NLML). Performance was evaluated on synthetic and in-vivo T₂ - and T₁ -weighted brain imaging data, and compared to the nonlocal-means (NLM) and its multispectral version, that is, MS-NLM, and the nonlocal maximum likelihood (NLML) filters. Visual inspection of filtered images and quantitative analyses showed that all filters provided substantial reduction of noise. Further, as expected, the use of multispectral information improves filtering quality. In addition, numerical and experimental analyses indicated that the new multispectral NLML filter, MS-NLML, demonstrated markedly less blurring and loss of image detail than seen with the other filters evaluated. In addition, since noise standard deviation (SD) is an important parameter for all of these nonlocal filters, a multispectral extension of the method of maximum likelihood estimation (MLE) of noise amplitude is presented and compared to both local and nonlocal MLE methods. Numerical and experimental analyses indicated the superior performance of this multispectral method for estimation of noise SD.

Exploring the complementarity of THz pulse imaging and DCE-MRIs: Toward a unified multi-channel classification and a deep learning framework

We provide a comprehensive account of recent advances in biomedical image analysis and classification from two complementary imaging modalities: terahertz (THz) pulse imaging and dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). The work aims to highlight underlining commonalities in both data structures so that a common multi-channel data fusion framework can be developed. Signal pre-processing in both datasets is discussed briefly taking into consideration advances in multi-resolution analysis and model based fractional order calculus system identification. Developments in statistical signal processing using principal component and independent component analysis are also considered. These algorithms have been developed independently by the THz-pulse imaging and DCE-MRI communities, and there is scope to place them in a common multi-channel framework to provide better software standardization at the pre-processing de-noising stage. A comprehensive discussion of feature selection strategies is also provided and the importance of preserving textural information is highlighted. Feature extraction and classification methods taking into consideration recent advances in support vector machine (SVM) and extreme learning machine (ELM) classifiers and their complex extensions are presented. An outlook on Clifford algebra classifiers and deep learning techniques suitable to both types of datasets is also provided. The work points toward the direction of developing a new unified multi-channel signal processing framework for biomedical image analysis that will explore synergies from both sensing modalities for inferring disease proliferation.

Wavelet-domain TI Wiener-like filtering for complex MR data denoising

Magnetic resonance (MR) images are affected by random noises, which degrade many image processing and analysis tasks. It has been shown that the noise in magnitude MR images follows a Rician distribution. Unlike additive Gaussian noise, the noise is signal-dependent, and consequently difficult to reduce, especially in low signal-to-noise ratio (SNR) images. Wirestam et al. in [20] proposed a Wiener-like filtering technique in wavelet-domain to reduce noise before construction of the magnitude MR image. Based on Wirestam's study, we propose a wavelet-domain translation-invariant (TI) Wiener-like filtering algorithm for noise reduction in complex MR data. The proposed denoising algorithm shows the following improvements compared with Wirestam's method: (1) we introduce TI property into the Wiener-like filtering in wavelet-domain to suppress artifacts caused by translations of the signal; (2) we integrate one Stein's Unbiased Risk Estimator (SURE) thresholding with two Wiener-like filters to make the hard-thresholding scale adaptive; and (3) the first Wiener-like filtering is used to filter the original noisy image in which the noise obeys Gaussian distribution and it provides more reasonable results. The proposed algorithm is applied to denoise the real and imaginary parts of complex MR images. To evaluate our proposed algorithm, we conduct extensive denoising experiments using T1-weighted simulated MR images, diffusion-weighted (DW) phantom and in vivo data. We compare our algorithm with other popular denoising methods. The results demonstrate that our algorithm outperforms others in term of both efficiency and robustness.

Analytical estimation of ultrasound properties, thermal diffusivity, and perfusion using magnetic resonance-guided focused ultrasound temperature data

For thermal modeling to play a significant role in treatment planning, monitoring, and control of magnetic resonance-guided focused ultrasound (MRgFUS) thermal therapies, accurate knowledge of ultrasound and thermal properties is essential. This study develops a new analytical solution for the temperature change observed in MRgFUS which can be used with experimental MR temperature data to provide estimates of the ultrasound initial heating rate, Gaussian beam variance, tissue thermal diffusivity, and Pennes perfusion parameter. Simulations demonstrate that this technique provides accurate and robust property estimates that are independent of the beam size, thermal diffusivity, and perfusion levels in the presence of realistic MR noise. The technique is also demonstrated in vivo using MRgFUS heating data in rabbit back muscle. Errors in property estimates are kept less than 5% by applying a third order Taylor series approximation of the perfusion term and ensuring the ratio of the fitting time (the duration of experimental data utilized for optimization) to the perfusion time constant remains less than one.

Monte Carlo-based noise compensation in coil intensity corrected endorectal MRI

Background

Prostate cancer is one of the most common forms of cancer found in males making early diagnosis important. Magnetic resonance imaging (MRI) has been useful in visualizing and localizing tumor candidates and with the use of endorectal coils (ERC), the signal-to-noise ratio (SNR) can be improved. The coils introduce intensity inhomogeneities and the surface coil intensity correction built into MRI scanners is used to reduce these inhomogeneities. However, the correction typically performed at the MRI scanner level leads to noise amplification and noise level variations.

Methods

In this study, we introduce a new Monte Carlo-based noise compensation approach for coil intensity corrected endorectal MRI which allows for effective noise compensation and preservation of details within the prostate. The approach accounts for the ERC SNR profile via a spatially-adaptive noise model for correcting non-stationary noise variations. Such a method is useful particularly for improving the image quality of coil intensity corrected endorectal MRI data performed at the MRI scanner level and when the original raw data is not available.

Results

SNR and contrast-to-noise ratio (CNR) analysis in patient experiments demonstrate an average improvement of 11.7 and 11.2 dB respectively over uncorrected endorectal MRI, and provides strong performance when compared to existing approaches.

Discussion

Experimental results using both phantom and patient data showed that ACER provided strong performance in terms of SNR, CNR, edge preservation, subjective scoring when compared to a number of existing approaches.

Conclusions

A new noise compensation method was developed for the purpose of improving the quality of coil intensity corrected endorectal MRI data performed at the MRI scanner level. We illustrate that promising noise compensation performance can be achieved for the proposed approach, which is particularly important for processing coil intensity corrected endorectal MRI data performed at the MRI scanner level and when the original raw data is not available.

An Optimized LMMSE Based Method for 3D MRI Denoising

Post-acquisition denoising of magnetic resonance (MR) images is an important step to improve any quantitative measurement of the acquired data. In this paper, assuming a Rician noise model, a new filtering method based on the linear minimum mean square error (LMMSE) estimation is introduced, which employs the self-similarity property of the MR data to restore the noise-less signal. This method takes into account the structural characteristics of images and the Bayesian mean square error (Bmse) of the estimator to address the denoising problem. In general, a twofold data processing approach is developed; first, the noisy MR data is processed using a patch-based L(2)-norm similarity measure to provide the primary set of samples required for the estimation process. Afterwards, the Bmse of the estimator is derived as the optimization function to analyze the pre-selected samples and minimize the error between the estimated and the underlying signal. Compared to the LMMSE method and also its recently proposed SNR-adapted realization (SNLMMSE), the optimized way of choosing the samples together with the automatic adjustment of the filtering parameters lead to a more robust estimation performance with our approach. Experimental results show the competitive performance of the proposed method in comparison with related state-of-the-art methods.

The EM Method in a Probabilistic Wavelet-Based MRI Denoising

Human body heat emission and others external causes can interfere in magnetic resonance image acquisition and produce noise. In this kind of images, the noise, when no signal is present, is Rayleigh distributed and its wavelet coefficients can be approximately modeled by a Gaussian distribution. Noiseless magnetic resonance images can be modeled by a Laplacian distribution in the wavelet domain. This paper proposes a new magnetic resonance image denoising method to solve this fact. This method performs shrinkage of wavelet coefficients based on the conditioned probability of being noise or detail. The parameters involved in this filtering approach are calculated by means of the expectation maximization (EM) method, which avoids the need to use an estimator of noise variance. The efficiency of the proposed filter is studied and compared with other important filtering techniques, such as Nowak's, Donoho-Johnstone's, Awate-Whitaker's, and nonlocal means filters, in different 2D and 3D images.

Improving the performance of the prony method using a wavelet domain filter for MRI denoising

The Prony methods are used for exponential fitting. We use a variant of the Prony method for abnormal brain tissue detection in sequences of T 2 weighted magnetic resonance images. Here, MR images are considered to be affected only by Rician noise, and a new wavelet domain bilateral filtering process is implemented to reduce the noise in the images. This filter is a modification of Kazubek's algorithm and we use synthetic images to show the ability of the new procedure to suppress noise and compare its performance with respect to the original filter, using quantitative and qualitative criteria. The tissue classification process is illustrated using a real sequence of T 2 MR images, and the filter is applied to each image before using the variant of the Prony method.

Improved guided image fusion for magnetic resonance and computed tomography imaging

Improved guided image fusion for magnetic resonance and computed tomography imaging is proposed. Existing guided filtering scheme uses Gaussian filter and two-level weight maps due to which the scheme has limited performance for images having noise. Different modifications in filter (based on linear minimum mean square error estimator) and weight maps (with different levels) are proposed to overcome these limitations. Simulation results based on visual and quantitative analysis show the significance of proposed scheme.