Quantitative platform for accurate and reproducible assessment of transverse (T2 ) relaxation time

Cartilage compositional MRI-a narrative review of technical development and clinical applications over the past three decades

Articular cartilage damage and degeneration are among hallmark manifestations of joint injuries and arthritis, classically osteoarthritis. Cartilage compositional MRI (Cart-C MRI), a quantitative technique, which aims to detect early-stage cartilage matrix changes that precede macroscopic alterations, began development in the 1990s. However, despite the significant advancements over the past three decades, Cart-C MRI remains predominantly a research tool, hindered by various technical and clinical hurdles. This paper will review the technical evolution of Cart-C MRI, delve into its clinical applications, and conclude by identifying the existing gaps and challenges that need to be addressed to enable even broader clinical application of Cart-C MRI.

Efficient standardization of clinical T2-weighted images: Phase-conjugacy e-CAMP with projected gradient descent

PURPOSE

To standardize T  2 $$ {}_2 $$ -weighted images from clinical Turbo Spin Echo (TSE) scans by generating corresponding T  2 $$ {}_2 $$ maps with the goal of removing scanner- and/or protocol-specific heterogeneity.

METHODS

The T  2 $$ {}_2 $$ map is estimated by minimizing an objective function containing a data fidelity term in a Virtual Conjugate Coils (VCC) framework, where the signal evolution model is expressed as a linear constraint. The objective function is minimized by Projected Gradient Descent (PGD).

RESULTS

The algorithm achieves accuracy comparable to methods with customized sampling schemes for accelerated T  2 $$ {}_2 $$ mapping. The results are insensitive to the tunable parameters, and the relaxed background phase prior produces better T  2 $$ {}_2 $$ maps compared to the strict real-value enforcement. It is worth noting that the algorithm works well with challenging T  2 $$ {}_2 $$ w-TSE data using typical clinical parameters. The observed normalized root mean square error ranges from 6.8% to 12.3% over grey and white matter, a clinically common level of quantitative map error.

CONCLUSION

The novel methodological development creates an efficient algorithm that allows for T  2 $$ {}_2 $$ map generated from TSE data with typical clinical parameters, such as high resolution, long echo train length, and low echo spacing. Reconstruction of T  2 $$ {}_2 $$ maps from TSE data with typical clinical parameters has not been previously reported.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

Improvements in Between-Vendor MRI Harmonization of Renal T2 Mapping using Stimulated Echo Compensation

BACKGROUND

T2 mapping is valuable to evaluate pathophysiology in kidney disease. However, variations in T2 relaxation time measurements across MR scanners and vendors may occur requiring additional correction.

PURPOSE

To harmonize renal T2 measurements between MR vendor platforms, and use an extended-phase-graph-based fitting method ("StimFit") to correct stimulated echoes and reduce between-vendor variations.

STUDY TYPE

Prospective.

SUBJECTS

8 healthy "travelling" volunteers (37.5% female, 32 ± 6 years) imaged on four MRI systems across three vendors at four sites, 10 healthy volunteers (50% female, 32 ± 8 years) scanned multiple times on a given MR scanner for repeatability evaluation. ISMRM/NIST system phantom scanned for evaluation of T2 accuracy.

FIELD STRENGTH/SEQUENCE

3T, multiecho spin-echo sequence.

ASSESSMENT

T2 images fit using conventional monoexponential fitting and "StimFit." Mean absolute percentage error (MAPE) of phantom measurements with reference T2 values. Average cortex and medulla T2 values compared between MR vendors, with masks obtained from T2 -weighted images and T1 maps. Full-width-at-half-maximum (FWHM) T2 distributions to evaluate local homogeneity of measurements.

STATISTICAL TESTS

Coefficient of variation (CV), linear mixed-effects model, analysis of variance, student's t-tests, Bland-Altman plots, P-value <0.05 considered statistically significant.

RESULTS

