Denoising High-Field Multi-Dimensional MRI With Local Complex PCA

Locally low-rank denoising in transform domains

Purpose

To develop an extension to locally low rank (LLR) denoising techniques based on transform domain processing that reduces the number of images required in the MR image series for high-quality denoising.

Theory and Methods

LLR methods with random matrix theory-based thresholds are successfully used in the denoising of MR image series in a number of applications. The performance of these methods depend on how well the LLR assumption is satisfied, which deteriorates with few numbers of images, as is commonly encountered in quantitative MRI applications. We propose a transform-domain approach for denoising of MR image series to represent the underlying signal with higher fidelity when using a locally low rank approximation. The efficacy of the method is demonstrated for fully-sampled k-space, undersampled k-space, DICOM images, and complex-valued SENSE-1 images in quantitative MRI applications with as few as 4 images.

Results

For both MSK and brain applications, the transform domain denoising preserves local subtle variability, whereas the quantitative maps based on image domain LLR methods tend to be locally more homogeneous.

Conclusion

A transform domain extension to LLR denoising produces high quality images and is compatible with both raw k-space data and vendor reconstructed data. This allows for improved imaging and more accurate quantitative analyses and parameters obtained therefrom.

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.

MRI Evaluation of Indomethacin Suppositories in the Prevention of Complications of Pancreatitis and Hyperamylasemia after Choledocholithiasis ERCP Based on Image Denoising Algorithm

Objective

To explore the value of MRI evaluation of indomethacin suppositories in the prevention of pancreatitis and hyperamylasemia in patients with common bile duct calculi after endoscopic retrograde cholangiopancreatography (ERCP) based on image denoising algorithm.

Methods

A retrospective analysis in August 2020 to December 2021. Because of the common bile duct calculi hospitalized parallel ERCP operation, 89 cases of patients, according to the different postoperative treatments, were divided into group A (n = 44) and group B (n = 45), in which A set of separate application inhibits the pancreatic enzyme secretion after surgery drug treatment, and B group on the basis of group A linked with indole beauty Xinshuan treatment. The incidence of postoperative pancreatitis and hyperamylasemia was compared between the two groups. The levels of serum amylase were compared between the two groups. Patients in group B were diagnosed with pancreatitis by conventional MRI and MRI with denoising algorithm, respectively, and the imaging characteristics and diagnosis rate differences of the two methods were observed. ROC curve was drawn to evaluate the diagnostic efficacy of MRI denoising algorithm for postoperative pancreatitis and serum amylase level detection for hyperamylasemia.

Results

The incidence of postoperative pancreatitis and hyperamylasemia in group B was significantly lower than that in group A (P < 0.05). There were 6 cases of postoperative pancreatitis in group B, 2 cases (33.33%) were diagnosed by conventional MRI, and 5 cases (83.33%) were diagnosed by MRI based on denoising algorithm. Although there was no significant difference in diagnosis rate between the two methods, the number of cases of pancreatitis diagnosed by MRI based on denoising algorithm was slightly higher than that by conventional MRI. Compared with conventional MRI images, MRI images with denoising algorithm showed that the number of cases with pancreatic swelling, the number of cases with pancreatic duct/bile duct dilation, and the number of cases with abdominal effusion were all high (all P < 0.05). ROC results showed that the area under the curve of MRI with denoising algorithm for the diagnosis of postoperative pancreatitis was 0.855, and the sensitivity was 89.40%. The specificity was 83.20%, and the area under the curve of serum amylase for postoperative hyperamylasemia was 0.893, the sensitivity was 89.80%, and the specificity was 85.20%, all of which had high diagnostic efficacy.

Conclusion

MRI results of denoising algorithm suggest that indomethacin suppositories can effectively reduce the incidence of postoperative pancreatitis and hyperamylasemia after ERCP, which is worthy of clinical application.

SDnDTI: Self-supervised deep learning-based denoising for diffusion tensor MRI

