Integrated image reconstruction and gradient nonlinearity correction

Open-source hardware and software for the measurement, characterization, reporting, and correction of geometric distortion in MRI

BACKGROUND

Geometric distortion is a serious problem in MRI, particularly in MRI guided therapy. A lack of affordable and adaptable tools in this area limits research progress and harmonized quality assurance.

PURPOSE

To develop and test a suite of open-source hardware and software tools for the measurement, characterization, reporting, and correction of geometric distortion in MRI.

METHODS

An open-source python library was developed, comprising modules for parametric phantom design, data processing, spherical harmonics, distortion correction, and interactive reporting. The code was used to design and manufacture a distortion phantom consisting of 618 oil filled markers covering a sphere of radius 150 mm. This phantom was imaged on a CT scanner and a novel split-bore 1.0 T MRI magnet. The CT images provide distortion-free dataset. These data were used to test all modules of the open-source software.

RESULTS

All markers were successfully extracted from all images. The distorted MRI markers were mapped to undistorted CT data using an iterative search approach. Spherical harmonics reconstructed the fitted gradient data to 1.0 ± 0.6% of the input data. High resolution data were reconstructed via spherical harmonics and used to generate an interactive report. Finally, distortion correction on an independent data set reduced distortion inside the DSV from 5.5 ± 3.1 to 1.6 ± 0.8 mm.

CONCLUSION

Open-source hardware and software for the measurement, characterization, reporting, and correction of geometric distortion in MRI have been developed. The utility of these tools has been demonstrated via their application on a novel 1.0 T split bore magnet.

Nonlinear Gradient Field Mapping Using a Spherical Grid Phantom for 3 and 7 Tesla MR Imaging Systems Equipped with High-performance Gradient Coils

Recent high-performance gradient coils are fabricated mainly at the expense of spatial linearity. In this study, we measured the spatial nonlinearity of the magnetic field generated by the gradient coils of two MRI systems with high-performance gradient coils. The nonlinearity of the gradient fields was measured using 3D gradient echo sequences and a spherical phantom with a built-in lattice structure. The spatial variation of the gradient field was approximated to the 3rd order polynomials. The coefficients of the polynomials were calculated using the steepest descent method. The geometric distortion of the acquired 3D MR images was corrected using the polynomials and compared with the 3D images corrected using the harmonic functions provided by the MRI venders. As a result, it was found that the nonlinearity correction formulae provided by the vendors were insufficient and needed to be verified or corrected using a geometric phantom such as used in this study.

Optimization of flip angle and radiofrequency pulse phase to maximize steady-state magnetization in three-dimensional missing pulse steady-state free precession

Missing pulse (MP) steady-state free precession (SSFP) is a magnetic resonance imaging (MRI) pulse sequence that is highly tolerant to the magnetic field inhomogeneity. In this study, optimal flip angle and radiofrequency (RF) phase scheduling in three-dimensional (3D) MP-SSFP is introduced to maximize the steady-state magnetization while keeping broadband excitation to cover widely distributed frequencies generated by inhomogeneous magnetic fields. Numerical optimization based on extended phase graph (EPG) simulation was performed to maximize the MP-SSFP steady-state magnetization. To limit the specific absorption rate (SAR) associated with the broadband excitation in 3D MP-SSFP, SAR constraint was introduced in the numerical optimization. Optimized flip angle and RF phase settings were experimentally tested by introducing a linear inhomogeneous magnetic field in a range of 10-20 mT/m and using a phantom with known T1/T2 relaxation and diffusion parameters at 3 T. The experimental results were validated through comparisons with EPG simulation. Image contrasts and molecular diffusion effects were investigated in in vivo human brain imaging with 3D MP-SSFP with the optimal flip angle and RF phase settings. In the phantom measurements, the optimal flip angle and RF phase settings improved the MP-SSFP steady-state magnetization/signal-to-noise ratio by up to 41% under the fixed SAR conditions, which matched well with EPG simulation results. In vivo brain imaging with the optimal RF pulse settings provided T2-like image contrasts. Diffusion effects were relatively minor with the linear inhomogeneous field of 10-20 mT/m for white and gray matter, but cerebrospinal fluid showed conspicuous signal intensity attenuation as the linear inhomogeneous field increased. Numerical optimization achieved significant improvement in the steady-state magnetization in MP-SSFP compared with the RF pulse settings used in previous studies. The proposed flip angle and RF phase optimization is promising to improve 3D MP-SSFP image quality for MRI in inhomogeneous magnetic fields.

Model-based iterative reconstruction for direct imaging with point spread function encoded echo planar MRI

BACKGROUND

Echo planar imaging (EPI) is a fast measurement technique commonly used in magnetic resonance imaging (MRI), but is highly sensitive to measurement non-idealities in reconstruction. Point spread function (PSF)-encoded EPI is a multi-shot strategy which alleviates distortion, but acquisition of encodings suitable for direct distortion-free imaging prolongs scan time. In this work, a model-based iterative reconstruction (MBIR) framework is introduced for direct imaging with PSF-EPI to improve image quality and acceleration potential.

METHODS

An MBIR platform was developed for accelerated PSF-EPI. The reconstruction utilizes a subspace representation, is regularized to promote local low-rankedness (LLR), and uses variable splitting for efficient iteration. Comparisons were made against standard reconstructions from prospectively accelerated PSF-EPI data and with retrospective subsampling. Exploring aggressive partial Fourier acceleration of the PSF-encoding dimension, additional comparisons were made against an extension of Homodyne to direct PSF-EPI in numerical experiments. A neuroradiologists' assessment was completed comparing images reconstructed with MBIR from retrospectively truncated data directly against images obtained with standard reconstructions from non-truncated datasets.

RESULTS

Image quality results were consistently superior for MBIR relative to standard and Homodyne reconstructions. As the MBIR signal model and reconstruction allow for arbitrary sampling of the PSF space, random sampling of the PSF-encoding dimension was also demonstrated, with quantitative assessments indicating best performance achieved through nonuniform PSF sampling combined with partial Fourier. With retrospective subsampling, MBIR reconstructs high-quality images from sub-minute scan datasets. MBIR was shown to be superior in a neuroradiologists' assessment with respect to three of five performance criteria, with equivalence for the remaining two.