In the ISMRM/NIST phantom, "StimFit" reduced the MAPE from 4.9%, 9.1%, 24.4%, and 18.1% for the four sites (three vendors) to 3.3%, 3.0%, 6.6%, and 4.1%, respectively. In vivo, there was a significant difference in kidney T2 measurements between vendors using a monoexponential fit, but not with "StimFit" (P = 0.86 and 0.92, cortex and medulla, respectively). The intervendor CVs of T2 measures were reduced from 8.0% to 2.6% (cortex) and 7.1% to 2.8% (medulla) with StimFit, resulting in no significant differences for the CVs of intravendor repeat acquisitions (P = 0.13 and 0.05). "StimFit" significantly reduced the FWHM of T2 distributions in the cortex and whole kidney.

DATA CONCLUSION

Stimulated-echo correction reduces renal T2 variation across MR vendor platforms.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Probing muscle recovery following downhill running using precise mapping of MRI T2 relaxation times

PURPOSE

Postexercise recovery rate is a vital component of designing personalized training protocols and rehabilitation plans. Tracking exercise-induced muscle damage and recovery requires sensitive tools that can probe the muscles' state and composition noninvasively.

METHODS

Twenty-four physically active males completed a running protocol consisting of a 60-min downhill run on a treadmill at -10% incline and 65% of maximal heart rate. Quantitative mapping of MRI T2 was performed using the echo-modulation-curve algorithm before exercise, and at two time points: 1 h and 48 h after exercise.

RESULTS

T2 values increased by 2%-4% following exercise in the primary mover muscles and exhibited further elevation of 1% after 48 h. For the antagonist muscles, T2 values increased only at the 48-h time point (2%-3%). Statistically significant decrease in the SD of T2 values was found following exercise for all tested muscles after 1 h (16%-21%), indicating a short-term decrease in the heterogeneity of the muscle tissue.

CONCLUSION

MRI T2 relaxation time constitutes a useful quantitative marker for microstructural muscle damage, enabling region-specific identification for short-term and long-term systemic processes, and sensitive assessment of muscle recovery following exercise-induced muscle damage. The variability in T2 changes across different muscle groups can be attributed to their different role during downhill running, with immediate T2 elevation occurring in primary movers, followed by delayed elevation in both primary and antagonist muscle groups, presumably due to secondary damage caused by systemic processes.

P2T2: A physically-primed deep-neural-network approach for robust T2 distribution estimation from quantitative T2-weighted MRI

Estimating T₂ relaxation time distributions from multi-echo T₂-weighted MRI (T₂W) data can provide valuable biomarkers for assessing inflammation, demyelination, edema, and cartilage composition in various pathologies, including neurodegenerative disorders, osteoarthritis, and tumors. Deep neural network (DNN) based methods have been proposed to address the complex inverse problem of estimating T₂ distributions from MRI data, but they are not yet robust enough for clinical data with low Signal-to-Noise ratio (SNR) and are highly sensitive to distribution shifts such as variations in echo-times (TE) used during acquisition. Consequently, their application is hindered in clinical practice and large-scale multi-institutional trials with heterogeneous acquisition protocols. We propose a physically-primed DNN approach, called P₂T₂, that incorporates the signal decay forward model in addition to the MRI signal into the DNN architecture to improve the accuracy and robustness of T₂ distribution estimation. We evaluated our P₂T₂ model in comparison to both DNN-based methods and classical methods for T₂ distribution estimation using 1D and 2D numerical simulations along with clinical data. Our model improved the baseline model's accuracy for low SNR levels (SNR<80) which are common in the clinical setting. Further, our model achieved a ∼35% improvement in robustness against distribution shifts in the acquisition process compared to previously proposed DNN models. Finally, Our P₂T₂ model produces the most detailed Myelin-Water fraction maps compared to baseline approaches when applied to real human MRI data. Our P₂T₂ model offers a reliable and precise means of estimating T₂ distributions from MRI data and shows promise for use in large-scale multi-institutional trials with heterogeneous acquisition protocols. Our source code is available at: https://github.com/Hben-atya/P2T2-Robust-T2-estimation.git.

Signal-to-noise ratio penalties from a loss of stimulated echoes when using slab-selective excitation in three-dimensional fast spin echo imaging with long echo trains

