Simultaneous variable flip angle-actual flip angle imaging method for improved accuracy and precision of three-dimensional T1 and B1 measurements

Accurate actual flip angle imaging (AFI) in the presence of fat

PURPOSE

To investigate the effect of incomplete fat spoiling on the accuracy of B1 mapping with actual flip angle imaging (AFI) and to propose a method to minimize the errors using the chemical shift properties of fat.

THEORY AND METHODS

Diffusion-based dephasing is the main spoiling mechanism exploited in AFI. However, a very low diffusion in fat may make the spoiling insufficient, leading to ghosts in the B1 maps. As the errors retain the chemical-shift signature of fat, their impact can be minimized using chemical-shift-based fat signal removal from AFI acquisition modified to include multi-echo readout. The source of the errors and the proposed correction were studied in simulations and phantom and in-vivo imaging experiments.

RESULTS

Our results support that AFI artifacts are caused by the incomplete fat spoiling present in clinically attractive short TR acquisition regimes. The correction eliminated the ghosting and significantly improved the B1 mapping accuracy as well as the accuracy of R1 mapping performed with AFI-derived B1 maps.

CONCLUSIONS

The incomplete fat signal spoiling may be a source of AFI B1 mapping errors, especially in subjects with high fat content. Achieving complete fat spoiling requires longer TR, which is undesirable in clinical applications. The proposed approach based on fat signal removal can reduce errors without significant prolongation of the AFI pulse sequence. We propose that, when attaining complete fat spoiling is not feasible, AFI mapping should be performed in a multi-echo regime with appropriate fat separation or suppression to minimize these errors.

A preparation pulse for fast steady state approach in Actual Flip angle Imaging

BACKGROUND

Actual Flip angle Imaging (AFI) is a sequence used for B1 mapping, also embedded in the Variable flip angle with AFI for simultaneous estimation of T1 , B1 and equilibrium magnetization.

PURPOSE

To investigate the design of a preparation module for AFI to allow a fast approach to steady state (SS) without requiring the use of dummy acquisitions.

METHODS

The features of a preparation module with a B1 insensitive adiabatic pulse, spoiler gradients, and a recovery timeT r e c $T_{rec}$ were studied with simulations and validated via experiments and acquired with different k-space traveling strategies. The robustness of the flip angle of the preparation pulse on the acquired signal is studied.

RESULTS

When a 90° adiabatic pulse is used, the forthcomingT r e c $T_{rec}$ can be expressed as a function of repetition times and AFI flip angle only asTR 1 ( n + cos α ) / ( 1 - cos 2 α ) $\mathrm{TR_1}(n+\cos \alpha )/(1-\cos ^2\alpha )$ , where n represents the ratio between the two repetition times of AFI. The robustness of the method is demonstrated by showing that using the values further away from 90° still allows for a faster approach to SS than the use of dummy pulses.

CONCLUSIONS

The preparation module is particularly advantageous for low flip angles, as well as for AFI sequences that sample the center of the k-space early in the sequence, such as centric ordering acquisitions, and for ultrafast EPI-based AFI methods, thus allowing to reduce scanner overhead time.

A relaxometry method that emphasizes practicality and availability

PURPOSE

Although many methods have been proposed to quantitatively map the main MRI parameters (e.g., T₁ , T₂ , C × M₀ ), these methods often involve special sequences not readily available on clinical scanners and/or may require long scan times. In contrast, the proposed method can readily run on most scanners, offer flexible tradeoffs between scan time and image quality, and map MRI parameters jointly to ensure spatial alignment.

METHODS

The approach is based on the multi-shot spin-echo (SE) EPI sequence. The corresponding signal equation was derived and strategies for solving it were developed. As usual with multi-shot EPI, scan time can readily be traded-off against image quality by adjusting the echo train length. Validation was performed against reference relaxometry methods, in gel phantoms with varying concentrations of gadobutrol and gadoterate meglumine contrast agents. In vivo examples are further presented, from 3 neuroradiology patients.

RESULTS

Bland-Altman analysis was performed: for T₂ , as compared to 2D SE, bias was 0.29 ms and the 95% limits of agreement ranged from -1.15 to +1.73 ms. For T₁ , compared to inversion-recovery SE (and MOLLI), bias was -20.2 ms (and -14.5 ms) and the limits of agreement ranged from -62.4 to +22.0 ms (and -53.8 to +24.9 ms). The mean relative T₁ error between the proposed method and each of the 2 reference methods was similar to that of the reference methods among themselves.

CONCLUSION

In the constellation of existing relaxometry methods, the proposed method is meant to stand out in terms of its practicality and availability.

Model-based multi-parameter mapping

Quantitative MR imaging is increasingly favoured for its richer information content and standardised measures. However, computing quantitative parameter maps, such as those encoding longitudinal relaxation rate (R1), apparent transverse relaxation rate (R2*) or magnetisation-transfer saturation (MTsat), involves inverting a highly non-linear function. Many methods for deriving parameter maps assume perfect measurements and do not consider how noise is propagated through the estimation procedure, resulting in needlessly noisy maps. Instead, we propose a probabilistic generative (forward) model of the entire dataset, which is formulated and inverted to jointly recover (log) parameter maps with a well-defined probabilistic interpretation (e.g., maximum likelihood or maximum a posteriori). The second order optimisation we propose for model fitting achieves rapid and stable convergence thanks to a novel approximate Hessian. We demonstrate the utility of our flexible framework in the context of recovering more accurate maps from data acquired using the popular multi-parameter mapping protocol. We also show how to incorporate a joint total variation prior to further decrease the noise in the maps, noting that the probabilistic formulation allows the uncertainty on the recovered parameter maps to be estimated. Our implementation uses a PyTorch backend and benefits from GPU acceleration. It is available at https://github.com/balbasty/nitorch.

Toward magnetic resonance fingerprinting for low-field MR-guided radiation therapy

PURPOSE

The acquisition of multiparametric quantitative magnetic resonance imaging (qMRI) is becoming increasingly important for functional characterization of cancer prior to- and throughout the course of radiation therapy. The feasibility of a qMRI method known as magnetic resonance fingerprinting (MRF) for rapid T₁ and T₂ mapping was assessed on a low-field MR-linac system.

METHODS

A three-dimensional MRF sequence was implemented on a 0.35T MR-guided radiotherapy system. MRF-derived measurements of T₁ and T₂ were compared to those obtained with gold standard single spin echo methods, and the impacts of the radiofrequency field homogeneity and scan times ranging between 6 and 48 min were analyzed by acquiring between 1 and 8 spokes per time point in a standard quantitative system phantom. The short-term repeatability of MRF was assessed over three measurements taken over a 10-h period. To evaluate transferability, MRF measurements were acquired on two additional MR-guided radiotherapy systems. Preliminary human volunteer studies were performed.

RESULTS