CONCLUSIONS

A novel image reconstruction framework is introduced for direct imaging with PSF-EPI, enabling arbitrary PSF space sampling and reconstruction of diagnostic-quality images from highly accelerated PSF-encoded EPI data.

Distortion-corrected image reconstruction with deep learning on an MRI-Linac

PURPOSE

MRI is increasingly utilized for image-guided radiotherapy due to its outstanding soft-tissue contrast and lack of ionizing radiation. However, geometric distortions caused by gradient nonlinearities (GNLs) limit anatomical accuracy, potentially compromising the quality of tumor treatments. In addition, slow MR acquisition and reconstruction limit the potential for effective image guidance. Here, we demonstrate a deep learning-based method that rapidly reconstructs distortion-corrected images from raw k-space data for MR-guided radiotherapy applications.

METHODS

We leverage recent advances in interpretable unrolling networks to develop a Distortion-Corrected Reconstruction Network (DCReconNet) that applies convolutional neural networks (CNNs) to learn effective regularizations and nonuniform fast Fourier transforms for GNL-encoding. DCReconNet was trained on a public MR brain dataset from 11 healthy volunteers for fully sampled and accelerated techniques, including parallel imaging (PI) and compressed sensing (CS). The performance of DCReconNet was tested on phantom, brain, pelvis, and lung images acquired on a 1.0T MRI-Linac. The DCReconNet, CS-, PI-and UNet-based reconstructed image quality was measured by structural similarity (SSIM) and RMS error (RMSE) for numerical comparisons. The computation time and residual distortion for each method were also reported.

RESULTS

Imaging results demonstrated that DCReconNet better preserves image structures compared to CS- and PI-based reconstruction methods. DCReconNet resulted in the highest SSIM (0.95 median value) and lowest RMSE (<0.04) on simulated brain images with four times acceleration. DCReconNet is over 10-times faster than iterative, regularized reconstruction methods.

CONCLUSIONS

DCReconNet provides fast and geometrically accurate image reconstruction and has the potential for MRI-guided radiotherapy applications.

Rapid distortion correction enables accurate magnetic resonance imaging-guided real-time adaptive radiotherapy

Background and purpose

Magnetic resonance imaging (MRI)-Linac systems combine simultaneous MRI with radiation delivery, allowing treatments to be guided by anatomically detailed, real-time images. However, MRI can be degraded by geometric distortions that cause uncertainty between imaged and actual anatomy. In this work, we develop and integrate a real-time distortion correction method that enables accurate real-time adaptive radiotherapy.

Materials and methods

The method was based on the pre-treatment calculation of distortion and the rapid correction of intrafraction images. A motion phantom was set up in an MRI-Linac at isocentre (P0 ), the edge (P 1) and just outside (P 2) the imaging volume. The target was irradiated and tracked during real-time adaptive radiotherapy with and without the distortion correction. The geometric tracking error and latency were derived from the measurements of the beam and target positions in the EPID images.

Results

Without distortion correction, the mean geometric tracking error was 1.3 mm at P 1 and 3.1 mm at P 2. When distortion correction was applied, the error was reduced to 1.0 mm at P 1 and 1.1 mm at P 2. The corrected error was similar to an error of 0.9 mm at P0 where the target was unaffected by distortion indicating that this method has accurately accounted for distortion during tracking. The latency was 319 ± 12 ms without distortion correction and 335 ± 34 ms with distortion correction.

Conclusions

We have demonstrated a real-time distortion correction method that maintains accurate radiation delivery to the target, even at treatment locations with large distortion.

Increased brain volumetric measurement precision from multi-site 3D T1-weighted 3 T magnetic resonance imaging by correcting geometric distortions

PURPOSE

Magnetic resonance imaging (MRI) scanner-specific geometric distortions may contribute to scanner induced variability and decrease volumetric measurement precision for multi-site studies. The purpose of this study was to determine whether geometric distortion correction increases the precision of brain volumetric measurements in a multi-site multi-scanner study.

METHODS

Geometric distortion variation was quantified over a one-year period at 10 sites using the distortion fields estimated from monthly 3D T1-weighted MRI geometrical phantom scans. The variability of volume and distance measurements were quantified using synthetic volumes and a standard quantitative MRI (qMRI) phantom. The effects of geometric distortion corrections on MRI derived volumetric measurements of the human brain were assessed in two subjects scanned on each of the 10 MRI scanners and in 150 subjects with cerebrovascaular disease (CVD) acquired across imaging sites.

RESULTS

Geometric distortions were found to vary substantially between different MRI scanners but were relatively stable on each scanner over a one-year interval. Geometric distortions varied spatially, increasing in severity with distance from the magnet isocenter. In measurements made with the qMRI phantom, the geometric distortion correction decreased the standard deviation of volumetric assessments by 35% and distance measurements by 42%. The average coefficient of variance decreased by 16% in gray matter and white matter volume estimates in the two subjects scanned on the 10 MRI scanners.

CONCLUSION

Geometric distortion correction using an up-to-date correction field is recommended to increase precision in volumetric measurements made from MRI images.

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

Purpose

To develop and evaluate an improved strategy for compensating concomitant field effects in non‐Cartesian MRI at the time of image reconstruction.

Theory

We present a higher‐order reconstruction method, denoted as MaxGIRF, for non‐Cartesian imaging that simultaneously corrects off‐resonance, concomitant fields, and trajectory errors without requiring specialized hardware. Gradient impulse response functions are used to predict actual gradient waveforms, which are in turn used to estimate the spatiotemporally varying concomitant fields based on analytic expressions. The result, in combination with a reference field map, is an encoding matrix that incorporates a correction for all three effects.

Methods