Three-dimensional fast spin echo imaging with long echo trains combines high resolution with reasonable acquisition times and reduced specific absorption rate due to low refocusing flip angles. Typically, an entire volume is encoded (nonselective excitation) or localization can be performed with slab select excitation, which uses a long 90° pulse for precise localization, followed by a preliminary nonselective 180° pulse bounded by spoiler gradients to destroy signal outside of the volume of interest. Subsequent flip angles in the train are nonselective and identical between the two methods. The inclusion of the initial selective pulse and spoiler gradients results in a signal-to-noise ratio (SNR) penalty for slab selection, beyond the slice-averaging dependence, arising from a loss of stimulated echoes. SNR differences are explored using Bloch equation simulations of a T2-weighted 96 echo train sequence with varying parameters including T2, T1, and B1⁺ and compared with phantom and in vivo brain, neck, and knee experiments. In vivo SNR measurements in the three regions showed a maximum decrease of selective SNR by 29% (gastrocnemius muscle), 25% (pons), and 22% (globus pallidus), despite similar experimental parameters to nonselective experiments. Decreased SNR was compounded by B1⁺ variation affecting prescribed flip angles with further smaller reductions with T2 and T1 times. In conclusion, the elimination of coherences via the preliminary nominal 180° pulse and spoiler gradients in addition to the extended echo timing from the long excitation pulse resulted in a reduction in SNR compared with the nonselective case. Consideration of the required SNR and chosen anatomy as well as sequence restrictions should be weighed before choosing slab-selective excitation.

Fast and accurate T2 mapping using Bloch simulations and low-rank plus sparse matrix decomposition

PURPOSE

MRI's T2 relaxation time is one of the key contrast mechanisms for clinical diagnosis and prognosis of pathologies. Mapping this relaxation time, however, involves extensive scan times, which are needed to collect quantitative data, thereby impeding its integration into clinical routine. This study employs a low-rank plus sparse (L + S) signal decomposition approach in order to reconstruct accurate T2-maps from highly undersampled multi-echo spin-echo (MESE) MRI data.

METHODS

Two new algorithms are presented: the first uses standard L + S approach, where both L and S are iteratively updated. The second technique, dubbed SPArse and fixed RanK (SPARK), uses a fixed-rank L, under the assumption that most MESE information is found in the L component and that this rank can be pre-calculated. The utility of these new techniques is demonstrated on in vivo brain and calf data at x2 to x6 acceleration factors.

RESULTS

Accelerated T2 maps showed improved accuracy compared to fully sampled ground truth maps, when using L + S and SPARK techniques vis-à-vis standard GRAPPA acceleration.

CONCLUSION

SPARK provides accurate T2 maps with increased robustness to the selection of reconstruction parameters making it suitable to a wide range of applications and facilitating the use of quantitative T2 information in clinical settings.

Estimation of subvoxel fat infiltration in neurodegenerative muscle disorders using quantitative multi-T2 analysis

MRI's T₂ relaxation time is a valuable biomarker for neuromuscular disorders and muscle dystrophies. One of the hallmarks of these pathologies is the infiltration of adipose tissue and a loss of muscle volume. This leads to a mixture of two signal components, from fat and from water, to appear in each imaged voxel, each having a specific T₂ relaxation time. In this proof-of-concept work, we present a technique that can separate the signals from water and from fat within each voxel, measure their separate T₂ values, and calculate their relative fractions. The echo modulation curve (EMC) algorithm is a dictionary-based technique that offers accurate and reproducible mapping of T₂ relaxation times. We present an extension of the EMC algorithm for estimating subvoxel fat and water fractions, alongside the T₂ and proton-density values of each component. To facilitate data processing, calf and thigh anatomy were automatically segmented using a fully convolutional neural network and FSLeyes software. The preprocessing included creating two signal dictionaries, for water and for fat, using Bloch simulations of the prospective protocol. Postprocessing included voxelwise fitting for two components, by matching the experimental decay curve to a linear combination of the two simulated dictionaries. Subvoxel fat and water fractions and relaxation times were generated and used to calculate a new quantitative biomarker, termed viable muscle index, and reflecting disease severity. This biomarker indicates the fraction of remaining muscle out of the entire muscle region. The results were compared with those using the conventional Dixon technique, showing high agreement (R = 0.98, p < 0.001). It was concluded that the new extension of the EMC algorithm can be used to quantify abnormal fat infiltration as well as identify early inflammatory processes corresponding to elevation in the T₂ value of the water (muscle) component. This new ability may improve the diagnostic accuracy of neuromuscular diseases, help stratification of patients according to disease severity, and offer an efficient tool for tracking disease progression.