The phantom benchmarking studies showed that MRF is capable of mapping T₁ and T₂ values within 8% and 10% of gold standard measures, respectively, at 0.35T. The coefficient of variation of T₁ and T₂ estimates over three repeated scans was < 5% over a broad range of relaxation times. The T₁ and T₂ times derived using a single-spoke MRF acquisition across three scanners were near unity and mean percent errors in T₁ and T₂ estimates using the same phantom were < 3%. The mean percent differences in T₁ and T₂ as a result of truncating the scan time to 6 min over the large range of relaxation times in the system phantom were 0.65% and 4.05%, respectively.

CONCLUSIONS

The technical feasibility and accuracy of MRF on a low-field MR-guided radiation therapy device has been demonstrated. MRF can be used to measure accurate T₁ and T₂ maps in three dimensions from a brief 6-min scan, offering strong potential for efficient and reproducible qMRI for future clinical trials in functional plan adaptation and tumor/normal tissue response assessment.

Fast T1 measurement of cortical bone using 3D UTE actual flip angle imaging and single-TR acquisition (3D UTE-AFI-STR)

PURPOSE

To describe a new method for accurate T₁ measurement of cortical bone that fits the data sets of both 3D UTE actual flip angle imaging (UTE-AFI) and UTE with a single TR (UTE-STR) simultaneously (UTE-AFI-STR).

THEORY AND METHODS

To make both the constant values and longitudinal mapping functions in the signal equations for UTE-AFI and UTE-STR identical, the same RF pulses and flip angles were used. Therefore, there were three unknowns in the three equations. This was sufficient to fit the data. Numerical simulation as well as ex vivo and in vivo cortical bone studies were performed to validate the T₁ measurement accuracy with the UTE-AFI-STR method. The original UTE-AFI variable TR (VTR) (ie, combined UTE-AFI and UTE with VTR) and simultaneous fitting (sf) of UTE-AFI and UTE-VTR (sf-UTE-AFI-VTR) methods were performed for comparison.

RESULTS

The numerical simulation study showed that the UTE-AFI-STR method provided accurate value of T₁ when the SNR of the UTE-STR image was higher than 40. The ex vivo study showed that the UTE-AFI-STR method measured the T₁ of cortical bone accurately, with difference ratios ranging from -5.0% to 0.4%. The in vivo study showed a mean T₁ of 246 ms with the UTE-AFI-STR method, and mean difference ratios of 2.4% and 5.0%, respectively, compared with the other two methods.

CONCLUSION

The 3D UTE-AFI-STR method provides accurate mapping of the T₁ of cortical bone with improved time efficiency compared with the UTE-AFI-VTR/sf-UTE-AFI-VTR methods.

A practical method for post-acquisition reduction of bias in fast, whole-brain B1-maps

Large consistent differences have been observed between maps of the flip angle correction factor (commonly called "B₁-maps") produced with different fast methods in the human brain. We present an empirical procedure for first-order multiplicative bias correction that can be applied when more than one B₁-mapping method is available. We use a B₁-map measurement in a calibration phantom as a reference and the voxel-wise histogram mode between ratios of B₁-maps produced from different methods to calculate determine the bias as a multiplicative correcting scale factor. Institutional implementations of four common methods of B₁-mapping were assessed: Method of Slopes, FSE and EPI double angle methods (DAM), and Bloch-Siegert. In human subjects, the multiplicative bias used to correct for each of the four methods was: Method of Slopes = 1.005, FSE-DAM = 0.956, EPI-DAM = 1.080, and Bloch-Siegert = 1.128. Scaling to remove this bias between methods produces more consistent B₁-maps which enable more consistent values for any computations requiring flip angle correction. In addition, we present evidence that the corrected B₁ maps, using our calibration method, are also more accurate.

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

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

Magnetization transfer in magnetic resonance fingerprinting

PURPOSE

To study the effects of magnetization transfer (MT, in which a semi-solid spin pool interacts with the free pool), in the context of magnetic resonance fingerprinting (MRF).

METHODS

Simulations and phantom experiments were performed to study the impact of MT on the MRF signal and its potential influence on T1 and T2 estimation. Subsequently, an MRF sequence implementing off-resonance MT pulses and a dictionary with an MT dimension, generated by incorporating a two-pool model, were used to estimate the fractional pool size in addition to the B 1 + , T1 , and T2 values. The proposed method was evaluated in the human brain.

RESULTS

Simulations and phantom experiments showed that an MRF signal obtained from a cross-linked bovine serum sample is influenced by MT. Using a dictionary based on an MT model, a better match between simulations and acquired MR signals can be obtained (NRMSE 1.3% vs. 4.7%). Adding off-resonance MT pulses can improve the differentiation of MT from T1 and T2 . In vivo results showed that MT affects the MRF signals from white matter (fractional pool-size ~16%) and gray matter (fractional pool-size ~10%). Furthermore, longer T1 (~1060 ms vs. ~860 ms) and T2 values (~47 ms vs. ~35 ms) can be observed in white matter if MT is accounted for.

CONCLUSION

Our experiments demonstrated a potential influence of MT on the quantification of T1 and T2 with MRF. A model that encompasses MT effects can improve the accuracy of estimated relaxation parameters and allows quantification of the fractional pool size.

Evidence That Cannabis Exposure, Abuse, and Dependence Are Related to Glutamate Metabolism and Glial Function in the Anterior Cingulate Cortex: A 1H-Magnetic Resonance Spectroscopy Study

There is evidence that long-term cannabis use is associated with alterations to glutamate neurotransmission and glial function. In this study, 26 long-term cannabis users (males=65.4%) and 47 non-cannabis using healthy controls (males=44.6%) underwent proton magnetic resonance spectroscopy (¹H-MRS) of the anterior cingulate cortex (ACC) in order to characterize neurometabolite alterations in cannabis users and to examine associations between neurometabolites, cannabis exposure, and cannabis use behaviors. Myo-inositol, a marker of glial function, and glutamate metabolites did not differ between healthy controls and cannabis users or cannabis users who met criteria for DSM5 cannabis use disorder (n=17). Lower myo-inositol, a putative marker of glial function, was related to greater problematic drug use (F_1,22 = 11.95, p=.002; Cohen's f=0.59, large effect; Drug Abuse Screening Test) and severity of cannabis dependence (F_1,22 = 6.61, p=.17; Cohen's f=0.44, large effect). Further, past-year cannabis exposure exerted different effects on glutamate and glutamate+glutamine in males and females (glutamate: F_1,21 = 6.31, p=.02; glutamate+glutamine: F_1,21 = 7.20, p=.014), such that greater past-year cannabis exposure was related to higher concentrations of glutamate metabolites in male cannabis users (glutamate: F_1,14 = 25.94, p=.00016; Cohen's f=1.32, large effect; glutamate+glutamine: F_1,14 = 23.24, p=.00027, Cohen's f=1.24, large effect) but not in female cannabis users (glutamate: F_1,6 = 1.37, p=0.78; glutamate+glutamine: F_1,6 = 0.001, p=.97). The present results extend existing evidence of altered glial function and glutamate metabolism with cannabis use by providing evidence linking problematic drug use behaviors with glial function as measured with myo-inositol and recent chronic cannabis exposure to alterations in glutamate metabolism. This provides novel directions for the interrogation of the impact of cannabis use on brain neurochemistry.

