Fast 3D magnetic resonance fingerprinting for a whole-brain coverage

Fast and reliable quantitative measures of white matter development with magnetic resonance fingerprinting

Developmental cognitive neuroscience aims to shed light on evolving relationships between brain structure and cognitive development. To this end, quantitative methods that reliably measure individual differences in brain tissue properties are fundamental. Standard qualitative MRI sequences are influenced by scan parameters and hardware-related biases, and also lack physical units, making the analysis of individual differences problematic. In contrast, quantitative MRI can measure physical properties of the tissue but with the cost of long scan durations and sensitivity to motion. This poses a critical limitation for studying young children. Here, we examine the reliability and validity of an efficient quantitative multiparameter mapping method - Magnetic Resonance Fingerprinting (MRF) - in children scanned longitudinally. We focus on T1 values in white matter, since quantitative T1 values are known to primarily reflect myelin content, a key factor in brain development. Forty-nine children aged 8-13y (mean 10.3y ±1.4) completed two scanning sessions 2-4 months apart. In each session, two 2-minute 3D-MRF scans at 1mm isotropic resolution were collected to evaluate the effect of scan duration on image quality and scan-rescan reliability. A separate calibration scan was used to measure B0 inhomogeneity and correct for bias. We examined the impact of scan time and B0 inhomogeneity correction on scan-rescan reliability of values in white matter, by comparing single 2-min and combined two 2-min scans, with and without B0-correction. Whole-brain voxel-based reliability analysis showed that combining two 2-min MRF scans improved reliability (pearson's r=0.87) compared with a single 2-min scan (r=0.84), while B0-correction had no effect on reliability in white matter (r=0.86 and 0.83 4-min vs 2-min). Using diffusion tractography, we delineated MRF-derived T1 profiles along major white matter fiber tracts and found similar or higher reliability for T1 from MRF compared to diffusion parameters (based on a 10-minute dMRI scan). Lastly, we found that T1 values in multiple white matter tracts were significantly correlated with age. In sum, MRF-derived T1 values were highly reliable in a longitudinal sample of children and replicated known age effects. Reliability in white matter was improved by longer scan duration but was not affected by B0-correction, making it a quick and straightforward scan to collect. We propose that MRF provides a promising avenue for acquiring quantitative brain metrics in children and patient populations where scan time and motion are of particular concern.

Magnetic resonance relaxometry in quantitative imaging of brain gliomas: A literature review

Magnetic resonance (MR) relaxometry is a quantitative imaging method that measures tissue relaxation properties. This review discusses the state of the art of clinical proton MR relaxometry for glial brain tumors. Current MR relaxometry technology also includes MR fingerprinting and synthetic MRI, which solve the inefficiencies and challenges of earlier techniques. Despite mixed results regarding its capability for brain tumor differential diagnosis, there is growing evidence that MR relaxometry can differentiate between gliomas and metastases and between glioma grades. Studies of the peritumoral zones have demonstrated their heterogeneity and possible directions of tumor infiltration. In addition, relaxometry offers T2* mapping that can define areas of tissue hypoxia not discriminated by perfusion assessment. Studies of tumor therapy response have demonstrated an association between survival and progression terms and dynamics of native and contrast-enhanced tumor relaxometric profiles. In conclusion, MR relaxometry is a promising technique for glial tumor diagnosis, particularly in association with neuropathological studies and other imaging techniques.

Quantifying 3D MR fingerprinting (3D-MRF) reproducibility across subjects, sessions, and scanners automatically using MNI atlases

PURPOSE

Quantitative MRI techniques such as MR fingerprinting (MRF) promise more objective and comparable measurements of tissue properties at the point-of-care than weighted imaging. However, few direct cross-modal comparisons of MRF's repeatability and reproducibility versus weighted acquisitions have been performed. This work proposes a novel fully automated pipeline for quantitatively comparing cross-modal imaging performance in vivo via atlas-based sampling.

METHODS

We acquire whole-brain 3D-MRF, turbo spin echo, and MPRAGE sequences three times each on two scanners across 10 subjects, for a total of 60 multimodal datasets. The proposed automated registration and analysis pipeline uses linear and nonlinear registration to align all qualitative and quantitative DICOM stacks to Montreal Neurological Institute (MNI) 152 space, then samples each dataset's native space through transformation inversion to compare performance within atlas regions across subjects, scanners, and repetitions.

RESULTS

Voxel values within MRF-derived maps were found to be more repeatable (σT1  = 1.90, σT2  = 3.20) across sessions than vendor-reconstructed MPRAGE (σT1w  = 6.04) or turbo spin echo (σT2w  = 5.66) images. Additionally, MRF was found to be more reproducible across scanners (σT1  = 2.21, σT2  = 3.89) than either qualitative modality (σT1w  = 7.84, σT2w  = 7.76). Notably, differences between repeatability and reproducibility of in vivo MRF were insignificant, unlike the weighted images.

CONCLUSION

MRF data from many sessions and scanners can potentially be treated as a single dataset for harmonized analysis or longitudinal comparisons without the additional regularization steps needed for qualitative modalities.

Motion Robust MR Fingerprinting Scan to Image Neonates With Prenatal Opioid Exposure

PURPOSE

To explore whether MR fingerprinting (MRF) scans provide motion-robust and quantitative brain tissue measurements for non-sedated infants with prenatal opioid exposure (POE).

STUDY TYPE

Prospective.

POPULATION

13 infants with POE (3 male; 12 newborns (age 7-65 days) and 1 infant aged 9-months).

FIELD STRENGTH/SEQUENCE

3T, 3D T1-weighted MPRAGE, 3D T2-weighted TSE and MRF sequences.

ASSESSMENT

The image quality of MRF and MRI was assessed in a fully crossed, multiple-reader, multiple-case study. Sixteen image quality features in three types-image artifacts, structure and myelination visualization-were ranked by four neuroradiologists (8, 7, 5, and 8 years of experience respectively), using a 3-point scale. MRF T1 and T2 values in 8 white matter brain regions were compared between babies younger than 1 month and babies between 1 and 2 months.

STATISTICAL TESTS

Generalized estimating equations model to test the significance of differences of regional T1 and T2 values of babies under 1 month and those older. MRI and MRF image quality was assessed using Gwet's second order auto-correlation coefficient (AC2) with confidence levels. The Cochran-Mantel-Haenszel test was used to assess the difference in proportions between MRF and MRI for all features and stratified by the type of features. A P value <0.05 was considered statistically significant.

RESULTS

The MRF of two infants were excluded in T1 and T2 value analysis due to severe motion artifact but were included in the image quality assessment. In infants under 1 month of age (N = 6), the T1 and T2 values were significantly higher compared to those between 1 and 2 months of age (N = 4). MRF images showed significantly higher image quality ratings in all three feature types compared to MRI images.

CONCLUSIONS

MR Fingerprinting scans have potential to be a motion-robust and efficient method for nonsedated infants.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Physics-Informed Discretization for Reproducible and Robust Radiomic Feature Extraction Using Quantitative MRI

OBJECTIVE

Given the limited repeatability and reproducibility of radiomic features derived from weighted magnetic resonance imaging (MRI), there may be significant advantages to using radiomics in conjunction with quantitative MRI. This study introduces a novel physics-informed discretization (PID) method for reproducible radiomic feature extraction and evaluates its performance using quantitative MRI sequences including magnetic resonance fingerprinting (MRF) and apparent diffusion coefficient (ADC) mapping.

MATERIALS AND METHODS

A multiscanner, scan-rescan dataset comprising whole-brain 3D quantitative (MRF T1, MRF T2, and ADC) and weighted MRI (T1w MPRAGE, T2w SPACE, and T2w FLAIR) from 5 healthy subjects was prospectively acquired. Subjects underwent 2 repeated acquisitions on 3 distinct 3 T scanners each, for a total of 6 scans per subject (30 total scans). First-order statistical (n = 23) and second-order texture (n = 74) radiomic features were extracted from 56 brain tissue regions of interest using the proposed PID method (for quantitative MRI) and conventional fixed bin number (FBN) discretization (for quantitative MRI and weighted MRI). Interscanner radiomic feature reproducibility was measured using the intraclass correlation coefficient (ICC), and the effect of image sequence (eg, MRF T1 vs T1w MPRAGE), as well as image discretization method (ie, PID vs FBN), on radiomic feature reproducibility was assessed using repeated measures analysis of variance. The robustness of PID and FBN discretization to segmentation error was evaluated by simulating segmentation differences in brainstem regions of interest. Radiomic features with ICCs greater than 0.75 following simulated segmentation were determined to be robust to segmentation.

RESULTS

First-order features demonstrated higher reproducibility in quantitative MRI than weighted MRI sequences, with 30% (n = 7/23) features being more reproducible in MRF T1 and MRF T2 than weighted MRI. Gray level co-occurrence matrix (GLCM) texture features extracted from MRF T1 and MRF T2 were significantly more reproducible using PID compared with FBN discretization; for all quantitative MRI sequences, PID yielded the highest number of texture features with excellent reproducibility (ICC > 0.9). Comparing texture reproducibility of quantitative and weighted MRI, a greater proportion of MRF T1 (n = 225/370, 61%) and MRF T2 (n = 150/370, 41%) texture features had excellent reproducibility (ICC > 0.9) compared with T1w MPRAGE (n = 148/370, 40%), ADC (n = 115/370, 32%), T2w SPACE (n = 98/370, 27%), and FLAIR (n = 102/370, 28%). Physics-informed discretization was also more robust than FBN discretization to segmentation error, as 46% (n = 103/222, 46%) of texture features extracted from quantitative MRI using PID were robust to simulated 6 mm segmentation shift compared with 19% (n = 42/222, 19%) of weighted MRI texture features extracted using FBN discretization.

CONCLUSIONS

The proposed PID method yields radiomic features extracted from quantitative MRI sequences that are more reproducible and robust than radiomic features extracted from weighted MRI using conventional (FBN) discretization approaches. Quantitative MRI sequences also demonstrated greater scan-rescan robustness and first-order feature reproducibility than weighted MRI.

Deep learning-assisted preclinical MR fingerprinting for sub-millimeter T1 and T2 mapping of entire macaque brain

PURPOSE

Preclinical MR fingerprinting (MRF) suffers from long acquisition time for organ-level coverage due to demanding image resolution and limited undersampling capacity. This study aims to develop a deep learning-assisted fast MRF framework for sub-millimeter T1 and T2 mapping of entire macaque brain on a preclinical 9.4 T MR system.

METHODS

Three dimensional MRF images were reconstructed by singular value decomposition (SVD) compressed reconstruction. T1 and T2 mapping for each axial slice exploited a self-attention assisted residual U-Net to suppress aliasing-induced quantification errors, and the transmit-field (B1 + ) measurements for robustness against B1 + inhomogeneity. Supervised network training used MRF images simulated via virtual parametric maps and a desired undersampling scheme. This strategy bypassed the difficulties of acquiring fully sampled preclinical MRF data to guide network training. The proposed fast MRF framework was tested on experimental data acquired from ex vivo and in vivo macaque brains.

RESULTS

The trained network showed reasonable adaptability to experimental MRF images, enabling robust delineation of various T1 and T2 distributions in the brain tissues. Further, the proposed MRF framework outperformed several existing fast MRF methods in handling the aliasing artifacts and capturing detailed cerebral structures in the mapping results. Parametric mapping of entire macaque brain at nominal resolution of 0.35× $$ \times $$ 0.35× $$ \times $$ 1 mm3 can be realized via a 20-min 3D MRF scan, which was sixfold faster than the baseline protocol.

CONCLUSION

Introducing deep learning to MRF framework paves the way for efficient organ-level high-resolution quantitative MRI in preclinical applications.

Emerging Trends in Magnetic Resonance Fingerprinting for Quantitative Biomedical Imaging Applications: A Review

Magnetic resonance imaging (MRI) stands as a vital medical imaging technique, renowned for its ability to offer high-resolution images of the human body with remarkable soft-tissue contrast. This enables healthcare professionals to gain valuable insights into various aspects of the human body, including morphology, structural integrity, and physiological processes. Quantitative imaging provides compositional measurements of the human body, but, currently, either it takes a long scan time or is limited to low spatial resolutions. Undersampled k-space data acquisitions have significantly helped to reduce MRI scan time, while compressed sensing (CS) and deep learning (DL) reconstructions have mitigated the associated undersampling artifacts. Alternatively, magnetic resonance fingerprinting (MRF) provides an efficient and versatile framework to acquire and quantify multiple tissue properties simultaneously from a single fast MRI scan. The MRF framework involves four key aspects: (1) pulse sequence design; (2) rapid (undersampled) data acquisition; (3) encoding of tissue properties in MR signal evolutions or fingerprints; and (4) simultaneous recovery of multiple quantitative spatial maps. This paper provides an extensive literature review of the MRF framework, addressing the trends associated with these four key aspects. There are specific challenges in MRF for all ranges of magnetic field strengths and all body parts, which can present opportunities for further investigation. We aim to review the best practices in each key aspect of MRF, as well as for different applications, such as cardiac, brain, and musculoskeletal imaging, among others. A comprehensive review of these applications will enable us to assess future trends and their implications for the translation of MRF into these biomedical imaging applications.

Correcting synthetic MRI contrast-weighted images using deep learning

Synthetic magnetic resonance imaging (MRI) offers a scanning paradigm where a fast multi-contrast sequence can be used to estimate underlying quantitative tissue parameter maps, which are then used to synthesize any desirable clinical contrast by retrospectively changing scan parameters in silico. Two benefits of this approach are the reduced exam time and the ability to generate arbitrary contrasts offline. However, synthetically generated contrasts are known to deviate from the contrast of experimental scans. The reason for contrast mismatch is the necessary exclusion of some unmodeled physical effects such as partial voluming, diffusion, flow, susceptibility, magnetization transfer, and more. The inclusion of these effects in signal encoding would improve the synthetic images, but would make the quantitative imaging protocol impractical due to long scan times. Therefore, in this work, we propose a novel deep learning approach that generates a multiplicative correction term to capture unmodeled effects and correct the synthetic contrast images to better match experimental contrasts for arbitrary scan parameters. The physics inspired deep learning model implicitly accounts for some unmodeled physical effects occurring during the scan. As a proof of principle, we validate our approach on synthesizing arbitrary inversion recovery fast spin-echo scans using a commercially available 2D multi-contrast sequence. We observe that the proposed correction visually and numerically reduces the mismatch with experimentally collected contrasts compared to conventional synthetic MRI. Finally, we show results of a preliminary reader study and find that the proposed method statistically significantly improves in contrast and SNR as compared to synthetic MR images.

High-fidelity direct contrast synthesis from magnetic resonance fingerprinting

PURPOSE

This work was aimed at proposing a supervised learning-based method that directly synthesizes contrast-weighted images from the Magnetic Resonance Fingerprinting (MRF) data without performing quantitative mapping and spin-dynamics simulations.

METHODS

To implement our direct contrast synthesis (DCS) method, we deploy a conditional generative adversarial network (GAN) framework with a multi-branch U-Net as the generator and a multilayer CNN (PatchGAN) as the discriminator. We refer to our proposed approach as N-DCSNet. The input MRF data are used to directly synthesize T1-weighted, T2-weighted, and fluid-attenuated inversion recovery (FLAIR) images through supervised training on paired MRF and target spin echo-based contrast-weighted scans. The performance of our proposed method is demonstrated on in vivo MRF scans from healthy volunteers. Quantitative metrics, including normalized root mean square error (nRMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), learned perceptual image patch similarity (LPIPS), and Fréchet inception distance (FID), were used to evaluate the performance of the proposed method and compare it with others.

RESULTS

In-vivo experiments demonstrated excellent image quality with respect to that of simulation-based contrast synthesis and previous DCS methods, both visually and according to quantitative metrics. We also demonstrate cases in which our trained model is able to mitigate the in-flow and spiral off-resonance artifacts typically seen in MRF reconstructions, and thus more faithfully represent conventional spin echo-based contrast-weighted images.

CONCLUSION

We present N-DCSNet to directly synthesize high-fidelity multicontrast MR images from a single MRF acquisition. This method can significantly decrease examination time. By directly training a network to generate contrast-weighted images, our method does not require any model-based simulation and therefore can avoid reconstruction errors due to dictionary matching and contrast simulation (code available at:https://github.com/mikgroup/DCSNet).

Magnetic resonance fingerprinting in multiple sclerosis

BACKGROUND

In this cross sectional study, we used MRF to investigate tissue properties of normal-appearing white matter, gray matter, and lesions in relapsing remitting MS (n = 21), secondary progressive MS (n = 16) and healthy controls (n = 9). A FISP-based MRF sequence was used for acquisition, imaging time 5 min 15 s. MRF T1 and T2 relaxation times were measured from lesional tissue, normal-appearing frontal white matter, corpus callous, thalamus, and caudate. Differences between healthy controls and MS were examined using ANCOVA adjusted for age and sex. Spearman rank correlations were assessed between T1 and T2 relaxation times and clinical measures.

OBJECTIVES

To examine brain T1 and T2 values using magnetic resonance fingerprinting (MRF) in healthy controls and MS.

METHODS

The subjects included 21 relapsing-remitting (RR) MS, 16 secondary progressive (SP) MS, and 9 age- and sex-matched HC without manifest neurological disease participating in a longitudinal MRI study. A 3T/ FISP-based MRF sequence was acquired. Regions of interest were drawn for lesions and normal appearing white matter. ANCOVA adjusted for age and sex were used to compare the groups with significance set at 0.05.

RESULTS

A step-wise increase in T1 and T2 relaxation times was found between healthy controls, relapsing remitting MS, and secondary progressive MS. Significant differences were found in T1 and T2 between MS and healthy controls in the frontal normal-appearing white matter, corpus callosum, and thalamus (p < 0.04 for all). Significant differences in T1 and T2 between RR and SPMS were found in the frontal normal-appearing white matter and T2 lesions (p < 0.02 for all). T1 relaxation from the frontal normal-appearing white matter correlated with the Expanded Disability Status Scale [ρ = 0.62, p < 0.001], timed 25 foot walk (ρ = 0.45, p = 0.01), 9 hole peg test (ρ = 0.62, p < 0.001), and paced auditory serial addition test (ρ = -0.4, p = 0.01).

CONCLUSION

These results suggest that MRF may be a clinically feasible quantitative approach for characterizing tissue damage in MS.

Improving motion robustness of 3D MR fingerprinting with a fat navigator

PURPOSE

To develop a 3D MR fingerprinting (MRF) method in combination with fat navigators to improve its motion robustness for neuroimaging.

METHODS

A rapid fat navigator was developed using the stack-of-spirals acquisition and non-Cartesian spiral GRAPPA. The fat navigator module was implemented in the 3D MRF sequence with high scan efficiency. The developed method was first validated in phantoms and five healthy subjects with intentional head motion. The method was further applied to infants with neonatal opioid withdrawal symptoms. The 3D MRF scans with fat navigators acquired with and without acceleration along the partition-encoding direction were both examined in the study.

RESULTS

Both phantom and in vivo results demonstrated that the added fat navigator modules did not influence the quantification accuracy in MRF. In combination with non-Cartesian spiral GRAPPA, a rapid fat navigator sampling with whole-brain coverage was achieved in ˜0.5 s at 3T, reducing its sensitivity to potential motion. Based on the motion waveforms extracted from fat navigators, the motion robustness of the 3D MRF was largely improved. With the proposed method, the motion-corrupted MRF datasets yielded T1 and T2 maps with significantly reduced artifacts and high correlations with measurements from the reference motion-free MRF scans.

CONCLUSION

We developed a 3D MRF method coupled with rapid fat navigators to improve its motion robustness for quantitative neuroimaging. Our results demonstrate that (1) accurate tissue quantification was preserved with the fat navigator modules and (2) the motion robustness for quantitative tissue mapping was largely improved with the developed method.

SSL-QALAS: Self-Supervised Learning for rapid multiparameter estimation in quantitative MRI using 3D-QALAS

PURPOSE

To develop and evaluate a method for rapid estimation of multiparametric T1 , T2 , proton density, and inversion efficiency maps from 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) measurements using self-supervised learning (SSL) without the need for an external dictionary.

METHODS

An SSL-based QALAS mapping method (SSL-QALAS) was developed for rapid and dictionary-free estimation of multiparametric maps from 3D-QALAS measurements. The accuracy of the reconstructed quantitative maps using dictionary matching and SSL-QALAS was evaluated by comparing the estimated T1 and T2 values with those obtained from the reference methods on an International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom. The SSL-QALAS and the dictionary-matching methods were also compared in vivo, and generalizability was evaluated by comparing the scan-specific, pre-trained, and transfer learning models.

RESULTS

Phantom experiments showed that both the dictionary-matching and SSL-QALAS methods produced T1 and T2 estimates that had a strong linear agreement with the reference values in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom. Further, SSL-QALAS showed similar performance with dictionary matching in reconstructing the T1 , T2 , proton density, and inversion efficiency maps on in vivo data. Rapid reconstruction of multiparametric maps was enabled by inferring the data using a pre-trained SSL-QALAS model within 10 s. Fast scan-specific tuning was also demonstrated by fine-tuning the pre-trained model with the target subject's data within 15 min.

CONCLUSION

The proposed SSL-QALAS method enabled rapid reconstruction of multiparametric maps from 3D-QALAS measurements without an external dictionary or labeled ground-truth training data.

Three-dimensional Multi-parameter Mapping of Relaxation Times and Susceptibility Using Partially RF-spoiled Gradient Echo

PURPOSE

MR parameter mapping is a technique that obtains distributions of parameters such as relaxation time and proton density (PD) and is starting to be used for disease quantification in clinical diagnoses. Quantitative susceptibility mapping is also promising for the early diagnosis of brain disorders such as degenerative neurological disorders. Therefore, we developed an MR quantitative parameter mapping (QPM) method to map four tissue-related parameters (T1, T2*, PD, and susceptibility) and B1 simultaneously by using a 3D partially RF-spoiled gradient echo (pRSGE). We verified the accuracy and repeatability of QPM in phantom and volunteer experiments.

METHODS

Tissue-related parameters are estimated by varying four scan parameters of the 3D pRSGE: flip angle, RF-pulse phase increment, TR and TE, performing multiple image scans, and finding a least-squares fit for an intensity function (which expresses the relationship between the scan parameters and intensity values). The intensity function is analytically complex, but by using a Bloch simulation to create it numerically, the least-squares fitting can be used to estimate the quantitative values. This has the advantage of shortening the image-reconstruction processing time needed to estimate the quantitative values than with methods using pattern matching.

RESULTS

A 1.1-mm isotropic resolution scan covering the whole brain was completed with a scan time of approximately 12 minutes, and the reconstruction time using a GPU was approximately 1 minute. The phantom experiments confirmed that both the accuracy and repeatability of the quantitative values were high. The volunteer scans also confirmed that the accuracy of the quantitative values was comparable to that of conventional methods.

CONCLUSION

The proposed QPM method can map T1, T2*, PD, susceptibility, and B1 simultaneously within a scan time that can be applied to human subjects.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

Magnetic Resonance Fingerprinting: A Review of Clinical Applications

ABSTRACT

Magnetic resonance fingerprinting (MRF) is an approach to quantitative magnetic resonance imaging that allows for efficient simultaneous measurements of multiple tissue properties, which are then used to create accurate and reproducible quantitative maps of these properties. As the technique has gained popularity, the extent of preclinical and clinical applications has vastly increased. The goal of this review is to provide an overview of currently investigated preclinical and clinical applications of MRF, as well as future directions. Topics covered include MRF in neuroimaging, neurovascular, prostate, liver, kidney, breast, abdominal quantitative imaging, cardiac, and musculoskeletal applications.

Comparing Magnetic Resonance Fingerprinting (MRF) and the MAGiC Sequence for Simultaneous T1 and T2 Quantitative Measurements in the Female Pelvis: A Prospective Study

The aim of this study was to explore the potential of magnetic resonance fingerprinting (MRF), an emerging quantitative MRI technique, in measuring relaxation values of female pelvic tissues compared to the conventional magnetic resonance image compilation (MAGiC) sequence. The study included 32 female patients who underwent routine pelvic MRI exams using anterior and posterior array coils on a 3T clinical scanner. Our findings demonstrated significant correlations between MRF and MAGiC measured T1 and T2 values (p < 0.0001) for various pelvic tissues, including ilium, femoral head, gluteus, obturator, iliopsoas, erector spinae, uterus, cervix, and cutaneous fat. The tissue contrasts generated from conventional MRI and synthetic MRF also showed agreement in bone, muscle, and uterus for both T1-weighted and T2-weighted images. This study highlights the strengths of MRF in providing simultaneous T1 and T2 mapping. MRF offers distinct tissue contrast and has the potential for accurate diagnosis of female pelvic diseases, including tumors, fibroids, endometriosis, and pelvic inflammatory disease. Additionally, MRF shows promise in monitoring disease progression or treatment response. Overall, the study demonstrates the potential of MRF in the field of female pelvic organ imaging and suggests that it could be a valuable addition to the clinical practice of pelvic MRI exams. Further research is needed to establish the clinical utility of MRF and to develop standardized protocols for its implementation in clinical practice.

Adaptive model-based Magnetic Resonance

PURPOSE

Conventional sequences are static in nature, fixing measurement parameters in advance in anticipation of a wide range of expected tissue parameter values. We set out to design and benchmark a new, personalized approach-termed adaptive MR-in which incoming subject data is used to update and fine-tune the pulse sequence parameters in real time.

METHODS

We implemented an adaptive, real-time multi-echo (MTE) experiment for estimating T₂ s. Our approach combined a Bayesian framework with model-based reconstruction. It maintained and continuously updated a prior distribution of the desired tissue parameters, including T₂ , which was used to guide the selection of sequence parameters in real time.

RESULTS

Computer simulations predicted accelerations between 1.7- and 3.3-fold for adaptive multi-echo sequences relative to static ones. These predictions were corroborated in phantom experiments. In healthy volunteers, our adaptive framework accelerated the measurement of T₂ for n-acetyl-aspartate by a factor of 2.5.

CONCLUSION

Adaptive pulse sequences that alter their excitations in real time could provide substantial reductions in acquisition times. Given the generality of our proposed framework, our results motivate further research into other adaptive model-based approaches to MRI and MRS.

A fast MR fingerprinting simulator for direct error estimation and sequence optimization

Magnetic resonance fingerprinting (MRF) is a novel quantitative MR technique that simultaneously provides multiple tissue property maps. When optimizing MRF scans, modeling undersampling errors and field imperfections in cost functions for direct measurement of quantitative errors will make the optimization results more practical and robust. However, optimizing such cost function is computationally expensive and impractical for MRF optimization with tens of thousands of iterations. Here, we introduce a fast MRF simulator to simulate aliased images from actual scan scenarios including undersampling and system imperfections, which substantially reduces computational time and allows for direct error estimation of the quantitative maps and efficient sequence optimization. We evaluate the performance and computational speed of the proposed approach by simulations and in vivo experiments. The simulations from the proposed method closely approximate the signals and MRF maps from in vivo scans, with 158 times shorter processing time than the conventional simulation method using Non-uniform Fourier transform. We also demonstrate the power of applying the fast MRF simulator in MRF sequence optimization. The optimized sequences are validated with in vivo scans to assess the image quality and accuracy. The optimized sequences produce artifact-free T1 and T2 maps in 2D and 3D scans with equivalent mapping accuracy as the human-designed sequence but at shorter scan times. Incorporating the proposed simulator in the MRF optimization framework makes direct estimation of undersampling errors during the optimization process feasible, and provide optimized MRF sequences that are robust against undersampling artifacts and field inhomogeneity.

Optimized three-dimensional ultrashort echo time: Magnetic resonance fingerprinting for myelin tissue fraction mapping

Quantitative assessment of brain myelination has gained attention for both research and diagnosis of neurological diseases. However, conventional pulse sequences cannot directly acquire the myelin-proton signals due to its extremely short T2 and T2* values. To obtain the myelin-proton signals, dedicated short T2 acquisition techniques, such as ultrashort echo time (UTE) imaging, have been introduced. However, it remains challenging to isolate the myelin-proton signals from tissues with longer T2. In this article, we extended our previous two-dimensional ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) with dual-echo acquisition to three dimensional (3D). Given a relatively low proton density (PD) of myelin-proton, we utilized Cramér-Rao Lower Bound to encode myelin-proton with the maximal SNR efficiency for optimizing the MR fingerprinting design, in order to improve the sensitivity of the sequence to myelin-proton. In addition, with a second echo of approximately 3 ms, myelin-water component can be also captured. A myelin-tissue (myelin-proton and myelin-water) fraction mapping can be thus calculated. The optimized 3D UTE-MRF with dual-echo acquisition is tested in simulations, physical phantom and in vivo studies of both healthy subjects and multiple sclerosis patients. The results suggest that the rapidly decayed myelin-proton and myelin-water signal can be depicted with UTE signals of our method at clinically relevant resolution (1.8 mm isotropic) in 15 min. With its good sensitivity to myelin loss in multiple sclerosis patients demonstrated, our method for the whole brain myelin-tissue fraction mapping in clinical friendly scan time has the potential for routine clinical imaging.

