MRI noise estimation and denoising using non-local PCA

DeepQA: A Unified Transcriptome-Based Aging Clock Using Deep Neural Networks

Understanding the complex biological process of aging is of great value, especially as it can help develop therapeutics to prolong healthy life. Predicting biological age from gene expression data has shown to be an effective means to quantify aging of a subject, and to identify molecular and cellular biomarkers of aging. A typical approach for estimating biological age, adopted by almost all existing aging clocks, is to train machine learning models only on healthy subjects, but to infer on both healthy and unhealthy subjects. However, the inherent bias in this approach results in inaccurate biological age as shown in this study. Moreover, almost all existing transcriptome-based aging clocks were built around an inefficient procedure of gene selection followed by conventional machine learning models such as elastic nets, linear discriminant analysis etc. To address these limitations, we proposed DeepQA, a unified aging clock based on mixture of experts. Unlike existing methods, DeepQA is equipped with a specially designed Hinge-Mean-Absolute-Error (Hinge-MAE) loss so that it can train on both healthy and unhealthy subjects of multiple cohorts to reduce the bias of inferring biological age of unhealthy subjects. Our experiments showed that DeepQA significantly outperformed existing methods for biological age estimation on both healthy and unhealthy subjects. In addition, our method avoids the inefficient exhaustive search of genes, and provides a novel means to identify genes activated in aging prediction, alternative to such as differential gene expression analysis.

High-resolution perfusion imaging in rodents using pCASL at 9.4 T

Adequate perfusion is critical for maintaining normal brain function and aberrations thereof are hallmarks of many diseases. Pseudo-Continuous Arterial Spin Labeling (pCASL) MRI enables noninvasive quantitative perfusion mapping without contrast agent injection and with a higher signal-to-noise ratio (SNR) than alternative methods. Despite its great potential, pCASL remains challenging, unstable, and relatively low-resolution in rodents - especially in mice - thereby limiting the investigation of perfusion properties in many transgenic or other relevant rodent models of disease. Here, we address this gap by developing a novel experimental setup for high-resolution pCASL imaging in mice and combining it with the enhanced SNR of cryogenic probes. We show that our new experimental setup allows for optimal positioning of the carotids within the cryogenic coil, rendering labeling reproducible. With the proposed methodology, we managed to increase the spatial resolution of pCASL perfusion images by a factor of 15 in mice; a factor of 6 in rats is gained compared to the state of the art just by virtue of the cryogenic coil. We also show that the improved pCASL perfusion imaging allows much better delineation of specific brain areas, both in healthy animals as well as in rat and mouse models of stroke. Our results bode well for future high-definition pCASL perfusion imaging in rodents.

Reducing thermal noise in high-resolution quantitative magnetic resonance imaging rotating frame relaxation mapping of the human brain at 3 T

Quantitative maps of rotating frame relaxation (RFR) time constants are sensitive and useful magnetic resonance imaging tools with which to evaluate tissue integrity in vivo. However, to date, only moderate image resolutions of 1.6 x 1.6 x 3.6 mm3 have been used for whole-brain coverage RFR mapping in humans at 3 T. For more precise morphometrical examinations, higher spatial resolutions are desirable. Towards achieving the long-term goal of increasing the spatial resolution of RFR mapping without increasing scan times, we explore the use of the recently introduced Transform domain NOise Reduction with DIstribution Corrected principal component analysis (T-NORDIC) algorithm for thermal noise reduction. RFR acquisitions at 3 T were obtained from eight healthy participants (seven males and one female) aged 52 ± 20 years, including adiabatic T1ρ, T2ρ, and nonadiabatic Relaxation Along a Fictitious Field (RAFF) in the rotating frame of rank n = 4 (RAFF4) with both 1.6 x 1.6 x 3.6 mm3 and 1.25 x 1.25 x 2 mm3 image resolutions. We compared RFR values and their confidence intervals (CIs) obtained from fitting the denoised versus nondenoised images, at both voxel and regional levels separately for each resolution and RFR metric. The comparison of metrics obtained from denoised versus nondenoised images was performed with a two-sample paired t-test and statistical significance was set at p less than 0.05 after Bonferroni correction for multiple comparisons. The use of T-NORDIC on the RFR images prior to the fitting procedure decreases the uncertainty of parameter estimation (lower CIs) at both spatial resolutions. The effect was particularly prominent at high-spatial resolution for RAFF4. Moreover, T-NORDIC did not degrade map quality, and it had minimal impact on the RFR values. Denoising RFR images with T-NORDIC improves parameter estimation while preserving the image quality and accuracy of all RFR maps, ultimately enabling high-resolution RFR mapping in scan times that are suitable for clinical settings.

Deep learning methods for 3D magnetic resonance image denoising, bias field and motion artifact correction: a comprehensive review

Magnetic resonance imaging (MRI) provides detailed structural information of the internal body organs and soft tissue regions of a patient in clinical diagnosis for disease detection, localization, and progress monitoring. MRI scanner hardware manufacturers incorporate various post-acquisition image-processing techniques into the scanner's computer software tools for different post-processing tasks. These tools provide a final image of adequate quality and essential features for accurate clinical reporting and predictive interpretation for better treatment planning. Different post-acquisition image-processing tasks for MRI quality enhancement include noise removal, motion artifact reduction, magnetic bias field correction, and eddy electric current effect removal. Recently, deep learning (DL) methods have shown great success in many research fields, including image and video applications. DL-based data-driven feature-learning approaches have great potential for MR image denoising and image-quality-degrading artifact correction. Recent studies have demonstrated significant improvements in image-analysis tasks using DL-based convolutional neural network techniques. The promising capabilities and performance of DL techniques in various problem-solving domains have motivated researchers to adapt DL methods to medical image analysis and quality enhancement tasks. This paper presents a comprehensive review of DL-based state-of-the-art MRI quality enhancement and artifact removal methods for regenerating high-quality images while preserving essential anatomical and physiological feature maps without destroying important image information. Existing research gaps and future directions have also been provided by highlighting potential research areas for future developments, along with their importance and advantages in medical imaging.

