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

Higher spatial resolution and sensitivity in whole brain functional MRI at 7T using 3D EPI accelerated with variable density CAIPI sampling and temporal random walk

PURPOSE

To develop an efficient 3D EPI encoding technique for high spatiotemporal resolution functional MRI.

METHODS

To exploit spatiotemporal fMRI data structure, we introduce a variable density 2D CAIPI sampling in the spatial domain combined with time-wise extra random encoding in the time domain, thus achieving pseudo-regular sampling with a regular blip while allowing incoherent sampling in a complementary manner across time. This enabled temporally regularized reconstruction of highly accelerated functional data acquisition. The encoding scheme was then validated against temporally invariant CAIPI encoding by applying to locally confined and whole-brain around the primary visual cortex, respectively, with increasing the spatial resolutions.

RESULTS

For partial brain imaging, our proposed method achieved higher reconstruction accuracy, resulting in a substantial increase of SSIM compared to an alternative method for 0.64 mm-isotropic resolution. When used for whole brain imaging at 0.56 mm-isotropic resolution, our method showed a decreased spatial extent of activation and produced high-quality images for a clear distinction between activated and non-activated regions around calcarine fissure with high spatial specificity.

CONCLUSION

The proposed 3D EPI encoding scheme, which exploits coherent and incoherent sampling properties, can significantly improve the image quality while providing a good balance between sensitivity and specificity in the activated regions.

3D Wasserstein Generative Adversarial Network with Dense U-Net-Based Discriminator for Preclinical fMRI Denoising

Functional magnetic resonance imaging (fMRI) is extensively used in clinical and preclinical settings to study brain function; however, fMRI data is inherently noisy due to physiological processes, hardware, and external noise. Denoising is one of the main preprocessing steps in any fMRI analysis pipeline. This process is challenging in preclinical data in comparison to clinical data due to variations in brain geometry, image resolution, and low signal-to-noise ratios. In this paper, we propose a structure-preserved algorithm based on a 3D Wasserstein generative adversarial network with a 3D dense U-net-based discriminator called 3D U-WGAN. We apply a 4D data configuration to effectively denoise temporal and spatial information in analyzing preclinical fMRI data. GAN-based denoising methods often utilize a discriminator to identify significant differences between denoised and noise-free images, focusing on global or local features. To refine the fMRI denoising model, our method employs a 3D dense U-Net discriminator to learn both global and local distinctions. To tackle potential oversmoothing, we introduce an adversarial loss and enhance perceptual similarity by measuring feature space distances. Experiments illustrate that 3D U-WGAN significantly improves image quality in resting-state and task preclinical fMRI data, enhancing signal-to-noise ratio without introducing excessive structural changes in existing methods. The proposed method outperforms state-of-the-art methods when applied to simulated and real data in a fMRI analysis pipeline.

NORDIC denoising on VASO data

The use of submillimeter resolution functional magnetic resonance imaging (fMRI) is increasing in popularity due to the prospect of studying human brain activation non-invasively at the scale of cortical layers and columns. This method, known as laminar fMRI, is inherently signal-to-noise ratio (SNR)-limited, especially at lower field strengths, with the dominant noise source being of thermal origin. Furthermore, laminar fMRI is challenged with signal displacements due to draining vein effects in conventional gradient-echo blood oxygen level-dependent (BOLD) imaging contrasts. fMRI contrasts such as cerebral blood volume (CBV)-sensitive vascular space occupancy (VASO) sequences have the potential to mitigate draining vein effects. However, VASO comes along with another reduction in detection sensitivity. NOise Reduction with DIstribution Corrected (NORDIC) PCA (principal component analysis) is a denoising technique specifically aimed at suppressing thermal noise, which has proven useful for increasing the SNR of high-resolution functional data. While NORDIC has been examined for BOLD acquisitions, its application to VASO data has been limited, which was the focus of the present study. We present a preliminary analysis to evaluate NORDIC's capability to suppress thermal noise while preserving the VASO signal across a wide parameter space at 3T. For the data presented here, with a proper set of parameters, NORDIC reduced thermal noise with minimal bias on the underlying signal and preserved spatial resolution. Denoising performance was found to vary with different implementation strategies and parameter choices, for which we provide recommendations. We conclude that when applied properly, NORDIC has the potential to overcome the sensitivity limitations of laminar-specific VASO fMRI. Since very few groups currently have 3T VASO data, by sharing our analysis and code, we can compile and compare the effects of NORDIC across a broader range of acquisition parameters and study designs. Such a communal effort will help develop robust recommendations that will increase the utility of laminar fMRI at lower field strengths.

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.