Diffusion tensor magnetic resonance imaging (DTI) is a widely adopted neuroimaging method for the in vivo mapping of brain tissue microstructure and white matter tracts. Nonetheless, the noise in the diffusion-weighted images (DWIs) decreases the accuracy and precision of DTI derived microstructural parameters and leads to prolonged acquisition time for achieving improved signal-to-noise ratio (SNR). Deep learning-based image denoising using convolutional neural networks (CNNs) has superior performance but often requires additional high-SNR data for supervising the training of CNNs, which reduces the feasibility of supervised learning-based denoising in practice. In this work, we develop a self-supervised deep learning-based method entitled "SDnDTI" for denoising DTI data, which does not require additional high-SNR data for training. Specifically, SDnDTI divides multi-directional DTI data into many subsets of six DWI volumes and transforms DWIs from each subset to along the same diffusion-encoding directions through the diffusion tensor model, generating multiple repetitions of DWIs with identical image contrasts but different noise observations. SDnDTI removes noise by first denoising each repetition of DWIs using a deep 3-dimensional CNN with the average of all repetitions with higher SNR as the training target, following the same approach as normal supervised learning based denoising methods, and then averaging CNN-denoised images for achieving higher SNR. The denoising efficacy of SDnDTI is demonstrated in terms of the similarity of output images and resultant DTI metrics compared to the ground truth generated using substantially more DWI volumes on two datasets with different spatial resolutions, b-values and numbers of input DWI volumes provided by the Human Connectome Project (HCP) and the Lifespan HCP in Aging. The SDnDTI results preserve image sharpness and textural details and substantially improve upon those from the raw data. The results of SDnDTI are comparable to those from supervised learning-based denoising and outperform those from state-of-the-art conventional denoising algorithms including BM4D, AONLM and MPPCA. By leveraging domain knowledge of diffusion MRI physics, SDnDTI makes it easier to use CNN-based denoising methods in practice and has the potential to benefit a wider range of research and clinical applications that require accelerated DTI acquisition and high-quality DTI data for mapping of tissue microstructure, fiber tracts and structural connectivity in the living human brain.

Functional Coupling of the Locus Coeruleus Is Linked to Successful Cognitive Control

The locus coeruleus (LC) is a brainstem structure that sends widespread efferent projections throughout the mammalian brain. The LC constitutes the major source of noradrenaline (NE), a modulatory neurotransmitter that is crucial for fundamental brain functions such as arousal, attention, and cognitive control. This role of the LC-NE is traditionally not believed to reflect functional influences on the frontoparietal network or the striatum, but recent advances in chemogenetic manipulations of the rodent brain have challenged this notion. However, demonstrations of LC-NE functional connectivity with these areas in the human brain are surprisingly sparse. Here, we close this gap. Using an established emotional stroop task, we directly compared trials requiring response conflict control with trials that did not require this, but were matched for visual stimulus properties, response modality, and controlled for pupil dilation differences across both trial types. We found that LC-NE functional coupling with the parietal cortex and regions of the striatum is substantially enhanced during trials requiring response conflict control. Crucially, the strength of this functional coupling was directly related to individual reaction time differences incurred by conflict resolution. Our data concur with recent rodent findings and highlight the importance of converging evidence between human and nonhuman neurophysiology to further understand the neural systems supporting adaptive and maladaptive behavior in health and disease.

A local multi-transmit coil combined with a high-density receive array for cerebellar fMRI at 7 T

The human cerebellum is involved in a wide array of functions, ranging from motor control to cognitive control, and as such is of great neuroscientific interest. However, its function is underexplored in vivo, due to its small size, its dense structure and its placement at the bottom of the brain, where transmit and receive fields are suboptimal. In this study, we combined two dense coil arrays of 16 small surface receive elements each with a transmit array of three antenna elements to improve BOLD sensitivity in the human cerebellum at 7 T. Our results showed improved B₁⁺ and SNR close to the surface as well as g-factor gains compared with a commercial coil designed for whole-head imaging. This resulted in improved signal stability and large gains in the spatial extent of the activation close to the surface (<3.5 cm), while good performance was retained deeper in the cerebellum. Modulating the phase of the transmit elements of the head coil to constructively interfere in the cerebellum improved the B₁⁺ , resulting in a temporal SNR gain. Overall, our results show that a dedicated transmit array along with the SNR gains of surface coil arrays can improve cerebellar imaging, at the cost of a decreased field of view and increased signal inhomogeneity.

Real-world stress resilience is associated with the responsivity of the locus coeruleus

Individuals may show different responses to stressful events. Here, we investigate the neurobiological basis of stress resilience, by showing that neural responsitivity of the noradrenergic locus coeruleus (LC-NE) and associated pupil responses are related to the subsequent change in measures of anxiety and depression in response to prolonged real-life stress. We acquired fMRI and pupillometry data during an emotional-conflict task in medical residents before they underwent stressful emergency-room internships known to be a risk factor for anxiety and depression. The LC-NE conflict response and its functional coupling with the amygdala was associated with stress-related symptom changes in response to the internship. A similar relationship was found for pupil-dilation, a potential marker of LC-NE firing. Our results provide insights into the noradrenergic basis of conflict generation, adaptation and stress resilience.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI

An efficient multi-slice inversion–recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T₁ mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T₁ maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T₁ measurements are assessed using both T₁-calibrated and T₂-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T₁ values. We model the effect of these fat suppression pulses on T₁ measurements and show that the model can remove their MT contribution from the measured T₁, thus providing accurate T₁ quantification. High spatial resolution T₁ maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T₁ maps than MP2RAGE, demonstrating the potential for ultra-high resolution T₁ mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

Improved cortical surface reconstruction using sub-millimeter resolution MPRAGE by image denoising

Automatic cerebral cortical surface reconstruction is a useful tool for cortical anatomy quantification, analysis and visualization. Recently, the Human Connectome Project and several studies have shown the advantages of using T₁-weighted magnetic resonance (MR) images with sub-millimeter isotropic spatial resolution instead of the standard 1-mm isotropic resolution for improved accuracy of cortical surface positioning and thickness estimation. Nonetheless, sub-millimeter resolution images are noisy by nature and require averaging multiple repetitions to increase the signal-to-noise ratio for precisely delineating the cortical boundary. The prolonged acquisition time and potential motion artifacts pose significant barriers to the wide adoption of cortical surface reconstruction at sub-millimeter resolution for a broad range of neuroscientific and clinical applications. We address this challenge by evaluating the cortical surface reconstruction resulting from denoised single-repetition sub-millimeter T₁-weighted images. We systematically characterized the effects of image denoising on empirical data acquired at 0.6 mm isotropic resolution using three classical denoising methods, including denoising convolutional neural network (DnCNN), block-matching and 4-dimensional filtering (BM4D) and adaptive optimized non-local means (AONLM). The denoised single-repetition images were found to be highly similar to 6-repetition averaged images, with a low whole-brain averaged mean absolute difference of ~0.016, high whole-brain averaged peak signal-to-noise ratio of ~33.5 dB and structural similarity index of ~0.92, and minimal gray matter–white matter contrast loss (2% to 9%). The whole-brain mean absolute discrepancies in gray matter–white matter surface placement, gray matter–cerebrospinal fluid surface placement and cortical thickness estimation were lower than 165 μm, 155 μm and 145 μm—sufficiently accurate for most applications. These discrepancies were approximately one third to half of those from 1-mm isotropic resolution data. The denoising performance was equivalent to averaging ~2.5 repetitions of the data in terms of image similarity, and 1.6–2.2 repetitions in terms of the cortical surface placement accuracy. The scan-rescan variability of the cortical surface positioning and thickness estimation was lower than 170 μm. Our unique dataset and systematic characterization support the use of denoising methods for improved cortical surface reconstruction at sub-millimeter resolution.

Multi-contrast anatomical subcortical structures parcellation

The human subcortex is comprised of more than 450 individual nuclei which lie deep in the brain. Due to their small size and close proximity, up until now only 7% have been depicted in standard MRI atlases. Thus, the human subcortex can largely be considered as terra incognita. Here, we present a new open-source parcellation algorithm to automatically map the subcortex. The new algorithm has been tested on 17 prominent subcortical structures based on a large quantitative MRI dataset at 7 Tesla. It has been carefully validated against expert human raters and previous methods, and can easily be extended to other subcortical structures and applied to any quantitative MRI dataset. In sum, we hope this novel parcellation algorithm will facilitate functional and structural neuroimaging research into small subcortical nuclei and help to chart terra incognita.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Transferring principles of solid-state and Laplace NMR to the field of in vivo brain MRI

Magnetic resonance imaging (MRI) is the primary method for noninvasive investigations of the human brain in health, disease, and development but yields data that are difficult to interpret whenever the millimeter-scale voxels contain multiple microscopic tissue environments with different chemical and structural properties. We propose a novel MRI framework to quantify the microscopic heterogeneity of the living human brain as spatially resolved five-dimensional relaxation-diffusion distributions by augmenting a conventional diffusion-weighted imaging sequence with signal encoding principles from multidimensional solid-state nuclear magnetic resonance (NMR) spectroscopy, relaxation-diffusion correlation methods from Laplace NMR of porous media, and Monte Carlo data inversion. The high dimensionality of the distribution space allows resolution of multiple microscopic environments within each heterogeneous voxel as well as their individual characterization with novel statistical measures that combine the chemical sensitivity of the relaxation rates with the link between microstructure and the anisotropic diffusivity of tissue water. The proposed framework is demonstrated on a healthy volunteer using both exhaustive and clinically viable acquisition protocols.