Measurements of cerebral blood volume using quantitative susceptibility mapping, R2 * relaxometry, and ferumoxytol-enhanced MRI

Ferumoxytol-enhanced MRI holds potential for the non-invasive assessment of vascular architecture using estimates of cerebral blood volume (CBV). Ferumoxytol specifically enables steady-state imaging with extended acquisition times, for substantial improvements in resolution and contrast-to-noise ratio. With such data, quantitative susceptibility mapping (QSM) can be used to obtain images of local tissue magnetic susceptibility and hence estimate the increase in blood susceptibility after administration of a contrast agent, which in turn can be correlated to tissue CBV. Here, we explore the use of QSM for CBV estimation and compare it with R₂ * (1/T₂ *)-based results. Institutional review board approval was obtained, and all subjects provided written informed consent. For this prospective study, MR images were acquired on a 3.0 T scanner in 19 healthy subjects using a multiple-echo T₂ *-weighted sequence. Scanning was performed before and after the administration of two doses of ferumoxytol (1 mg FE/kg and 4 mg FE/kg). Different QSM approaches were tested on numerical phantom simulations. Results showed that the accuracy of magnetic susceptibility measurements improved with increasing image resolution and decreasing vascular density. In vivo changes in magnetic susceptibility were measured after the administration of ferumoxytol utilizing QSM, and significantly higher QSM-based CBV was measured in gray matter compared with white matter. QSM- and R₂ *-based CBV estimates correlated well, with similar average values, but a larger variance was found in QSM-based estimates.

Multi-pathway multi-echo acquisition and neural contrast translation to generate a variety of quantitative and qualitative image contrasts

PURPOSE

Clinical exams typically involve acquiring many different image contrasts to help discriminate healthy from diseased states. Ideally, 3D quantitative maps of all of the main MR parameters would be obtained for improved tissue characterization. Using data from a 7-min whole-brain multi-pathway multi-echo (MPME) scan, we aimed to synthesize several 3D quantitative maps (T₁ and T₂ ) and qualitative contrasts (MPRAGE, FLAIR, T₁ -weighted, T₂ -weighted, and proton density [PD]-weighted). The ability of MPME acquisitions to capture large amounts of information in a relatively short amount of time suggests it may help reduce the duration of neuro MR exams.

METHODS

Eight healthy volunteers were imaged at 3.0T using a 3D isotropic (1.2 mm) MPME sequence. Spin-echo, MPRAGE, and FLAIR scans were performed for training and validation. MPME signals were interpreted through neural networks for predictions of different quantitative and qualitative contrasts. Predictions were compared to reference values at voxel and region-of-interest levels.

RESULTS

Mean absolute errors (MAEs) for T₁ and T₂ maps were 216 ms and 11 ms, respectively. In ROIs containing white matter (WM) and thalamus tissues, the mean T₁ /T₂ predicted values were 899/62 ms and 1139/58 ms, consistent with reference values of 850/66 ms and 1126/58 ms, respectively. For qualitative contrasts, signals were normalized to those of WM, and MAEs for MPRAGE, FLAIR, T₁ -weighted, T₂ -weighted, and PD-weighted contrasts were 0.14, 0.15, 0.13, 0.16, and 0.05, respectively.

CONCLUSIONS

Using an MPME sequence and neural-network contrast translation, whole-brain results were obtained with a variety of quantitative and qualitative contrast in ~6.8 min.

Four-angle method for practical ultra-high-resolution magnetic resonance mapping of brain longitudinal relaxation time and apparent proton density

Changes in longitudinal relaxation time (T₁) and proton density (PD) are sensitive indicators of microstructural alterations associated with various central nervous system diseases as well as brain maturation and aging. In this work, we introduce a new approach for rapid and accurate high-resolution (HR) or ultra HR (UHR) mapping of T₁ and apparent PD (APD) of the brain with correction of radiofrequency field, B₁, inhomogeneities. The four-angle method (FAM) uses four spoiled-gradient recalled-echo (SPGR) images acquired at different flip angles (FA) and short repetition times (TRs). The first two SPGR images are acquired at low-spatial resolution and used to accurately map the active B₁⁺ field with the recently introduced steady-state double angle method (SS-DAM). The estimated B₁⁺ map is used in conjunction with the two other SPGR images, acquired at HR or UHR, to map T₁ and APD. The method is evaluated with numerical, phantom, and in-vivo imaging measurements. Furthermore, we investigated imaging acceleration methods to further shorten the acquisition time. Our results indicate that FAM provides an accurate method for simultaneous HR or UHR mapping of T₁ and APD in human brain in clinical high-field MRI. Derived parameter maps without B₁⁺correction suffer from large inaccuracies, but this issue is well-corrected through use of the SS-DAM. Furthermore, the use of SPGR imaging with short TR and phased-array coil acquisition permits substantial imaging acceleration and enables robust HR or UHR T₁ and APD mapping in a clinically acceptable time frame, with whole brain coverage obtained in less than 2 min or 5 min, respectively. The method exhibits high reproducibility and benefits from the use of the conventional SPGR sequence, available in all preclinical and clinical MRI machines, and very simple modeling to address a critical outstanding issue in neuroimaging.

More than BOLD: Dual-spin populations create functional contrast

PURPOSE

Functional MRI contrast has generally been associated with changes in transverse relaxivity caused by blood oxygen concentration, the so-called blood oxygen level dependent contrast. However, this interpretation of fMRI contrast has been called into question by several recent experiments at high spatial resolution. Experiments were conducted to examine contrast dependencies that cannot be explained only by differences in relaxivity in a single-spin population.

METHODS

Measurements of functional signal and contrast were obtained in human early visual cortex during a high-contrast visual stimulation over a large range of TEs and for several flip angles. Small voxels (1.5 mm) were used to restrict the measurements to cortical gray matter in early visual areas identified using retinotopic mapping procedures.

RESULTS

Measurements were consistent with models that include 2 spin populations. The dominant population has a relatively short transverse lifetime that is strongly modulated by activation. However, functional contrast is also affected by volume changes between this short-lived population and the longer-lived population.

CONCLUSION

Some of the previously observed "nonclassical" behaviors of functional contrast can be explained by these interacting dual-spin populations.

A compact solution for estimation of physiological parameters from ultrafast prostate dynamic contrast enhanced MRI