NOise Reduction with DIstribution Corrected (NORDIC) principal component analysis improves brain activity detection across rodent and human functional MRI contexts

NOise Reduction with DIstribution Corrected (NORDIC) principal component analysis (PCA) has been shown to selectively suppress thermal noise and improve the temporal signal-to-noise ratio (tSNR) in human functional magnetic resonance imaging (fMRI). However, the feasibility to improve data quality for rodent fMRI using NORDIC PCA remains uncertain. NORDIC PCA may also be particularly beneficial for improving topological brain mapping, as conventional mapping requires precise spatiotemporal signals from large datasets (ideally ~1 hour acquisition) for individual representations. In this study, we evaluated the effects of NORDIC PCA compared with "Standard" processing in various rodent fMRI contexts that range from task-evoked optogenetic fMRI to resting-state fMRI. We also evaluated the effects of NORDIC PCA on human resting-state and retinotopic mapping fMRI via population receptive field (pRF) modeling. In rodent optogenetic fMRI, apart from doubling the tSNR, NORDIC PCA resulted in a larger number of activated voxels and a significant decrease in the variance of evoked brain responses without altering brain morphology. In rodent resting-state fMRI, we found that NORDIC PCA induced a nearly threefold increase in tSNR and preserved task-free relative cerebrovascular reactivity (rCVR) across cortical depth. NORDIC PCA further improved the detection of TGN020-induced aquaporin-4 inhibition on rCVR compared with Standard processing without NORDIC PCA. NORDIC PCA also increased the tSNR for both human resting-state and pRF fMRI, and for the latter also increased activation cluster sizes while retaining retinotopic organization. This suggests that NORDIC PCA preserves the spatiotemporal precision of fMRI signals needed for pRF analysis, and effectively captures small activity changes with high sensitivity. Taken together, these results broadly demonstrate the value of NORDIC PCA for the enhanced detection of neural dynamics across various rodent and human fMRI contexts. This can in turn play an important role in improving fMRI image quality and sensitivity for translational and preclinical neuroimaging research.

Enhancing clinical diagnostics: novel denoising methodology for brain MRI with adaptive masking and modified non-local block

Medical image denoising has been a subject of extensive research, with various techniques employed to enhance image quality and facilitate more accurate diagnostics. The evolution of denoising methods has highlighted impressive results but struggled to strike equilibrium between noise reduction and edge preservation which limits its applicability in various domains. This paper manifests the novel methodology that integrates an adaptive masking strategy, transformer-based U-Net Prior generator, edge enhancement module, and modified non-local block (MNLB) for denoising brain MRI clinical images. The adaptive masking strategy maintains the vital information through dynamic mask generation while the prior generator by capturing hierarchical features regenerates the high-quality prior MRI images. Finally, these images are fed to the edge enhancement module to boost structural information by maintaining crucial edge details, and the MNLB produces the denoised output by deriving non-local contextual information. The comprehensive experimental assessment is performed by employing two datasets namely the brain tumor MRI dataset and Alzheimer's dataset for diverse metrics and compared with conventional denoising approaches. The proposed denoising methodology achieves a PSNR of 40.965 and SSIM of 0.938 on the Alzheimer's dataset and also achieves a PSNR of 40.002 and SSIM of 0.926 on the brain tumor MRI dataset at a noise level of 50% revealing its supremacy in noise minimization. Furthermore, the impact of different masking ratios on denoising performance is analyzed which reveals that the proposed method showed PSNR of 40.965, SSIM of 0.938, MAE of 5.847, and MSE of 3.672 at the masking ratio of 60%. Moreover, the findings pave the way for the advancement of clinical image processing, facilitating precise detection of tumors in clinical MRI images.

High-resolution prostate diffusion MRI using eddy current-nulled convex optimized diffusion encoding and random matrix theory-based denoising

OBJECTIVE