The MaxGIRF reconstruction is applied to noiseless phantom simulations, spiral gradient‐echo imaging of an International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom, and axial and sagittal multislice spiral spin‐echo imaging of a healthy volunteer at 0.55 T. The MaxGIRF reconstruction was compared against previously established concomitant field‐compensation and image‐correction methods. Reconstructed images are evaluated qualitatively and quantitatively using normalized RMS error. Finally, a low‐rank approximation of MaxGIRF is used to reduce computational burden. The accuracy of the low‐rank approximation is studied as a function of minimum rank.

Results

The MaxGIRF reconstruction successfully mitigated blurring artifacts both in phantoms and in vivo and was effective in regions where concomitant fields counteract static off‐resonance, superior to the comparator method. A minimum rank of 8 and 30 for axial and sagittal scans, respectively, gave less than 2% error compared with the full‐rank reconstruction.

Conclusions

The MaxGIRF reconstruction simultaneously corrects off‐resonance, trajectory errors, and concomitant field effects. The impact of this method is greatest when imaging with longer readouts and/or at lower field strength.

Click here for author‐reader discussions

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Left–Right Intensity Asymmetries Vary Depending on Scanner Model for FLAIR and T 1 Weighted MRI Images

Background

Localized regions of left–right image intensity asymmetry (LRIA) were incidentally observed on T₂‐weighted (T₂‐w) and T₁‐weighted (T₁‐w) diagnostic magnetic resonance imaging (MRI) images. Suspicion of herpes encephalitis resulted in unnecessary follow‐up imaging. A nonbiological imaging artifact that can lead to diagnostic uncertainty was identified.

Purpose

To investigate whether systematic LRIA exist for a range of scanner models and to determine if LRIA can introduce diagnostic uncertainty.

Study Type

A retrospective study using the Alzheimer's Disease Neuroimaging Initiative (ADNI) data base.

Subjects

One thousand seven hundred fifty‐three (median age: 72, males/females: 878/875) unique participants with longitudinal data were included.

Field Strength

3T.

Sequences

T₁‐w three‐dimensional inversion‐recovery spoiled gradient‐echo (IR‐SPGR) or magnetization‐prepared rapid gradient‐echo (MP‐RAGE) and T₂‐w fluid‐attenuated inversion recovery (FLAIR) long tau fast spin echo inversion recovery (LT‐FSE‐IR). Only General Electric, Philips, and Siemens' product sequences were used.

Assessment

LRIA was calculated as the left–right percent difference with respect to the mean intensity from automated anatomical atlas segmented regions. Three neuroradiologists with 37 (**), 32 (**), and 3 (**) years of experience rated the clinical impact of 30 T₂‐w three‐dimensional FLAIR exams with LRIA to determine the diagnostic uncertainty. Statistical comparisons between retrospective intensity normalized T₁m and original T₁‐w images were made.

Statistical Tests

For each image type, a linear mixed effects model was fit using LRIA scores from all scanners, regions, and participants as the outcome and age and sex as predictors. Statistical significance was defined as having a P‐value <0.05.

Results

LRIA scores were significantly different from zero on most scanners. All clinicians were uncertain or recommended definite diagnostic follow‐up in 62.5% of cases with LRIA >10%. Individuals with acute brain pathology or focal neurologic deficits are not enrolled in ADNI; therefore, focal signal abnormalities were considered false positives.

Data Conclusion

LRIA is system specific, systematic, creates diagnostic uncertainty, and impacts IR‐SPGR, MP‐RAGE, and LT‐FSE‐IR product sequences.

Level of Evidence: 2

Technical Efficacy Stage: 3

Task group 284 report: magnetic resonance imaging simulation in radiotherapy: considerations for clinical implementation, optimization, and quality assurance

The use of dedicated magnetic resonance simulation (MR-SIM) platforms in Radiation Oncology has expanded rapidly, introducing new equipment and functionality with the overall goal of improving the accuracy of radiation treatment planning. However, this emerging technology presents a new set of challenges that need to be addressed for safe and effective MR-SIM implementation. The major objectives of this report are to provide recommendations for commercially available MR simulators, including initial equipment selection, siting, acceptance testing, quality assurance, optimization of dedicated radiation therapy specific MR-SIM workflows, patient-specific considerations, safety, and staffing. Major contributions include guidance on motion and distortion management as well as MRI coil configurations to accommodate patients immobilized in the treatment position. Examples of optimized protocols and checklists for QA programs are provided. While the recommendations provided here are minimum requirements, emerging areas and unmet needs are also highlighted for future development.

ReUINet: A fast GNL distortion correction approach on a 1.0 T MRI-Linac scanner

PURPOSE

The hybrid system combining a magnetic resonance imaging (MRI) scanner with a linear accelerator (Linac) has become increasingly desirable for tumor treatment because of excellent soft tissue contrast and nonionizing radiation. However, image distortions caused by gradient nonlinearity (GNL) can have detrimental impacts on real-time radiotherapy using MRI-Linac systems, where accurate geometric information of tumors is essential.

METHODS

In this work, we proposed a deep convolutional neural network-based method to efficiently recover undistorted images (ReUINet) for real-time image guidance. The ReUINet, based on the encoder-decoder structure, was created to learn the relationship between the undistorted images and distorted images. The ReUINet was pretrained and tested on a publically available brain MR image dataset acquired from 23 volunteers. Then, transfer learning was adopted to implement the pretrained model (i.e., network with optimal weights) on the experimental three-dimensional (3D) grid phantom and in-vivo pelvis image datasets acquired from the 1.0 T Australian MRI-Linac system.

RESULTS

Evaluations on the phantom (768 slices) and pelvis data (88 slices) showed that the ReUINet achieved improvement over 15 times and 45 times on computational efficiency in comparison with standard interpolation and GNL-encoding methods, respectively. Moreover, qualitative and quantitative results demonstrated that the ReUINet provided better correction results than the standard interpolation method, and comparable performance compared to the GNL-encoding approach.

CONCLUSIONS

Validated by simulation and experimental results, the proposed ReUINet showed promise in obtaining accurate MR images for the implementation of real-time MRI-guided radiotherapy.

Empirical field mapping for gradient nonlinearity correction of multi-site diffusion weighted MRI

BACKGROUND