The Tofts pharmacokinetic model requires multiple calculations for analysis of dynamic contrast enhanced (DCE) MRI. In addition, the Tofts model may not be appropriate for the prostate. This can result in error propagation that reduces the accuracy of pharmacokinetic measurements. In this study, we present a compact solution allowing estimation of physiological parameters K ^trans and v _e from ultrafast DCE acquisitions, without fitting DCE-MRI data to the standard Tofts pharmacokinetic model. Since the standard Tofts model can be simplified to the Patlak model at early times when contrast efflux from the extravascular extracellular space back to plasma is negligible, K ^trans can be solved explicitly for a specific time. Further, v _e can be estimated directly from the late steady-state signal using the derivative form of Tofts model. Ultrafast DCE-MRI data were acquired from 18 prostate cancer patients on a Philips Achieva 3T-TX scanner. Regions-of-interest (ROIs) for prostate cancer, normal tissue, gluteal muscle, and iliac artery were manually traced. The contrast media concentration as function of time was calculated over each ROI using gradient echo signal equation with pre-contrast tissue T1 values, and using the 'reference tissue' model with a linear approximation. There was strong correlation (r  =  0.88-0.91, p   <  0.0001) between K ^trans extracted from the Tofts model and K ^trans estimated from the compact solution for prostate cancer and normal tissue. Additionally, there was moderate correlation (r  =  0.65-0.73, p   <  0.0001) between extracted versus estimated v _e. Bland-Altman analysis showed moderate to good agreement between physiological parameters extracted from the Tofts model and those estimated from the compact solution with absolute bias less than 0.20 min^-1 and 0.10 for K ^trans and v _e, respectively. The compact solution may decrease systematic errors and error propagation, and could increase the efficiency of clinical workflow. The compact solution requires high temporal resolution DCE-MRI due to the need to adequately sample the early phase of contrast media uptake.

Differential levels of prefrontal cortex glutamate+glutamine in adults with antisocial personality disorder and bipolar disorder: A proton magnetic resonance spectroscopy study

As the main excitatory neurotransmitter in the central nervous system, glutamate, as measured in combination with glutamine (Glx), is implicated in several psychopathologies when levels are aberrant. One illness that shows heightened Glx levels is bipolar disorder (BD), an illness characterized by high impulsivity. In addition, although animal studies have reported elevated levels of Glx in aggressive and impulsive phenotypes, no study, to our knowledge, has reported Glx in the human cortex in relation to aggression. Here, we addressed the question of whether elevated levels of Glx would be present in patients with BD and antisocial personality disorder (ASPD), a condition associated with aggression and, like BD, also presents high impulsivity. We recruited individuals with ASPD (n = 18), individuals with BD (n = 16), and a healthy control group (n = 24). We used proton magnetic resonance spectroscopy to measure relative neurometabolite concentrations in the left dorsolateral prefrontal cortex (dlPFC) and supra-genual anterior cingulate cortex (ACC), two brain regions associated with impulsivity and behavior control. We found significantly elevated levels of Glx in the ASPD group relative to the BD and healthy control groups in the dlPFC (p = .014), and a positive correlation between Glx levels and aggression in the dlPFC in the ASPD group alone (r = .59, p = .026). These findings suggest a link between aggression in ASPD and Glx levels.

Whole knee joint T1 values measured in vivo at 3T by combined 3D ultrashort echo time cones actual flip angle and variable flip angle methods

PURPOSE

To measure T1 relaxations for the major tissues in whole knee joints on a clinical 3T scanner.

METHODS

The 3D UTE-Cones actual flip angle imaging (AFI) method was used to map the transmission radiofrequency field (B1 ) in both short and long T2 tissues, which was then used to correct the 3D UTE-Cones variable flip angle (VFA) fitting to generate accurate T1 maps. Numerical simulation was carried out to investigate the accuracy of T1 measurement for a range of T2 values, excitation pulse durations, and B1 errors. Then, the 3D UTE-Cones AFI-VFA method was applied to healthy volunteers (N = 16) to quantify the T1 of knee tissues including cartilage, meniscus, quadriceps tendon, patellar tendon, anterior cruciate ligament (ACL), posterior cruciate ligament (PCL), marrow, and muscles at 3T.

RESULTS

Numerical simulation showed that the 3D UTE-Cones AFI-VFA technique can provide accurate T1 measurements (error <1%) when the tissue T2 is longer than 1 ms and a 150 μs excitation RF pulse is used and therefore is suitable for most knee joint tissues. The proposed 3D UTE-Cones AFI-VFA method showed an average T1 of 1098 ± 67 ms for cartilage, 833 ± 47 ms for meniscus, 800 ± 66 ms for quadriceps tendon, 656 ± 43 ms for patellar tendon, 873 ± 38 ms for ACL, 832 ± 49 ms for PCL, 379 ± 18 ms for marrow, and 1393 ± 46 ms for muscles.

CONCLUSION

The 3D UTE-Cones AFI-VFA method allows volumetric T1 measurement of the major tissues in whole knee joints on a clinical 3T scanner.

Measurement of microvascular cerebral blood volume changes over the cardiac cycle with ferumoxytol-enhanced T2 * MRI

PURPOSE

This feasibility study investigates the non-invasive measurement of microvascular cerebral blood volume (BV) changes over the cardiac cycle using cardiac-gated, ferumoxytol-enhanced T 2 ∗ MRI.

METHODS

Institutional review board approval was obtained and all subjects provided written informed consent. Cardiac gated MR scans were prospectively acquired on a 3.0T scanner in 22 healthy subjects using T 2 ∗ -weighted sequences with 2D-EPI and 3D spiral trajectories. Images were collected before and after the intravenous administration of 2 doses of ferumoxytol (1 mg FE/kg and 4 mg FE/kg). Cardiac cycle-induced R 2 ∗ (1/ T 2 ∗ ) changes (Δ R 2 ∗ ) and BV changes (ΔBV) throughout the cardiac cycle in gray matter (GM) and white matter (WM) were quantified and differences assessed using ANOVA followed by post hoc analysis.

RESULTS

Δ R 2 ∗ was found to increase in a dose-dependent fashion. A significantly larger increase was observed in GM compared to WM in both 2D and 3D acquisitions (P < 0.050). In addition, Δ R 2 ∗ increased significantly (P < 0.001) post versus pre-contrast injection in GM in both T 2 ∗ MRI acquisitions. Mean GM Δ R 2 ∗ derived from 2D-EPI images was 0.14 ± 0.06 s^-1 pre-contrast and 0.33 ± 0.13 s^-1 after 5 mg FE/kg. In WM, Δ R 2 ∗ was 0.19 ± 0.06 s^-1 pre-contrast, and 0.23 ± 0.06 s^-1 after 5 mg FE/kg. The fractional changes in BV throughout the cardiac cycle were 0.031 ± 0.019% in GM and 0.011 ± 0.008% in WM (P < 0.001) after 5 mg FE/kg.