Accuracy, repeatability, and reproducibility of T1 and T2 relaxation times measurement by 3D magnetic resonance fingerprinting with different dictionary resolutions

OBJECTIVES

To assess the accuracy, repeatability, and reproducibility of T₁ and T₂ relaxation time measurements by three-dimensional magnetic resonance fingerprinting (3D MRF) using various dictionary resolutions.

METHODS

The ISMRM/NIST phantom was scanned daily for 10 days in two 3 T MR scanners using a 3D MRF sequence reconstructed using four dictionaries with varying step sizes and one dictionary with wider ranges. Thirty-nine healthy volunteers were enrolled: 20 subjects underwent whole-brain MRF scans in both scanners and the rest in one scanner. ROI/VOI analyses were performed on phantom and brain MRF maps. Accuracy, repeatability, and reproducibility metrics were calculated.

RESULTS

In the phantom study, all dictionaries showed high T₁ linearity to the reference values (R² > 0.99), repeatability (CV < 3%), and reproducibility (CV < 3%) with lower linearity (R² > 0.98), repeatability (CV < 6%), and reproducibility (CV ≤ 4%) for T₂ measurement. The volunteer study demonstrated high T₁ reproducibility of within-subject CV (wCV) < 4% by all dictionaries with the same ranges, both in the brain parenchyma and CSF. Yet, reproducibility was moderate for T₂ measurement (wCV < 8%). In CSF measurement, dictionaries with a smaller range showed a seemingly better reproducibility (T₁, wCV 3%; T₂, wCV 8%) than the much wider range dictionary (T₁, wCV 5%; T₂, wCV 13%). Truncated CSF relaxometry values were evident in smaller range dictionaries.