Inline dual-echo T2 quantification in brain using a fast mapping reconstruction technique

T2 mapping from 2D proton density and T2-weighted images (PD-T2) using Bloch equation simulations can be time consuming and introduces a latency between image acquisition and T2 map production. A fast T2 mapping reconstruction method is investigated and compared with a previous modeling approach to reduce computation time and allow inline T2 maps on the MRI console. Brain PD-T2 images from five multiple sclerosis patients were used to compare T2 map reconstruction times between the new subtraction method and the Euclidean norm minimization technique. Bloch equation simulations were used to create the lookup table for decay curve matching in both cases. Agreement of the two techniques used Bland-Altman analysis for investigating individual subsets of data and all image points in the five volumes (meta-analysis). The subtraction method resulted in an average reduction of computation time for single slices from 134 s (minimization method) to 0.44 s. Comparing T2 values between the subtraction and minimization methods resulted in a confidence interval ranging from -0.06 to 0.06 ms (95% of values were within ± 0.06 ms between the techniques). Using identical reconstruction code based on the subtraction method, inline T2 maps were produced from PD-T2 images directly on the scanner console. The excellent agreement between the two methods permits the subtraction technique to be interchanged with the previous method, reducing computation time and allowing inline T2 map reconstruction based on Bloch simulations directly on the scanner.

Mapping of magnetic resonance imaging's transverse relaxation time at low signal-to-noise ratio using Bloch simulations and principal component analysis image denoising

High-resolution mapping of magnetic resonance imaging (MRI)'s transverse relaxation time (T₂ ) can benefit many clinical applications by offering improved anatomic details, enhancing the ability to probe tissues' microarchitecture, and facilitating the identification of early pathology. Increasing spatial resolutions, however, decreases data's signal-to-noise ratio (SNR), particularly at clinical scan times. This impairs imaging quality, and the accuracy of subsequent radiological interpretation. Recently, principal component analysis (PCA) was employed for denoising diffusion-weighted MR images and was shown to be effective for improving parameter estimation in multiexponential relaxometry. This study combines the Marchenko-Pastur PCA (MP-PCA) signal model with the echo modulation curve (EMC) algorithm for denoising multiecho spin-echo (MESE) MRI data and improving the precision of EMC-generated single T₂ relaxation maps. The denoising technique was validated on simulations, phantom scans, and in vivo brain and knee data. MESE scans were performed on a 3-T Siemens scanner. The acquired images were denoised using the MP-PCA algorithm and were then provided as input for the EMC T₂ -fitting algorithm. Quantitative analysis of the denoising quality included comparing the standard deviation and coefficient of variation of T₂ values, along with gold standard SNR estimation of the phantom scans. The presented denoising technique shows an increase in T₂ maps' precision and SNR, while successfully preserving the morphological features of the tissue. Employing MP-PCA denoising as a preprocessing step decreases the noise-related variability of T₂ maps produced by the EMC algorithm and thus increases their precision. The proposed method can be useful for a wide range of clinical applications by facilitating earlier detection of pathologies and improving the accuracy of patients' follow-up.

Determining the internal orientation, degree of ordering, and volume of elongated nanocavities by NMR: Application to studies of plant stem

