Noise estimation in single- and multiple-coil magnetic resonance data based on statistical models

Noise estimation is a challenging task in magnetic resonance imaging (MRI), with applications in quality assessment, filtering or diffusion tensor estimation. Main noise estimators based on the Rician model are revisited and classified in this article, and new useful methods are proposed. Additionally, all the surveyed estimators are extended to the noncentral chi model, which applies to multiple-coil MRI and some important parallel imaging algorithms for accelerated acquisitions. The proposed new noise estimation procedures, based on the distribution of local moments, show better performance in terms of smaller variance and unbiased estimation over a wide range of experiments, with the additional advantage of not needing to explicitly segment the background of the image.     

Denoising 3D MR images by the enhanced non-local means filter for Rician noise

The non-local means (NLM) filter removes noise by calculating the weighted average of the pixels in the global area and shows superiority over existing local filter methods that only consider local neighbor pixels. This filter has been successfully extended from 2D images to 3D images and has been applied to denoising 3D magnetic resonance (MR) images. In this article, a novel filter based on the NLM filter is proposed to improve the denoising effect. Considering the characteristics of Rician noise in the MR images, denoising by the NLM filter is first performed on the squared magnitude images. Then, unbiased correcting is carried out to eliminate the biased deviation. When performing the NLM filter, the weight is calculated based on the Gaussian-filtered image to reduce the disturbance of the noise. The performance of this filter is evaluated by carrying out a qualitative and quantitative comparison of this method with three other filters, namely, the original NLM filter, the unbiased NLM (UNLM) filter and the Rician NLM (RNLM) filter. Experimental results demonstrate that the proposed filter achieves better denoising performance over the other filters being compared.     

Wavelet domain non-linear filtering for MRI denoising

Feature-preserved denoising is of great interest in medical image processing. This article presents a wavelet-based bilateral filtering scheme for noise reduction in magnetic resonance images. Undecimated wavelet transform is employed to provide effective representation of the noisy coefficients. Bilateral filtering of the approximate coefficients improves the denoising efficiency and effectively preserves the edge features. Denoising is done in the square magnitude domain, where the noise tends to be signal independent and is additive. The proposed method has been adapted specifically to Rician noise. The visual and the diagnostic quality of the denoised image is well preserved. The quantitative and the qualitative measures used as the quality metrics demonstrate the ability of the proposed method for noise suppression.     

去除乘性噪音的主成分分析算法

雷达成像系统的进一步应用依赖于对图像中噪音的有效抑制.在目前现有消除噪音方法的基础上,基于图像的局部相似性,结合主成分分析法,提出一种新的有效去除乘性噪音的滤波算法.乘性噪音经对数变换后可转化为加性噪音处理.分析了对数域中噪音的类型.首先在图像的对数域,通过非局部方法选取局部相似块作为训练样本,利用主成分分析法提取出信号的主要特征.然后基于统计理论中最小均方误差估计法给出了一种适用于图像信息的阈值原则.最后分析了变换过程引起的偏差,由对数域的偏估计得到滤波图像.数值实验验证了新算法的有效性.对比于目前提出的变分方法,新算法处理后的图像有更高的信噪比和更好的视觉效果,且具有一定的实用性.     

An MRI denoising method using image data redundancy and local SNR estimation

This paper presents an LMMSE-based method for the three-dimensional (3D) denoising of MR images assuming a Rician noise model. Conventionally, the LMMSE method estimates the noise-less signal values using the observed MR data samples within local neighborhoods. This is not an efficient procedure to deal with this issue while the 3D MR data intrinsically includes many similar samples that can be used to improve the estimation results. To overcome this problem, we model MR data as random fields and establish a principled way which is capable of choosing the samples not only from a local neighborhood but also from a large portion of the given data. To follow the similar samples within the MR data, an effective similarity measure based on the local statistical moments of images is presented. The parameters of the proposed filter are automatically chosen from the estimated local signal-to-noise ratio. To further enhance the denoising performance, a recursive version of the introduced approach is also addressed. The proposed filter is compared with related state-of-the-art filters using both synthetic and real MR datasets. The experimental results demonstrate the superior performance of our proposal in removing the noise and preserving the anatomical structures of MR images.     