Denoising complex-valued diffusion MR images using a two-step non-local principal component analysis approach

Purpose

to propose a two-step non-local principal component analysis (PCA) method and demonstrate its utility for denoising diffusion tensor MRI (DTI) with a few diffusion directions.

Methods

A two-step denoising pipeline was implemented to ensure accurate patch selection even with high noise levels and was coupled with data preprocessing for g-factor normalization and phase stabilization before data denoising with a non-local PCA algorithm. At the heart of our proposed pipeline was the use of a data-driven optimal shrinkage algorithm to manipulate the singular values in a way that would optimally estimate the noise-free signal. Our approach's denoising performances were evaluated using simulation and in-vivo human data experiments. The results were compared to those obtained with existing local-PCA-based methods.

Results

In both simulation and human data experiments, our approach substantially enhanced image quality relative to the noisy counterpart, yielding improved performances for estimation of relevant DTI metrics. It also outperformed existing local-PCA-based methods in reducing noise while preserving anatomic details. It also led to improved whole-brain tractography relative to the noisy counterpart.

Conclusion

The proposed denoising method has the utility for improving image quality for DTI with reduced diffusion directions and is believed to benefit many applications especially those aiming to achieve quality parametric mapping using only a few image volumes.

MRI Denoising Using Pixel-Wise Threshold Selection

Magnetic resonance imaging (MRI) has emerged as a promising technique for non-invasive medical imaging. The primary challenge in MRI is the trade-off between image visual quality and acquisition time. Current MRI image denoising algorithms employ global thresholding to denoise the whole image, which leads to inadequate denoising or image distortion. This study introduces a novel pixel-wise (localized) thresholding approach of singular vectors, obtained from singular value decomposition, to denoise magnetic resonance (MR) images. The pixel-wise thresholding of singular vectors is performed using separate singular values as thresholds at each pixel, which is advantageous given the spatial noise variation throughout the image. The method presented is validated on MR images of a standard phantom approved by the magnetic resonance accreditation program (MRAP). The denoised images display superior visual quality and recover minute structural information otherwise suppressed in the noisy image. The increase in peak-signal-to-noise-ratio (PSNR) and contrast-to-noise-ratio (CNR) values of ≥ 18% and ≥ 200% of the denoised images, respectively, imply efficient noise removal and visual quality enhancement. The structural similarity index (SSIM) of ≥ 0.95 for denoised images indicates that the crucial structural information is recovered through the presented method. A comparison with the standard filtering methods widely used for MRI denoising establishes the superior performance of the presented method. The presented pixel-wise denoising technique reduces the scan time by 2-3 times and has the potential to be integrated into any MRI system to obtain faster and better quality images.

Fast reconstruction of SMS bSSFP myocardial perfusion images using noise map estimation network (NoiseMapNet): a head-to-head comparison with parallel imaging and iterative reconstruction

Background

Simultaneous multi-slice (SMS) bSSFP imaging enables stress myocardial perfusion imaging with high spatial resolution and increased spatial coverage. Standard parallel imaging techniques (e.g., TGRAPPA) can be used for image reconstruction but result in high noise level. Alternatively, iterative reconstruction techniques based on temporal regularization (ITER) improve image quality but are associated with reduced temporal signal fidelity and long computation time limiting their online use. The aim is to develop an image reconstruction technique for SMS-bSSFP myocardial perfusion imaging combining parallel imaging and image-based denoising using a novel noise map estimation network (NoiseMapNet), which preserves both sharpness and temporal signal profiles and that has low computational cost.

Methods

The proposed reconstruction of SMS images consists of a standard temporal parallel imaging reconstruction (TGRAPPA) with motion correction (MOCO) followed by image denoising using NoiseMapNet. NoiseMapNet is a deep learning network based on a 2D Unet architecture and aims to predict a noise map from an input noisy image, which is then subtracted from the noisy image to generate the denoised image. This approach was evaluated in 17 patients who underwent stress perfusion imaging using a SMS-bSSFP sequence. Images were reconstructed with (a) TGRAPPA with MOCO (thereafter referred to as TGRAPPA), (b) iterative reconstruction with integrated motion compensation (ITER), and (c) proposed NoiseMapNet-based reconstruction. Normalized mean squared error (NMSE) with respect to TGRAPPA, myocardial sharpness, image quality, perceived SNR (pSNR), and number of diagnostic segments were evaluated.

Results

NMSE of NoiseMapNet was lower than using ITER for both myocardium (0.045 ± 0.021 vs. 0.172 ± 0.041, p < 0.001) and left ventricular blood pool (0.025 ± 0.014 vs. 0.069 ± 0.020, p < 0.001). There were no significant differences between all methods for myocardial sharpness (p = 0.77) and number of diagnostic segments (p = 0.36). ITER led to higher image quality than NoiseMapNet/TGRAPPA (2.7 ± 0.4 vs. 1.8 ± 0.4/1.3 ± 0.6, p < 0.001) and higher pSNR than NoiseMapNet/TGRAPPA (3.0 ± 0.0 vs. 2.0 ± 0.0/1.3 ± 0.6, p < 0.001). Importantly, NoiseMapNet yielded higher pSNR (p < 0.001) and image quality (p < 0.008) than TGRAPPA. Computation time of NoiseMapNet was only 20s for one entire dataset.

Conclusion

NoiseMapNet-based reconstruction enables fast SMS image reconstruction for stress myocardial perfusion imaging while preserving sharpness and temporal signal profiles.

