Comparison of uniform-density, variable-density, and dual-density spiral samplings for multi-shot DWI

PURPOSE

To compare the performances of uniform-density spiral (UDS), variable-density spiral (VDS), and dual-density spiral (DDS) samplings in multi-shot diffusion imaging, and determine a sampling strategy that balances reliability of shot navigator and overall DWI image quality.

THEORY AND METHODS

UDS, VDS, and DDS trajectories were implemented to achieve four-shot diffusion-weighted spiral imaging. First, the static B0 off-resonance effects in UDS, VDS, and DDS acquisitions were analyzed based on a signal model. Then, in vivo experiments were performed to verify the theoretical analyses, and fractional anisotropy (FA) fitting residuals were used to quantitatively assess the quality of spiral diffusion data for tensor estimation. Finally, the SNR performances and g-factor behavior of the three spiral samplings were evaluated using a Monte Carlo-based pseudo multiple replica method.

RESULTS

Among the three spiral trajectories with the same readout duration, UDS sampling exhibited the least off-resonance artifacts. This was most evident when the static B0 off-resonance effect was severe. The UDS diffusion images had higher anatomical fidelity and lower FA fitting residuals than the other two counterparts. Furthermore, the four-shot UDS acquisition achieved the best SNR performance in diffusion imaging with 12.11% and 40.85% improvements over the VDS and DDS acquisitions with the same readout duration, respectively.

CONCLUSION

UDS sampling is an efficient spiral acquisition scheme for high-resolution diffusion imaging with reliable navigator information. It provides superior off-resonance performance and SNR efficiency over the VDS and DDS samplings for the tested scenarios.

Quantitative effects of off-resonance related distortion on brain mechanical property estimation with magnetic resonance elastography

Off-resonance related geometric distortion can impact quantitative MRI techniques, such as magnetic resonance elastography (MRE), and result in errors to these otherwise sensitive metrics of brain health. MRE is a phase contrast technique to determine the mechanical properties of tissue by imaging shear wave displacements and estimating tissue stiffness through inverse solution of Navier's equation. In this study, we systematically examined the quantitative effects of distortion and corresponding correction approaches on MRE measurements through a series of simulations, phantom models, and in vivo brain experiments. We studied two different k-space trajectories, echo-planar imaging and spiral, and we determined that readout time, off-resonance gradient strength, and the combination of readout direction and off-resonance gradient direction, impact the estimated mechanical properties. Images were also processed through traditional distortion correction pipelines, and we found that each of the correction mechanisms works well for reducing stiffness errors, but are limited in cases of very large distortion. The ability of MRE to detect subtle changes to neural tissue health relies on accurate, artifact-free imaging, and thus off-resonance related geometric distortion must be considered when designing sequences and protocols by limiting readout time and applying correction where appropriate.

Fast 3D MR elastography of the whole brain using spiral staircase: Data acquisition, image reconstruction, and joint deblurring

PURPOSE

To address the need for a method to acquire 3D data for MR elastography (MRE) of the whole brain with substantially improved spatial resolution, high SNR, and reduced acquisition time compared with conventional methods.

METHODS

We combined a novel 3D spiral staircase data-acquisition method with a spoiled gradient-echo pulse sequence and MRE motion-encoding gradients (MEGs). The spiral-out acquisition permitted use of longer-duration motion-encoding gradients (ie, over two full oscillatory cycles) to enhance displacement SNR, while still maintaining a reasonably short TE for good phase-SNR. Through-plane parallel imaging with low noise penalties was implemented to accelerate acquisition along the slice direction. Shared anatomical information was exploited in the deblurring procedure to further boost SNR for stiffness inversion.

RESULTS

In vivo and phantom experiments demonstrated the feasibility of the proposed method in producing brain MRE results comparable to the spin-echo-based approaches, both qualitatively and quantitatively. High-resolution (2-mm isotropic) brain MRE data were acquired in 5 minutes using our method with good SNR. Joint deblurring with shared anatomical information produced SNR-enhanced images, leading to upward stiffness estimation.

CONCLUSION

A novel 3D gradient-echo-based approach has been designed and implemented, and shown to have promising potential for fast and high-resolution in vivo MRE of the whole brain.

Technical Note: Effects of rotating gantry on magnetic field and eddy currents in 0.35 T MRI-guided radiotherapy (MR-IGRT) system

PURPOSE

The purpose of this study was to identify the cause of severe image artifacts that occurred during gantry rotation in a 0.35 T MRI-Linac by comparing measurements of eddy currents, center frequency, and field inhomogeneities made with the gantry in motion and stationary.

METHODS

Gradient and B₀ eddy currents were calculated from the free induction decays (FIDs) resulting from selective excitation at a temporal resolution of 200 ms/measurement. B₀ eddy currents were also calculated from FIDs acquired with nonselective excitation at a temporal resolution of 100 ms/measurement. Center frequencies and B₀ inhomogeneities were measured by acquiring FIDs with a repetition time (TR) of 290 ms. Cartesian and radial 2D true fast imaging with steady-state precession (TrueFISP) pulse sequences used in real-time MRI-guided radiation therapy (MR-IGRT) were acquired. To assess artifact severity, the normalized root mean square error (nRMSE) was calculated between a reference MRI (static gantry) and MRIs acquired during gantry rotation for each serial acquisition. Image artifacts were qualitatively graded as nominal, minor, or severe. Measurements were conducted while the gantry was rotated through its entire range for both clockwise and counterclockwise. Measurements during gantry rotation were compared to measurements with a stationary gantry (every 30°).

RESULTS

Severe image artifacts were observed 22-35% of the time while the gantry was rotating. Short time constant eddy currents were not affected by gantry rotation. The peak to peak center frequency and FWHM rose by factors of 13.2-14.5 and 1.1-1.6, respectively, for the rotating versus stationary gantry. The magnitude of the center frequency offset and field inhomogeneities depended on the direction of the gantry rotation.