Estimation and application of spatially variable noise fields in diffusion tensor imaging

Optimal interpretation of magnetic resonance image content often requires an estimate of the underlying image noise, which is typically realized as a spatially invariant estimate of the noise distribution. This is not an ideal practice in diffusion tensor imaging because the noise distribution is usually spatially varying due to the use of fast imaging and noise suppression techniques. A new estimation approach for spatially varying noise fields (NFs) is proposed in this article. The approach is based on a noise invariance property in scenarios in which more than one image, each with potentially different signal levels, is acquired on each slice, as in diffusion-weighted MRI. This technique leads to improved NF estimates in simulations, phantom experiments and in vivo studies when compared to traditional NF estimators that use regional variability or background intensity histograms. The proposed method reduces the NF estimation error by a factor of 100 in simulations, shows a strong linear correlation (R²=0.99) between theoretical and estimated noise changes in phantoms and demonstrates consistent (<5% variability) NF estimates in vivo. The advantages of spatially varying NF estimation are demonstrated for power analysis, outlier detection and tensor estimation.     

Application of perceptual difference model on regularization techniques of parallel MR imaging

Parallel magnetic resonance imaging through sensitivity encoding using multiple receiver coils has emerged as an effective tool to reduce imaging time or to improve image SNR. The quality of reconstructed images is limited by the inaccurate estimation of the sensitivity map, noise in the acquired k-space data and the ill-conditioned nature of the coefficient matrix. Tikhonov regularization is a popular method to reduce or eliminate the ill-conditioned nature of the problem. In this approach, selection of the regularization map and the regularization parameter is very important. Perceptual difference model (PDM) is a quantitative image quality evaluation tool that has been successfully applied to varieties of MR applications. High correlation between the human rating and PDM score shows that PDM should be suitable to evaluate image quality in parallel MR imaging. By applying PDM, we compared four methods of selecting the regularization map and four methods of selecting the regularization parameter. We found that a regularization map obtained using generalized series (GS) together with a spatially adaptive regularization parameter gave the best reconstructions. PDM was also used as an objective function for optimizing two important parameters in the spatially adaptive method. We conclude that PDM enables one to do comprehensive experiments and that it is an effective tool for designing and optimizing reconstruction methods in parallel MR imaging.     

Laplacian based non-local means denoising of MR images with Rician noise

Magnetic Resonance (MR) image is often corrupted with a complex white Gaussian noise (Rician noise) which is signal dependent. Considering the special characteristics of Rician noise, we carry out nonlocal means denoising on squared magnitude images and compensate the introduced bias. In this paper, we propose an algorithm which not only preserves the edges and fine structures but also performs efficient denoising. For this purpose we have used a Laplacian of Gaussian (LoG) filter in conjunction with a nonlocal means filter (NLM). Further, to enhance the edges and to accelerate the filtering process, only a few similar patches have been preselected on the basis of closeness in edge and inverted mean values. Experiments have been conducted on both simulated and clinical data sets. The qualitative and quantitative measures demonstrate the efficacy of the proposed method.     

用于MR图像非局域均值去噪的边缘增强策略

非局域均值(NLM)滤波有很好的去噪效果并已成功地应用于磁共振图像的去噪中,但与所有去噪方法相同,总是会在一定程度上模糊图像细节. 该文提出将从原始图像中提取出来的高频信息与NLM去噪图像相融合,来还原在去噪过程中丢失的细节. 首先利用一种基于拉普拉斯金字塔的多分辨率方法,从原始图像中提取出包含丰富的边缘信息的高频组分. 然后利用作者提出的一种新的基于SUSAN算子的边缘检测算子产生一幅连续的边缘图,并利用该边缘图将高频组分与NLM方法去噪的图像相融合. 该方法在图像的平滑区域取得了良好的去噪效果,同时可以保留甚至增强图像的细节. 同时,该方法对图像的增强不会导致增强图像中常见的伪影.     