New non-local mean methods for MRI denoising based on global self-similarity between values

Magnetic resonance imaging (MRI) is a non-invasive medical imaging technique that provides high-resolution 3D images and valuable insights into human tissue conditions. Even at present, the refinement of denoising methods for MRI remains a crucial concern for improving the quality of the images. This study aims to improve the prefiltered rotationally invariant non-local principal component analysis (PRI-NL-PCA) algorithm. We relaxed the original restrictions using particle swarm optimization to determine optimal parameters for the PCA part of the original algorithm. In addition, we adjusted the prefiltered rotationally invariant non-local mean (PRI-NLM) part by traversing the signal intensities of voxels instead of their spatial positions to reduce duplicate calculations and expand the search volume to the whole image when estimating voxels' signal intensities. The new method demonstrated superior denoising performance compared to the original approach. Moreover, in most cases, the new algorithm ran faster. Furthermore, our proposed method can also be applied to process Gaussian noise in natural images and has the potential to enhance other NLM-based denoising algorithms.

Brain fiber structure estimation based on principal component analysis and RINLM filter

Diffusion magnetic resonance imaging is a technique for non-invasive detection of microstructure in the white matter of the human brain, which is widely used in neuroscience research of the brain. However, diffusion-weighted images(DWI) are sensitive to noise, which affects the subsequent reconstruction of fiber orientation direction, microstructural parameter estimation and fiber tracking. In order to better eliminate the noise in diffusion-weighted images, this study proposes a noise reduction method combining Marchenko-Pastur principal component analysis(MPPCA) and rotation-invariant non-local means filter(RINLM) to further remove residual noise and preserve the image texture detail information. In this study, the algorithm is applied to the fiber structure and the prevailing microstructural models within the human brain voxels based on simulated and real human brain datasets. Experimental comparisons between the proposed method and the state-of-the-art methods are performed in single-fiber, multi-fiber, crossed and curved-fiber regions as well as in different microstructure estimation models. Results demonstrated the superior performance of the proposed method in denoising DWI data, which can reduce the angular error in fiber orientation reconstruction to obtain more valid fiber structure estimation and enable more complete fiber tracking trajectories with higher coverage. Meanwhile, the method reduces the estimation errors of various white matter microstructural parameters and verifies the performance of the method in white matter microstructure estimation.

Survey of Denoising, Segmentation and Classification of Pancreatic Cancer Imaging

BACKGROUND

Pancreatic cancer is one of the most serious problems that has taken many lives worldwide. The diagnostic procedure using the traditional approaches was manual by visually analyzing the large volumes of the dataset, making it time-consuming and prone to subjective errors. Hence the need for the computer-aided diagnosis system (CADs) emerged that comprises the machine and deep learning approaches for denoising, segmentation and classification of pancreatic cancer.

INTRODUCTION

There are different modalities used for the diagnosis of pancreatic cancer, such as Positron Emission Tomography/Computed Tomography (PET/CT), Magnetic Resonance Imaging (MRI), Multiparametric-MRI (Mp-MRI), Radiomics and Radio-genomics. Although these modalities gave remarkable results in diagnosis on the basis of different criteria. CT is the most commonly used modality that produces detailed and fine contrast images of internal organs of the body. However, it may also contain a certain amount of gaussian and rician noise that is necessary to be preprocessed before segmentation of the required region of interest (ROI) from the images and classification of cancer.

METHOD

This paper analyzes different methodologies used for the complete diagnosis of pancreatic cancer, including the denoising, segmentation and classification, along with the challenges and future scope for the diagnosis of pancreatic cancer.

RESULT

Various filters are used for denoising and image smoothening and filters as gaussian scale mixture process, non-local means, median filter, adaptive filter and average filter have been used more for better results.

CONCLUSION

In terms of segmentation, atlas based region-growing method proved to give better results as compared to the state of the art whereas, for the classification, deep learning approaches outperformed other methodologies to classify the images as cancerous and non- cancerous. These methodologies have proved that CAD systems have become a better solution to the ongoing research proposals for the detection of pancreatic cancer worldwide.

Boosting quantification accuracy of chemical exchange saturation transfer MRI with a spatial-spectral redundancy-based denoising method

Chemical exchange saturation transfer (CEST) is a versatile technique that enables noninvasive detections of endogenous metabolites present in low concentrations in living tissue. However, CEST imaging suffers from an inherently low signal-to-noise ratio (SNR) due to the decreased water signal caused by the transfer of saturated spins. This limitation challenges the accuracy and reliability of quantification in CEST imaging. In this study, a novel spatial-spectral denoising method, called BOOST (suBspace denoising with nOnlocal lOw-rank constraint and Spectral local-smooThness regularization), was proposed to enhance the SNR of CEST images and boost quantification accuracy. More precisely, our method initially decomposes the noisy CEST images into a low-dimensional subspace by leveraging the global spectral low-rank prior. Subsequently, a spatial nonlocal self-similarity prior is applied to the subspace-based images. Simultaneously, the spectral local-smoothness property of Z-spectra is incorporated by imposing a weighted spectral total variation constraint. The efficiency and robustness of BOOST were validated in various scenarios, including numerical simulations and preclinical and clinical conditions, spanning magnetic field strengths from 3.0 to 11.7 T. The results demonstrated that BOOST outperforms state-of-the-art algorithms in terms of noise elimination. As a cost-effective and widely available post-processing method, BOOST can be easily integrated into existing CEST protocols, consequently promoting accuracy and reliability in detecting subtle CEST effects.

MP-PCA denoising of fMRI time-series data can lead to artificial activation "spreading"