CONCLUSIONS

Image artifacts during gantry rotation were primarily caused by center frequency variations and field inhomogeneities. Therefore, dynamic B₀ compensation techniques should be able to reduce artifacts during gantry rotation.

Removal of off-resonance xenon gas artifacts in pulmonary gas-transfer MRI

PURPOSE

Hyperpolarized xenon (¹²⁹ Xe) gas-transfer imaging allows different components of pulmonary gas transfer-alveolar air space, lung interstitium/blood plasma (barrier), and red blood cells (RBCs)-to be assessed separately in a single breath. However, quantitative analysis is challenging because dissolved-phase ¹²⁹ Xe images are often contaminated by off-resonant gas-phase signal generated via imperfectly selective excitation. Although previous methods required additional data for gas-phase removal, the method reported here requires no/minimal sequence modifications/data acquisitions, allowing many previously acquired images to be corrected retroactively.

METHODS

¹²⁹ Xe imaging was implemented at 3.0T via an interleaved three-dimensional radial acquisition of the gaseous and dissolved phases (using one-point Dixon reconstruction for the dissolved phase) in 46 human subjects and a phantom. Gas-phase contamination (9.5% ± 4.8%) was removed from gas-transfer data using a modified gas-phase image. The signal-to-noise ratio (SNR) and signal distributions were compared before and after contamination removal. Additionally, theoretical gaseous contaminations were simulated at different magnetic field strengths for comparison.

RESULTS

Gas-phase contamination at 3.0T was more diffuse and located predominantly outside the lungs, relative to simulated 1.5T contamination caused by the larger frequency offset. Phantom experiments illustrated a 91% removal efficiency. In human subjects, contamination removal produced significant changes in dissolved signal SNR (+7.8%), mean (-1.4%), and standard deviation (-2.3%) despite low contamination. Repeat measurements showed reduced variance (dissolved mean, -1.0%; standard deviation, -8.4%).

CONCLUSION

Off-resonance gas-phase contamination can be removed robustly with no/minimal sequence modifications. Contamination removal permits more accurate quantification, reduces radiofrequency stringency requirements, and increases data consistency, providing improved sensitivity needed for multicenter trials.

Deblurring for spiral real-time MRI using convolutional neural networks

PURPOSE

To develop and evaluate a fast and effective method for deblurring spiral real-time MRI (RT-MRI) using convolutional neural networks.

METHODS

We demonstrate a 3-layer residual convolutional neural networks to correct image domain off-resonance artifacts in speech production spiral RT-MRI without the knowledge of field maps. The architecture is motivated by the traditional deblurring approaches. Spatially varying off-resonance blur is synthetically generated by using discrete object approximation and field maps with data augmentation from a large database of 2D human speech production RT-MRI. The effect of off-resonance range, shift-invariance of blur, and readout durations on deblurring performance are investigated. The proposed method is validated using synthetic and real data with longer readouts, quantitatively using image quality metrics and qualitatively via visual inspection, and with a comparison to conventional deblurring methods.

RESULTS

Deblurring performance was found superior to a current autocalibrated method for in vivo data and only slightly worse than an ideal reconstruction with perfect knowledge of the field map for synthetic test data. Convolutional neural networks deblurring made it possible to visualize articulator boundaries with readouts up to 8 ms at 1.5 T, which is 3-fold longer than the current standard practice. The computation time was 12.3 ± 2.2 ms per frame, enabling low-latency processing for RT-MRI applications.

CONCLUSION

Convolutional neural networks deblurring is a practical, efficient, and field map-free approach for the deblurring of spiral RT-MRI. In the context of speech production imaging, this can enable 1.7-fold improvement in scan efficiency and the use of spiral readouts at higher field strengths such as 3 T.

Off-Resonant Absorption Enhancement in Single Nanowires via Graded Dual-Shell Design

Single nanowires (NWs) are of great importance for optoelectronic applications, especially solar cells serving as powering nanoscale devices. However, weak off-resonant absorption can limit its light-harvesting capability. Here, we propose a single NW coated with the graded-index dual shells (DSNW). We demonstrate that, with appropriate thickness and refractive index of the inner shell, the DSNW exhibits significantly enhanced light trapping compared with the bare NW (BNW) and the NW only coated with the outer shell (OSNW) and the inner shell (ISNW), which can be attributed to the optimal off-resonant absorption mode profiles due to the improved coupling between the reemitted light of the transition modes of the leak mode resonances of the Si core and the nanofocusing light from the dual shells with the graded refractive index. We found that the light absorption can be engineered via tuning the thickness and the refractive index of the inner shell, the photocurrent density is significantly enhanced by 134% (56%, 12%) in comparison with that of the BNW (OSNW, ISNW). This work advances our understanding of how to improve off-resonant absorption by applying graded dual-shell design and provides a new choice for designing high-efficiency single NW photovoltaic devices.

Slice-selective extended phase graphs in gradient-crushed, transient-state free precession sequences: An application to MR fingerprinting

PURPOSE

Slice-selective, gradient-crushed, transient-state sequences such as those used in MR fingerprinting (MRF) relaxometry are sensitive to slice profile effects. Whereas balanced steady-state free precession MRF profile effects have been studied, less attention has been given to gradient-crushed MRF forms. Extensions of the extended phase graph (EPG) formalism, called slice-selective EPG (ssEPG), are proposed that model slice profile effects.

THEORY AND METHODS

The hard-pulse approximation to slice-selective excitation in the spatial domain is reformulated in k-space. Excitation is modeled by standard EPG shift and transition operators. This ssEPG modeling is validated against Bloch simulations and phantom slice profile measurements. ssEPG relaxometry accuracy and variability are compared with other EPG methods in phantoms and human leg in vivo. The role of ∆B₀ interactions with slice profile and gradient crushers is investigated.