CONCLUSIONS

The accuracy, repeatability, and reproducibility of 3D MRF across various dictionary resolutions were high for T₁ and moderate for T₂ measurements. A lower-resolution dictionary with a well-defined range may be adequate, thus significantly reducing the computational load.

KEY POINTS

• A lower-resolution dictionary with a well-defined range may be sufficient for 3D MRF reconstruction. • CSF relaxation times might be underestimated due to truncation by the upper dictionary range. • Dictionary with a higher upper range might be advisable, especially for CSF evaluation and elderly subjects whose perivascular spaces are more prominent.

Normative quantitative relaxation atlases for characterization of cortical regions using magnetic resonance fingerprinting

Quantitative magnetic resonance (MR) has been used to study cyto- and myelo-architecture of the human brain non-invasively. However, analyzing brain cortex using high-resolution quantitative MR acquisition can be challenging to perform using 3T clinical scanners. MR fingerprinting (MRF) is a highly efficient and clinically feasible quantitative MR technique that simultaneously provides T1 and T2 relaxation maps. Using 3D MRF from 40 healthy subjects (mean age = 25.6 ± 4.3 years) scanned on 3T magnetic resonance imaging, we generated whole-brain gyral-based normative MR relaxation atlases and investigated cortical-region-based T1 and T2 variations. Gender and age dependency of T1 and T2 variations were additionally analyzed. The coefficient of variation of T1 and T2 for each cortical-region was 3.5% and 7.3%, respectively, supporting low variability of MRF measurements across subjects. Significant differences in T1 and T2 were identified among 34 brain regions (P < 0.001), lower in the precentral, postcentral, paracentral lobule, transverse temporal, lateral occipital, and cingulate areas, which contain sensorimotor, auditory, visual, and limbic functions. Significant correlations were identified between age and T1 and T2 values. This study established whole-brain MRF T1 and T2 atlases of healthy subjects using a clinical 3T scanner, which can provide a quantitative and region-specific baseline for future brain studies and pathology detection.