MP-PCA denoising has become the method of choice for denoising MRI data since it provides an objective threshold to separate the signal components from unwanted thermal noise components. In rodents, thermal noise in the coils is an important source of noise that can reduce the accuracy of activation mapping in fMRI. Further confounding this problem, vendor data often contains zero-filling and other post-processing steps that may violate MP-PCA assumptions. Here, we develop an approach to denoise vendor data and assess activation "spreading" caused by MP-PCA denoising in rodent task-based fMRI data. Data was obtained from N = 3 mice using conventional multislice and ultrafast acquisitions (1 s and 50 ms temporal resolution, respectively), during visual stimulation. MP-PCA denoising produced SNR gains of 64% and 39% and Fourier Spectral Amplitude (FSA) increases in BOLD maps of 9% and 7% for multislice and ultrafast data, respectively, when using a small [2 2] denoising window. Larger windows provided higher SNR and FSA gains with increased spatial extent of activation that may or may not represent real activation. Simulations showed that MP-PCA denoising can incur activation "spreading" with increased false positive rate and smoother functional maps due to local "bleeding" of principal components, and that the optimal denoising window for improved specificity of functional mapping, based on Dice score calculations, depends on the data's tSNR and functional CNR. This "spreading" effect applies also to another recently proposed low-rank denoising method (NORDIC), although to a lesser degree. Our results bode well for enhancing spatial and/or temporal resolution in future fMRI work, while taking into account the sensitivity/specificity trade-offs of low-rank denoising methods.

Parametric MRI Detects Aristolochic Acid Induced Acute Kidney Injury

Exposure to aristolochic acid (AA) is of increased concern due to carcinogenic and nephrotoxic effects, and incidence of aristolochic acid nephropathy (AAN) is increasing. This study characterizes renal alterations during the acute phase of AAN using parametric magnetic resonance imaging (MRI). An AAN and a control group of male Wistar rats received administration of aristolochic acid I (AAI) and polyethylene glycol (PEG), respectively, for six days. Both groups underwent MRI before and 2, 4 and 6 days after AAI or PEG administration. T₂ relaxation times and apparent diffusion coefficients (ADCs) were determined for four renal layers. Serum creatinine levels (sCr) and blood urea nitrogen (BUN) were measured. Tubular injury scores (TIS) were evaluated based on histologic findings. Increased T₂ values were detected since day 2 in the AAN group, but decreased ADCs and increased sCr levels and BUN were not detected until day 4. Significant linear correlations were observed between T₂ of the cortex and the outer stripe of outer medulla and TIS. Our results demonstrate that parametric MRI facilitates early detection of renal injury induced by AAI in a rat model. T₂ mapping may be a valuable tool for assessing kidney injury during the acute phase of AAN.

Blind Additive Gaussian White Noise Level Estimation from a Single Image by Employing Chi-Square Distribution

The additive Gaussian white noise (AGWN) level in real-life images is usually unknown, for which the empirical setting will make the denoising methods over-smooth fine structures or remove noise incompletely. The previous noise level estimation methods are easily lost in accurately estimating them from images with complicated structures. To cope with this issue, we propose a novel noise level estimation scheme based on Chi-square distribution, including the following key points: First, a degraded image is divided into many image patches through a sliding window. Then, flat patches are selected by using a patch selection strategy on the gradient maps of those image patches. Next, the initial noise level is calculated by employing Chi-square distribution on the selected flat patches. Finally, the stable noise level is optimized by an iterative strategy. Quantitative, with association, to qualitative results of experiments on synthetic real-life images validate that the proposed noise level estimation method is effective and even superior to the state-of-the-art methods. Extensive experiments on noise removal using BM3D further illustrate that the proposed noise level estimation method is more beneficial for achieving favorable denoising performance with detail preservation.

Lossless Medical Image Compression by Using Difference Transform

This paper introduces a new method of compressing digital images by using the Difference Transform applied in medical imaging. The Difference Transform algorithm performs the decorrelation process of image data, and in this way improves the encoding process, achieving a file with a smaller size than the original. The proposed method proves to be competitive and in many cases better than the standards used for medical images such as TIFF or PNG. In addition, the Difference Transform can replace other transforms like Cosine or Wavelet.

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Sensitivity schemes for dynamic glucose-enhanced magnetic resonance imaging to detect glucose uptake and clearance in mouse brain at 3 T

We investigated three dynamic glucose-enhanced (DGE) MRI methods for sensitively monitoring glucose uptake and clearance in both brain parenchyma and cerebrospinal fluid (CSF) at clinical field strength (3 T). By comparing three sequences, namely, Carr-Purcell-Meiboom-Gill (CPMG), on-resonance variable delay multipulse (onVDMP), and on-resonance spin-lock (onSL), a high-sensitivity DGE MRI scheme with truncated multilinear singular value decomposition (MLSVD) denoising was proposed. The CPMG method showed the highest sensitivity in detecting the parenchymal DGE signal among the three methods, while both onVDMP and onSL were more robust for CSF DGE imaging. Here, onVDMP was applied for CSF imaging, as it displayed the best stability of the DGE results in this study. The truncated MLSVD denoising method was incorporated to further improve the sensitivity. The proposed DGE MRI scheme was examined in mouse brain with 50%/25%/12.5% w/w D-glucose injections. The results showed that this combination could detect DGE signal changes from the brain parenchyma and CSF with as low as a 12.5% w/w D-glucose injection. The proposed DGE MRI schemes could sensitively detect the glucose signal change from brain parenchyma and CSF after D-glucose injection at a clinically relevant concentration, demonstrating high potential for clinical translation.

Artifact- and content-specific quality assessment for MRI with image rulers