RESULTS

Simulations and slice profile measurements show that ssEPG can model highly dynamic slice profile effects of gradient-crushed sequences. The MRF ssEPG T₂ estimates over 0 < T₂ < 100 ms improve accuracy by > 10 ms at some values relative to other modeling approaches. Small deviations in B₀ can produce substantial bias in T₂ estimations from a range of MRF sequence types, and these effects can be modeled and understood by ssEPG.

CONCLUSION

Transient-state, gradient-crushed sequences such as those used in MRF are sensitive to slice profile effects, and these effects depend on RF pulse choice, gradient crusher strength, and ∆B₀ . It was found ssEPG was the most accurate EPG-based means to model these effects.

Improved reconstruction stability for chemical shift encoded hyperpolarized 13 C magnetic resonance spectroscopic imaging using k-t spiral acquisitions

PURPOSE

A multiecho, field of view (FOV)-oversampled k-t spiral acquisition and direct iterative decomposition of water and fat with echo asymmetry and least-squares estimation reconstruction is demonstrated to improve the stability of hyperpolarized ¹³ C magnetic resonance spectroscopic imaging (MRSI) in the presence of signal ambiguities attributed to low-SNR (signal-to-noise-ratio) species, local uncertainties in metabolite peaks, and echo-to-echo signal inconsistencies.

THEORY

k-t spiral acquisitions redistribute readout points to be more densely spaced radially in k-space by acquiring an FOV and matrix that are oversampled by η. These more densely spaced spiral turns constitute effective intraspiral echoes and can supplement conventional interspiral echoes to improve spectral separation and reduce spectral cross-talk to better resolve ¹³ C-labeled species for spectroscopic imaging.

METHODS

Digital simulations and imaging phantom experiments were performed for a range of interspiral echo spacings and η using multiecho, k-t spiral acquisitions. Image spectral cross-talk artifacts were evaluated both qualitatively and quantitatively as the percent error in measured metabolite ratios. In vivo murine experiments evaluated the feasibility of multiecho, k-t spiral [1-¹³ C]pyruvate MRSI to reduce spectral cross-talk for 3 scenarios of different expected reconstruction stability.

RESULTS

Digital simulations and imaging phantom experiments both demonstrated reduced or comparable image spectral cross-talk and percent errors in measured metabolite ratios with increasing η and better choices of echo spacings. In vivo images displayed markedly reduced spectral cross-talk in lactate images acquired with η = 7 versus η = 1.

CONCLUSION

The precision of hyperpolarized ¹³ C metabolic imaging and quantification in the presence of low-SNR species, local uncertainties in metabolite resonances, and echo-to-echo signal inconsistencies can be improved with the use of FOV-oversampled k-t spiral acquisitions.

Improved off-resonance phase behavior using a phase-inverted adiabatic half-passage pulse for 13 C MRS in humans at 7 T

In vivo ¹³ C MRS at high field benefits from an improved SNR and spectral resolution especially when using surface coils in combination with adiabatic pulses, such as the adiabatic half-passage (AHP) pulse for ¹³ C excitation. However, the excitation profile of the AHP pulse is asymmetric relative to the carrier frequency, which could lead to asymmetric excitation of the spectral lines relative to the center of the spectrum. In this study, a pulse-acquire sequence was designed for adiabatic ¹³ C excitation with a symmetric bandwidth, utilizing a combination of two AHP pulses with inverted phases in alternate scans. Magnetization and phase behavior as a function of frequency offset and RF amplitude of the B₁ field, as well as the steady-state transverse magnetization response to off-resonance, were simulated. Excitation properties of the combined pulse sequence were studied by ²³ Na imaging and ¹³ C spectroscopy in vitro on a phantom and in vivo on the human calf at 7 T. Simulations demonstrated symmetric transverse magnetization and phase with respect to positive and negative frequency offsets when using two AHP pulses with inverted phases in alternate scans, thereby minimizing baseline distortion and achieving symmetric T₁ weighting, as confirmed by in vitro measurements. The intensities of the lipid peaks at 15, 30, 62, 73, and 130 ppm were in agreement with those theoretically predicted using two AHP pulses with inverted phases in alternate scans. We conclude that using two phase-inverted AHP pulses improves the symmetry of the ¹³ C excitation profile and phase response to off-resonance effects at 7 T in comparison with using a single AHP pulse.

Practical considerations for territorial perfusion mapping in the cerebral circulation using super-selective pseudo-continuous arterial spin labeling

PURPOSE

This paper discusses several challenges faced by super-selective pseudo-continuous arterial spin labeling, which is used to quantify territorial perfusion in the cerebral circulation. The effects of off-resonance, pulsatility, vessel movement, and label rotation scheme are investigated, and methods to maximize labeling efficiency and overall image quality are evaluated. A strategy to calculate the territorial perfusion fractions of individual vessels is proposed.

METHODS

The effects of off-resonance, label rotation scheme, and vessel movement on labeling efficiency were simulated. Two off-resonance compensation strategies (multiphase prescan, field map), cardiac triggering, and vessel movement were studied in vivo in a group of 10 subjects. Subsequently, a territorial perfusion fraction map was acquired in 2 subjects based on the mean vessel labeling efficiency.

RESULTS

Multiphase calibration provided the highest labeling efficiency (P = .002) followed by the field map compensation (P = .037) compared with the uncompensated acquisition. Cardiac triggering resulted in a qualitative improvement of the image and an increase in signal contrast between the perfusion territory and the surrounding tissue (P = .010) but failed to show a significant change in temporal and spatial SNR. The constant clockwise label rotation scheme yielded the highest labeling efficiency. Significant vessel movement (>2 mm according to simulations) was observed in 50% of subjects. The measured territorial perfusion fractions showed good agreement with anatomical data.