High-Resolution 3D Magnetic Resonance Fingerprinting With a Graph Convolutional Network

Magnetic resonance fingerprinting (MRF) is a novel quantitative imaging framework for rapid and simultaneous quantification of multiple tissue properties. 3D MRF allows higher through-plane resolution, but the acquisition process is slow when whole-brain coverage is needed. Existing methods for acceleration mainly rely on GRAPPA for k-space interpolation in the partition-encoding direction, limiting the acceleration factor to 2 or 3. In this work, we replace GRAPPA with a deep learning approach for accurate tissue quantification with greater acceleration. Specifically, a graph convolution network (GCN) is developed to cater to the non-Cartesian spiral sampling trajectories typical in MRF acquisition. The GCN maintains high quantification accuracy with up to 6-fold acceleration and allows 1mm isotropic resolution whole-brain 3D MRF data to be acquired in 3min and submillimeter 3D MRF (0.8mm) in 5min, greatly improving the feasibility of MRF in clinical settings.

Quantitative Relaxometry Metrics for Brain Metastases Compared to Normal Tissues: A Pilot MR Fingerprinting Study

The purpose of the present pilot study was to estimate T1 and T2 metric values derived simultaneously from a new, rapid Magnetic Resonance Fingerprinting (MRF) technique, as well as to assess their ability to characterize-brain metastases (BM) and normal-appearing brain tissues. Fourteen patients with BM underwent MRI, including prototype MRF, on a 3T scanner. In total, 108 measurements were analyzed: 42 from solid parts of BM's (21 each on T1 and T2 maps) and 66 from normal-appearing brain tissue (11 ROIs each on T1 and T2 maps for gray matter [GM], white matter [WM], and cerebrospinal fluid [CSF]). The BM's mean T1 and T2 values differed significantly from normal-appearing WM (p < 0.05). The mean T1 values from normal-appearing GM, WM, and CSF regions were 1205 ms, 840 ms, and 4233 ms, respectively. The mean T2 values were 108 ms, 78 ms, and 442 ms, respectively. The mean T1 and T2 values for untreated BM (n = 4) were 2035 ms and 168 ms, respectively. For treated BM (n = 17) the T1 and T2 values were 2163 ms and 141 ms, respectively. MRF technique appears to be a promising and rapid quantitative method for the characterization of free water content and tumor morphology in BMs.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Multicenter Repeatability and Reproducibility of MR Fingerprinting in Phantoms and in Prostatic Tissue