In clinical practice MR images are often first seen by radiologists long after the scan. If image quality is inadequate either patients have to return for an additional scan, or a suboptimal interpretation is rendered. An automatic image quality assessment (IQA) would enable real-time remediation. Existing IQA works for MRI give only a general quality score, agnostic to the cause of and solution to low-quality scans. Furthermore, radiologists' image quality requirements vary with the scan type and diagnostic task. Therefore, the same score may have different implications for different scans. We propose a framework with multi-task CNN model trained with calibrated labels and inferenced with image rulers. Labels calibrated by human inputs follow a well-defined and efficient labeling task. Image rulers address varying quality standards and provide a concrete way of interpreting raw scores from the CNN. The model supports assessments of two of the most common artifacts in MRI: noise and motion. It achieves accuracies of around 90%, 6% better than the best previous method examined, and 3% better than human experts on noise assessment. Our experiments show that label calibration, image rulers, and multi-task training improve the model's performance and generalizability.

Super Resolution of Magnetic Resonance Images

In this work, novel denoising and super resolution (SR) approaches for magnetic resonance (MR) images are addressed, and are integrated in a unified framework, which do not require example low resolution (LR)/high resolution (HR)/cross-modality/noise-free images and prior information of noise-noise variance. The proposed method categorizes the patches as either smooth or textured and then denoises them by deploying different denoising strategies for efficient denoising. The denoising algorithm is integrated into the SR approach, which uses a gradient profile-based constraint in a sparse representation-based framework to improve the resolution of MR images with reduced smearing of image details. This constraint regularizes the estimation of HR images such that the estimated HR image has gradient profiles similar to the gradient profiles of the original HR image. For this, the gradient profile sharpness (GPS) values of an unknown HR image are estimated using an approximated piece-wise linear relation among GPS values of LR and upsampled LR images. The experiments are performed on three different publicly available datasets. The proposed SR approach outperforms the existing unsupervised SR approach addressed for real MR images that exploits low rank and total variation (LRTV) regularization, by an average peak signal to noise ratio (PSNR) of 0.73 dB and 0.38 dB for upsampling factors 2 and 3, respectively. For the super resolution of noisy real MR images (degraded with 2% noise), the proposed approach outperforms the LRTV approach by an average PSNR of 0.54 dB and 0.46 dB for upsampling factors 2 and 3, respectively. The qualitative analysis is shown for real MR images from healthy subjects and subjects with Alzheimer's disease and structural deformity, i.e., cavernoma.

Joint low-rank prior and difference of Gaussian filter for magnetic resonance image denoising

The low-rank matrix approximation (LRMA) is an efficient image denoising method to reduce additive Gaussian noise. However, the existing low-rank matrix approximation does not perform well in terms of Rician noise removal for magnetic resonance imaging (MRI). To this end, we propose a novel MR image denoising approach based on the extended difference of Gaussian (DoG) filter and nonlocal low-rank regularization. In the proposed method, a novel nonlocal self-similarity evaluation with the tight frame is exploited to improve the patch matching. To remove the Rician noise and preserve the edge details, the extended DoG filter is exploited to the nonlocal low-rank regularization model. The experimental results demonstrate that the proposed method can preserve more edge and fine structures while removing noise in MR image as compared with some of the existing methods.

Lowering latency and processing burden in computational imaging through dimensionality reduction of the sensing matrix

Recent demonstrations have shown that frequency-diverse computational imaging systems can greatly simplify conventional architectures developed for imaging by transferring constraints into the digital layer. Here, in order to limit the latency and processing burden involved in image reconstruction, we propose to truncate insignificant principal components of the sensing matrix that links the measurements to the scene to be imaged. In contrast to recent work using principle component analysis to synthesize scene illuminations, our generic approach is fully unsupervised and is applied directly to the sensing matrix. We impose no restrictions on the type of imageable scene, no training data is required, and no actively reconfigurable radiating apertures are employed. This paper paves the way to the constitution of a new degree of freedom in image reconstructions, allowing one to place the performance emphasis either on image quality or latency and computational burden. The application of such relaxations will be essential for widespread deployment of computational microwave and millimeter wave imagers in scenarios such as security screening. We show in this specific context that it is possible to reduce both the processing time and memory consumption with a minor impact on the quality of the reconstructed images.

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.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

Classifying MRI motion severity using a stacked ensemble approach

Motion artifacts are a common occurrence in Magnetic Resonance Imaging exam. Motion during acquisition has a profound impact on workflow efficiency, often requiring a repeat of sequences. Furthermore, motion artifacts may escape notice by technologists, only to be revealed at the time of reading by the radiologists, affecting their diagnostic quality. There is a paucity of clinical tools to identify and quantitatively assess the severity of motion artifacts in MRI. An image with subtle motion may still have diagnostic value, while severe motion may be uninterpretable by radiologists and requires the exam to be repeated. Therefore, a tool for the automatic identification of motion artifacts would aid in maintaining diagnostic quality, while potentially driving workflow efficiencies. Here we aim to quantify the severity of motion artifacts from MRI images using deep learning. Impact of subject movement parameters like displacement and rotation on image quality is also studied. A state-of-the-art, stacked ensemble model was developed to classify motion artifacts into five levels (no motion, slight, mild, moderate and severe) in brain scans. The stacked ensemble model is able to robustly predict rigid-body motion severity across different acquisition parameters, including T1-weighted and T2-weighted slices acquired in different anatomical planes. The ensemble model with XGBoost metalearner achieves 91.6% accuracy, 94.8% area under the curve, 90% Cohen's Kappa, and is observed to be more accurate and robust than the individual base learners.

Image Denoising Using Non-Local Means (NLM) Approach in Magnetic Resonance (MR) Imaging: A Systematic Review

