Mapping subcomponents of numerical cognition in relation to functional and anatomical landmarks of human parietal cortex

n/

An Empirical Study of the Maximum Degree of Undersampling in Compressed Sensing for T2*-weighted MRI

International audienceMagnetic Resonance Imaging (MRI) is one of the most dynamic and safe imaging modalities used in clinical routine today. Yet, one major limitation to this technique resides in its long acquisition times. Over the last decade, Compressed Sensing (CS) has been increasingly used to address this issue and offers to shorten MR scans by reconstructing images from undersampled Fourier data. Nevertheless, a quantitative guide on the degree of acceleration applicable to a given acquisition scenario is still lacking today, leading in practice to a trial-and-error approach in the selection of the appropriate undersampling factor. In this study, we shortly point out the existing theoretical sampling results in CS and their limitations which motivate the focus of this work: an empirical and quantitative analysis of the maximum degree of undersampling allowed by CS in the specific context of T2*-weighted MRI. We make use of a generic method based on retrospective undersampling to quantitatively deduce the maximum acceleration factor R max which preserves a desired image quality as a function of the image resolution and the available signal-to-noise ratio (SNR). Our results quantify how larger acceleration factors can be applied to higher resolution images as long as a minimum SNR is guaranteed. In practice however, the maximum acceleration factor for a given resolution appears to be constrained by the available SNR inherent to the considered acquisition. Our analysis enables to take this a priori knowledge into account, allowing to derive a sequence-specific maximum acceleration factor adapted to the intrinsic SNR of any MR pipeline. These results obtained on an analytical T2*-weighted phantom image were corroborated by prospective experiments performed on MR data collected with radial trajectories on a 7 Tesla scanner with the same contrast. The proposed framework allows to study other sequence weightings and therefore better optimize sequences when accelerated using CS

Optimization of sample preparation for MRI of formaldehyde-fixed brains

International audienceMagnetic resonance imaging of post-mortem brains allows long acquisition times up to several days and can be used to obtain high-resolution images at high field (7 T) which can be readily correlated with histological examination of the tissue. However, death and formaldehyde fixation are known to modify severely the relaxivity and diffusion properties of brain tissue. In particular, formaldehyde is known to shorten T2, which drastically reduces SNR.In order to counteract this effect and recover better SNR, free fixative can be washed out by soaking the sample in isotonic saline solution. This has been demonstrated in small biopsy-sized tissue samples, but little data is available concerning whole brain specimens.This study was designed to describe the kinetics of the change of relaxivity and diffusion properties of whole brain specimen at 7 T, during fixation, and during soaking in saline solution, in order to determine optimal soaking times.In the ewe brain, the fixation was found to stabilize after approximately 8 weeks, and the optimal duration of saline soaking is found to be around 3 weeks. These durations can be expected to be longer for larger specimen, such as human brains, which require longer penetration times

Is good old GRAPPA dead?

Presented at ISMRM 2021International audienceWe perform a qualitative analysis of performance of XPDNet, a state-of-the-art deep learning approach for MRI reconstruction, compared to GRAPPA, a classical approach. We do this in multiple settings, in particular testing the robustness of the XPDNet to unseen settings, and show that the XPDNet can to some degree generalize well

In Vivo High-Resolution 7 Tesla MRI Shows Early and Diffuse Cortical Alterations in CADASIL