PURPOSE

To evaluate multicenter repeatability and reproducibility of T1 and T2 maps generated using MR fingerprinting (MRF) in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology MRI system phantom and in prostatic tissues.

METHODS

MRF experiments were performed on 5 different 3 Tesla MRI scanners at 3 different institutions: University Hospitals Cleveland Medical Center (Cleveland, OH), Brigham and Women's Hospital (Boston, MA) in the United States, and Diagnosticos da America (Rio de Janeiro, RJ) in Brazil. Raw MRF data were reconstructed using a Gadgetron-based MRF online reconstruction pipeline to yield quantitative T1 and T2 maps. The repeatability of T1 and T2 values over 6 measurements in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology MRI system phantom was assessed to demonstrate intrascanner variation. The reproducibility between the 4 clinical scanners was assessed to demonstrate interscanner variation. The same-day test-retest normal prostate mean T1 and T2 values from peripheral zone and transitional zone were also compared using the intraclass correlation coefficient and Bland-Altman analysis.

RESULTS

The intrascanner variation of values measured using MRF was less than 2% for T1 and 4.7% for T2 for relaxation values, within the range of 307.7 to 2360 ms for T1 and 19.1 to 248.5 ms for T2 . Interscanner measurements showed that the T1 variation was less than 4.9%, and T2 variation was less than 8.1% between multicenter scanners. Both T1 and T2 values in in vivo prostatic tissue demonstrated high test-retest reliability (intraclass correlation coefficient > 0.92) and strong linear correlation (R2  > 0.840).

CONCLUSION

Prostate MRF measurements of T1 and T2 are repeatable and reproducible between MRI scanners at different centers on different continents for the above measurement ranges.

Characterizing Thalamic and Basal Ganglia Nuclei in Medically Intractable Focal Epilepsy by MR Fingerprinting

OBJECTIVES

Magnetic resonance fingerprinting (MRF) is a novel, quantitative and noninvasive technique to measure brain tissue properties. We aim to use MRF for characterizing normal-appearing thalamic and basal ganglia nuclei in the epileptic brain.

METHODS

A 3D MRF protocol (1mm³ isotropic resolution) was acquired from 48 patients with unilateral medically refractory focal epilepsy and 39 healthy controls (HCs). Whole-brain T1 and T2 maps (containing T1 and T2 relaxation times) were reconstructed for each subject. Ten subcortical nuclei in the thalamus and basal ganglia were segmented as regions of interest (ROIs), within which the mean T1 and T2 values, as well as their coefficient of variation (CV) were compared between the patients and HCs at group level. Subgroup and correlation analyses were performed to examine the relationship between significant MRF measures and various clinical characteristics. Using significantly abnormal MRF measures from the group-level analyses, support vector machine (SVM) and logistic regression machine learning models were built and tested with 5-fold and 10-fold cross-validations, to separate patients from HCs, and to separate patients with left-sided and right-sided epilepsy, at individual level.

RESULTS

MRF revealed increased T1 mean value in the ipsilateral thalamus and nucleus accumbens; increased T1 CV in bilateral thalamus, bilateral pallidum, and ipsilateral caudate; and increased T2 CV in the ipsilateral thalamus in patients compared to HCs (P<0.05, FDR corrected). The SVM classifier produced 78.2% average accuracy to separate individual patients from HCs, with AUC of 0.83. The logistic regression classifier produced 67.4% average accuracy to separate patients with left-sided and right-sided epilepsy, with AUC of 0.72.

SIGNIFICANCE

MRF revealed bilateral tissue-property changes in the normal-appearing thalamus and basal ganglia, with ipsilateral predominance and thalamic preference, suggesting subcortical involvement/impairment in patients with medically intractable focal epilepsy. The individual-level performance of the MRF-based machine-learning models suggests potential opportunities for predicting lateralization.

Simultaneous 3D acquisition of 1 H MRF and 23 Na MRI

PURPOSE

To develop a 3D MR technique to simultaneously acquire proton multiparametric maps (T₁ , T₂ , and proton density) and sodium density weighted images over the whole brain.

METHODS

We implemented a 3D stack-of-stars MR pulse sequence which consists of interleaved proton (¹ H) and sodium (²³ Na) excitations, tailored slice encoding gradients that can encode the same slice for both nuclei, and simultaneous readout with different radial trajectories (¹ H, full-radial; ²³ Na, center-out radial). The receive chain of our 7T scanner was modified to enable simultaneous acquisition of ¹ H and ²³ Na signal. A heuristically optimized flip angle train was implemented for proton MR fingerprinting (MRF). The SNR and the accuracy of proton T₁ and T₂ were evaluated in phantoms. Finally, in vivo application of the method was demonstrated in five healthy subjects.

RESULTS

The SNR for the simultaneous measurement was almost identical to that for the single-nucleus measurements (<2% change). The proton T₁ and T₂ maps remained similar to the results from a reference 2D MRF technique (normalized RMS error in T₁ ≈ 4.2% and T₂ ≈ 11.3%). Measurements in healthy subjects corroborated these results and demonstrated the feasibility of our method for in vivo application. The in vivo T₁ values measured using our method were lower than the results measured by other conventional techniques.

CONCLUSIONS

With the 3D simultaneous implementation, we were able to acquire sodium and proton density weighted images in addition to proton T₁ , T₂ , and B1+ from ¹ H MRF that covers the whole brain volume within 21 min.

Low‐rank inversion reconstruction for through‐plane accelerated radial MR fingerprinting applied to relaxometry at 0. 35 T

Purpose

To reduce scan time, methods to accelerate phase‐encoded/non‐Cartesian MR fingerprinting (MRF) acquisitions for variable density spiral acquisitions have recently been developed. These methods are not applicable to MRF acquisitions, wherein a single k‐space spoke is acquired per frame. Therefore, we propose a low‐rank inversion method to resolve MRF contrast dynamics from through‐plane accelerated Cartesian/radial measurements applied to quantitative relaxation‐time mapping on a 0.35T system.

Methods

An algorithm was implemented to reconstruct through‐plane aliased low‐rank images describing the contrast dynamics occurring because of the transient‐state MRF acquisition. T₁ and T₂ times from accelerated acquisitions were compared with those from unaccelerated linear reconstructions in a standardized system phantom and within in vivo brain and prostate experiments on a hybrid 0.35T MRI/linear accelerator.

Results

No significant differences between T₁ and T₂ times for the accelerated reconstructions were observed compared to fully sampled acquisitions (p = 0.41 and p = 0.36, respectively). The mean absolute errors in T₁ and T₂ were 5.6% and 2.9%, respectively, between the full and accelerated acquisitions. The SDs in T₁ and T₂ decreased with the advanced accelerated reconstruction compared with the unaccelerated reconstruction (p = 0.02 and p = 0.03, respectively). The quality of the T₁ and T₂ maps generated with the proposed approach are comparable to those obtained using the unaccelerated data sets.

Conclusions

Through‐plane accelerated MRF with radial k‐space coverage was demonstrated at a low field strength of 0.35 T. This method enabled 3D T₁ and T₂ mapping at 0.35 T with a 3‐min scan.

Rigid real-time prospective motion-corrected three-dimensional multiparametric mapping of the human brain

PURPOSE

To develop a rigid real-time prospective motion-corrected multiparametric mapping technique and to test the performance of quantitative estimates.

METHODS

Motion tracking and correction were performed by integrating single-shot spiral navigators into a multiparametric imaging technique, three-dimensional quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS). The spiral navigator was optimized, and quantitative measurements were validated using a standard system phantom. The effect of motion correction on whole-brain T1 and T2 mapping under different types of head motion during the scan was evaluated in 10 healthy volunteers. Finally, six patients with Parkinson's disease, which is known to be associated with a high prevalence of motion artifacts, were scanned to evaluate the effectiveness of our method in the real world.

RESULTS

The phantom study demonstrated that the proposed motion correction method did not introduce quantitative bias. Improved parametric map quality and repeatability were shown in volunteer experiments with both in-plane and through-plane motions, comparable to the no-motion ground truth. In real-life validation in patients, the approach showed improved parametric map quality compared to images obtained without motion correction.

CONCLUSIONS

Real-time prospective motion-corrected multiparametric relaxometry based on 3D-QALAS provided robust and repeatable whole-brain multiparametric mapping.

Using magnetic resonance fingerprinting to characterize periventricular nodular heterotopias in pharmacoresistant epilepsy

OBJECTIVE