CONCLUSION

Cardiac-gated, ferumoxytol-enhanced T 2 ∗ MRI enables characterization of microvascular BV changes throughout the cardiac cycle in GM and WM tissue of healthy subjects.

NOVIFAST: A Fast Algorithm for Accurate and Precise VFA MRI Mapping

In quantitative magnetic resonance mapping, the variable flip angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution weighted images in a clinically feasible time. Fast, linear methods that estimate maps from these weighted images have been proposed, such as DESPOT1 and iterative re-weighted linear least squares. More accurate, non-linear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this paper, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared with conventional gradient-based NLLS estimators while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than the efficient implementations of the variable projection (VARPRO) method. Furthermore, NOVIFAST is shown to be robust against initialization.

Multipathway multi-echo (MPME) imaging: all main MR parameters mapped based on a single 3D scan

PURPOSE

Quantitative parameter maps, as opposed to qualitative grayscale images, may represent the future of diagnostic MRI. A new quantitative MRI method is introduced here that requires a single 3D acquisition, allowing good spatial coverage to be achieved in relatively short scan times.

METHODS

A multipathway multi-echo sequence was developed, and at least 3 pathways with 2 TEs were needed to generate T₁ , T₂ , T₂ ^* , B₁ ⁺ , and B₀ maps. The method required the central k-space region to be sampled twice, with the same sequence but with 2 very different nominal flip angle settings. Consequently, scan time was only slightly longer than that of a single scan. The multipathway multi-echo data were reconstructed into parameter maps, for phantom as well as brain acquisitions, in 5 healthy volunteers at 3 T. Spatial resolution, matrix size, and FOV were 1.2 × 1.0 × 1.2 mm³ , 160 × 192 × 160, and 19.2 × 19.2 × 19.2 cm³ (whole brain), acquired in 11.5 minutes with minimal acceleration. Validation was performed against T₁ , T₂ , and T₂ ^* maps calculated from gradient-echo and spin-echo data.

RESULTS

In Bland-Altman plots, bias and limits of agreement for T₁ and T₂ results in vivo and in phantom were -2.9/±125.5 ms (T₁ in vivo), -4.8/±20.8 ms (T₂ in vivo), -1.5/±18.1 ms (T₁ in phantom), and -5.3/±7.4 ms (T₂ in phantom), for regions of interest including given brain structures or phantom compartments. Due to relatively high noise levels, the current implementation of the approach may prove more useful for region of interest-based as opposed to pixel-based interpretation.

CONCLUSIONS

We proposed a novel approach to quantitatively map MR parameters based on a multipathway multi-echo acquisition.

Fast quantitative MRI using controlled saturation magnetization transfer

Purpose

This study demonstrates magnetization transfer (MT) effects directly affect relaxometry measurements and develops a framework that allows single‐pool models to be valid in 2‐pool MT systems.

Methods

A theoretical framework is developed in which a 2‐pool MT system effectively behaves as a single‐pool if the RMS RF magnetic field (B1rms{\text{B}}_{1}^{{{\text{rms}}}}) is kept fixed across all measurements. A practical method for achieving controlled saturation magnetization transfer (CSMT) using multiband RF pulses is proposed. Numerical, Phantom, and in vivo validations were performed directly comparing steady state (SS) estimation approaches that under correct single‐pool assumptions would be expected to vary in precision but not accuracy.

Results

Numerical simulations predict single‐pool estimates obtained from MT model generated data are not consistent for different SS estimation methods, and a systematic underestimation of T₂ is expected. Neither effect occurs under the proposed CSMT approach. Both phantom and in vivo experiments corroborate the numerical predictions. Experimental data highlights that even when using the same relaxometry method, different estimates are obtained depending on which combination of flip angles (FAs) and TRs are used if the CSMT approach is not used. Using CSMT, stable measurements of both T₁ and T₂ are obtained. The measured T₁(T1CSMT)) depends on B1rms{\text{B}}_{1}^{{{\text{rms}}}}, which is therefore an important parameter to specify.

Conclusion

This work demonstrates that conventional single pool relaxometry, which is highly efficient for human studies, results in unreliable parameter estimates in biological tissues because of MT effects. The proposed CSMT framework is shown to allow single‐pool assumptions to be valid, enabling reliable and efficient quantitative imaging to be performed.

Accurate T1 mapping of short T2 tissues using a three-dimensional ultrashort echo time cones actual flip angle imaging-variable repetition time (3D UTE-Cones AFI-VTR) method

PURPOSE

To develop an accurate T1 measurement method for short T2 tissues using a combination of a 3-dimensional ultrashort echo time cones actual flip angle imaging technique and a variable repetition time technique (3D UTE-Cones AFI-VTR) on a clinical 3T scanner.

METHODS

First, the longitudinal magnetization mapping function of the excitation pulse was obtained with the 3D UTE-Cones AFI method, which provided information about excitation efficiency and B1 inhomogeneity. Then, the derived mapping function was substituted into the VTR fitting to generate accurate T1 maps. Numerical simulation and phantom studies were carried out to compare the AFI-VTR method with a B1 -uncorrected VTR method, a B1 -uncorrected variable flip angle (VFA) method, and a B1 -corrected VFA method. Finally, the 3D UTE-Cones AFI-VTR method was applied to bovine bone samples (N = 6) and healthy volunteers (N = 3) to quantify the T1 of cortical bone.

RESULTS

Numerical simulation and phantom studies showed that the 3D UTE-Cones AFI-VTR technique provides more accurate measurement of the T1 of short T2 tissues than the B1 -uncorrected VTR and VFA methods or the B1 -corrected VFA method. The proposed 3D UTE-Cones AFI-VTR method showed a mean T1 of 240 ± 25 ms for bovine cortical bone and 218 ± 10 ms for the tibial midshaft of human volunteers, respectively, at 3 T.

CONCLUSION

The 3D UTE-Cones AFI-VTR method can provide accurate T1 measurements of short T2 tissues such as cortical bone. Magn Reson Med 80:598-608, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

Low-rank magnetic resonance fingerprinting

PURPOSE

Magnetic resonance fingerprinting (MRF) is a relatively new approach that provides quantitative MRI measures using randomized acquisition. Extraction of physical quantitative tissue parameters is performed offline, without the need of patient presence, based on acquisition with varying parameters and a dictionary generated according to the Bloch equation simulations. MRF uses hundreds of radio frequency (RF) excitation pulses for acquisition, and therefore, a high undersampling ratio in the sampling domain (k-space) is required for reasonable scanning time. This undersampling causes spatial artifacts that hamper the ability to accurately estimate the tissue's quantitative values. In this work, we introduce a new approach for quantitative MRI using MRF, called magnetic resonance fingerprinting with low rank (FLOR).