CONCLUSION

Optimized labeling efficiency resulted in increased image quality and accuracy of territorial perfusion fraction maps. Labeling efficiency depends critically on off-resonance calibration, cardiac triggering, optimal label rotation scheme, and vessel location tracking.

Cross-relaxation parameters in cortical bone assessed with different MR sequences

This study aimed to show evidence of MR cross-relaxation effects in cortical bone and to compare different MR sequences for the quantification of cross-relaxation parameters. Measurements were performed on bovine diaphysis samples with spectroscopic methods (inversion-recovery, off-resonance saturation) and with a variable flip angle (VFA) UTE imaging method on a 4.7 T laboratory-assembled scanner. Cross-relaxation parameter assessment was carried out via a two-pool model simulation with a matrix algebra approach. A proton signal amplitude of 28 Mol/L was observed (equivalent water fraction of 25%). It was attributed to collagen-bound water, with T2* values of ~ 0.3 ms, a "long-T₂ " proton pool, in exchange with protons from the collagen macromolecules ( T2* of 10-20 μs). Magnetization transfer (MT) effects were detected with all sequences. The best precision of model parameters was obtained with off-resonance saturation; the fraction of collagen methylene protons was found in the range of 22-28% and the transverse relaxation time for collagen methylene protons was 11 μs (1% precision). The model parameters obtained were compatible with VFA-UTE results but could not be assessed with acceptable accuracy and precision using this method. In vivo MT quantification using off-resonance saturation with a single B₁ amplitude and offset frequency may provide information about the relative amount of collagen per unit volume in cortical bone.

Dynamic water/fat separation and B 0 inhomogeneity mapping-joint estimation using undersampled triple-echo multi-spoke radial FLASH

PURPOSE

To achieve dynamic water/fat separation and B 0 field inhomogeneity mapping via model-based reconstructions of undersampled triple-echo multi-spoke radial FLASH acquisitions.

METHODS

This work introduces an undersampled triple-echo multi-spoke radial FLASH sequence, which uses (i) complementary radial spokes per echo train for faster spatial encoding, (ii) asymmetric echoes for flexible and nonuniform echo spacing, and (iii) a golden angle increment across frames for optimal k-space coverage. Joint estimation of water, fat, B 0 inhomogeneity, and coil sensitivity maps from undersampled triple-echo data poses a nonlinear and non-convex inverse problem which is solved by a model-based reconstruction with suitable regularization. The developed methods are validated using phantom experiments with different degrees of undersampling. Real-time MRI studies of the knee, liver, and heart are conducted without prospective gating or retrospective data sorting at temporal resolutions of 70, 158, and 40 ms, respectively.

RESULTS

Up to 18-fold undersampling is achieved in this work. Even in the presence of rapid physiological motion, large B 0 field inhomogeneities, and phase wrapping, the model-based reconstruction yields reliably separated water/fat maps in conjunction with spatially smooth inhomogeneity maps.

CONCLUSIONS

The combination of a triple-echo acquisition and joint reconstruction technique provides a practical solution to time-resolved and motion robust water/fat separation at high spatial and temporal resolution.

Reducing distortions in echo-planar breast imaging at ultrahigh field with high-resolution off-resonance maps

PURPOSE

DWI is a promising modality in breast MRI, but its clinical acceptance is slow. Analysis of DWI is hampered by geometric distortion artifacts, which are caused by off-resonant spins in combination with the low phase-encoding bandwidth of the EPI sequence used. Existing correction methods assume smooth off-resonance fields, which we show to be invalid in the human breast, where high discontinuities arise at tissue interfaces.

METHODS

We developed a distortion correction method that incorporates high-resolution off-resonance maps to better solve for severe distortions at tissue interfaces. The method was evaluated quantitatively both ex vivo in a porcine tissue phantom and in vivo in 5 healthy volunteers. The added value of high-resolution off-resonance maps was tested using a Wilcoxon signed rank test comparing the quantitative results obtained with a low-resolution off-resonance map with those obtained with a high-resolution map.

RESULTS

Distortion correction using low-resolution off-resonance maps corrected most of the distortions, as expected. Still, all quantitative comparison metrics showed increased conformity between the corrected EPI images and a high-bandwidth reference scan for both the ex vivo and in vivo experiments. All metrics showed a significant improvement when a high-resolution off-resonance map was used (P < 0.05), in particular at tissue boundaries.

CONCLUSION

The use of off-resonance maps of a resolution higher than EPI scans significantly improves upon existing distortion correction techniques, specifically by superior correction at glandular tissue boundaries.

Mapping tissue water T1 in the liver using the MOLLI T1 method in the presence of fat, iron and B0 inhomogeneity

Modified Look-Locker inversion recovery (MOLLI) T1 mapping sequences can be useful in cardiac and liver tissue characterization, but determining underlying water T1 is confounded by iron, fat and frequency offsets. This article proposes an algorithm that provides an independent water MOLLI T1 (referred to as on-resonance water T1 ) that would have been measured if a subject had no fat and normal iron, and imaging had been done on resonance. Fifteen NiCl2 -doped agar phantoms with different peanut oil concentrations and 30 adults with various liver diseases, nineteen (63.3%) with liver steatosis, were scanned at 3 T using the shortened MOLLI (shMOLLI) T1 mapping, multiple-echo spoiled gradient-recalled echo and 1 H MR spectroscopy sequences. An algorithm based on Bloch equations was built in MATLAB, and water shMOLLI T1 values of both phantoms and human participants were determined. The quality of the algorithm's result was assessed by Pearson's correlation coefficient between shMOLLI T1 values and spectroscopically determined T1 values of the water, and by linear regression analysis. Correlation between shMOLLI and spectroscopy-based T1 values increased, from r = 0.910 (P < 0.001) to r = 0.998 (P < 0.001) in phantoms and from r = 0.493 (for iron-only correction; P = 0.005) to r = 0.771 (for iron, fat and off-resonance correction; P < 0.001) in patients. Linear regression analysis revealed that the determined water shMOLLI T1 values in patients were independent of fat and iron. It can be concluded that determination of on-resonance water (sh)MOLLI T1 independent of fat, iron and macroscopic field inhomogeneities was possible in phantoms and human subjects.