We aimed to use a novel magnetic resonance fingerprinting (MRF) technique to examine in vivo tissue property characteristics of periventricular nodular heterotopia (PVNH). These characteristics were further correlated with stereotactic-electroencephalographic (SEEG) ictal onset findings.

METHODS

We included five patients with PVNH who had SEEG-guided surgery and at least 1 year of seizure freedom or substantial seizure reduction. High-resolution MRF scans were acquired at 3 T, generating three-dimensional quantitative T₁ and T₂  maps. We assessed the differences between T₁ and T₂  values from the voxels in the nodules located in the SEEG-defined seizure onset zone (SOZ) and non-SOZ, on -individual and group levels. Receiver operating characteristic analyses were performed to obtain the optimal classification performance. Quantification of SEEG ictal onset signals from the nodules was performed by calculating power spectrum density (PSD). The association between PSD and T₁ /T₂  values was further assessed at different frequency bands.

RESULTS

Individual-level analysis showed T₁ was significantly higher in SOZ voxels than non-SOZ voxels (p < .05), with an average 73% classification accuracy. Group-level analysis also showed higher T₁ was significantly associated with SOZ voxels (p < .001). At the optimal cutoff (normalized T₁ of 1.1), a 76% accuracy for classifying SOZ nodules from non-SOZ nodules was achieved. T₁  values were significantly associated with ictal onset PSD at the ultraslow, θ, β, γ, and ripple bands (p < .05). T₂  values were significantly associated with PSD only at the ultraslow band (p < .05).

SIGNIFICANCE

Quantitative MRF measures, especially T₁ , can provide additional noninvasive information to separate nodules in SOZ and non-SOZ. The T₁ and T₂ tissue property changes carry electrophysiological underpinnings relevant to the epilepsy, as shown by their significant positive associations with power changes during the SEEG seizure onset. The use of MRF as a supplementary noninvasive tool may improve presurgical evaluation for patients with PVNH and pharmacoresistant epilepsy.

3D MR fingerprinting using Seiffert spirals

PURPOSE

Seiffert spirals were recently explored as an efficient way to traverse 3D k-space compared to traditional 3D techniques. Several studies have shown the ability of 3D MR fingerprinting (MRF) techniques to acquire T₁ and T₂ relaxation maps in a short period of time. However, these sequences do not sample across a large region of 3D k-space every TR, especially in the way that Seiffert trajectories can.

METHODS

A 3D MRF sequence was designed using 8 Seiffert spirals rotated in 3D k-space, with flip angle modulation for T₁ and T₂ sensitivity. The sequence was compared to an MRF sequence using a 2D spiral rotated in 3D k-space using the tiny golden angle acquisition with similar resolution/readout duration. Both sequences were evaluated using simulations, phantom validation, and in vivo imaging.

RESULTS

In all experiments, the Seiffert spiral MRF sequence performed similar to if not better than the multi-axis 2D spiral MRF sequence. Strong intraclass correlation coefficients (> 0.9) were found between conventional and MRF sequences in phantoms, whereas the in vivo results showed slightly less aliasing artifact with the Seiffert trajectory.

CONCLUSION

In this study, Seiffert spirals were used within the MRF framework to acquire high-resolution T₁ and T₂ relaxation time maps in less than 2.5 min. The reduced aliasing artifacts seen with the Seiffert sequence suggests that sampling over 3D k-space evenly each TR can improve quantification or shorten scan times.

Radiomics with 3-dimensional magnetic resonance fingerprinting: influence of dictionary design on repeatability and reproducibility of radiomic features

Objectives

We aimed to investigate the influence of magnetic resonance fingerprinting (MRF) dictionary design on radiomic features using in vivo human brain scans.

Methods

Scan-rescans of three-dimensional MRF and conventional T1-weighted imaging were performed on 21 healthy volunteers (9 males and 12 females; mean age, 41.3 ± 14.6 years; age range, 22–72 years). Five patients with multiple sclerosis (3 males and 2 females; mean age, 41.2 ± 7.3 years; age range, 32–53 years) were also included. MRF data were reconstructed using various dictionaries with different step sizes. First- and second-order radiomic features were extracted from each dataset. Intra-dictionary repeatability and inter-dictionary reproducibility were evaluated using intraclass correlation coefficients (ICCs). Features with ICCs > 0.90 were considered acceptable. Relative changes were calculated to assess inter-dictionary biases.

Results

The overall scan-rescan ICCs of MRF-based radiomics ranged from 0.86 to 0.95, depending on dictionary step size. No significant differences were observed in the overall scan-rescan repeatability of MRF-based radiomic features and conventional T1-weighted imaging (p = 1.00). Intra-dictionary repeatability was insensitive to dictionary step size differences. MRF-based radiomic features varied among dictionaries (overall ICC for inter-dictionary reproducibility, 0.62–0.99), especially when step sizes were large. First-order and gray level co-occurrence matrix features were the most reproducible feature classes among different step size dictionaries. T1 map-derived radiomic features provided higher repeatability and reproducibility among dictionaries than those obtained with T2 maps.

Conclusion

MRF-based radiomic features are highly repeatable in various dictionary step sizes. Caution is warranted when performing MRF-based radiomics using datasets containing maps generated from different dictionaries.

Key Points

• MRF-based radiomic features are highly repeatable in various dictionary step sizes.

• Use of different MRF dictionaries may result in variable radiomic features, even when the same MRF acquisition data are used.

• Caution is needed when performing radiomic analysis using data reconstructed from different dictionaries.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08555-3.

Simultaneous relaxometry and morphometry of human brain structures with 3D magnetic resonance fingerprinting: a multicenter, multiplatform, multifield-strength study

Relaxation times and morphological information are fundamental magnetic resonance imaging-derived metrics of the human brain that reflect the status of the underlying tissue. Magnetic resonance fingerprinting (MRF) enables simultaneous acquisition of T1 and T2 maps inherently aligned to the anatomy, allowing whole-brain relaxometry and morphometry in a single scan. In this study, we revealed the feasibility of 3D MRF for simultaneous brain structure-wise morphometry and relaxometry. Comprehensive test-retest scan analyses using five 1.5-T and three 3.0-T systems from a single vendor including different scanner types across 3 institutions demonstrated that 3D MRF-derived morphological information and relaxation times are highly repeatable at both 1.5 T and 3.0 T. Regional cortical thickness and subcortical volume values showed high agreement and low bias across different field strengths. The ability to acquire a set of regional T1, T2, thickness, and volume measurements of neuroanatomical structures with high repeatability and reproducibility facilitates the ability of longitudinal multicenter imaging studies to quantitatively monitor changes associated with underlying pathologies, disease progression, and treatments.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging

PURPOSE

To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF).

METHODS

Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B₀ inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect.

RESULTS

The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T₁ and T₂ mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures.

CONCLUSIONS

The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

MR fingerprinting for semisolid magnetization transfer and chemical exchange saturation transfer quantification

Chemical exchange saturation transfer (CEST) MRI has positioned itself as a promising contrast mechanism, capable of providing molecular information at sufficient resolution and amplified sensitivity. However, it has not yet become a routinely employed clinical technique, due to a variety of confounding factors affecting its contrast-weighted image interpretation and the inherently long scan time. CEST MR fingerprinting (MRF) is a novel approach for addressing these challenges, allowing simultaneous quantitation of several proton exchange parameters using rapid acquisition schemes. Recently, a number of deep-learning algorithms have been developed to further boost the performance and speed of CEST and semi-solid macromolecule magnetization transfer (MT) MRF. This review article describes the fundamental theory behind semisolid MT/CEST-MRF and its main applications. It then details supervised and unsupervised learning approaches for MRF image reconstruction and describes artificial intelligence (AI)-based pipelines for protocol optimization. Finally, practical considerations are discussed, and future perspectives are given, accompanied by basic demonstration code and data.

Three-dimensional high-resolution T1 and T2 mapping of whole macaque brain at 9.4 T using magnetic resonance fingerprinting

PURPOSE

Quantitative T₁ and T₂ mapping in non-human primates with whole-brain coverage is challenged by the requirement of sub-millimeter resolution and the inhomogeneity of the transmit magnetic field (B₁ ⁺ ) covering a large field of view. The goal of the current study is to develop a magnetic resonance fingerprinting (MRF) method for simultaneous T₁ and T₂ mapping of the entire macaque brain within feasible scan time.

METHODS

A three-dimensional (3D) MRF sequence with both inversion- and T₂ -preparation modules was developed and evaluated on a 9.4 T preclinical scanner. Data acquisition used a 3D stack-of-spirals trajectory, with undersampling along both the in-plane and the through-plane directions. The effect of B₁ ⁺ inhomogeneity was accounted for by matching the acquired fingerprint to a dictionary simulated with the B₁ ⁺ factors measured from a separate scan. In vitro and ex vivo studies were performed to evaluate the accuracy and the undersampling capacity of the MRF method. The application of the MRF method for in vivo, brain-wide T₁ and T₂ mapping was demonstrated on macaques at 4, 6, and 12 years of age.