METHODS

We exploit the low-rank property of the concatenated temporal imaging contrasts, on top of the fact that the MRF signal is sparsely represented in the generated dictionary domain. We present an iterative recovery scheme that consists of a gradient step followed by a low-rank projection using the singular value decomposition.

RESULTS

Experimental results consist of retrospective sampling that allows comparison to a well defined reference, and prospective sampling that shows the performance of FLOR for a real-data sampling scenario. Both experiments demonstrate improved parameter accuracy compared to other compressed-sensing and low-rank based methods for MRF at 5% and 9% sampling ratios for the retrospective and prospective experiments, respectively.

CONCLUSIONS

We have shown through retrospective and prospective experiments that by exploiting the low-rank nature of the MRF signal, FLOR recovers the MRF temporal undersampled images and provides more accurate parameter maps compared to previous iterative approaches.

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

PURPOSE

The purpose of this study was to accelerate the acquisition and reconstruction time of 3D magnetic resonance fingerprinting scans.

METHODS

A 3D magnetic resonance fingerprinting scan was accelerated by using a single-shot spiral trajectory with an undersampling factor of 48 in the x-y plane, and an interleaved sampling pattern with an undersampling factor of 3 through plane. Further acceleration came from reducing the waiting time between neighboring partitions. The reconstruction time was accelerated by applying singular value decomposition compression in k-space. Finally, a 3D premeasured B1 map was used to correct for the B1 inhomogeneity.

RESULTS

The T1 and T2 values of the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology MRI phantom showed a good agreement with the standard values, with an average concordance correlation coefficient of 0.99, and coefficient of variation of 7% in the repeatability scans. The results from in vivo scans also showed high image quality in both transverse and coronal views.

CONCLUSIONS

This study applied a fast acquisition scheme for a fully quantitative 3D magnetic resonance fingerprinting scan with a total acceleration factor of 144 as compared with the Nyquist rate, such that 3D T1 , T2 , and proton density maps can be acquired with whole-brain coverage at clinical resolution in less than 5 min. Magn Reson Med 79:2190-2197, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Comparison of ferumoxytol-based cerebral blood volume estimates using quantitative R1 and R2* relaxometry

PURPOSE

Cerebral perfusion is commonly assessed clinically with dynamic susceptibility contrast MRI using a bolus injection of gadolinium-based contrast agents, resulting in semi-quantitative values of cerebral blood volume (CBV). Steady-state imaging with ferumoxytol allows estimation of CBV with the potential for higher precision and accuracy. Prior CBV studies have focused on the signal disrupting T2* effects, but ferumoxytol also has high signal-enhancing T₁ relaxivity. The purpose of this study was to investigate and compare CBV estimation using T₁ and T2*, with the goal of understanding the contrast mechanisms and quantitative differences.

METHODS

Changes in R₁ (1/T₁ ) and R2* (1/ T2*) were measured after the administration of ferumoxytol using high-resolution quantitative approaches. Images were acquired at 3.0T and R₁ was estimated from an ultrashort echo time variable flip angle approach, while R2* was estimated from a multiple gradient echo sequence. Twenty healthy volunteers were imaged at two doses. CBV was derived and compared from relaxometry in gray and white matter using different approaches.

RESULTS

R₁ measurements showed a linear dependence of blood R₁ with respect to dose in large vessels, in contrast to the nonlinear dose-dependence of blood R2* estimates. In the brain parenchyma, R2* showed linear dose-dependency whereas R₁ showed nonlinearity. CBV calculations based on R2* changes in tissue and ferumoxytol blood concentration estimates based on R₁ relaxivity showed the lowest variability in our cohort.

CONCLUSIONS

CBV measurements were successfully derived using a combined approach of R₁ and R2* relaxometry. Magn Reson Med 79:3072-3081, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

STrategically Acquired Gradient Echo (STAGE) imaging, part II: Correcting for RF inhomogeneities in estimating T1 and proton density

PURPOSE

To develop a method for mapping the B₁ transmit (B_1t) and B₁ receive (B_1r) fields from two gradient echo datasets each with a different flip angle and from these two images obtain accurate T₁ and proton density (PD) maps of the brain.

METHODS

A strategically acquired gradient echo (STAGE) data set is collected using two flip angles each with multiple echoes. The B_1t field extraction was based on forcing cortical gray matter and white matter to have specific T₁ values and fitting the resulting B_1t field to a quadratic function. The B_1r field extraction was based on synthesizing isointense images despite there being two or three tissue types present in the brain. This method was tested on 10 healthy volunteers and 20 stroke patients from data acquired at 3.0Tesla.

RESULTS

With the knowledge of the B_1t and B_1r fields, the uniformity of tissue T₁ and PD maps was considerably improved. T₁ values were measured for both the midbrain and basal ganglia and found to be in good agreement with the literature.

DISCUSSION AND CONCLUSIONS

STAGE provides a practical way to assess the B_1t and the B_1r fields which can then be used to correct for spatial variations in the images.

Propagation of error from parameter constraints in quantitative MRI: Example application of multiple spin echo T2 mapping

PURPOSE

Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T₂ mapping is explored.

METHODS

An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T₂ mapping protocols.

RESULTS

Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%.

CONCLUSION

The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T₂ mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Bias and SNR of T1 estimates derived from joint fitting of actual flip-angle and FLASH imaging data with variable flip angles

Previous work has suggested fitting joint AFI/FLASH data for T1 and B1+ by minimizing the 2-norm of the difference between the signal model and measurements will produce unbiased estimates of T1. We demonstrate that, contrary to previous results, the estimator has a substantial bias that varies with both the true T1 and B1+, and the receive channel count. We also demonstrate that the correct ML estimator removes the effect of channel count, and that the choice of AFI protocol has a larger impact of the quality of estimates than the addition of an extra FLASH scan.

Slice profile and B1 corrections in 2D magnetic resonance fingerprinting

PURPOSE

The goal of this study is to characterize and improve the accuracy of 2D magnetic resonance fingerprinting (MRF) scans in the presence of slice profile (SP) and B₁ imperfections, which are two main factors that affect quantitative results in MRF.

METHODS

The SP and B₁ imperfections are characterized and corrected separately. The SP effect is corrected by simulating the radiofrequency pulse in the dictionary, and the B₁ is corrected by acquiring a B₁ map using the Bloch-Siegert method before each scan. The accuracy, precision, and repeatability of the proposed method are evaluated in phantom studies. The effects of both SP and B₁ imperfections are also illustrated and corrected in the in vivo studies.

RESULTS

The SP and B₁ corrections improve the accuracy of the T₁ and T₂ values, independent of the shape of the radiofrequency pulse. The T₁ and T₂ values obtained from different excitation patterns become more consistent after corrections, which leads to an improvement of the robustness of the MRF design.