Achieving inter-site / inter-scanner reproducibility of diffusion weighted magnetic resonance imaging (DW-MRI) metrics has been challenging given differences in acquisition protocols, analysis models, and hardware factors.

PURPOSE

Magnetic field gradients impart scanner-dependent spatial variations in the applied diffusion weighting that can be corrected if the gradient nonlinearities are known. However, retrieving manufacturer nonlinearity specifications is not well supported and may introduce errors in interpretation of units or coordinate systems. We propose an empirical approach to mapping the gradient nonlinearities with sequences that are supported across the major scanner vendors.

STUDY TYPE

Prospective observational study.

SUBJECTS

A spherical isotropic diffusion phantom, and a single human control volunteer.

FIELD STRENGTH/SEQUENCE

3 T (two scanners). Stejskal-Tanner spin echo sequence with b-values of 1000, 2000 s/mm² with 12, 32, and 384 diffusion gradient directions per shell.

ASSESSMENT

We compare the proposed correction with the prior approach using manufacturer specifications against typical diffusion pre-processing pipelines (i.e., ignoring spatial gradient nonlinearities). In phantom data, we evaluate metrics against the ground truth. In human and phantom data, we evaluate reproducibility across scans, sessions, and hardware.

STATISTICAL TESTS

Wilcoxon rank-sum test between uncorrected and corrected data.

RESULTS

In phantom data, our correction method reduces variation in mean diffusivity across sessions over uncorrected data (p < 0.05). In human data, we show that this method can also reduce variation in mean diffusivity across scanners (p < 0.05).

CONCLUSION

Our method is relatively simple, fast, and can be applied retroactively. We advocate incorporating voxel-specific b-value and b-vector maps should be incorporated in DW-MRI harmonization preprocessing pipelines to improve quantitative accuracy of measured diffusion parameters.

Development and validation of 3D MP-SSFP to enable MRI in inhomogeneous magnetic fields

PURPOSE

We demonstrate the feasibility of MRI with missing-pulse steady-state free precession (MP-SSFP) in a 4T magnet with artificially degraded homogeneity.

METHODS

T₁ , T₂ , and diffusion contrast of MP-SSFP was simulated with constant and alternate radiofrequency (RF) phase using an extended phase graph. To validate MP-SSFP performance in human brain imaging, MP-SSFP was tested with two types of artificially introduced inhomogeneous magnetic fields: (1) a pure linear gradient field, and (2) a pseudo-linear gradient field introduced by mounting a head-gradient set at 36 cm from the magnet isocenter. Image distortion induced by the nonlinear inhomogeneous field was corrected using B₀ mapping measured with MP-SSFP.

RESULTS

The maximum flip angle in MP-SSFP was limited to ≤10° because of the large range of resonance frequencies in the inhomogeneous magnetic fields tested in this study. Under this flip-angle limitation, MP-SSFP with constant RF phase provided advantages of higher signal-to-noise ratio and insensitivity to B₁ ⁺ field inhomogeneity as compared with an alternate RF phase. In diffusion simulation, the steady-state magnetization in constant RF phase MP-SSFP increased with an increase of static field gradient up to 8 to 21 mT/m depending on simulation parameters. Experimental results at 4T validated these findings. In human brain imaging, MP-SSFP preserved sufficient signal intensities, but images showed severe image distortion from the pseudo-linear inhomogeneous field. However, following distortion correction, good-quality brain images were achieved.

CONCLUSION

MP-SSFP appears to be a feasible MRI technique for brain imaging in an inhomogeneous magnetic field.

Parallel implementation of L + S signal recovery in dynamic MRI

Dynamic MRI is useful to diagnose different diseases, e.g. cardiac ailments, by monitoring the structure and function of the heart and blood flow through the valves. Faster data acquisition is highly desirable in dynamic MRI, but this may lead to aliasing artifacts due to under-sampling. Advanced image reconstruction algorithms are required to obtain aliasing-free MR images from the acquired under-sampled data. One major limitation of using the advanced reconstruction algorithms is their computationally expensive and time-consuming nature, which make them infeasible for clinical use, especially for applications like cardiac MRI. L + S decomposition model is an approach provided in literature which separates the sparse and low-rank information in dynamic MRI. However, L + S decomposition model is a computationally complex process demanding significant computation time. In this paper, a parallel framework is proposed to accelerate the image reconstruction process of L + S decomposition model using GPU. Experiments are performed on cardiac perfusion dataset ([Formula: see text]) and cardiac cine dataset ([Formula: see text]) using NVIDIA's GeForce GTX780 GPU and Core-i7 CPU. The results show that the proposed method provides up to 18 × speed-up including the memory transfer time (i.e. data transfer between the CPU and GPU) and ~ 46 × speed-up without memory transfer for the cardiac perfusion dataset in our experiments. This level of improvement in the reconstruction time will increase the usefulness of L + S reconstruction by making it feasible for clinical applications.

Ultimate MRI

The basic principles of Magnetic Resonance have been understood for over 70 years and a mainstay of medical imaging for over 40. At this point, it's no longer about simply porting these principles to medical imaging. But we are by no means confined to simply polishing either. Significant innovation and even revolution can come to old technologies. The recent revolution in optical microscopy shattered the resolution constraint imposed by a seemingly fundamental physical law (the diffraction limit) and reinvigorated a 500-year-old modality. Progress comes from re-examining old-ways and sidestepping underlying assumptions. This is already underway for MRI; and is fueled by advances in image reconstruction. Reconstruction increasingly employs sophisticated general models often using subtle and hopefully innocuous prior knowledge about the object. This allows a careful re-examination of some basic prerequisites for MRI such as uniform static fields, linear encoding fields, full Nyquist sampling, or even a stationary object. These powerful reconstruction tools are driving changes in acquisition strategy and basic hardware. The scanner of the future will know more about itself and its patient and his/her biology than ever before. This strategy emboldens relaxed hardware constraints and more specialized scanners, hopefully expanding the reach and value offered by MR imaging.

Distortion-free imaging: A double encoding method (DIADEM) combined with multiband imaging for rapid distortion-free high-resolution diffusion imaging on a compact 3T with high-performance gradients