The non-local means (NLM) noise reduction algorithm is well known as an excellent technique for removing noise from a magnetic resonance (MR) image to improve the diagnostic accuracy. In this study, we undertook a systematic review to determine the effectiveness of the NLM noise reduction algorithm in MR imaging. A systematic literature search was conducted of three databases of publications dating from January 2000 to March 2020; of the 82 publications reviewed, 25 were included in this study. The subjects were categorized into four major frameworks and analyzed for each research result. Research in NLM noise reduction for MR images has been increasing worldwide; however, it was found to have slightly decreased since 2016. It was found that the NLM technique was most frequently used on brain images taken using the general MR imaging technique; these were most frequently performed during simultaneous real and simulated experimental studies. In particular, comparison parameters were frequently used to evaluate the effectiveness of the algorithm on MR images. The ultimate goal is to provide an accurate method for the diagnosis of disease, and our conclusion is that the NLM noise reduction algorithm is a promising method of achieving this goal.

Deep learning of brain magnetic resonance images: A brief review

Magnetic resonance imaging (MRI) is one of the most popular techniques in brain science and is important for understanding brain function and neuropsychiatric disorders. However, the processing and analysis of MRI is not a trivial task with lots of challenges. Recently, deep learning has shown superior performance over traditional machine learning approaches in image analysis. In this survey, we give a brief review of the recent popular deep learning approaches and their applications in brain MRI analysis. Furthermore, popular brain MRI databases and deep learning tools are also introduced. The strength and weaknesses of different approaches are addressed, and challenges as well as future directions are also discussed.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

A Non-Local Low-Rank Algorithm for Sub-Bottom Profile Sonar Image Denoising

Due to the influence of equipment instability and surveying environment, scattering echoes and other factors, it is sometimes difficult to obtain high-quality sub-bottom profile (SBP) images by traditional denoising methods. In this paper, a novel SBP image denoising method is developed for obtaining underlying clean images based on a non-local low-rank framework. Firstly, to take advantage of the inherent layering structures of the SBP image, a direction image is obtained and used as a guidance image. Secondly, the robust guidance weight for accurately selecting the similar patches is given. A novel denoising method combining the weight and a non-local low-rank filtering framework is proposed. Thirdly, after discussing the filtering parameter settings, the proposed method is tested in actual measurements of sub-bottom, both in deep water and shallow water. Experimental results validate the excellent performance of the proposed method. Finally, the proposed method is verified and compared with other methods quantificationally based on the synthetic images and has achieved the total average peak signal-to-noise ratio (PSNR) of 21.77 and structural similarity index (SSIM) of 0.573, which is far better than other methods.

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.

Denoise magnitude diffusion magnetic resonance images via variance-stabilizing transformation and optimal singular-value manipulation

Although shown to have a great utility for a wide range of neuroscientific and clinical applications, diffusion-weighted magnetic resonance imaging (dMRI) faces a major challenge of low signal-to-noise ratio (SNR), especially when pushing the spatial resolution for improved delineation of brain's fine structure or increasing the diffusion weighting for increased angular contrast or both. Here, we introduce a comprehensive denoising framework for denoising magnitude dMRI. The framework synergistically combines the variance stabilizing transform (VST) with optimal singular value manipulation. The purpose of VST is to transform the Rician data to Gaussian-like data so that an asymptotically optimal singular value manipulation strategy tailored for Gaussian data can be used. The output of the framework is the estimated underlying diffusion signal for each voxel in the image domain. The usefulness of the proposed framework for denoising magnitude dMRI is demonstrated using both simulation and real-data experiments. Our results show that the proposed denoising framework can significantly improve SNR across the entire brain, leading to substantially enhanced performances for estimating diffusion tensor related indices and for resolving crossing fibers when compared to another competing method. More encouragingly, the proposed method when used to denoise a single average of 7 ​Tesla Human Connectome Project-style diffusion acquisition provided comparable performances relative to those achievable with ten averages for resolving multiple fiber populations across the brain. As such, the proposed denoising method is expected to have a great utility for high-quality, high-resolution whole-brain dMRI, desirable for many neuroscientific and clinical applications.

Non-Local SVD Denoising of MRI Based on Sparse Representations

Magnetic Resonance (MR) Imaging is a diagnostic technique that produces noisy images, which must be filtered before processing to prevent diagnostic errors. However, filtering the noise while keeping fine details is a difficult task. This paper presents a method, based on sparse representations and singular value decomposition (SVD), for non-locally denoising MR images. The proposed method prevents blurring, artifacts, and residual noise. Our method is composed of three stages. The first stage divides the image into sub-volumes, to obtain its sparse representation, by using the KSVD algorithm. Then, the global influence of the dictionary atoms is computed to upgrade the dictionary and obtain a better reconstruction of the sub-volumes. In the second stage, based on the sparse representation, the noise-free sub-volume is estimated using a non-local approach and SVD. The noise-free voxel is reconstructed by aggregating the overlapped voxels according to the rarity of the sub-volumes it belongs, which is computed from the global influence of the atoms. The third stage repeats the process using a different sub-volume size for producing a new filtered image, which is averaged with the previously filtered images. The results provided show that our method outperforms several state-of-the-art methods in both simulated and real data.

Analysis of the Quantization Noise in Discrete Wavelet Transform Filters for 3D Medical Imaging

Denoising and compression of 2D and 3D images are important problems in modern medical imaging systems. Discrete wavelet transform (DWT) is used to solve them in practice. We analyze the quantization noise effect in coefficients of DWT filters for 3D medical imaging in this paper. The method for wavelet filters coefficients quantizing is proposed, which allows minimizing resources in hardware implementation by simplifying rounding operations. We develop the method for estimating the maximum error of 3D grayscale and color images DWT with various bits per color (BPC). The dependence of the peak signal-to-noise ratio (PSNR) of the images processing result on wavelet used, the effective bit-width of filters coefficients and BPC is revealed. We derive formulas for determining the minimum bit-width of wavelet filters coefficients that provide a high (PSNR ≥ 40 dB for images with 8 BPC, for example) and maximum (PSNR = ∞ dB) quality of 3D medical imaging by DWT depending on wavelet used. The experiments of 3D tomographic images processing confirmed the accuracy of theoretical analysis. All data are presented in the fixed-point format in the proposed method of 3D medical images DWT. It is making possible efficient, from the point of view of hardware and time resources, the implementation for image denoising and compression on modern devices such as field-programmable gate arrays and application-specific integrated circuits.