This study investigates the fibril nanostructure of fresh celery samples by modeling the anisotropic behavior of the transverse relaxation time (T₂) in nuclear magnetic resonance (NMR). Experimental results are interpreted within the framework of a previously developed theory, which was successfully used to model the nanostructures of several biological tissues as a set of water filled nanocavities, hence explaining the anisotropy the T₂ relaxation time in vivo. An important feature of this theory is to determine the degree of orientational ordering of the nanocavities, their characteristic volume, and their average direction with respect to the macroscopic sample. Results exhibit good agreement between theory and experimental data, which are, moreover, supported by optical microscopic resolution. The quantitative NMR approach presented herein can be potentially used to determine the internal ordering of biological tissues noninvasively.

Quantification of Intra-Muscular Adipose Infiltration in Calf/Thigh MRI Using Fully and Weakly Supervised Semantic Segmentation

Purpose: Infiltration of fat into lower limb muscles is one of the key markers for the severity of muscle pathologies. The level of fat infiltration varies in its severity across and within patients, and it is traditionally estimated using visual radiologic inspection. Precise quantification of the severity and spatial distribution of this pathological process requires accurate segmentation of lower limb anatomy into muscle and fat. Methods: Quantitative magnetic resonance imaging (qMRI) of the calf and thigh muscles is one of the most effective techniques for estimating pathological accumulation of intra-muscular adipose tissue (IMAT) in muscular dystrophies. In this work, we present a new deep learning (DL) network tool for automated and robust segmentation of lower limb anatomy that is based on the quantification of MRI’s transverse (T2) relaxation time. The network was used to segment calf and thigh anatomies into viable muscle areas and IMAT using a weakly supervised learning process. A new disease biomarker was calculated, reflecting the level of abnormal fat infiltration and disease state. A biomarker was then applied on two patient populations suffering from dysferlinopathy and Charcot–Marie–Tooth (CMT) diseases. Results: Comparison of manual vs. automated segmentation of muscle anatomy, viable muscle areas, and intermuscular adipose tissue (IMAT) produced high Dice similarity coefficients (DSCs) of 96.4%, 91.7%, and 93.3%, respectively. Linear regression between the biomarker value calculated based on the ground truth segmentation and based on automatic segmentation produced high correlation coefficients of 97.7% and 95.9% for the dysferlinopathy and CMT patients, respectively. Conclusions: Using a combination of qMRI and DL-based segmentation, we present a new quantitative biomarker of disease severity. This biomarker is automatically calculated and, most importantly, provides a spatially global indication for the state of the disease across the entire thigh or calf.

Correcting for imaging gradients-related bias of T2 relaxation times at high-resolution MRI

Purpose

High‐resolution animal imaging is an integral part of preclinical drug development and the investigation of diseases' pathophysiology. Quantitative mapping of T₂ relaxation times (qT₂) is a valuable tool for both preclinical and research applications, providing high sensitivity to subtle tissue pathologies. High‐resolution T₂ mapping, however, suffers from severe underestimation of T₂ values due to molecular diffusion. This affects both single‐echo and multi‐echo spin echo (SSE and MESE), on top of the well‐known contamination of MESE signals by stimulated echoes, and especially on high‐field and preclinical scanners in which high imaging gradients are used in comparison to clinical scanners.

Methods

Diffusion bias due to imaging gradients was analyzed by quantifying the effective b‐value for each coherence pathway in SSE and MESE protocols, and incorporating this information in a joint T₂‐diffusion reconstruction algorithm. Validation was done on phantoms and in vivo mouse brain using a 9.4T and a 7T MRI scanner.

Results

Underestimation of T₂ values due to strong imaging gradients can reach up to 70%, depending on scan parameters and on the sample's diffusion coefficient. The algorithm presented here produced T₂ values that agreed with reference spectroscopic measurements, were reproducible across scan settings, and reduced the average bias of T₂ values from −33.5 ± 20.5% to −0.1 ± 3.6%.

Conclusions

A new joint T₂‐diffusion reconstruction algorithm is able to negate imaging gradient–related underestimation of T₂ values, leading to reliable mapping of T₂ values at high resolutions.

Click here for author‐reader discussions

Advanced Magnetic Resonance Imaging (MRI) Techniques: Technical Principles and Applications in Nanomedicine

Simple Summary