BACKGROUND

Distortion-free, high-resolution diffusion imaging using DIADEM (Distortion-free Imaging: A Double Encoding Method), proposed recently, has great potential for clinical applications. However, it can suffer from prolonged scan times and its reliability for quantitative diffusion imaging has not been evaluated.

PURPOSE

To investigate the clinical feasibility of DIADEM-based high-resolution diffusion imaging on a novel compact 3T (C3T) by evaluating the reliability of quantitative diffusion measurements and utilizing both the high-performance gradients (80 mT/m, 700 T/m/s) and the sequence optimization with the navigator acquisition window reduction and simultaneous multislice (multiband) imaging.

STUDY TYPE

Prospective feasibility study.

PHANTOM/SUBJECTS

Diffusion quality control phantom scans to evaluate the reliability of quantitative diffusion measurements; 36 normal control scans for B₀ -field mapping; six healthy and two patient subject scans with a brain tumor for comparisons of diffusion and anatomical imaging.

FIELD STRENGTH/SEQUENCE

3T; the standard single-shot echo-planar-imaging (EPI), multishot DIADEM diffusion, and anatomical (2D-FSE [fast-spin-echo], 2D-FLAIR [fluid-attenuated-inversion-recovery], and 3D-MPRAGE [magnetization prepared rapid acquisition gradient echo]) imaging.

ASSESSMENT

The scan time reduction, the reliability of quantitative diffusion measurements, and the clinical efficacy for high-resolution diffusion imaging in healthy control and brain tumor volunteers.

STATISTICAL TEST

Bland-Altman analysis.

RESULTS

The scan time for high in-plane (0.86 mm² ) resolution, distortion-free, and whole brain diffusion imaging were reduced from 10 to 5 minutes with the sequence optimizations. All of the mean apparent diffusion coefficient (ADC) values in phantom were within the 95% confidence interval in the Bland-Altman plot. The proposed acquisition with a total off-resonance coverage of 597.2 Hz wider than the expected bandwidth of 500 Hz in human brain could yield a distortion-free image without foldover artifacts. Compared with EPI, therefore, this approach allowed direct image matching with the anatomical images and enabled improved delineation of the tumor boundaries.

DATA CONCLUSION

The proposed high-resolution diffusion imaging approach is clinically feasible on C3T due to a combination of hardware and sequence improvements.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1 J. Magn. Reson. Imaging 2020;51:296-310.

Large field of view distortion assessment in a low-field MR-linac

PURPOSE

MR-guided radiation therapy (RT) offers unparalleled soft tissue contrast for localization and target tracking. However, MRI distortions may be detrimental to high precision RT. This work characterizes the gradient nonlinearity (GNL) and total distortions over the first year of clinical operation of a 0.35T MR-linac.

METHODS

For GNL characterization, an in-house large field of view (FOV) phantom (60 × 42.5 × 55 cm³ , >6000 spherical landmarks) was configured and scanned at four timepoints with forward/reverse read polarities (Gradient Echo sequence, FA/TR/TE = 28°/30 ms/6 ms). GNL was measured in Anterior-Posterior (AP), Left-Right (LR), and Superior-Inferior (SI) frequency-encoding directions based on deviation of the auto-segmented landmark centroids between rigidly registered MR and CT images and assessed based on radial distance from magnet isocenter. Total distortion was assessed using a 30 × 30 cm² grid phantom oriented along the cardinal axes over >1 year of operation.

RESULTS

The scanner's spatial integrity within the first ~10 months was stable (maximum total distortion variation = 10/6/8%, maximum distortion = 1.41/0.99/1.56 mm in Axial/Coronal/Sagittal planes, respectively). GNL distortions measured during this time period <10 cm from isocenter were (-0.74, 0.45), (-0.67, 0.53), and (-0.86, 0.70) mm in AP/LR/SI directions. In the 10-20 cm range, <1.5% of the distortions exceeded 2 mm in the AP and LR axes while <4% of the distortions exceeded 2 mm for SI. After major repairs and magnet re-shim, detectable changes were observed in total and GNL distortions (20% reduction in AP and 36% increase in SI direction in the 20-25 cm range). Across all timepoints and axes, 38-53% of landmarks in the 20-25 cm range were displaced by >1 mm.

CONCLUSIONS

GNL distortions were negligible within a 10 cm radius from isocenter. However, in the periphery, non-negligible distortions of up to ~7 mm were observed, which may necessitate GNL corrections for MR-IGRT for treatment sites distant from magnet isocenter.

Partial fourier shells trajectory for non-cartesian MRI

Non-Cartesian MRI acquisition has demonstrated various advantages in many clinical applications. The shells trajectory is a 3D non-Cartesian MRI acquisition technique that samples the k-space using a series of concentric shells to achieve efficient 3D isotropic acquisition. Partial Fourier acquisition is an acceleration technique that is widely used in Cartesian MRI. It exploits the conjugate symmetry of k-space measurement to reduce the number of k-space samples compared to full-k-space acquisition, without loss of spatial resolution. For a Cartesian MRI acquisition, the direction of partial Fourier acceleration is aligned either with the phase encoded or frequency encoded direction. In those cases, the underlying image matrix can be reconstructed from the undersampled k-space data using a non-iterative, homodyne reconstruction framework. However, designing a non-Cartesian acquisition trajectory that is compatible with non-iterative homodyne reconstruction is not nearly as straightforward as in the Cartesian case. One reason is the non-iterative homodyne reconstruction requires (slightly over) half of the k-space to be fully sampled. Since the direction of partial Fourier acceleration varies throughout the acquisition in the non-Cartesian trajectory, directly applying the same partial Fourier acquisition pattern (as in Cartesian acquisitions) to a non-Cartesian trajectory does not necessarily yield a continuous, physically-achievable trajectory. In this work, we develop an asymmetric shells trajectory with fully-automated trajectory and gradient waveform design to achieve partial Fourier acquisition for the shells trajectory. We then demonstrate a non-iterative image reconstruction framework for the proposed trajectory. Phantom and in vivo brain scans based on spoiled gradient echo (SPGR) shells and magnetization-prepared shells (MP-shells) were performed to test the proposed trajectory design and reconstruction method. Our phantom and in vivo results demonstrate that the proposed partial Fourier shells trajectory maintains the desirable image contrast and high sampling efficiency from the fully sampled shells, while further reducing data acquisition time.