CONCLUSION

This study demonstrates that MRF is sensitive to both SP and B₁ effects, and that corrections can be made to improve the accuracy of MRF with only a 2-s increase in acquisition time. Magn Reson Med 78:1781-1789, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Data-driven optimized flip angle selection for T1 estimation from spoiled gradient echo acquisitions

PURPOSE

Define criteria for selection of optimal flip angle sets for T1 estimation and evaluate effects on T1 mapping.

THEORY AND METHODS

Flip angle sets for spoiled gradient echo-based T1 mapping were selected by minimizing T1 estimate variance weighted by the joint density of M0 and T1 in an initial acquisition. The effect of optimized flip angle selection on T1 estimate error was measured using simulations and experimental data in the human and rat brain.

RESULTS

For two-point acquisitions, optimized angle sets were similar to those proposed by other groups and, therefore, performed similarly. For multipoint acquisitions, optimal angle sets for T1 mapping in the brain consisted of a repetition of two angles. Implementation of optimal angles reduced T1 estimate variance by 30-40% compared with a multipoint acquisition using a range of angles. Performance of the optimal angle set was equivalent to that of a repetition of the two-angle set selected using criteria proposed by other researchers.

CONCLUSION

Repetition of two carefully selected flip angles notably improves the precision of resulting T1 estimates compared with acquisitions using a range of flip angles. This work provides a flexible and widely applicable optimization method of particular use for those who repeatedly perform T1 estimation. Magn Reson Med 76:792-802, 2016. © 2015 Wiley Periodicals, Inc.

Motion-insensitive determination of B1+ amplitudes based on the bloch-siegert shift in single voxels of moving organs including the human heart

PURPOSE

To reliably determine the amplitude of the transmit radiofrequency ( B1+) field in moving organs like the liver and heart, where most current techniques are usually not feasible.

METHODS

B1+ field measurement based on the Bloch-Siegert shift induced by a pair of Fermi pulses in a double-triggered modified Point RESolved Spectroscopy (PRESS) sequence with motion-compensated crusher gradients has been developed. Performance of the sequence was tested in moving phantoms and in muscle, liver, and heart of six healthy volunteers each, using different arrangements of transmit/receive coils.

RESULTS

B1+ determination in a moving phantom was almost independent of type and amplitude of the motion and agreed well with theory. In vivo, repeated measurements led to very small coefficients of variance (CV) if the amplitude of the Fermi pulse was chosen above an appropriate level (CV in muscle 0.6%, liver 1.6%, heart 2.3% with moderate amplitude of the Fermi pulses and 1.2% with stronger Fermi pulses).

CONCLUSION

The proposed sequence shows a very robust determination of B1+ in a single voxel even under challenging conditions (transmission with a surface coil or measurements in the heart without breath-hold).

B1 and T1 mapping of the breast with a reference tissue method

PURPOSE

To develop a method for mapping the B1 field using a reference signal from a tissue with known T1.

METHODS

Flip angle correction factors were calculated in a region with a known "gold standard" T1; by comparing T1 values from a variable flip angle (VFA) sequence to the "gold standard" and correcting the value of the Ernst angle. The resulting partial B1 map was interpolated for all other regions. In the breast, fat is an ideal reference tissue because its T1 is spatially homogeneous and interpatient variability is low. This method was tested with scans of phantoms and patients (n = 4) on a 3T magnet. The performance of the method was evaluated by comparing the results of VFA T1 mapping with and without B1 correction to inversion recovery (IR) T1 maps.

RESULTS

Phantom data determined that a linear inverse distance weighted interpolation accurately recovered the full B1 map. Use of interpolated maps to correct the VFA data in vivo, reduced the average difference in the T1 of parenchyma between VFA and IR results from 58% to 8%.

CONCLUSION

This proof-of-principle study showed that it is possible to recover a full and accurate map of the B1 field in the breast by using a reference tissue (fat) with an accurately measured T1.

Variable flip angle T1 mapping in the human brain with reduced T2 sensitivity using fast radiofrequency-spoiled gradient echo imaging

PURPOSE

Variable flip angle (VFA) T1 quantification using three-dimensional (3D) radiofrequency (RF) spoiled gradient echo imaging offers the acquisition of whole-brain T1 maps in clinically acceptable times. However, conventional VFA T1 relaxometry is biased by incomplete spoiling (i.e., residual T2 dependency). A new postprocessing approach is proposed to overcome this T2-related bias.

METHODS

T1 is quantified from the signal ratio of two spoiled gradient echo (SPGR) images acquired at different flip angles using an analytical solution for the RF-spoiled steady-state signal in combination with a global T2 guess. T1 accuracy is evaluated from simulations and in vivo 3D SPGR imaging of the human brain at 3 Tesla.

RESULTS

The simulations demonstrated that the sensitivity of VFA T1 mapping to T2 can considerably be reduced using a global T2 guess. The method proved to deliver reliable and accurate T1 values in vivo for white and gray matter in good agreement with inversion recovery reference measurements.

CONCLUSION

Based on a global T2 estimate, the accuracy of VFA T1 relaxometry in the human brain can substantially be improved compared with conventional approaches which rely on the generally wrong assumption of ideal spoiling.

MPnRAGE: A technique to simultaneously acquire hundreds of differently contrasted MPRAGE images with applications to quantitative T1 mapping

PURPOSE

To introduce a new technique called MPnRAGE, which produces hundreds of images with different T1 contrasts and a B1 corrected T1 map.

THEORY AND METHODS

An interleaved three-dimensional radial k-space trajectory with a sliding window reconstruction is used in conjunction with magnetization preparation pulses. This work modifies the SNAPSHOT-FLASH T1 fitting equations for radial imaging with view-sharing and develops a new rapid B1 correction procedure. MPnRAGE is demonstrated in phantoms and volunteers, including two volunteers with eight scans each and eight volunteers with two scans each. T1 values from MPnRAGE were compared with those from fast spin echo inversion recovery (FSE-IR) in phantoms and a healthy human brain at 3 Tesla (T).

RESULTS

The T1 fit for human white and gray matter was T1MPnRAGE = 1.00 · T1FSE-IR + 24 ms, r(2) = 0.990. Voxel-wise coefficient of variation in T1 measurements across eight time points was between 0.02 and 0.08. Region of interest-based T1 values were reproducible to within 2% and agree well with literature values.

CONCLUSION

In the same amount of time as a traditional MPRAGE exam (7.5 min), MPnRAGE was shown to produce hundreds of images with alternate T1 contrasts as well as an accurate and reproducible T1 map that is robust to B1 errors.

Iterative multi-channel radio frequency pulse calibration with improving B1 field uniformity in high field MRI

Background

