MaxGIRF: Image reconstruction incorporating concomitant field and gradient impulse response function effects

Low-Field, Low-Cost, Point-of-Care Magnetic Resonance Imaging

Low-field magnetic resonance imaging (MRI) has recently experienced a renaissance that is largely attributable to the numerous technological advancements made in MRI, including optimized pulse sequences, parallel receive and compressed sensing, improved calibrations and reconstruction algorithms, and the adoption of machine learning for image postprocessing. This new attention on low-field MRI originates from a lack of accessibility to traditional MRI and the need for affordable imaging. Low-field MRI provides a viable option due to its lack of reliance on radio-frequency shielding rooms, expensive liquid helium, and cryogen quench pipes. Moreover, its relatively small size and weight allow for easy and affordable installation in most settings. Rather than replacing conventional MRI, low-field MRI will provide new opportunities for imaging both in developing and developed countries. This article discusses the history of low-field MRI, low-field MRI hardware and software, current devices on the market, advantages and disadvantages, and low-field MRI's global potential.

Ultrafast T2 and T2* mapping using single-shot spatiotemporally encoded MRI with reduced field of view and spiral out-in-out-in trajectory

BACKGROUND

T2 and T2* mapping are crucial components of quantitative magnetic resonance imaging, offering valuable insights into tissue characteristics and pathology. Single-shot methods can achieve ultrafast T2 or T2* mapping by acquiring multiple readout echo trains. However, the extended echo trains pose challenges, such as compromised image quality and diminished quantification accuracy.

PURPOSE

In this study, we develop a single-shot method for ultrafast T2 and T2* mapping with reduced echo train length.

METHODS

The proposed method is based on ultrafast single-shot spatiotemporally encoded (SPEN) MRI combined with reduced field of view (FOV) and spiral out-in-out-in (OIOI) trajectory. Specifically, a biaxial SPEN excitation scheme was employed to excite the spin signal into the spatiotemporal encoding domain. The OIOI trajectory with high acquisition efficiency was employed to acquire signals within targeted reduced FOV. Through non-Cartesian super-resolved (SR) reconstruction, 12 aliasing-free images with different echo times were obtained within 150 ms. These images were subsequently fitted to generate T2 or T2* mapping simultaneously using a derived model.

RESULTS

Accurate and co-registered T2 and T2* maps were generated, closely resembling the reference maps. Numerical simulations demonstrated substantial consistency (R2 > 0.99) with the ground truth values. A mean difference of 0.6% and 1.7% was observed in T2 and T2*, respectively, in in vivo rat brain experiments compared to the reference. Moreover, the proposed method successfully obtained T2 and T2* mappings of rat kidney in free-breathing mode, demonstrating its superiority over multishot methods lacking respiratory navigation.

CONCLUSIONS

The results suggest that the proposed method can achieve ultrafast and accurate T2 and T2* mapping, potentially facilitating the application of T2 and T2* mapping in scenarios requiring high temporal resolution.

Trajectory correction enables free-running chemical shift encoded imaging for accurate cardiac proton-density fat fraction quantification at 3T

BACKGROUND

Metabolic diseases can negatively alter epicardial fat accumulation and composition, which can be probed using quantitative cardiac chemical shift encoded (CSE) cardiovascular magnetic resonance (CMR) by mapping proton-density fat fraction (PDFF). To obtain motion-resolved high-resolution PDFF maps, we proposed a free-running cardiac CSE-CMR framework at 3T. To employ faster bipolar readout gradients, a correction for gradient imperfections was added using the gradient impulse response function (GIRF) and evaluated on intermediate images and PDFF quantification.

METHODS

Ten minutes free-running cardiac 3D radial CSE-CMR acquisitions were compared in vitro and in vivo at 3T. Monopolar and bipolar readout gradient schemes provided 8 echoes (TE1/ΔTE = 1.16/1.96 ms) and 13 echoes (TE1/ΔTE = 1.12/1.07 ms), respectively. Bipolar-gradient free-running cardiac fat and water images and PDFF maps were reconstructed with or without GIRF correction. PDFF values were evaluated in silico, in vitro on a fat/water phantom, and in vivo in 10 healthy volunteers and 3 diabetic patients.