The effect of concomitant fields in fast spin echo acquisition on asymmetric MRI gradient systems

PURPOSE

To investigate the effect of the asymmetric gradient concomitant fields (CF) with zeroth and first-order spatial dependence on fast/turbo spin-echo acquisitions, and to demonstrate the effectiveness of their real-time compensation.

METHODS

After briefly reviewing the CF produced by asymmetric gradients, the effects of the additional zeroth and first-order CFs on these systems are investigated using extended-phase graph simulations. Phantom and in vivo experiments are performed to corroborate the simulation. Experiments are performed before and after the real-time compensations using frequency tracking and gradient pre-emphasis to demonstrate their effectiveness in correcting the additional CFs. The interaction between the CFs and prescan-based correction to compensate for eddy currents is also investigated.

RESULTS

It is demonstrated that, unlike the second-order CFs on conventional gradients, the additional zeroth/first-order CFs on asymmetric gradients cause substantial signal loss and dark banding in fast spin-echo acquisitions within a typical brain-scan field of view. They can confound the prescan correction for eddy currents and degrade image quality. Performing real-time compensation successfully eliminates the artifacts.

CONCLUSIONS

We demonstrate that the zeroth/first-order CFs specific to asymmetric gradients can cause substantial artifacts, including signal loss and dark bands for brain imaging. These effects can be corrected using real-time compensation. Magn Reson Med 79:1354-1364, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

B0 concomitant field compensation for MRI systems employing asymmetric transverse gradient coils

PURPOSE

Imaging gradients result in the generation of concomitant fields, or Maxwell fields, which are of increasing importance at higher gradient amplitudes. These time-varying fields cause additional phase accumulation, which must be compensated for to avoid image artifacts. In the case of gradient systems employing symmetric design, the concomitant fields are well described with second-order spatial variation. Gradient systems employing asymmetric design additionally generate concomitant fields with global (zeroth-order or B₀ ) and linear (first-order) spatial dependence.

METHODS

This work demonstrates a general solution to eliminate the zeroth-order concomitant field by applying the correct B₀ frequency shift in real time to counteract the concomitant fields. Results are demonstrated for phase contrast, spiral, echo-planar imaging (EPI), and fast spin-echo imaging.

RESULTS

A global phase offset is reduced in the phase-contrast exam, and blurring is virtually eliminated in spiral images. The bulk image shift in the phase-encode direction is compensated for in EPI, whereas signal loss, ghosting, and blurring are corrected in the fast-spin echo images.

CONCLUSION

A user-transparent method to compensate the zeroth-order concomitant field term by center frequency shifting is proposed and implemented. This solution allows all the existing pulse sequences-both product and research-to be retained without any modifications. Magn Reson Med 79:1538-1544, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Gradient nonlinearity calibration and correction for a compact, asymmetric magnetic resonance imaging gradient system