Improving the robustness of pseudo-continuous arterial spin labeling to off-resonance and pulsatile flow velocity

PURPOSE

To improve pseudo-continuous arterial spin labeling (PCASL) robustness to off-resonance and pulsatile blood flow velocity.

METHODS

Bloch equations were solved to evaluate the effect of labeling parameters in a pulsatile flow model for a range of off-resonance. Experimental confirmation was achieved in volunteers using linear phase increase between labeling pulses to approximate off-resonance errors. We first assessed the location of the labeling plane in four volunteers. Next, we explored a range of parameters-including balanced and unbalanced gradients-in five more volunteers at an optimal labeling plane location.

RESULTS

Simulations demonstrated that 1) high velocities are vulnerable to off-resonance, 2) unbalanced PCASL outperforms balanced PCASL, 3) increased B1 and low average gradient improve the labeling efficiency for high-velocity flow, and 4) a low ratio of selective to average gradient improves off-resonance robustness. A good agreement between theory and experiment was observed.

CONCLUSION

The robustness of PCASL can be increased by selecting an unbalanced scheme with a low average gradient (0.5 mT/m), a low ratio (7×) of selective to average gradients, and the highest feasible B1 (1.8 μT). Placing the labeling plane above the carotid bifurcation and below the V3 segment, usually between the second and third vertebrae, yielded robust results. Magn Reson Med 78:1342-1351, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

PLANET: An ellipse fitting approach for simultaneous T1 and T2 mapping using phase-cycled balanced steady-state free precession

Purpose

To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T₁ and T₂ from a single 3D phase‐cycled balanced steady‐state free precession (bSSFP) sequence.

Methods

A method is presented in which the elliptical signal model is used to describe the phase‐cycled bSSFP steady‐state signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T₁ and T₂, the banding free magnitude, and the off‐resonance are calculated from the fitting results.

Results

Maps of T₁ and T₂, as well as an off‐resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase‐cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T₁ was observed in the brain.

Conclusion

The presented method allows for accurate mapping of relaxation times and off‐resonance, and for the reconstruction of banding free magnitude images at realistic signal‐to‐noise ratios. Magn Reson Med 79:711–722, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

Externally calibrated parallel imaging for 3D multispectral imaging near metallic implants using broadband ultrashort echo time imaging

PURPOSE

To develop an externally calibrated parallel imaging technique for three-dimensional multispectral imaging (3D-MSI) in the presence of metallic implants.

THEORY AND METHODS

A fast, ultrashort echo time (UTE) calibration acquisition is proposed to enable externally calibrated parallel imaging techniques near metallic implants. The proposed calibration acquisition uses a broadband radiofrequency (RF) pulse to excite the off-resonance induced by the metallic implant, fully phase-encoded imaging to prevent in-plane distortions, and UTE to capture rapidly decaying signal. The performance of the externally calibrated parallel imaging reconstructions was assessed using phantoms and in vivo examples.

RESULTS

Phantom and in vivo comparisons to self-calibrated parallel imaging acquisitions show that significant reductions in acquisition times can be achieved using externally calibrated parallel imaging with comparable image quality. Acquisition time reductions are particularly large for fully phase-encoded methods such as spectrally resolved fully phase-encoded three-dimensional (3D) fast spin-echo (SR-FPE), in which scan time reductions of up to 8 min were obtained.

CONCLUSION

A fully phase-encoded acquisition with broadband excitation and UTE enabled externally calibrated parallel imaging for 3D-MSI, eliminating the need for repeated calibration regions at each frequency offset. Significant reductions in acquisition time can be achieved, particularly for fully phase-encoded methods like SR-FPE. Magn Reson Med 77:2303-2309, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

On the accuracy and precision of cardiac magnetic resonance T2 mapping: A high-resolution radial study using adiabatic T2 preparation at 3 T

PURPOSE

The goal of this study was to characterize the accuracy and precision of cardiac T₂ mapping as a function of different factors including low signal-to-noise ratio (SNR), imaging in systole, and off-resonance frequencies.

METHODS

Bloch equation and Monte Carlo simulations were used to determine the influence of SNR and the choice of T₂ preparation echo time (TE_T2prep ) increments on the accuracy and precision of high-resolution radial cardiac T₂ mapping at 3.0 T. Healthy volunteers were scanned to establish the difference in precision and inter- and intraobserver variability between T₂ mapping in diastole and systole, as well as the effect of SNR and off-resonance frequencies on the accuracy of T₂ maps.

RESULTS

The simulations demonstrated that a TE_T2prep increment of ∼0.75 times the T₂ value of interest optimally increases the precision of the T₂ fit. Systolic T₂ maps were found to have a higher precision (P = 0.002), but similar inter- and intraobserver variability compared with diastolic T₂ maps, whereas off-resonance frequencies beyond ± 100 Hz cause a significant decrease in both accuracy and precision (P < 0.05).

CONCLUSION

This evaluation of the accuracy and precision of cardiac T₂ mapping characterizes the major vulnerabilities of the technique and will help guide protocol definition of studies that include T₂ mapping. Magn Reson Med 77:159-169, 2017. © 2016 Wiley Periodicals, Inc.

Gradient-Modulated PETRA MRI