To develop and evaluate a technique combining eddy current-nulled convex optimized diffusion encoding (ENCODE) with random matrix theory (RMT)-based denoising to accelerate and improve the apparent signal-to-noise ratio (aSNR) and apparent diffusion coefficient (ADC) mapping in high-resolution prostate diffusion-weighted MRI (DWI). MATERIALS AND METHODS: Eleven subjects with clinical suspicion of prostate cancer were scanned at 3T with high-resolution (HR) (in-plane: 1.0 × 1.0 mm2) ENCODE and standard-resolution (1.6 × 2.2 mm2) bipolar DWI sequences (both had 7 repetitions for averaging, acquisition time [TA] of 5 min 50 s). HR-ENCODE was retrospectively analyzed using three repetitions (accelerated effective TA of 2 min 30 s). The RMT-based denoising pipeline utilized complex DWI signals and Marchenko-Pastur distribution-based principal component analysis to remove additive Gaussian noise in images from multiple coils, b-values, diffusion encoding directions, and repetitions. HR-ENCODE with RMT-based denoising (HR-ENCODE-RMT) was compared with HR-ENCODE in terms of aSNR in prostate peripheral zone (PZ) and transition zone (TZ). Precision and accuracy of ADC were evaluated by the coefficient of variation (CoV) between repeated measurements and mean difference (MD) compared to the bipolar ADC reference, respectively. Differences were compared using two-sided Wilcoxon signed-rank tests (P < 0.05 considered significant).

RESULTS

HR-ENCODE-RMT yielded 62% and 56% higher median aSNR than HR-ENCODE (b = 800 s/mm2) in PZ and TZ, respectively (P < 0.001). HR-ENCODE-RMT achieved 63% and 70% lower ADC-CoV than HR-ENCODE in PZ and TZ, respectively (P < 0.001). HR-ENCODE-RMT ADC and bipolar ADC had low MD of 22.7 × 10-6 mm2/s in PZ and low MD of 90.5 × 10-6 mm2/s in TZ.

CONCLUSIONS

HR-ENCODE-RMT can shorten the acquisition time and improve the aSNR of high-resolution prostate DWI and achieve accurate and precise ADC measurements in the prostate.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise-the dominant contributing noise component in high-resolution fMRI. NOise Reduction with DIstribution Corrected Principal Component Analysis (NORDIC PCA) is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. While investigating auditory functional responses poses unique challenges, we anticipated NORDIC to have a positive impact on the data on the basis of previous applications. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we did observe a reduction in the average response amplitude (percent signal change) within regions of interest, which may suggest that a portion of the signal of interest, which could not be distinguished from general i.i.d. noise, was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high-resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

Tensor denoising of quantitative multi-parameter mapping

PURPOSE

Quantitative multi-parameter mapping (MPM) provides maps of physical quantities representing physiologically meaningful tissue characteristics, which allows to investigate microstructure-function relationships reflecting normal or pathologic processes in the brain. However, the achievable spatial resolution and stability of MPM for basic research or clinical applications is severely constrained by SNR limits of the MR acquisition process, resulting in relatively long acquisition times. To increase SNR, we denoise MPM acquisitions using principal component analysis along tensors exploiting the Marchenko-Pastur law (tMPPCA).

METHODS

tMPPCA denoising was applied to three sets of MPM raw data before the quantification of maps of proton density, magnetization transfer saturation, R1, and R2*. The regional SNR gain for high-resolution MPM was investigated as well as reproducibility gains for clinically optimized protocols with moderate and high acceleration factors at different image resolutions.

RESULTS

Substantial noise reduction in raw data was achieved, resulting in reduced noise for quantitative mapping up to sixfold without introducing bias of mean values (below 1%). Scan-rescan fluctuations were reduced up to threefold. Denoising allowed to decrease the voxel volume fourfold at the same scan time or reduce the scan time twofold at same voxel volume without loss of sensitivity.

CONCLUSIONS

tMPPCA denoising can (a) improve of fine spatial and temporal patterns, (b) considerably reduce scan time for clinical applications, or (c) increase resolution to potentially push cutting-edge MPM protocols from the upper to the lower limit of the mesoscopic scale.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise - the dominant contributing noise component in high resolution fMRI. NORDIC PCA is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. As investigating auditory functional responses poses unique challenges, we anticipated that the benefit of this technique would be especially pronounced. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we also observed a reduction in the average response amplitude (percent signal), which may suggest that a small amount of signal was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.