RESULTS

The MRF method enabled highly repeatable T₁ and T₂ mapping at high spatial resolution (0.35 × 0.35 × 1 mm³ ) with an acceleration factor of 24. In vivo studies showed significant age-related T₂ reduction in deep gray nuclei including the globus pallidus, the putamen, and the caudate nucleus.

CONCLUSIONS

This study demonstrates the first MRF study for brain-wide, multi-parametric quantification in non-human primates with sub-millimeter resolution.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T₁, T₂, T₂*, PD and B₁⁺ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

Three-Dimensional Magnetic Resonance Fingerprinting in Neonates: Quantifying Regional Difference and Maturation in the Brain

OBJECTIVES

Magnetic resonance fingerprinting (MRF) allows the simultaneous measurement of multiple tissue properties in a single acquisition. Three-dimensional (3D) MRF with high spatial resolution can be used for neonatal brain imaging. The aim of this study is to apply 3D MRF to neonates and show regional differences and maturation in the brain.

MATERIALS AND METHODS

In this prospective study, 3D MRF using hybrid radial-interleaved acquisition was performed on phantoms and neonates from December 2019 to October 2020. For the reconstruction of 3D MRF, singular value decomposition was applied to reduce reconstruction time, and the iterative reconstruction technique was applied to improve image quality. The accuracies of T1 and T2 values derived from 3D MRF were evaluated in a phantom experiment. Regional T1 and T2 values were obtained from neonates' brain T1 and T2 maps derived from 3D MRF. Regional T1 and T2 values were compared, and their changes according to corrected gestational age were evaluated.

RESULTS

The acquisition time for 3D MRF with a spatial resolution of 0.7 × 0.7 × 2 mm3 was less than 5 minutes. The phantom study showed high correlation between T1 and T2 values derived from 3D MRF and those from conventional spin echo sequences (T1, R2 = 0.998, P < 0.001; T2, R2 = 0.998, P < 0.001). Three-dimensional MRF was performed in 25 neonates (15 boys, 10 girls; median corrected gestational age, 263 days; interquartile range, 10 days). In neonates, T1 and T2 values differed in the frontal (median [interquartile range], 2785 [2684-2888] milliseconds and 189.8 [176.7-222.9] milliseconds), parietal (2849 [2741-2950] milliseconds and 191.6 [167.5-232.9] milliseconds), and occipital white matter (2621 [2513-2722] milliseconds and 162.9 [143.5-186.1] milliseconds), showing lower values in occipital white matter (P < 0.001). Regional T1 values showed a negative relationship with corrected gestational age (coefficient, -0.775 to -0.480; P < 0.05).

CONCLUSIONS

Fast and high spatial resolution 3D MRF was applied to neonates. T1 and T2 maps derived from 3D MRF enabled the quantification of regional differences and maturation in the neonatal brain.

Rapid In Vivo Quantification of Creatine Kinase Activity by Phosphorous-31 Magnetic Resonance Spectroscopic Fingerprinting (31P-MRSF)

Creatine kinase (CK) plays an important role in tissue metabolism by providing a buffering mechanism for maintaining a constant supply of adenosine triphosphate (ATP) during metabolic perturbations. Phosphorous-31 magnetic resonance spectroscopy (31P-MRS) employing magnetization transfer techniques is the only noninvasive method for measuring the rate of ATP synthesis via creatine kinase. However, due to the low concentrations of phosphate metabolites, current 31P-MRS methods require long acquisition time to achieve adequate measurement accuracy. In this chapter, we present a new framework of data acquisition and parameter estimation, the 31P magnetic resonance spectroscopic fingerprinting (31P-MRSF) method, for rapid quantification of CK reaction rate constant in the hindlimb of small laboratory animals.

Magnetic Resonance Fingerprinting of the Pediatric Brain

Magnetic resonance fingerprinting (MRF) is increasingly being used to evaluate brain development and differentiate normal and pathologic tissues in children. MRF can provide reliable and accurate intrinsic tissue properties, such as T1 and T2 relaxation times. MRF is a powerful tool in evaluating brain disease in pediatric population. MRF is a new quantitative MR imaging technique for rapid and simultaneous quantification of multiple tissue properties.

Technical overview of magnetic resonance fingerprinting and its applications in radiation therapy

Magnetic resonance fingerprinting (MRF) is an emerging imaging technique for rapid and simultaneous quantification of multiple tissue properties. The technique has been developed for quantitative imaging of different organs. The obtained quantitative measures have the potential to improve multiple steps of a typical radiotherapy workflow and potentially further improve integration of magnetic resonance imaging guided clinical decision making. In this review paper, we first provide a technical overview of the MRF method from data acquisition to postprocessing, along with recent development in advanced reconstruction methods. We further discuss critical aspects that could influence its usage in radiation therapy, such as accuracy and precision, repeatability and reproducibility, geometric distortion, and motion robustness. Finally, future directions for MRF application in radiation therapy are discussed.

Automated design of pulse sequences for magnetic resonance fingerprinting using physics-inspired optimization

Magnetic resonance fingerprinting (MRF) is a method to extract quantitative tissue properties such as [Formula: see text] and [Formula: see text] relaxation rates from arbitrary pulse sequences using conventional MRI hardware. MRF pulse sequences have thousands of tunable parameters, which can be chosen to maximize precision and minimize scan time. Here, we perform de novo automated design of MRF pulse sequences by applying physics-inspired optimization heuristics. Our experimental data suggest that systematic errors dominate over random errors in MRF scans under clinically relevant conditions of high undersampling. Thus, in contrast to prior optimization efforts, which focused on statistical error models, we use a cost function based on explicit first-principles simulation of systematic errors arising from Fourier undersampling and phase variation. The resulting pulse sequences display features qualitatively different from previously used MRF pulse sequences and achieve fourfold shorter scan time than prior human-designed sequences of equivalent precision in [Formula: see text] and [Formula: see text] Furthermore, the optimization algorithm has discovered the existence of MRF pulse sequences with intrinsic robustness against shading artifacts due to phase variation.

Current Applications and Future Development of Magnetic Resonance Fingerprinting in Diagnosis, Characterization, and Response Monitoring in Cancer

Magnetic resonance imaging (MRI) has enabled non-invasive cancer diagnosis, monitoring, and management in common clinical settings. However, inadequate quantitative analyses in MRI continue to limit its full potential and these often have an impact on clinicians' judgments. Magnetic resonance fingerprinting (MRF) has recently been introduced to acquire multiple quantitative parameters simultaneously in a reasonable timeframe. Initial retrospective studies have demonstrated the feasibility of using MRF for various cancer characterizations. Further trials with larger cohorts are still needed to explore the repeatability and reproducibility of the data acquired by MRF. At the moment, technical difficulties such as undesirable processing time or lack of motion robustness are limiting further implementations of MRF in clinical oncology. This review summarises the latest findings and technology developments for the use of MRF in cancer management and suggests possible future implications of MRF in characterizing tumour heterogeneity and response assessment.

Quantitative T1 and T2 mapping by magnetic resonance fingerprinting (MRF) of the placenta before and after maternal hyperoxia

INTRODUCTION

MR relaxometry has been used to assess placental exchange function, but methods to date are not sufficiently fast to be robust to placental motion. Magnetic resonance fingerprinting (MRF) permits rapid, voxel-wise, intrinsically co-registered T₁ and T₂ mapping. After characterizing measurement error, we scanned pregnant women during air and oxygen breathing to demonstrate MRF's ability to detect placental oxygenation changes.

METHODS

The accuracy of FISP-based, sliding-window reconstructed MRF was tested on phantoms. MRF scans in 9-s breath holds were acquired at 3T in 31 pregnant women during air and oxygen breathing. A mixed effects model was used to test for changes in placenta relaxation times between physiological states, to assess the dependency on gestational age (GA), and the impact of placental motion.

RESULTS

MRF estimates of known phantom relaxation times resulted in mean absolute errors for T₁ of 92 ms (4.8%), but T₂ was less accurate at 16 ms (13.6%). During normoxia, placental T₁ = 1825 ± 141 ms (avg ± standard deviation) and T₂ = 60 ± 16 ms (gestational age range 24.3-36.7, median 32.6 weeks). In the statistical model, placental T₂ rose and T₁ remained contant after hyperoxia, and no GA dependency was observed for T₁ or T₂.

DISCUSSION

Well-characterized, motion-robust MRF was used to acquire T₁ and T₂ maps of the placenta. Changes with hyperoxia are consistent with a net increase in oxygen saturation. Toward the goal of whole-placenta quantitative oxygenation imaging over time, we aim to implement 3D MRF with integrated motion correction to improve T₂ accuracy.

Efficiency analysis for quantitative MRI of T1 and T2 relaxometry methods