Due to engineering limitations, the spatial encoding gradient fields in conventional magnetic resonance imaging cannot be perfectly linear and always contain higher-order, nonlinear components. If ignored during image reconstruction, gradient nonlinearity (GNL) manifests as image geometric distortion. Given an estimate of the GNL field, this distortion can be corrected to a degree proportional to the accuracy of the field estimate. The GNL of a gradient system is typically characterized using a spherical harmonic polynomial model with model coefficients obtained from electromagnetic simulation. Conventional whole-body gradient systems are symmetric in design; typically, only odd-order terms up to the 5th-order are required for GNL modeling. Recently, a high-performance, asymmetric gradient system was developed, which exhibits more complex GNL that requires higher-order terms including both odd- and even-orders for accurate modeling. This work characterizes the GNL of this system using an iterative calibration method and a fiducial phantom used in ADNI (Alzheimer's Disease Neuroimaging Initiative). The phantom was scanned at different locations inside the 26 cm diameter-spherical-volume of this gradient, and the positions of fiducials in the phantom were estimated. An iterative calibration procedure was utilized to identify the model coefficients that minimize the mean-squared-error between the true fiducial positions and the positions estimated from images corrected using these coefficients. To examine the effect of higher-order and even-order terms, this calibration was performed using spherical harmonic polynomial of different orders up to the 10th-order including even- and odd-order terms, or odd-order only. The results showed that the model coefficients of this gradient can be successfully estimated. The residual root-mean-squared-error after correction using up to the 10th-order coefficients was reduced to 0.36 mm, yielding spatial accuracy comparable to conventional whole-body gradients. The even-order terms were necessary for accurate GNL modeling. In addition, the calibrated coefficients improved image geometric accuracy compared with the simulation-based coefficients.

Image-based gradient non-linearity characterization to determine higher-order spherical harmonic coefficients for improved spatial position accuracy in magnetic resonance imaging

PURPOSE

Spatial position accuracy in magnetic resonance imaging (MRI) is an important concern for a variety of applications, including radiation therapy planning, surgical planning, and longitudinal studies of morphologic changes to study neurodegenerative diseases. Spatial accuracy is strongly influenced by gradient linearity. This work presents a method for characterizing the gradient non-linearity fields on a per-system basis, and using this information to provide improved and higher-order (9th vs. 5th) spherical harmonic coefficients for better spatial accuracy in MRI.

METHODS

A large fiducial phantom containing 5229 water-filled spheres in a grid pattern is scanned with the MR system, and the positions all the fiducials are measured and compared to the corresponding ground truth fiducial positions as reported from a computed tomography (CT) scan of the object. Systematic errors from off-resonance (i.e., B0) effects are minimized with the use of increased receiver bandwidth (±125kHz) and two acquisitions with reversed readout gradient polarity. The spherical harmonic coefficients are estimated using an iterative process, and can be subsequently used to correct for gradient non-linearity. Test-retest stability was assessed with five repeated measurements on a single scanner, and cross-scanner variation on four different, identically-configured 3T wide-bore systems.

RESULTS

A decrease in the root-mean-square error (RMSE) over a 50cm diameter spherical volume from 1.80mm to 0.77mm is reported here in the case of replacing the vendor's standard 5th order spherical harmonic coefficients with custom fitted 9th order coefficients, and from 1.5mm to 1mm by extending custom fitted 5th order correction to the 9th order. Minimum RMSE varied between scanners, but was stable with repeated measurements in the same scanner.

CONCLUSIONS

The results suggest that the proposed methods may be used on a per-system basis to more accurately calibrate MR gradient non-linearity coefficients when compared to vendor standard corrections.

Technical Note: Compact three-tesla magnetic resonance imager with high-performance gradients passes ACR image quality and acoustic noise tests

PURPOSE

A compact, three-tesla magnetic resonance imaging (MRI) system has been developed. It features a 37 cm patient aperture, allowing the use of commercial receiver coils. Its design allows simultaneously for gradient amplitudes of 85 millitesla per meter (mT/m) sustained and 700 tesla per meter per second (T/m/s) slew rates. The size of the gradient system allows for these simultaneous performance targets to be achieved with little or no peripheral nerve stimulation, but also raises a concern about the geometric distortion as much of the imaging will be done near the system's maximum 26 cm field-of-view. Additionally, the fast switching capability raises acoustic noise concerns. This work evaluates the system for both the American College of Radiology's (ACR) MRI image quality protocol and the Food and Drug Administration's (FDA) nonsignificant risk (NSR) acoustic noise limits for MR. Passing these two tests is critical for clinical acceptance.

METHODS

In this work, the gradient system was operated at the maximum amplitude and slew rate of 80 mT/m and 500 T/m/s, respectively. The geometric distortion correction was accomplished by iteratively determining up to the tenth order spherical harmonic coefficients using a fiducial phantom and position-tracking software, with seventh order correction utilized in the ACR test. Acoustic noise was measured with several standard clinical pulse sequences.

RESULTS

The system passes all the ACR image quality tests. The acoustic noise as measured when the gradient coil was inserted into a whole-body MRI system conforms to the FDA NSR limits.

CONCLUSIONS

The compact system simultaneously allows for high gradient amplitude and high slew rate. Geometric distortion concerns have been mitigated by extending the spherical harmonic correction to higher orders. Acoustic noise is within the FDA limits.

NonCartesian MR image reconstruction with integrated gradient nonlinearity correction

PURPOSE

To derive a noniterative gridding-type reconstruction framework for nonCartesian magnetic resonance imaging (MRI) that prospectively accounts for gradient nonlinearity (GNL)-induced image geometrical distortion during MR image reconstruction, as opposed to the standard, image-domain based GNL correction that is applied after reconstruction; to demonstrate that such framework is able to reduce the image blurring introduced by the conventional GNL correction, while still offering effective correction of GNL-induced geometrical distortion and compatibility with off-resonance correction.

METHODS

After introducing the nonCartesian MRI signal model that explicitly accounts for the effects of GNL and off-resonance, a noniterative gridding-type reconstruction framework with integrated GNL correction based on the type-III nonuniform fast Fourier transform (NUFFT) is derived. A novel type-III NUFFT implementation is then proposed as a numerically efficient solution to the proposed framework. The incorporation of simultaneous B0 off-resonance correction to the proposed framework is then discussed. Several phantom and in vivo data acquired via various 2D and 3D nonCartesian acquisitions, including 2D Archimedean spiral, 3D shells with integrated radial and spiral, and 3D radial sampling, are used to compare the results of the proposed and the standard GNL correction methods.

RESULTS

Various phantom and in vivo data demonstrate that both the proposed and the standard GNL correction methods are able to correct the coarse-scale geometric distortion and blurring induced by GNL and off-resonance. However, the standard GNL correction method also introduces blurring effects to corrected images, causing blurring of resolution inserts in the phantom images and loss of small vessel clarity in the angiography examples. On the other hand, the results after the proposed GNL correction show better depiction of resolution inserts and higher clarity of small vessel.

CONCLUSIONS

The proposed GNL-integrated nonCartesian reconstruction method can mitigate the resolution loss that occurs during standard image-domain GNL correction, while still providing effective correction of coarse-scale geometric distortion and blurring induced by GNL and off-resonance.

Gradient pre-emphasis to counteract first-order concomitant fields on asymmetric MRI gradient systems

PURPOSE

To develop a gradient pre-emphasis scheme that prospectively counteracts the effects of the first-order concomitant fields for any arbitrary gradient waveform played on asymmetric gradient systems, and to demonstrate the effectiveness of this approach using a real-time implementation on a compact gradient system.

METHODS

After reviewing the first-order concomitant fields that are present on asymmetric gradients, we developed a generalized gradient pre-emphasis model assuming arbitrary gradient waveforms to counteract their effects. A numerically straightforward, easily implemented approximate solution to this pre-emphasis problem was derived that was compatible with the current hardware infrastructure of conventional MRI scanners for eddy current compensation. The proposed method was implemented on the gradient driver subsystem, and its real-time use was tested using a series of phantom and in vivo data acquired from two-dimensional Cartesian phase-difference, echo-planar imaging, and spiral acquisitions.

RESULTS

The phantom and in vivo results demonstrated that unless accounted for, first-order concomitant fields introduce considerable phase estimation error into the measured data and result in images with spatially dependent blurring/distortion. The resulting artifacts were effectively prevented using the proposed gradient pre-emphasis.

CONCLUSION

We have developed an efficient and effective gradient pre-emphasis framework to counteract the effects of first-order concomitant fields of asymmetric gradient systems. Magn Reson Med 77:2250-2262, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion

Because magnetic resonance imaging-guided radiation therapy (MRIgRT) offers exquisite soft tissue contrast and the ability to image tissues in arbitrary planes, the interest in this technology has increased dramatically in recent years. However, intrinsic geometric distortion stemming from both the system hardware and the magnetic properties of the patient affects MR images and compromises the spatial integrity of MRI-based radiation treatment planning, given that for real-time MRIgRT, precision within 2 mm is desired. In this article, we discuss the causes of geometric distortion, describe some well-known distortion correction algorithms, and review geometric distortion measurements from 12 studies, while taking into account relevant imaging parameters. Eleven of the studies reported phantom measurements quantifying system-dependent geometric distortion, while 2 studies reported simulation data quantifying magnetic susceptibility-induced geometric distortion. Of the 11 studies investigating system-dependent geometric distortion, 5 reported maximum measurements less than 2 mm. The simulation studies demonstrated that magnetic susceptibility-induced distortion is typically smaller than system-dependent distortion but still nonnegligible, with maximum distortion ranging from 2.1 to 2.6 mm at a field strength of 1.5 T. As expected, anatomic landmarks containing interfaces between air and soft tissue had the largest distortions. The evidence indicates that geometric distortion reduces the spatial integrity of MRI-based radiation treatment planning and likely diminishes the efficacy of MRIgRT. Better phantom measurement techniques and more effective distortion correction algorithms are needed to achieve the desired spatial precision.

High slew-rate head-only gradient for improving distortion in echo planar imaging: Preliminary experience

PURPOSE

To investigate the effects on echo planar imaging (EPI) distortion of using high gradient slew rates (SR) of up to 700 T/m/s for in vivo human brain imaging, with a dedicated, head-only gradient coil.

MATERIALS AND METHODS

Simulation studies were first performed to determine the expected echo spacing and distortion reduction in EPI. A head gradient of 42-cm inner diameter and with asymmetric transverse coils was then installed in a whole-body, conventional 3T magnetic resonance imaging (MRI) system. Human subject imaging was performed on five subjects to determine the effects of EPI on echo spacing and signal dropout at various gradient slew rates. The feasibility of whole-brain imaging at 1.5 mm-isotropic spatial resolution was demonstrated with gradient-echo and spin-echo diffusion-weighted EPI.

RESULTS

As compared to a whole-body gradient coil, the EPI echo spacing in the head-only gradient coil was reduced by 48%. Simulation and in vivo results, respectively, showed up to 25-26% and 19% improvement in signal dropout. Whole-brain imaging with EPI at 1.5 mm spatial resolution provided good whole-brain coverage, spatial linearity, and low spatial distortion effects.

CONCLUSION

Our results of human brain imaging with EPI using the compact head gradient coil at slew rates higher than in conventional whole-body MR systems demonstrate substantially improved image distortion, and point to a potential for benefits to non-EPI pulse sequences. J. Magn. Reson. Imaging 2016;44:653-664.

Correction of B 0-induced geometric distortion variations in prospective motion correction for 7T MRI

OBJECTIVE

Prospective motion correction can effectively fix the imaging volume of interest. For large motion, this can lead to relative motion of coil sensitivities, distortions associated with imaging gradients and B 0 field variations. This work accounts for the B 0 field change due to subject movement, and proposes a method for correcting tissue magnetic susceptibility-related distortion in prospective motion correction.

MATERIALS AND METHODS

The B 0 field shifts at the different head orientations were characterized. A volunteer performed large motion with prospective motion correction enabled. The acquired data were divided into multiple groups according to the object positions. The correction of B 0-related distortion was applied to each group of data individually via augmented sensitivity encoding with additionally integrated gradient nonlinearity correction.

RESULTS

The relative motion of the gradients, B 0 field and coil sensitivities in prospective motion correction results in residual spatial distortion, blurring, and coil artifacts. These errors can be mitigated by the proposed method. Moreover, iterative conjugate gradient optimization with regularization provided superior results with smaller RMSE in comparison to standard conjugate gradient.

CONCLUSION

The combined correction of B 0-related distortion and gradient nonlinearity leads to a reduction of residual motion artifacts in prospective motion correction data.

A majorize-minimize framework for Rician and non-central chi MR images

The statistics of many MR magnitude images are described by the non-central chi (NCC) family of probability distributions, which includes the Rician distribution as a special case. These distributions have complicated negative log-likelihoods that are nontrivial to optimize. This paper describes a novel majorize-minimize framework for NCC data that allows penalized maximum likelihood estimates to be obtained by solving a series of much simpler regularized least-squares surrogate problems. The proposed framework is general and can be useful in a range of applications. We illustrate the potential advantages of the framework with real and simulated data in two contexts: 1) MR image denoising and 2) diffusion profile estimation in high angular resolution diffusion MRI. The proposed approach is shown to yield improved results compared to methods that model the noise statistics inaccurately and faster computation relative to commonly-used nonlinear optimization techniques.