Image blurring due to off-resonance and fast T 2* signal decay is a common issue in radial ultrashort echo time MRI sequences. One solution is to use a higher readout bandwidth, but this may be impractical for some techniques like pointwise encoding time reduction with radial acquisition (PETRA), which is a hybrid method of zero echo time and single point imaging techniques. Specifically, PETRA has severe specific absorption rate (SAR) and radiofrequency (RF) pulse peak power limitations when using higher bandwidths in human measurements. In this study, we introduce gradient modulation (GM) to PETRA to reduce image blurring artifacts while keeping SAR and RF peak power low. Tolerance of GM-PETRA to image blurring was evaluated in simulations and experiments by comparing with the conventional PETRA technique. We performed inner ear imaging of a healthy subject at 7T. GM-PETRA showed significantly less image blurring due to off-resonance and fast T2* signal decay compared to PETRA. In in vivo imaging, GM-PETRA nicely captured complex structures of the inner ear such as the cochlea and semicircular canals. Gradient modulation can improve the PETRA image quality and mitigate SAR and RF peak power limitations without special hardware modification in clinical scanners.

Characterizing the limits of MRI near metallic prostheses

PURPOSE

To characterize the fundamental limits of MRI near metallic implants on RF excitation, frequency encoding, and chemical shift-encoding water-fat separation.

METHODS

Multicomponent three-dimensional (3D) digital models of a total hip and a total knee replacement were used to construct material-specific susceptibility maps. The fundamental limits and spatial relationship of imaging near metallic prostheses were investigated as a function of distance from the prosthetic surface by calculating 3D field map perturbations using a well-validated k-space based dipole kernel.

RESULTS

Regions limited by the bandwidth of RF excitation overlap substantially with those fundamentally limited by frequency encoding. Rapid breakdown of water-fat separation occurs once the intravoxel off-resonance exceeds ∼6 ppm over a full range of fat fractions (0%-100%) and SNR (5-100).

CONCLUSION

Current 3D multispectral imaging methods would not benefit greatly from exciting spins beyond ±12 kHz despite the presence of signal that lies outside of this range from tissue directly adjacent to the metallic implants. Methods such as phase encoding in all three spatial dimensions are necessary to spatially resolve spins beyond an excitation bandwidth of ±12 kHz. The approach described in this study provides a benchmark for the capabilities of current imaging techniques to guide development of new MRI methods for imaging near metal.

Accelerating sequences in the presence of metal by exploiting the spatial distribution of off-resonance

PURPOSE

To demonstrate feasibility of exploiting the spatial distribution of off-resonance surrounding metallic implants for accelerating multispectral imaging techniques.

THEORY AND METHODS

Multispectral imaging (MSI) techniques perform time-consuming independent three-dimensional acquisitions with varying radio frequency offsets to address the extreme off-resonance from metallic implants. Each off-resonance bin provides a unique spatial sensitivity that is analogous to the sensitivity of a receiver coil and, therefore, provides a unique opportunity for acceleration. Fully sampled MSI was performed to demonstrate retrospective acceleration. A uniform sampling pattern across off-resonance bins was compared with several adaptive sampling strategies using a total hip replacement phantom. Monte Carlo simulations were performed to compare noise propagation of two of these strategies. With a total knee replacement phantom, positive and negative off-resonance bins were strategically sampled with respect to the B0 field to minimize aliasing. Reconstructions were performed with a parallel imaging framework to demonstrate retrospective acceleration.

RESULTS

An adaptive sampling scheme dramatically improved reconstruction quality, which was supported by the noise propagation analysis. Independent acceleration of negative and positive off-resonance bins demonstrated reduced overlapping of aliased signal to improve the reconstruction.

CONCLUSION

This work presents the feasibility of acceleration in the presence of metal by exploiting the spatial sensitivities of off-resonance bins.

T1-mapping in the heart: accuracy and precision

The longitudinal relaxation time constant (T1) of the myocardium is altered in various disease states due to increased water content or other changes to the local molecular environment. Changes in both native T1 and T1 following administration of gadolinium (Gd) based contrast agents are considered important biomarkers and multiple methods have been suggested for quantifying myocardial T1 in vivo. Characterization of the native T1 of myocardial tissue may be used to detect and assess various cardiomyopathies while measurement of T1 with extracellular Gd based contrast agents provides additional information about the extracellular volume (ECV) fraction. The latter is particularly valuable for more diffuse diseases that are more challenging to detect using conventional late gadolinium enhancement (LGE). Both T1 and ECV measures have been shown to have important prognostic significance. T1-mapping has the potential to detect and quantify diffuse fibrosis at an early stage provided that the measurements have adequate reproducibility. Inversion recovery methods such as MOLLI have excellent precision and are highly reproducible when using tightly controlled protocols. The MOLLI method is widely available and is relatively mature. The accuracy of inversion recovery techniques is affected significantly by magnetization transfer (MT). Despite this, the estimate of apparent T1 using inversion recovery is a sensitive measure, which has been demonstrated to be a useful tool in characterizing tissue and discriminating disease. Saturation recovery methods have the potential to provide a more accurate measurement of T1 that is less sensitive to MT as well as other factors. Saturation recovery techniques are, however, noisier and somewhat more artifact prone and have not demonstrated the same level of reproducibility at this point in time.This review article focuses on the technical aspects of key T1-mapping methods and imaging protocols and describes their limitations including the factors that influence their accuracy, precision, and reproducibility.

Parameter estimation approach to banding artifact reduction in balanced steady-state free precession

PURPOSE

The balanced steady-state free precession (bSSFP) pulse sequence has shown to be of great interest due to its high signal-to-noise ratio efficiency. However, bSSFP images often suffer from banding artifacts due to off-resonance effects, which we aim to minimize in this article.

METHODS