This study presents a comparison of quantitative MRI methods based on an efficiency metric that quantifies their intrinsic ability to extract information about tissue parameters. Under a regime of unbiased parameter estimates, an intrinsic efficiency metricηwas derived for fully-sampled experiments which can be used to both optimize and compare sequences. Here we optimize and compare several steady-state and transient gradient-echo based qMRI methods, such as magnetic resonance fingerprinting (MRF), for jointT1andT2mapping. The impact of undersampling was also evaluated, assuming incoherent aliasing that is treated as noise by parameter estimation.In vivovalidation of the efficiency metric was also performed. Transient methods such as MRF can be up to 3.5 times more efficient than steady-state methods, when spatial undersampling is ignored. If incoherent aliasing is treated as noise during least-squares parameter estimation, the efficiency is reduced in proportion to the SNR of the data, with reduction factors of 5 often seen for practical SNR levels.In vivovalidation showed a very good agreement between the theoretical and experimentally predicted efficiency. This work presents and validates an efficiency metric to optimize and compare the performance of qMRI methods. Transient methods were found to be intrinsically more efficient than steady-state methods, however the effect of spatial undersampling can significantly erode this advantage.

Water-Fat Separation in MR Fingerprinting for Quantitative Monitoring of the Skeletal Muscle in Neuromuscular Disorders

Background Quantitative MRI is increasingly proposed in clinical trials related to neuromuscular disorders (NMDs). Purpose To investigate the potential of an MR fingerprinting sequence for water and fat fraction (FF) quantification (MRF T1-FF) for providing markers of fatty replacement and disease activity in patients with NMDs and to establish the sensitivity of water T1 as a marker of disease activity compared with water T2 mapping. Materials and Methods Data acquired between March 2018 and March 2020 from the legs of patients with NMDs were retrospectively analyzed. The MRI examination comprised fat-suppressed T2-weighted imaging, mapping of the FF measured with the three-point Dixon technique (FF_Dixon), water T2 mapping, and MRF T1-FF, from which the FF measured with MRF T1-FF (FF_MRF) and water T1 were derived. Data from the legs of healthy volunteers were prospectively acquired between January and July 2020 to derive abnormality thresholds for FF, water T2, and water T1 values. Kruskal-Wallis tests and receiver operating characteristic curve analysis were performed, and linear models were used. Results A total of 73 patients (mean age ± standard deviation, 47 years ± 12; 45 women) and 15 healthy volunteers (mean age, 33 years ± 8; three women) were evaluated. A linear correlation was observed between FF_MRF and FF_Dixon (R² = 0.97, P < .001). Water T1 values were higher in muscles with high signal intensity at fat-suppressed T2-weighted imaging than in muscles with low signal intensity (mean value, 1281 msec [95% CI: 1165, 1604] vs 1198 msec [95% CI: 1099, 1312], respectively; P < .001), and a correlation was found between water T1 and water T2 distribution metrics (R² = 0.66 and 0.79 for the median and 90th percentile values, respectively; P < .001). Water T1 classified the patients' muscles as abnormal based on quantitative water T2, with high sensitivity (93%; 68 of 73 patients) and specificity (80%; 53 of 73 patients) (area under the receiver operating characteristic curve, 0.92 [95% CI: 0.83, 0.97]; P < .001). Conclusion Water-fat separation in MR fingerprinting is robust for deriving quantitative imaging markers of intramuscular fatty replacement and disease activity in patients with neuromuscular disorders. © RSNA, 2021 Online supplemental material is available for this article.

Feasibility of MR fingerprinting using a high-performance 0.55 T MRI system

BACKGROUND

MR fingerprinting (MRF) is a versatile method for rapid multi-parametric quantification. The application of MRF for lower MRI field could enable multi-contrast imaging and improve exam efficiency on these systems. The purpose of this work is to demonstrate the feasibility of 3D whole-brain T1 and T2 mapping using MR fingerprinting on a contemporary 0.55 T MRI system.

MATERIALS AND METHODS

A 3D whole brain stack-of-spirals FISP MRF sequence was implemented for 0.55 T. Quantification was validated using the NIST/ISMRM Quantitative MRI phantom, and T1 and T2 values of white matter, gray matter, and cerebrospinal fluid were measured in 19 healthy subjects. To assess MRF performance in the lower SNR regime of 0.55 T, measurement precision was calculated from 100 simulated pseudo-replicas of in vivo data and within-session measurement repeatability was evaluated.

RESULTS

T1 and T2 values calculated by MRF were strongly correlated to standard measurements in the ISMRM/NIST MRI system phantom (R² > 0.99), with a small constant bias of approximately 5 ms in T2 values. 3D stack-of-spirals MRF was successfully applied for whole brain quantitative T1 and T2 at 0.55 T, with spatial resolution of 1.2 mm × 1.2 mm × 5 mm, and acquisition time of 8.5 min. Moreover, the T1 and T2 quantifications had precision <5%, despite the lower SNR of 0.55 T.

CONCLUSION

A 3D whole-brain stack-of-spirals FISP MRF sequence is feasible for T1 and T2 mapping at 0.55 T.

Rapid B1-Insensitive MR Fingerprinting for Quantitative Kidney Imaging

Background MR fingerprinting (MRF) provides rapid and simultaneous quantification of multiple tissue parameters in a single scan. Purpose To evaluate a rapid kidney MRF technique at 3.0 T in phantoms, healthy volunteers, and patients. Materials and Methods A 15-second kidney MRF acquisition was designed with 12 acquisition segments, a range of low flip angles (5°-12°), multiple magnetization preparation schema (T1, T2, and fat suppression), and an undersampled spiral trajectory. This technique was first validated in vitro using standardized T1 and T2 phantoms. Kidney T1 and T2 maps were then obtained for 10 healthy adult volunteers (mean age ± standard deviation, 35 years ± 13; six men) and three pediatric patients with autosomal recessive polycystic kidney disease (ARPKD) (mean age, 10 years ± 3; two boys) between August 2019 and October 2020 to evaluate the method in vivo. Results Results in nine phantoms showed good agreement with spin-echo-based T1 and T2 values (R² > 0.99). In vivo MRF kidney T1 and T2 assessments in healthy adult volunteers (cortex: T1, 1362 msec ± 5; T2, 64 msec ± 5; medulla: T1, 1827 msec ± 94; T2, 69 msec ± 3) were consistent with values in the literature but with improved precision in comparison with prior MRF implementations. In vivo MRF-based kidney T1 and T2 values with and without B₁ correction were in good agreement (R² > 0.96, P < .001), demonstrating limited sensitivity to B₁ field inhomogeneities. Additional MRF reconstructions using the first nine segments of the MRF profiles (11-second acquisition time) were in good agreement with the reconstructions using 12 segments (15-second acquisition time) (R² > 0.87, P < .001). Repeat kidney MRF scans for the three patients with ARPKD on successive days also demonstrated good reproducibility (T1 and T2: <3% difference). Conclusion A kidney MR fingerprinting method provided in vivo kidney T1 and T2 maps at 3.0 T in a single breath hold with improved precision and no need for B₁ correction. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Laustsen in this issue.

Rapid high-resolution volumetric T1 mapping using a highly accelerated stack-of-stars Look Locker technique

PURPOSE

To develop a fast volumetric T1 mapping technique.

MATERIALS AND METHODS

A stack-of-stars (SOS) Look Locker technique based on the acquisition of undersampled radial data (>30× relative to Nyquist) and an efficient multi-slab excitation scheme is presented. A principal-component based reconstruction is used to reconstruct T1 maps. Computer simulations were performed to determine the best choice of partitions per slab and degree of undersampling. The technique was validated in phantoms against reference T1 values measured with a 2D Cartesian inversion-recovery spin-echo technique. The SOS Look Locker technique was tested in brain (n = 4) and prostate (n = 5). Brain T1 mapping was carried out with and without kz acceleration and results between the two approaches were compared. Prostate T1 mapping was compared to standard techniques. A reproducibility study was conducted in brain and prostate. Statistical analyses were performed using linear regression and Bland Altman analysis.

RESULTS

Phantom T1 values showed excellent correlations between SOS Look Locker and the inversion-recovery spin-echo reference (r2 = 0.9965; p < 0.0001) and between SOS Look Locker with slab-selective and non-slab selective inversion pulses (r2 = 0.9999; p < 0.0001). In vivo results showed that full brain T1 mapping (1 mm3) with kz acceleration is achieved in 4 min 21 s. Full prostate T1 mapping (0.9 × 0.9 × 4 mm3) is achieved in 2 min 43 s. T1 values for brain and prostate were in agreement with literature values. A reproducibility study showed coefficients of variation in the range of 0.18-0.2% (brain) and 0.15-0.18% (prostate).

CONCLUSION

A rapid volumetric T1 mapping technique was developed. The technique enables high-resolution T1 mapping with adequate anatomical coverage in a clinically acceptable time.