Computational modeling of MR flow imaging by the lattice Boltzmann method and Bloch equation

In this work, a computational model of magnetic resonance (MR) flow imaging is proposed. The first model component provides fluid dynamics maps by applying the lattice Boltzmann method. The second one uses the flow maps and couples MR imaging (MRI) modeling with a new magnetization transport algorithm based on the Eulerian coordinate approach. MRI modeling is based on the discrete time solution of the Bloch equation by analytical local magnetization transformations (exponential scaling and rotations).     

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.     

Nonlocal maximum likelihood estimation method for denoising multiple-coil magnetic resonance images

Effective denoising is vital for proper analysis and accurate quantitative measurements from magnetic resonance (MR) images. Even though many methods were proposed to denoise MR images, only few deal with the estimation of true signal from MR images acquired with phased-array coils. If the magnitude data from phased array coils are reconstructed as the root sum of squares, in the absence of noise correlations and subsampling, the data is assumed to follow a non central-χ distribution. However, when the k-space is subsampled to increase the acquisition speed (as in GRAPPA like methods), noise becomes spatially varying. In this note, we propose a method to denoise multiple-coil acquired MR images. Both the non central-χ distribution and the spatially varying nature of the noise is taken into account in the proposed method. Experiments were conducted on both simulated and real data sets to validate and to demonstrate the effectiveness of the proposed method.     

Brain MRI tissue classification based on local Markov random fields

A new method for tissue classification of brain magnetic resonance images (MRI) of the brain is proposed. The method is based on local image models where each models the image content in a subset of the image domain. With this local modeling approach, the assumption that tissue types have the same characteristics over the brain needs not to be evoked. This is important because tissue type characteristics, such as T1 and T2 relaxation times and proton density, vary across the individual brain and the proposed method offers improved protection against intensity non-uniformity artifacts that can hamper automatic tissue classification methods in brain MRI. A framework in which local models for tissue intensities and Markov Random Field (MRF) priors are combined into a global probabilistic image model is introduced. This global model will be an inhomogeneous MRF and it can be solved by standard algorithms such as iterative conditional modes. The division of the whole image domain into local brain regions possibly having different intensity statistics is realized via sub-volume probabilistic atlases. Finally, the parameters for the local intensity models are obtained without supervision by maximizing the weighted likelihood of a certain finite mixture model. For the maximization task, a novel genetic algorithm almost free of initialization dependency is applied. The algorithm is tested on both simulated and real brain MR images. The experiments confirm that the new method offers a useful improvement of the tissue classification accuracy when the basic tissue characteristics vary across the brain and the noise level of the images is reasonable. The method also offers better protection against intensity non-uniformity artifact than the corresponding method based on a global (whole image) modeling scheme.     

An automatic cerebellum extraction method in T1-weighted brain MR images using an active contour model with a shape prior

Purpose

The objective of this paper was to automatically segment the cerebellum from T1-weighted human brain magnetic resonance (MR) images.

Materials and Methods

The proposed method constructs a cerebellum template using five sets of 3-T MR imaging (MRI) data, which are used to determine the initial position and the shape prior of the cerebellum for the active contour model. Our formulation includes the active contour model with shape prior, which thereby maintains the shape of the template. The proposed active contour model is sequentially applied to sagittal-, coronal- and transverse-view images. To evaluate the proposed method, it is applied to BrainWeb data and a 3-T MRI data set and compared with FreeSurfer with respect to performance assessment metrics.

Results

The segmented cerebellum was compared with the results from FreeSurfer. Using the manually segmented cerebellum as reference, we measured the average Jaccard coefficients of the proposed method, which were 0.882 and 0.885 for the BrainWeb data and 3-T MRI data set, respectively.

Conclusion

We presented the active contour model with shape prior for extracting the cerebellum from T1-weighted brain MR images. The proposed method yielded a robust and accurate segmentation result.     