RESULTS

In monopolar mode, fat-water swaps were demonstrated in silico and confirmed in vitro. Using bipolar readout gradients, PDFF quantification was reliable and accurate with GIRF correction with a mean bias of 0.03% in silico and 0.36% in vitro while it suffered from artifacts without correction, leading to a PDFF bias of 4.9% in vitro and swaps in vivo. Using bipolar readout gradients, in vivo PDFF of epicardial adipose tissue was significantly lower compared to subcutaneous fat (80.4 ± 7.1% vs 92.5 ± 4.3%, P < 0.0001).

CONCLUSIONS

Aiming for an accurate PDFF quantification, high-resolution free-running cardiac CSE-MRI imaging proved to benefit from bipolar echoes with k-space trajectory correction at 3T. This free-breathing acquisition framework enables to investigate epicardial adipose tissue PDFF in metabolic diseases.

Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging

PURPOSE

To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques.

METHODS

Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method.

RESULTS

The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m.

CONCLUSION

A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

Beat phenomena of oscillating readouts

PURPOSE

To demonstrate slowly varying, erroneous magnetic field gradients for oscillating readouts due to the mechanically resonant behavior of gradient systems.

METHODS

Projections of a static phantom were acquired using a one-dimensional (1D) EPI sequence with varying EPI frequencies ranging from 1121 to 1580 Hz on clinical 3T systems (30 mT/m, 200 T/m/s). Phase due to static B0 inhomogeneities was eliminated by a complex division of two separate scans with different polarities of the EPI readout. The temporal evolution of phase was evaluated and related to the mechanical resonances of the gradient systems derived from the gradient modulation transfer function. Additionally, the impact of temporally varying mechanical resonance effects on EPI was evaluated using an echo-planar spectroscopic imaging sequence.

RESULTS

A beat phenomenon resulting in a slowly varying phase was observed. Its temporal frequency was given by the difference between the EPI frequency and the mechanical resonance frequency of the activated gradient axis. The maximum erroneous, oscillating phase during phase encoding was ±0.5 rad for an EPI frequency of 1281 Hz. Echo-planar spectroscopic imaging images showed the resulting time-dependent stretching/compression of the FOV.

CONCLUSION

Oscillating readouts such as those used in EPI can result in low-frequency, erroneous phase contributions, which are explained by the beat phenomenon. Therefore, EPI phase-correction approaches may need to include beat effects for accurate image reconstruction.

Real-time water/fat imaging at 0.55T with spiral out-in-out-in sampling

PURPOSE

To develop an efficient and flexible water/fat separated real-time MRI (RT-MRI) method using spiral out-in-out-in (OIOI) sampling and balanced SSFP (bSSFP) at 0.55T.

METHODS

A bSSFP sequence with golden-angle spiral OIOI readout was developed, capturing three echoes to allow water/fat separation. A low-latency reconstruction that combines all echoes was available for online visualization. An offline reconstruction provided water and fat RT-MRI in two steps: (1) image reconstruction with spatiotemporally constrained reconstruction (STCR) and (2) water/fat separation with hierarchical iterative decomposition of water and fat with echo asymmetry and least-squares estimation (HIDEAL). In healthy volunteers, spiral OIOI was acquired in the wrist during a radial-to-ulnar deviation maneuver, in the heart without breath-hold and cardiac gating, and in the lower abdomen during free-breathing for visualizing small bowel motility.

RESULTS

We demonstrate successful water/fat separated RT-MRI for all tested applications. In the wrist, resulting images provided clear depiction of ligament gaps and their interactions during the radial-to-ulnar deviation maneuver. In the heart, water/fat RT-MRI depicted epicardial fat, provided improved delineation of epicardial coronary arteries, and provided high blood-myocardial contrast for ventricular function assessment. In the abdomen, water-only RT-MRI captured small bowel mobility clearly with improved water-fat contrast.

CONCLUSIONS

We have demonstrated a novel and flexible bSSFP spiral OIOI sequence at 0.55T that can provide water/fat separated RT-MRI with a variety of application-specific temporal resolution and spatial resolution requirements.