Magnetic Resonance Imaging (MRI) is a consolidated imaging tool for the multiparametric assessment of tissues in various pathologies from degenerative and inflammatory diseases to cancer. In recent years, the continuous technological evolution of the equipment has led to the development of sequences that provide not only anatomical but also functional and metabolic information. In addition, there is a growing and emerging field of research in clinical applications using MRI to exploit the diagnostic and therapeutic capabilities of nanocompounds. This review illustrates the application of the most advanced magnetic resonance techniques in the field of nanomedicine.

Abstract

In the last decades, nanotechnology has been used in a wide range of biomedical applications, both diagnostic and therapeutic. In this scenario, imaging techniques represent a fundamental tool to obtain information about the properties of nanoconstructs and their interactions with the biological environment in preclinical and clinical settings. This paper reviews the state of the art of the application of magnetic resonance imaging in the field of nanomedicine, as well as the use of nanoparticles as diagnostic and therapeutic tools, especially in cancer, including the characteristics that hinder the use of nanoparticles in clinical practice.

Accelerated reconstruction of dictionary-based T2 relaxation maps based on dictionary compression and gradient descent search algorithms

Background Quantitative T₂-relaxation-based contrast maps have shown to be highly beneficial for clinical diagnosis and follow-up. The generation of quantitative maps, however, is impaired by long acquisition times, and time-consuming post-processing schemes. The EMC platform is a dictionary-based technique, which involves simulating theoretical signal curves for different physical and experimental values, followed by matching the experimentally acquired signals to the set simulated ones. Purpose Although the EMC technique has shown to produce accurate T₂ maps, it involves computationally intensive post-processing procedures. In this work we present an approach for accelerating the reconstruction of T₂ relaxation maps. Methods This work presents two alternative post-processing approaches for accelerating the reconstruction of EMC-based T₂ relaxation maps. These are (a) Dictionary compression using principal component analysis (PCA) and (b) gradient-descent search algorithm. Additional acceleration was achieved by finding the optimal MATLAB C++ compiler. The utility of the two suggested approaches was examined by calculating the relative error, produced by each technique. Results Gradient descent method was in perfect agreement with the ground truth exhaustive search matching process. PCA based acceleration produced root mean square error (RMSE) of up to 4% compared to exhaustive matching process. Overall acceleration of x16 was achieved using gradient descent in addition to x7 acceleration by choosing the optimal MATLAB C++ compiler. Conclusions Postprocessing of EMC-based T₂ relaxation maps can be accelerated without impairing the accuracy of the ensuing T₂ values.

Data-driven algorithm for myelin water imaging: Probing subvoxel compartmentation based on identification of spatially global tissue features

PURPOSE

Multicomponent analysis of MRI T₂ relaxation time (mcT₂ ) is commonly used for estimating myelin content by separating the signal at each voxel into its underlying distribution of T₂ values. This voxel-based approach is challenging due to the large ambiguity in the multi-T₂ space and the low SNR of MRI signals. Herein, we present a data-driven mcT₂ analysis, which utilizes the statistical strength of identifying spatially global mcT₂ motifs in white matter segments before deconvolving the local signal at each voxel.

METHODS

Deconvolution is done using a tailored optimization scheme, which incorporates the global mcT₂ motifs without additional prior assumptions regarding the number of microscopic components. The end results of this process are voxel-wise myelin water fraction maps.

RESULTS

Validations are shown for computer-generated signals, uniquely designed subvoxel mcT₂ phantoms, and in vivo human brain. Results demonstrated excellent fitting accuracy, both for the numerical and the physical mcT₂ phantoms, exhibiting excellent agreement between calculated myelin water fraction and ground truth. Proof-of-concept in vivo validation is done by calculating myelin water fraction maps for white matter segments of the human brain. Interscan stability of myelin water fraction values was also estimated, showing good correlation between scans.

CONCLUSION

We conclude that studying global tissue motifs prior to performing voxel-wise mcT₂ analysis stabilizes the optimization scheme and efficiently overcomes the ambiguity in the T₂ space. This new approach can improve myelin water imaging and the investigation of microstructural compartmentation in general.

A unifying view on extended phase graphs and Bloch simulations for quantitative MRI

Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.