Background and Purpose: Recent data suggest that early symptoms may be related to cortex alterations in CADASIL (Cerebral Autosomal-Dominant Arteriopathy with Subcortical Infarcts and Leukoencephalopathy), a monogenic model of cerebral small vessel disease (SVD). The aim of this study was to investigate cortical alterations using both high-resolution T2* acquisitions obtained with 7 Tesla MRI and structural T1 images with 3 Tesla MRI in CADASIL patients with no or only mild symptomatology (modified Rankin's scale = 24). Methods: Complete reconstructions of the cortex using 7 Tesla T2* acquisitions with 0.7 mm isotropic resolution were obtained in 11 patients (52.1 +/- 13.2 years, 36% male) and 24 controls (54.8 +/- 11.0 years, 42% male). Seven Tesla T2* within the cortex and cortical thickness and morphology obtained from 3 Tesla images were compared between CADASIL and control subjects using general linear models. Results: MMSE, brain volume, cortical thickness and global sulcal morphology did not differ between groups. By contrast, T2* measured by 7 Tesla MRI was significantly increased in frontal, parietal, occipital and cingulate cortices in patients after correction for multiple testing. These changes were not related to white matter lesions, lacunes or microhemorrhages in patients having no brain atrophy compared to controls. Conclusions: Seven Tesla MRI, by contrast to state of the art post-processing of 3 Tesla acquisitions, shows diffuse T2* alterations within the cortical mantle in CADASIL whose origin remains to be determined

Quantitative investigation of the effect of the extra-cerebral vasculature in diffuse optical imaging: a simulation study

Diffuse optical imaging (DOI) is a non invasive technique allowing the recovery of hemodynamic changes in the brain. Due to the diffusive nature of photon propagation in turbid media and the fact that cerebral tissues are located around 1.5 cm under the adult human scalp, DOI measurements are subject to partial volume errors. DOI measurements are also sensitive to large pial vessels because oxygenated and deoxygenated hemoglobin are the dominant chromophores in the near infrared window. In this study, the effect of the extra-cerebral vasculature in proximity of the sagittal sinus was investigated for its impact on DOI measurements simulated over the human adult visual cortex. Numerical Monte Carlo simulations were performed on two specific models of the human head derived from magnetic resonance imaging (MRI) scans. The first model included the extra-cerebral vasculature in which constant hemoglobin concentrations were assumed while the second did not. The screening effect of the vasculature was quantified by comparing recovered hemoglobin changes from each model for different optical arrays and regions of activation. A correction factor accounting for the difference between the recovered and the simulated hemoglobin changes was computed in each case. The results show that changes in hemoglobin concentration are better estimated when the extra-cerebral vasculature is modeled and the correction factors obtained in this case were at least 1.4-fold lower. The effect of the vasculature was also examined in a high-density diffuse optical tomography configuration. In this case, the difference between changes in hemoglobin concentration recovered with each model was reduced down to 10%

Simultaneous PET/MRI with 13C magnetic resonance spectroscopic imaging (hyperPET): phantom-based evaluation of PET quantification

BACKGROUND: Integrated PET/MRI with hyperpolarized (13)C magnetic resonance spectroscopic imaging ((13)C-MRSI) offers simultaneous, dual-modality metabolic imaging. A prerequisite for the use of simultaneous imaging is the absence of interference between the two modalities. This has been documented for a clinical whole-body system using simultaneous (1)H-MRI and PET but never for (13)C-MRSI and PET. Here, the feasibility of simultaneous PET and (13)C-MRSI as well as hyperpolarized (13)C-MRSI in an integrated whole-body PET/MRI hybrid scanner is evaluated using phantom experiments. METHODS: Combined PET and (13)C-MRSI phantoms including a NEMA [(18)F]-FDG phantom, (13)C-acetate and (13)C-urea sources, and hyperpolarized (13)C-pyruvate were imaged repeatedly with PET and/or (13)C-MRSI. Measurements evaluated for interference effects included PET activity values in the largest sphere and a background region; total number of PET trues; and (13)C-MRSI signal-to-noise ratio (SNR) for urea and acetate phantoms. Differences between measurement conditions were evaluated using t tests. RESULTS: PET and (13)C-MRSI data acquisition could be performed simultaneously without any discernible artifacts. The average difference in PET activity between acquisitions with and without simultaneous (13)C-MRSI was 0.83 (largest sphere) and −0.76 % (background). The average difference in net trues was −0.01 %. The average difference in (13)C-MRSI SNR between acquisitions with and without simultaneous PET ranged from −2.28 to 1.21 % for all phantoms and measurement conditions. No differences were significant. The system was capable of (13)C-MRSI of hyperpolarized (13)C-pyruvate. CONCLUSIONS: Simultaneous PET and (13)C-MRSI in an integrated whole-body PET/MRI hybrid scanner is feasible. Phantom experiments showed that possible interference effects introduced by acquiring data from the two modalities simultaneously are small and non-significant. Further experiments can now investigate the benefits of simultaneous PET and hyperpolarized (13)C-MRI in vivo studies

Iterative static ∆B0 field map estimation for off-resonance correction in non-Cartesian susceptibility weighted imaging

International audiencePatient-induced inhomogeneities in the magnetic field cause distortions and blurring during acquisitions with long readouts such as in susceptibility-weighted imaging (SWI). Most correction methods require collecting an additional ∆B 0 field map to remove these artifacts. The static ∆B0 field map can be approximated with an acceptable error directly from a single echo acquisition in SWI. The main component of the observed phase is linearly related to ∆B0 and the echo time TE , and the relative impact of non-∆B0 terms becomes insignificant with TE > 20ms at 3T for a well-tuned system. The main step is to combine and unfold the multi-channel phase maps wrapped many times, and several competing algorithms are compared for this purpose. Four in vivo brain data sets collected using the recently proposed 3D SPARKLING readouts are used to assess the proposed method. The estimated 3D field maps generated with a 0.6 mm isotropic spatial resolution provide overall similar off-resonance corrections compared to reference corrections based on an external ∆B0 acquisitions, and even improved for two out of four individuals. Although a small estimation error is expected, no aftermath was observed in the proposed corrections, while degradations were observed in the references. A static ∆B0 field map estimation method was proposed to take advantage of acquisitions with long echo times, and outperformed the reference technique based on an external field map. The difference can be attributed to an inherent robustness to mismatches between volumes and external ∆B0 maps, and diverse other sources investigated

Optimizing full 3D SPARKLING trajectories for high-resolution T2*-weighted Magnetic Resonance Imaging

International audienceThe Spreading Projection Algorithm for Rapid K-space samplING, or SPARKLING, is an optimization-driven method that has been recently introduced for accelerated 2D T2*-w MRI using compressed sensing. It has then been extended to address 3D imaging using either stacks of 2D sampling patterns or a local 3D strategy that optimizes a single sampling trajectory at a time. 2D SPARKLING actually performs variable density sampling (VDS) along a prescribed target density while maximizing sampling efficiency and meeting the gradient-based hardware constraints. However, 3D SPARKLING has remained limited in terms of acceleration factors along the third dimension if one wants to preserve a peaky point spread function (PSF) and thus good image quality.In this paper, in order to achieve higher acceleration factors in 3D imaging while preserving image quality, we propose a new efficient algorithm that performs optimization on full 3D SPARKLING. The proposed implementation based on fast multipole methods (FMM) allows us to design sampling patterns with up to 10^7 k-space samples, thus opening the door to 3D VDS. We compare multi-CPU and GPU implementations and demonstrate that the latter is optimal for 3D imaging in the high-resolution acquisition regime (600µm isotropic). Finally, we show that this novel optimization for full 3D SPARKLING outperforms stacking strategies or 3D twisted projection imaging through retrospective and prospective studies on NIST phantom and in vivo brain scans at 3 Tesla. Overall the proposed method allows for 2.5-3.75x shorter scan times compared to GRAPPA-4 parallel imaging acquisition at 3 Tesla without compromising image quality

Correction d'inhomogénéités de champs pour la SWI non-cartésienne par estimation des cartes de champs

International audiencePatient-induced inhomogeneities in the magnetic field cause distortions and blurring during acquisitions with long echo times, as in susceptibility-weighted imaging. Most correction methods require collecting an additional ΔB0 field map. To avoid that, we propose a method to approximate this field map using the single echo acquisition only. The main component of the observed phase is linearly related to ΔB0 and TE, and the relative impact of non-ΔB0 terms becomes insignificant with TE>20ms at 3T. The estimated 3D field maps, produced at 0.6 mm isotropic under 3 minutes, provide equivalent corrections to acquired ones.Les inhomogénéités de champs induites par les patients sont à l'origine de distorsions et de floutages durant les acquisitions à temps d'écho longs, comme pour l'imagerie pondérée en susceptibilité. La plupart des méthodes de correction nécessitent d'acquérir une carte de champ ΔB0 additionnelle. Pour éviter cela, nous proposons une méthode pour approximer cette carte de champs en utilisant seulement l'acquisition à écho unique. La composante principale de la phase observée est linéairement liée au ΔB0 et au TE, et l'impact relatif des termes indépendants du ΔB0 deviennent négligeables pour TE>20ms à 3T. Les cartes 3D estimées, produites à 0.6 mm isotrope en moins de 3 minutes, permettent d'obtenir une correction équivalente aux cartes acquises