3D MR image denoising using rough set and kernel PCA method

In this paper, we have presented a two stage method, using kernel principal component analysis (KPCA) and rough set theory (RST), for denoising volumetric MRI data. A rough set theory (RST) based clustering technique has been used for voxel based processing. The method groups similar voxels (3D cubes) using class and edge information derived from noisy input. Each clusters thus formed now represented via basis vector. These vectors now projected into kernel space and PCA is performed in the feature space. This work is motivated by idea that under Rician noise MRI data may be non-linear and kernel mapping will help to define linear separator between these clusters/basis vectors thus used for image denoising. We have further investigated various kernels for Rician noise for different noise levels. The best kernel is then selected on the performance basis over PSNR and structure similarity (SSIM) measures. The work has been compared with state-of-the-art methods under various measures for synthetic and real databases.     

A general strategy for anisotropic diffusion in MR image denoising and enhancement

Anisotropic diffusion (AD) has proven to be very effective in the denoising of magnetic resonance (MR) images. The result of AD filtering is highly dependent on several parameters, especially the conductance parameter. However, there is no automatic method to select the optimal parameter values. This paper presents a general strategy for AD filtering of MR images using an automatic parameter selection method. The basic idea is to estimate the parameters through an optimization step on a synthetic image model, which is different from traditional analytical methods. This approach can be easily applied to more sophisticated diffusion models for better denoising results. We conducted a systematic study of parameter selection for the AD filter, including the dynamic parameter decreasing rate, the parameter selection range for different noise levels and the influence of the image contrast on parameter selection. The proposed approach was validated using both simulated and real MR images. The model image generated using our approach was shown to be highly suitable for the purpose of parameter optimization. The results confirm that our method outperforms most state-of-the-art methods in both quantitative measurement and visual evaluation. By testing on real images with different noise levels, we demonstrated that our method is sufficiently general to be applied to a variety of MR images.     

Multispectral MR images segmentation based on fuzzy knowledge and modified seeded region growing

Magnetic resonance imaging (MRI) is a valuable diagnostic tool in medical science due to its capability for soft-tissue characterization and three-dimensional visualization. One potential application of MRI in clinical practice is brain parenchyma classification and segmentation. Based on fuzzy knowledge and modified seeded region growing, this work proposes a novel image segmentation method, called Fuzzy Knowledge-Based Seeded Region Growing (FKSRG), for multispectral MR images. In this work, fuzzy knowledge includes the fuzzy edge, fuzzy similarity and fuzzy distance, which are obtained from relationships between pixels in multispectral MR images and are applied to the modified seeded regions growing process. In conventional regions merging, the final number of regions is unknown. Therefore, a Target Generation Process is proposed and applied to support conventional regions merging, such that the FKSRG method does not over- or undersegment images. Finally, two image sets, namely, computer-generated phantom images and real MR images, are used in experiments to assess the effectiveness of the proposed FKSRG method. Experimental results demonstrate that the FKSRG method segments multispectral MR images much more effectively than the Functional MRI of the Brain Automated Segmentation Tool, K-means and Support Vector Machine methods.     

Low-field magnetic resonance imaging using multiplicative regularization

In this paper we present a magnetic resonance imaging (MRI) technique that is based on multiplicative regularization. Instead of adding a regularizing objective function to a data fidelity term, we multiply by such a regularizing function. By following this approach, no regularization parameter needs to be determined for each new data set that is acquired. Reconstructions are obtained by iteratively updating the images using short-term conjugate gradient-type update formulas and Polak-Ribière update directions. We show that the algorithm can be used as an image reconstruction algorithm and as a denoising algorithm. We illustrate the performance of the algorithm on two-dimensional simulated low-field MR data that is corrupted by noise and on three-dimensional measured data obtained from a low-field MR scanner. Our reconstruction results show that the algorithm effectively suppresses noise and produces accurate reconstructions even for low-field MR signals with a low signal-to-noise ratio.