Cardiac MRI at Low Field Strengths

Cardiac MR imaging is well established for assessment of cardiovascular structure and function, myocardial scar, quantitative flow, parametric mapping, and myocardial perfusion. Despite the clear evidence supporting the use of cardiac MRI for a wide range of indications, it is underutilized clinically. Recent developments in low-field MRI technology, including modern data acquisition and image reconstruction methods, are enabling high-quality low-field imaging that may improve the cost-benefit ratio for cardiac MRI. Studies to-date confirm that low-field MRI offers high measurement concordance and consistent interpretation with clinical imaging for several routine sequences. Moreover, low-field MRI may enable specific new clinical opportunities for cardiac imaging such as imaging near metal implants, MRI-guided interventions, combined cardiopulmonary assessment, and imaging of patients with severe obesity. In this review, we discuss the recent progress in low-field cardiac MRI with a focus on technical developments and early clinical validation studies. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

New clinical opportunities of low-field MRI: heart, lung, body, and musculoskeletal

Contemporary whole-body low-field MRI scanners (< 1 T) present new and exciting opportunities for improved body imaging. The fundamental reason is that the reduced off-resonance and reduced SAR provide substantially increased flexibility in the design of MRI pulse sequences. Promising body applications include lung parenchyma imaging, imaging adjacent to metallic implants, cardiac imaging, and dynamic imaging in general. The lower cost of such systems may make MRI favorable for screening high-risk populations and population health research, and the more open configurations allowed may prove favorable for obese subjects and for pregnant women. This article summarizes promising body applications for contemporary whole-body low-field MRI systems, with a focus on new platforms developed within the past 5 years. This is an active area of research, and one can expect many improvements as MRI physicists fully explore the landscape of pulse sequences that are feasible, and as clinicians apply these to patient populations.

Low-field MRI: A report on the 2022 ISMRM workshop

In March 2022, the first ISMRM Workshop on Low-Field MRI was held virtually. The goals of this workshop were to discuss recent low field MRI technology including hardware and software developments, novel methodology, new contrast mechanisms, as well as the clinical translation and dissemination of these systems. The virtual Workshop was attended by 368 registrants from 24 countries, and included 34 invited talks, 100 abstract presentations, 2 panel discussions, and 2 live scanner demonstrations. Here, we report on the scientific content of the Workshop and identify the key themes that emerged. The subject matter of the Workshop reflected the ongoing developments of low-field MRI as an accessible imaging modality that may expand the usage of MRI through cost reduction, portability, and ease of installation. Many talks in this Workshop addressed the use of computational power, efficient acquisitions, and contemporary hardware to overcome the SNR limitations associated with low field strength. Participants discussed the selection of appropriate clinical applications that leverage the unique capabilities of low-field MRI within traditional radiology practices, other point-of-care settings, and the broader community. The notion of "image quality" versus "information content" was also discussed, as images from low-field portable systems that are purpose-built for clinical decision-making may not replicate the current standard of clinical imaging. Speakers also described technical challenges and infrastructure challenges related to portability and widespread dissemination, and speculated about future directions for the field to improve the technology and establish clinical value.

Efficient imaging using spiral acquisitions on a portable 50-mT MR head scanner

Ultralow-field (ULF) magnetic resonance imaging (MRI) can suffer from inferior image quality because of low signal-to-noise ratio (SNR). As an efficient way to cover the k-space, the spiral acquisition technique has shown great potential in improving imaging SNR efficiency at ULF. The current study aimed to address the problems of noise and blurring cancelation in the ULF case with spiral trajectory, and we proposed a spiral-out sequence for brain imaging using a portable 50-mT MRI system. The proposed sequence consisted of three modules: noise calibration, field map acquisition, and imaging. In the calibration step, transfer coefficients were obtained between signals from primary and noise-pick-up coils to perform electromagnetic interference (EMI) cancelation. Embedded field map acquisition was performed to correct accumulated phase error due to main field inhomogeneity. Considering imaging SNR, a lower bandwidth for data sampling was adopted in the sequence design because the 50-mT scanner is in a low SNR regime. Image reconstruction proceeded with sampled data by leveraging system imperfections, such as gradient delays and concomitant fields. The proposed method can provide images with higher SNR efficiency compared with its Cartesian counterparts. An improvement in temporal SNR of approximately 23%-44% was measured via phantom and in vivo experiments. Distortion-free images with a noise suppression rate of nearly 80% were obtained by the proposed technique. A comparison was also made with a state-of-the-art EMI cancelation algorithm used in the ULF-MRI system. SNR efficiency-enhanced spiral acquisitions were investigated for ULF-MR scanners and future studies could focus on various image contrasts based on our proposed approach to widen ULF applications.