We present a general and fast two-step algorithm for 1) estimating the unknowns in the bSSFP signal model from multiple phase-cycled acquisitions, and 2) reconstructing band-free images. The first step, linearization for off-resonance estimation (LORE), solves the nonlinear problem approximately by a robust linear approach. The second step applies a Gauss-Newton algorithm, initialized by LORE, to minimize the nonlinear least squares criterion. We name the full algorithm LORE-GN.

RESULTS

We derive the Cramér-Rao bound, a theoretical lower bound of the variance for any unbiased estimator, and show that LORE-GN is statistically efficient. Furthermore, we show that simultaneous estimation of T1 and T2 from phase-cycled bSSFP is difficult, since the Cramér-Rao bound is high at common signal-to-noise ratio. Using simulated, phantom, and in vivo data, we illustrate the band-reduction capabilities of LORE-GN compared to other techniques, such as sum-of-squares.

CONCLUSION

Using LORE-GN we can successfully minimize banding artifacts in bSSFP.

Influence of Off-resonance in myocardial T1-mapping using SSFP based MOLLI method

BACKGROUND

Myocardial T1-mapping methods such as MOLLI use SSFP readout and are prone to frequency-dependent error in T1-measurement. A significant error in T1 may result at relatively small off-resonance frequencies that are well within the region without banding artifacts.

METHODS

The sensitivity of T1-estimates based on the SSFP based MOLLI sequence to errors in center frequency are calculated by means of a Bloch simulation and validated by phantom measurements. Typical off-resonance errors following local cardiac shimming are determined by field mapping at both 1.5 and 3.0T. In vivo examples demonstrate the artifactual appearance of T1-maps in the presence of off-resonance variation.

RESULTS

Off-resonance varied 61.8 ± 15.5 Hz (mean ± SD, n = 18) across the heart at 1.5T and 125.0 ± 40.6 Hz (mean ± SD, n = 18) at 3.0T. For T1 = 1000 ms, the variation in T1 due to off-resonance variation was approximately 20 ms at 62 Hz, and > 50 ms at 125 Hz.

CONCLUSIONS

Regional variations due to the inability to completely shim the B0-field variation around the heart appear as regional variation in T1, which is artifactual.

B0-informed variable density trajectory design for enhanced correction of off-resonance effects in parallel transmission

PURPOSE

To improve B1 and B0 inhomogeneity mitigation performance of spatially selective radio-frequency (RF) pulses in parallel transmission while decreasing RF pulse power. Further enhancement of off-resonance correction for rectilinear spoke-trajectory-based RF pulses with known residual geometric distortions after optimization.

METHODS

The appropriate definition of the target magnetization pattern is discussed regarding the maximum physical excitation resolution. Furthermore, a novel variable-density trajectory design is introduced, which subsamples accrued B0 phase error elevations in k-space. A simulation study (echo-planar and spiral 2DRF) at different off-resonance levels and pulse acceleration factors was pursued using data from a whole-body 2-channel parallel transmit 3T MRI system. The new trajectory design for echo-planar 2DRF was validated in human in-vivo experiments.

RESULTS

Proper target pattern definition can require spatial filtering, such that RF pulse optimization is prevented from lower excitation performance with significant higher RF power level. The new trajectory design proposed can considerably improve off-resonance compensation, while further reducing the RF power, e.g., 43% less RMSE with 79% less RF power for spoke based pulses.

CONCLUSION

The proposed methods offer significant improvements of the excitation performance (homogeneity and acceleration), while significantly decreasing the RF power. Furthermore, single-channel transmit RF pulse performance can be similarly improved.

Off-resonance magnetoresistance spike in irradiated ultraclean 2D electron systems

: We report on the theoretical studies of a recently discovered strong radiation-induced magnetoresistance spike obtained in ultraclean two-dimensional electron systems at low temperatures. The most striking feature of this spike is that it shows up on the second harmonic of the cyclotron resonance. We apply the radiation-driven electron orbits model in the ultraclean scenario. Accordingly, we calculate the new average advanced distance by the electron in a scattering event which will define the unexpected resonance spike position. Calculated results are in good agreement with experiments.

Off-resonance insensitive complementary SPAtial Modulation of Magnetization (ORI-CSPAMM) for quantification of left ventricular twist

PURPOSE

To evaluate Off Resonance Insensitive Complementary SPAtial Modulation of Magnetization (ORI-CSPAMM) and Fourier Analysis of STimulated echoes (FAST) for the quantification of left ventricular (LV) systolic and diastolic function and compare it with the previously validated FAST+SPAMM technique.

MATERIALS AND METHODS

LV short-axis tagged images were acquired with ORI-CSPAMM and SPAMM in healthy volunteers (n = 13). The FAST method was used to automatically estimate LV systolic and diastolic twist parameters from rotation of the stimulated echo and stimulated anti-echo about the middle of k-space subsequent to ∼3 min of user interaction.

RESULTS

There was no significant difference between measures obtained for FAST+ORI-CSPAMM and FAST+SPAMM for mean peak twist (12.9 ± 3.4° versus 11.9 ± 4.0°; P = 0.4), torsion (3.3 ± 0.9°/cm versus 2.9 ± 1.0°/cm, P = 0.3), circumferential-longitudinal shear angle (9.1 ± 3.0° versus 8.2 ± 3.4°, P = 0.3), twisting rate (79.6 ± 20.2°/s versus 68.2 ± 23.4°/s, P = 0.1), untwisting rate (-117.5 ± 31.4°/s versus -106.6 ± 32.4°/s, P = 0.3), normalized untwisting rate (-9.3 ± 2.0/s versus -9.9 ± 4.4/s, P = 0.7), and time of peak twist (281 ± 18 ms versus 293 ± 25 ms, P = 0.04). FAST+ORI-CSPAMM also provided measures of duration of untwisting (148 ± 21 ms) and the ratio of rapid untwisting to peak twist (0.9 ± 0.3). Bland-Altman analysis of FAST+ORI-CSPAMM and FAST+SPAMM twist data demonstrates excellent agreement with a bias of -0.1° and 95% confidence intervals of (-1.0°, 3.2°).