Partial fourier and parallel MR image reconstruction with integrated gradient nonlinearity correction

PURPOSE

To describe how integrated gradient nonlinearity (GNL) correction can be used within noniterative partial Fourier (homodyne) and parallel (SENSE and GRAPPA) MR image reconstruction strategies, and demonstrate that performing GNL correction during, rather than after, these routines mitigates the image blurring and resolution loss caused by postreconstruction image domain based GNL correction.

METHODS

Starting from partial Fourier and parallel magnetic resonance imaging signal models that explicitly account for GNL, noniterative image reconstruction strategies for each accelerated acquisition technique are derived under the same core mathematical assumptions as their standard counterparts. A series of phantom and in vivo experiments on retrospectively undersampled data were performed to investigate the spatial resolution benefit of integrated GNL correction over conventional postreconstruction correction.

RESULTS

Phantom and in vivo results demonstrate that the integrated GNL correction reduces the image blurring introduced by the conventional GNL correction, while still correcting GNL-induced coarse-scale geometrical distortion. Images generated from undersampled data using the proposed integrated GNL strategies offer superior depiction of fine image detail, for example, phantom resolution inserts and anatomical tissue boundaries.

CONCLUSION

Noniterative partial Fourier and parallel imaging reconstruction methods with integrated GNL correction reduce the resolution loss that occurs during conventional postreconstruction GNL correction while preserving the computational efficiency of standard reconstruction techniques. Magn Reson Med 75:2534-2544, 2016. © 2015 Wiley Periodicals, Inc.