An evolution of low-field strength MRI

The paper describes the evolution of low-field MRI from the very early pioneering days in the late 70 s until today. It is not meant to give a comprehensive historical account of the development of MRI, but rather to highlight the different research environments then and now. In the early 90 s, when low-field systems below 1.5 T essentially vanished, there were just no reasonable means available to make up for the factor of roughly three in signal-to-noise-ratio (SNR) between 0.5 and 1.5 T. This has drastically changed. Improvements in hardware-closed Helium-free magnets, RF receiver systems and especially much faster gradients, much more flexible sampling schemes including parallel imaging and compressed sensing and especially the use of AI at all stages of the imaging process have made low-field MRI a clinically viable supplement to conventional MRI. Ultralow-field MRI with magnets around 0.05 T are also back and constitute a bold and courageous endeavor to bring MRI to communities, which have neither the means nor the infrastructure to sustain a current standard of care MRI.

Concomitant field compensation of spiral turbo spin-echo at 0.55 T

OBJECTIVE

Diagnostic-quality neuroimaging methods are vital for widespread clinical adoption of low field MRI. Spiral imaging is an efficient acquisition method that can mitigate the reduced signal-to-noise ratio at lower field strengths. As concomitant field artifacts are worse at lower field, we propose a generalizable quadratic gradient-field nulling as an echo-to-echo compensation and apply it to spiral TSE at 0.55 T.

MATERIALS AND METHODS

A spiral in-out TSE acquisition was developed with a compensation for concomitant field variation between spiral interleaves, by adding bipolar gradients around each readout to minimize phase differences at each refocusing pulse. Simulations were performed to characterize concomitant field compensation approaches. We demonstrate our proposed compensation method in phantoms and (n = 8) healthy volunteers at 0.55 T.

RESULTS

Spiral read-outs with integrated spoiling demonstrated strong concomitant field artifacts but were mitigated using the echo-to-echo compensation. Simulations predicted a decrease of concomitant field phase RMSE between echoes of 42% using the proposed compensation. Spiral TSE improved SNR by 17.2 ± 2.3% compared to reference Cartesian acquisition.

DISCUSSION

We demonstrated a generalizable approach to mitigate concomitant field artifacts for spiral TSE acquisitions via the addition of quadratic-nulling gradients, which can potentially improve neuroimaging at low-field through increased acquisition efficiency.

Low-field 0.55 T MRI evaluation of the fetus

Fetal magnetic resonance imaging (MRI) is an important adjunct modality for the evaluation of fetal abnormalities. Recently, low-field MRI systems at 0.55 Tesla have become available which can produce images on par with 1.5 Tesla systems but with lower power deposition, acoustic noise, and artifact. In this article, we describe a technical innovation using low-field MRI to perform diagnostic quality fetal MRI.

Body composition profiling at 0.55T: Feasibility and precision

PURPOSE

Body composition MRI captures the distribution of fat and lean tissues throughout the body, and provides valuable biomarkers of obesity, metabolic disease, and muscle disorders, as well as risk assessment. Highly reproducible protocols have been developed for 1.5T and 3T MRI. The purpose of this work was to demonstrate the feasibility and test-retest repeatability of MRI body composition profiling on a 0.55T whole-body system.

METHODS

Healthy adult volunteers were scanned on a whole-body 0.55T MRI system using the integrated body RF coil. Experiments were performed to refine parameter settings such as TEs, resolution, flip angle, bandwidth, acceleration, and oversampling factors. The final protocol was evaluated using a test-retest study with subject removal and replacement in 10 adult volunteers (5 M/5F, age 25-60, body mass index 20-30).