CONCLUSION

FAST+ORI-CSPAMM is a semi-automated method that provides a quick and quantitative assessment of LV systolic and diastolic twist and torsion. ORI-CSPAMM corrects off-resonance accrued during tagging preparation and readout and visibly removes chemical shift from the tagging pattern, which confers greater robustness to the derived quantitative measures.

Robust selective signal suppression using binomial off-resonant rectangular (BORR) pulses

PURPOSE

To study the selective signal suppression capability of a binomial off-resonant rectangular (BORR) radiofrequency pulse method.

MATERIALS AND METHODS

The BORR pulse consists of two consecutive rectangular pulses with a phase difference of π. The exact solution of the Bloch equations was used to simulate its frequency response. The BORR pulse was implemented in a gradient echo sequence and tested on phantoms, the knee, and the breast.

RESULTS

The frequency response of the BORR pulse acquired on the phantom confirmed the theory. Broad suppression bands ensured high suppression efficiency and robustness in both in vitro and in vivo scans compared with other saturation pulses.

CONCLUSION

The BORR pulse method provides a simple, efficient, and robust selective signal suppression alternative for three-dimensional short TR (repetition time) imaging.

The Importance of k-Space Trajectory on Off-Resonance Artifact in Segmented Echo-Planar Imaging

Segmented interleaved echo planar imaging (EPI) is a highly efficient data acquisition technique; however, EPI is sensitive to artifacts from off-resonance spins. The choice of k-space trajectories is important in determining how off-resonance spins contribute to image artifacts. Top-down and center-out trajectories are theoretically analyzed, simulated, implemented, and tested in phantom and volunteer experiments. Theoretical results show off-resonance artifact manifests as a simple positional shift for the top-down trajectory, while for the center-out trajectory off-resonance artifact manifests as a splitting of the object, which entails both shift and blurring. These results were validated using simulation and phantom scan data where a frequency-offset was introduced ranging from -300 Hz to +300 Hz. As predicted by the theoretical results, inferior image quality was observed for the center-out trajectory in a single volunteer. Off-resonance produces more severe and complex artifacts with the center-out trajectory than the top-down trajectory.

Improved slice-selective adiabatic excitation

PURPOSE

The purpose of this work is to design an improved Slice-selective Tunable-flip AdiaBatic Low peak-power Excitation (STABLE) pulse with shorter duration and increased off-resonance immunity to make it suitable for use in a greater range of applications and at higher field strengths. An additional aim is to design a variant of this pulse to achieve B1 -insensitive, fat-suppressed excitation.

METHODS

The adiabatic SLR algorithm was used to generate a more uniform spectral pulse envelope for this improved radiofrequency pulse for adiabatic slice-selective excitation, called STABLE-2. Pulse parameters were adjusted to design a version of STABLE-2 with a spectral null centered on lipids.

RESULTS

In vivo images obtained of the human brain at 3 and 7 T demonstrate that STABLE-2 provides robust, uniform, slice-selective excitation over a range of B1 values. Phantom and in vivo knee images obtained at 3 T demonstrate the effectiveness of STABLE-2 for fat suppression.

CONCLUSIONS

STABLE-2 achieves B1 -insensitive slice-selective excitation while providing greater off-resonance immunity and a shorter pulse duration, when compared to the original STABLE pulse. In particular, the 9.8-ms STABLE-2 pulse provides slice selectivity over 120 Hz whereas the 21-ms STABLE pulse is limited to 80 Hz off-resonance. B1 -Insensitive fat-suppressed excitation may also be achieved by using a variant of this pulse.

Reducing artifacts in one-dimensional Fourier velocity encoding for fast and pulsatile flow

When evaluating the severity of valvular stenosis, the peak velocity of the blood flow is routinely used to estimate the transvalvular pressure gradient. One-dimensional Fourier velocity encoding effectively detects the peak velocity with an ungated time series of spatially resolved velocity spectra in real time. However, measurement accuracy can be degraded by the pulsatile and turbulent nature of stenotic flow and the existence of spatially varying off-resonance. In this work, we investigate the feasibility of improving the peak velocity detection capability of one-dimensional Fourier velocity encoding for stenotic flow using a novel echo-shifted interleaved readout combined with a variable-density circular k-space trajectory. The shorter echo and readout times of the echo-shifted interleaved acquisitions are designed to reduce sensitivity to off-resonance. Preliminary results from limited phantom and in vivo results also indicate that some artifacts from pulsatile flow appear to be suppressed when using this trajectory compared to conventional single-shot readouts, suggesting that peak velocity detection may be improved. The efficiency of the new trajectory improves the temporal and spatial resolutions. To realize the proposed readout, a novel multipoint-traversing algorithm is introduced for flexible and automated gradient-waveform design.

Turboprop+: enhanced Turboprop diffusion-weighted imaging with a new phase correction

Faster periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) diffusion-weighted imaging acquisitions, such as Turboprop and X-prop, remain subject to phase errors inherent to a gradient echo readout, which ultimately limits the applied turbo factor (number of gradient echoes between each pair of radiofrequency refocusing pulses) and, thus, scan time reductions. This study introduces a new phase correction to Turboprop, called Turboprop+. This technique employs calibration blades, which generate 2-D phase error maps and are rotated in accordance with the data blades, to correct phase errors arising from off-resonance and system imperfections. The results demonstrate that with a small increase in scan time for collecting calibration blades, Turboprop+ had a superior immunity to the off-resonance-related artifacts when compared to standard Turboprop and recently proposed X-prop with the high turbo factor (turbo factor = 7). Thus, low specific absorption rate and short scan time can be achieved in Turboprop+ using a high turbo factor, whereas off-resonance related artifacts are minimized.