In high field MRI capable of multi-channel radio frequency (RF) transmission, B₁ shimming is a time-consuming job because conventional B₁ shimming techniques require B₁ mapping for each channel. After acquiring the complex-numbered B₁ field maps, the optimal amplitude and phase of the driving RF pulse are determined for each channel to maximize the B₁ field uniformity in conventional B₁ shimming. However, time-consuming B₁ shimming procedures at the pre-scan may not be tolerated in the clinical imaging in which patient throughput is one of the important factors.

Methods

To avoid the time-consuming B₁ mapping, the first spin echo and the stimulated echo were repeatedly acquired in the slice-selective stimulated echo sequence without imaging gradients. A cost function of the amplitudes and phases of the driving RF pulse for every channel was defined in a way that the ratio between the spin echo and stimulated echo amplitudes rapidly converged to √ 2. The amplitude and phase of the driving RF pulse were iteratively modified over the repeating RF pulse sequence so that the cost function was minimized.

Results

From the finite-difference-time-domain (FDTD) electromagnetic field simulations with a human body model placed in a birdcage coil operating at 3 T, it was observed that the RF pulse calibration with iterative cost function minimization can give improvement of B₁ field uniformity as well as flip-angle calibration. The experiments at 3 T also showed improvement of RF field uniformity in the phantom imaging studies.

Conclusions

Since the proposed RF pulse calibration is not based on B₁ mapping, the RF pulse calibration time could be much shorter than the B₁-mapping based methods. The proposed method is expected to be a practical substitute for the B₁-mapping-based B₁ shimming methods when long pre-scan time is not tolerable.

Minimum acquisition methods for simultaneously imaging T(1), T(2), and proton density with B(1) correction and no spin-echoes

The spin lattice (T(1)) and spin-spin (T(2)) relaxation times, along with the proton density (PD) contain almost all of the information that (1)H MRI routinely uses in clinical diagnosis and research, but are seldom imaged directly. Here, three methods for directly imaging T(1), T(2), and PD with the least possible number of acquisitions - three, are presented. All methods utilize long 0° self-refocusing adiabatic pre-pulses instead of spin-echoes to encode the T(2) information prior to a conventional gradient-echo MRI sequence. T(1) information is encoded by varying the flip-angle (FA) in the 'Dual-τ Dual-FA' and 'Four-FA' methods, or the sequence repetition period, TR, in the 'Dual-τ Dual-TR' method. Inhomogeneity in the FA distribution and slice-selection profile are recognized as the main error sources for T(1) measurements. The former is remedied by integrating an extra FA-dependent acquisition into the 'Four-FA' method to provide self-corrected T(1), T(2), PD, and FA in just four acquisitions - again, the minimum possible. Slice profile errors - which manifest as differences between 2D and 3D T(1) measurements, can be addressed by Bloch equation analysis and experimental calibration. All three methods are validated in phantom studies, and the 'Dual-τ Dual-FA' and 'Four-FA' methods are validated in human brain studies using standard partial saturation and spin-echo methods for reference. The new methods offer a minimum-acquisition option for imaging single-component T(1), T(2), and PD. 'Four-FA' performs best overall in accuracy, with high efficiency per unit accuracy vs. existing methods when B(1)-inhomogeneity is appropriately addressed.

Fast B1 mapping based on interleaved-three-flip-angle (ITFA) excitation

PURPOSE

Fast B1 mapping based on short-TR sequences is prone to T1-induced errors. The purpose of this study is to develop a novel fast B1 mapping method that is less prone to T1-induced errors.

METHODS

The authors acquired three gradient echoes by applying three RF pulses of different flip angles in an interleaved manner. The new method, named interleaved-three-flip-angle (ITFA) method, employs a short TR for fast scan. Since the pixel intensity of the gradient echo images is dependent on both B1 and T1, the authors could compute a B1 map from the three gradient echo images with excluding the T1-effects. The authors simulated the proposed B1 mapping method for various T1 values, and they found optimal flip angles for ITFA experiments for a given repetition time. To evaluate the B1 mapping performance, the authors made a human-brain-mimicking phantom that had six compartments with different T1 and T2. The authors performed B1 mapping experiments at 3T on the phantom and a volunteer using the ITFA method, the actual flip angle imaging (AFI) method, and the double angle method (DAM), and then, the authors compared the B1 mapping results.

RESULTS

Using a birdcage coil as a transceiver at 3T, the authors performed ITFA scans of the phantom and a volunteer with TR of 60 ms and the nominal flips angles of (25°, 70°, 80°). The authors also performed AFI scans with TR1/TR2 of 30/150 ms and the nominal flip angle of 60°. In both the phantom and human head imaging performed with the same scan times for ITFA and AFI, ITFA showed smaller average B1 errors than AFI when they were compared to DAM.

CONCLUSIONS

ITFA excitations made it possible to reduce the T1-effects on B1 mapping of the human-brain-mimicking phantom and the human brain at 3T. The authors expect the ITFA method can be used for B1 shimming once the optimal flip angles have been predetermined for the target imaging region and for the preferred TR.

Improved encoding pulses for Bloch-Siegert B1(+) mapping

A new family of optimized encoding pulses for Bloch-Siegert (BS) |B(1)(+)| mapping is introduced, as well as an algorithm to design them. The pulses are designed by numerical maximization of BS sequence sensitivity, subject to constraints that ensure low on-resonance excitation. The pulses are in all cases characterized by a constant envelope and U-shaped frequency sweep. They are validated in simulations, 7T in vivo experiments, and an experiment to measure their on-resonance excitation, and are compared to a Fermi pulse conventionally used in the BS method. The pulses are shown to produce larger phase shifts in a shorter time and with lower on-resonance excitation than the Fermi pulse, which results in lower SAR and improved |B(1)(+)| accuracy in areas of the body with large main field inhomogeneities.

Accelerating MR parameter mapping using sparsity-promoting regularization in parametric dimension

MR parameter mapping requires sampling along additional (parametric) dimension, which often limits its clinical appeal due to a several-fold increase in scan times compared to conventional anatomic imaging. Data undersampling combined with parallel imaging is an attractive way to reduce scan time in such applications. However, inherent SNR penalties of parallel MRI due to noise amplification often limit its utility even at moderate acceleration factors, requiring regularization by prior knowledge. In this work, we propose a novel regularization strategy, which uses smoothness of signal evolution in the parametric dimension within compressed sensing framework (p-CS) to provide accurate and precise estimation of parametric maps from undersampled data. The performance of the method was demonstrated with variable flip angle T1 mapping and compared favorably to two representative reconstruction approaches, image space-based total variation regularization and an analytical model-based reconstruction. The proposed p-CS regularization was found to provide efficient suppression of noise amplification and preservation of parameter mapping accuracy without explicit utilization of analytical signal models. The developed method may facilitate acceleration of quantitative MRI techniques that are not suitable to model-based reconstruction because of complex signal models or when signal deviations from the expected analytical model exist.