RESULTS

Compared to 1.5T and 3T, the optimal flip angle at 0.55T was higher (15°), due to the shorter T1 times, and the optimal echo spacing was larger, due to smaller chemical shift between water and fat. Overall image quality was comparable to conventional field strengths, with no significant issues with fat/water swapping or inadequate SNR. Repeatability coefficient of visceral fat, subcutaneous fat, total thigh muscle volume, muscle fat infiltration, and liver fat were 11.8 cL (2.2%), 46.9 cL (1.9%), 14.6 cL (0.5%), 0.1 pp (2%), and 0.2 pp (5%), respectively (coefficient of variation in parenthesis).

CONCLUSIONS

We demonstrate that 0.55T body composition MRI is feasible and present optimized scan parameters. The resulting images provide satisfactory quality for automated post-processing and produce repeatable results.

Contrast-optimal simultaneous multi-slice bSSFP cine cardiac imaging at 0.55 T

PURPOSE

To determine if contemporary 0.55 T MRI supports the use of contrast-optimal flip angles (FA) for simultaneous multi-slice (SMS) balanced SSFP (bSSFP) cardiac function assessment, which is impractical at conventional field strengths because of excessive SAR and/or banding artifacts.

METHODS

Blipped-CAIPI bSSFP was combined with spiral sampling for ventricular function assessment at 0.55 T. Cine movies with single band and SMS factors of 2 and 3 (SMS 2 and 3), and FA ranging from 60° to 160°, were acquired in seven healthy volunteers. Left ventricular blood and myocardial signal intensity (SI) normalized by background noise and blood-myocardium contrast were measured and compared across acquisition settings.

RESULTS

Myocardial SI was slightly higher in single band than in SMS and decreased with an increasing FA. Blood SI increased as the FA increased for single band, and increment was small for FA ≥120°. Blood SI for SMS 2 and 3 increased with an increasing FA up to ∼100°. Blood-myocardium contrast increased with an increasing FA for single band, peaked at FA = 160° (systole: 28.43, diastole: 29.15), attributed mainly to reduced myocardial SI when FA ≥120°. For SMS 2, contrast peaked at 120° (systole: 21.43, diastole: 19.85). For SMS 3, contrast peaked at 120° in systole (16.62) and 100° in diastole (19.04).

CONCLUSIONS

Contemporary 0.55 T MR scanners equipped with high-performance gradient systems allow the use of contrast-optimal FA for SMS accelerated bSSFP cine examinations without compromising image quality. The contrast-optimal FA was found to be 140° to 160° for single band and 100° to 120° for SMS 2 and 3.

Calibration of concomitant field offsets using phase contrast MRI for asymmetric gradient coils

PURPOSE

Asymmetric gradient coils introduce zeroth- and first-order concomitant field terms, in addition to higher-order terms common to both asymmetric and symmetric gradients. Salient to compensation strategies is the accurate calibration of the concomitant field spatial offset parameters for asymmetric coils. A method that allows for one-time calibration of the offset parameters is described.

THEORY AND METHODS

A modified phase contrast pulse sequence with single-sided bipolar flow encoding is proposed to calibrate the offsets for asymmetric, transverse gradient coils. By fitting the measured phase offsets to different gradient amplitudes, the spatial offsets were calculated by fitting the phase variation. This was used for calibrating real-time pre-emphasis compensation of the zeroth- and first-order concomitant fields.

RESULTS

Image quality improvement with the proposed corrections was demonstrated in phantom and healthy volunteers with non-Cartesian and Cartesian trajectory acquisitions. Concomitant field compensation using the calibrated offsets resulted in a residual phase error <3% at the highest gradient amplitude and demonstrated substantial reduction of image blur and slice position/selection artifacts.

CONCLUSIONS

The proposed implementation provides an accurate method for calibrating spatial offsets that can be used for real-time concomitant field compensation of zeroth and first-order terms, substantially reducing artifacts without retrospective correction or sequence specific waveform modifications.