Adaptive denoising for chemical exchange saturation transfer MR imaging

High image signal-to-noise ratio (SNR) is required to reliably detect the inherently small chemical exchange saturation transfer (CEST) effects in vivo. In this study, it was demonstrated that identifying spectral redundancies of CEST data by principal component analysis (PCA) in combination with an appropriate data-driven extraction of relevant information can be used for an effective and robust denoising of CEST spectra. The relationship between the number of relevant principal components and SNR was studied on fitted in vivo Z-spectra with artificially introduced noise. Three different data-driven criteria to automatically determine the optimal number of necessary components were investigated. In addition, these criteria facilitate straightforward assessment of data quality that could provide guidance for CEST MR protocols in terms of SNR. Insights were applied to achieve a robust denoising of highly sampled low power Z-spectra of the human brain at 3 and 7 T. The median criterion provided the best estimation for the optimal number of components consistently for all three investigated artificial noise levels. Application of the denoising technique to in vivo data revealed a considerable increase in image quality for the amide and rNOE contrast with a considerable SNR gain. At 7 T the denoising capability was quantified to be comparable or even superior to an averaging of six measurements. The proposed denoising algorithm enables an efficient and robust denoising of CEST data by combining PCA with appropriate data-driven truncation criteria. With this generally applicable technique at hand, small CEST effects can be reliably detected without the need for repeated measurements.

Multi-shot diffusion-weighted MRI reconstruction with magnitude-based spatial-angular locally low-rank regularization (SPA-LLR)

PURPOSE

To resolve the motion-induced phase variations in multi-shot multi-direction diffusion-weighted imaging (DWI) by applying regularization to magnitude images.

THEORY AND METHODS

A nonlinear model was developed to estimate phase and magnitude images separately. A locally low-rank regularization (LLR) term was applied to the magnitude images from all diffusion-encoding directions to exploit the spatial and angular correlation. In vivo experiments with different resolutions and b-values were performed to validate the proposed method.

RESULTS

The proposed method significantly reduces the noise level compared to the conventional reconstruction method and achieves submillimeter (0.8mm and 0.9mm isotropic resolutions) DWI with a b-value of 1,000  s / mm 2 and 1-mm isotropic DWI with a b-value of 2,000  s / mm 2 without modification of the sequence.

CONCLUSIONS

A joint reconstruction method with spatial-angular LLR regularization on magnitude images substantially improves multi-direction DWI reconstruction, simultaneously removes motion-induced phase artifacts, and denoises images.

Towards an optimised processing pipeline for diffusion magnetic resonance imaging data: Effects of artefact corrections on diffusion metrics and their age associations in UK Biobank

Increasing interest in the structural and functional organisation of the human brain encourages the acquisition of big data sets comprising multiple neuroimaging modalities, often accompanied by additional information obtained from health records, cognitive tests, biomarkers and genotypes. Diffusion weighted magnetic resonance imaging data enables a range of promising imaging phenotypes probing structural connections as well as macroanatomical and microstructural properties of the brain. The reliability and biological sensitivity and specificity of diffusion data depend on processing pipeline. A state-of-the-art framework for data processing facilitates cross-study harmonisation and reduces pipeline-related variability. Using diffusion magnetic resonance imaging (MRI) data from 218 individuals in the UK Biobank, we evaluate the effects of different processing steps that have been suggested to reduce imaging artefacts and improve reliability of diffusion metrics. In lack of a ground truth, we compared diffusion metric sensitivity to age between pipelines. By comparing distributions and age sensitivity of the resulting diffusion metrics based on different approaches (diffusion tensor imaging, diffusion kurtosis imaging and white matter tract integrity), we evaluate a general pipeline comprising seven postprocessing blocks: noise correction; Gibbs ringing correction; evaluation of field distortions; susceptibility, eddy-current and motion-induced distortion corrections; bias field correction; spatial smoothing and final diffusion metric estimations. Based on this evaluation, we suggest an optimised processing pipeline for diffusion weighted MRI data.

Laplacian Eigenmaps Network-Based Nonlocal Means Method for MR Image Denoising

Magnetic resonance (MR) images are often corrupted by Rician noise which degrades the accuracy of image-based diagnosis tasks. The nonlocal means (NLM) method is a representative filter in denoising MR images due to its competitive denoising performance. However, the existing NLM methods usually exploit the gray-level information or hand-crafted features to evaluate the similarity between image patches, which is disadvantageous for preserving the image details while smoothing out noise. In this paper, an improved nonlocal means method is proposed for removing Rician noise in MR images by using the refined similarity measures. The proposed method firstly extracts the intrinsic features from the pre-denoised image using a shallow convolutional neural network named Laplacian eigenmaps network (LEPNet). Then, the extracted features are used for computing the similarity in the NLM method to produce the denoised image. Finally, the method noise of the denoised image is utilized to further improve the denoising performance. Specifically, the LEPNet model is composed of two cascaded convolutional layers and a nonlinear output layer, in which the Laplacian eigenmaps are employed to learn the filter bank in the convolutional layers and the Leaky Rectified Linear Unit activation function is used in the final output layer to output the nonlinear features. Due to the advantage of LEPNet in recovering the geometric structure of the manifold in the low-dimension space, the features extracted by this network can facilitate characterizing the self-similarity better than the existing NLM methods. Experiments have been performed on the BrainWeb phantom and the real images. Experimental results demonstrate that among several compared denoising methods, the proposed method can provide more effective noise removal and better details preservation in terms of human vision and such objective indexes as peak signal-to-noise ratio (PSNR) and structural similarity index measure (SSIM).

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry

PURPOSE

Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures.

METHODS

PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared.

RESULTS

Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of ≈3×, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution.

CONCLUSION

Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

Principal component analysis fosr fast and model-free denoising of multi b-value diffusion-weighted MR images

Despite the utility of tumour characterisation using quantitative parameter maps from multi-b-value diffusion-weighted MRI (DWI), clinicians often prefer the use of the image with highest diffusion-weighting (b-value), for instance for defining regions of interest (ROIs). However, these images are typically degraded by noise, as they do not utilize the information from the full acquisition. We present a principal component analysis (PCA) approach for model-free denoising of DWI data. PCA-denoising was compared to synthetic MRI, where a diffusion model is fitted for each voxel and a denoised image at a given b-value is generated from the model fit. A quantitative comparison of systematic and random errors was performed on data simulated using several diffusion models (mono-exponential, bi-exponential, stretched-exponential and kurtosis). A qualitative visual comparison was also performed for in vivo images in six healthy volunteers and three pancreatic cancer patients. In simulations, the reduction in random errors from PCA-denoising was substantial (up to 55%) and similar to synthetic MRI (up to 53%). Model-based synthetic MRI denoising resulted in substantial (up to 29% of signal) systematic errors, whereas PCA-denoising was able to denoise without introducing systematic errors (less than 2%). In vivo, the signal-to-noise ratio (SNR) and sharpness of PCA-denoised images were superior to synthetic MRI, resulting in clearer tumour boundaries. In the presence of motion, PCA-denoising did not cause image blurring, unlike image averaging or synthetic MRI. Multi-b-value MRI can be denoised model-free with our PCA-denoising strategy that reduces noise to a level similar to synthetic MRI, but without introducing systematic errors associated with the synthetic MRI method.

Registration-based image enhancement improves multi-atlas segmentation of the thalamic nuclei and hippocampal subfields

Magnetic resonance imaging (MRI) is an important tool for analysis of deep brain grey matter structures. However, analysis of these structures is limited due to low intensity contrast typically found in whole brain imaging protocols. Herein, we propose a big data registration-enhancement (BDRE) technique to augment the contrast of deep brain structures using an efficient large-scale non-rigid registration strategy. Direct validation is problematic given a lack of ground truth data. Rather, we validate the usefulness and impact of BDRE for multi-atlas (MA) segmentation on two sets of structures of clinical interest: the thalamic nuclei and hippocampal subfields. The experimental design compares algorithms using T1-weighted 3 T MRI for both structures (and additional 7 T MRI for the thalamic nuclei) with an algorithm using BDRE. As baseline comparisons, a recent denoising (DN) technique and a super-resolution (SR) method are used to preprocess the original 3 T MRI. The performance of each MA segmentation is evaluated by the Dice similarity coefficient (DSC). BDRE significantly improves mean segmentation accuracy over all methods tested for both thalamic nuclei (3 T imaging: 9.1%; 7 T imaging: 15.6%; DN: 6.9%; SR: 16.2%) and hippocampal subfields (3 T T1 only: 8.7%; DN: 8.4%; SR: 8.6%). We also present DSC performance for each thalamic nucleus and hippocampal subfield and show that BDRE can help MA segmentation for individual thalamic nuclei and hippocampal subfields. This work will enable large-scale analysis of clinically relevant deep brain structures from commonly acquired T1 images.

A post-processing method based on interphase motion correction and averaging to improve image quality of 4D magnetic resonance imaging: a clinical feasibility study

METHODS:

Nine patients (seven pancreas, one liver, and one lung) were recruited. 4D-MRI was performed using two prototype k-space sorted techniques, stack-of-stars (SOS) and koosh-ball (KB) acquisitions. Post-processing using MoCoAve was implemented for both methods. Image quality score, apparent SNR (aSNR), sharpness, motion trajectory and standard deviation (σ_GTV) of the gross tumor volumes were compared between original and MoCoAve image sets.

RESULTS:

All subjects successfully underwent 4D-MRI scans and MoCoAve was performed on all data sets. Significantly higher image quality scores (2.64 ± 0.39 vs 1.18 ± 0.34, p = 0.001) and aSNR (37.6 ± 15.3 vs 18.1 ± 5.7, p = 0.001) was observed in the MoCoAve images when compared to the original images. High correlation in tumor motion trajectories in the superoinferior direction (SI: 0.91 ± 0.08) and weaker in the anteroposterior (AP: 0.51 ± 0.44) and mediolateral (ML: 0.37 ± 0.23) directions, similar image sharpness (0.367 ± 0.068 vs 0.369 ± 0.072, p = 0.805), and minimal average absolute difference (0.47 ± 0.34  mm) of the motion trajectory profiles was found between the two image sets. The σ_GTV in pancreas patients was significantly (p = 0.039) lower in MoCoAve images (1.48 ± 1.35  cm3) than in the original images (2.17 ± 1.31  cm3).

CONCLUSION:

MoCoAve using interphase motion correction and averaging has shown promise as a post-processing method for improving k-space sorted (SOS and KB) 4D-MRI image quality in thoracic and abdominal cancer patients.

ADVANCES IN KNOWLEDGE:

The proposed method is an image based post-processing method that could be applied to many k-space sorted 4D-MRI methods for improved image quality and signal-to-noise ratio while preserving image sharpness and respiratory motion fidelity. It is a useful technique for the radiotherapy planning community who are interested in using 4D-MRI but aren't satisfied with their current MR image quality.