Trajectory correction based on the gradient impulse response function improves high-resolution UTE imaging of the musculoskeletal system

Comparing CT-Like Images Based on Ultra-Short Echo Time and Gradient Echo T1-Weighted MRI Sequences for the Assessment of Vertebral Disorders Using Histology and True CT as the Reference Standard

BACKGROUND

Several magnetic resonance (MR) techniques have been suggested for radiation-free imaging of osseous structures.

PURPOSE

To compare the diagnostic value of ultra-short echo time and gradient echo T1-weighted MRI for the assessment of vertebral pathologies using histology and computed tomography (CT) as the reference standard.

STUDY TYPE

Prospective.

SUBJECTS

Fifty-nine lumbar vertebral bodies harvested from 20 human cadavers (donor age 73 ± 13 years; 9 male).

FIELD STRENGTH/SEQUENCE

Ultra-short echo time sequence optimized for both bone (UTEb) and cartilage (UTEc) imaging and 3D T1-weighted gradient-echo sequence (T1GRE) at 3 T; susceptibility-weighted imaging (SWI) gradient echo sequence at 1.5 T. CT was performed on a dual-layer dual-energy CT scanner using a routine clinical protocol.

ASSESSMENT

Histopathology and conventional CT were acquired as standard of reference. Semi-quantitative and quantitative morphological features of degenerative changes of the spines were evaluated by four radiologists independently on CT and MR images independently and blinded to all other information. Features assessed were osteophytes, endplate sclerosis, visualization of cartilaginous endplate, facet joint degeneration, presence of Schmorl's nodes, and vertebral dimensions. Vertebral disorders were assessed by a pathologist on histology.

STATISTICAL TESTS

Agreement between T1GRE, SWI, UTEc, and UTEb sequences and CT imaging and histology as standard of reference were assessed using Fleiss' κ and intra-class correlation coefficients, respectively.

RESULTS

For the morphological assessment of osteophytes and endplate sclerosis, the overall agreement between SWI, T1GRE, UTEb, and UTEc with the reference standard (histology combined with CT) was moderate to almost perfect for all readers (osteophytes: SWI, κ range: 0.68-0.76; T1GRE: 0.92-1.00; UTEb: 0.92-1.00; UTEc: 0.77-0.85; sclerosis: SWI, κ range: 0.60-0.70; T1GRE: 0.77-0.82; UTEb: 0.81-0.92; UTEc: 0.61-0.71). For the visualization of the cartilaginous endplate, UTEc showed the overall best agreement with the reference standard (histology) for all readers (κ range: 0.85-0.93).

DATA CONCLUSIONS

Morphological assessment of vertebral pathologies was feasible and accurate using the MR-based bone imaging sequences compared to CT and histopathology. T1GRE showed the overall best performance for osseous changes and UTEc for the visualization of the cartilaginous endplate.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 2.

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.

3D balanced SSFP UTE MRI for multiple contrasts whole brain imaging

PURPOSE

This study aimed to develop a new high-resolution MRI sequence for the imaging of the ultra-short transverse relaxation time (uT2) components in the brain, while simultaneously providing proton density (PD) contrast for reference and quantification.

THEORY

The sequence combines low flip angle balanced SSFP (bSSFP) and UTE techniques, together with a 3D dual-echo rosette k-space trajectory for readout.

METHODS

The expected image contrast was evaluated by simulations. A study cohort of six healthy volunteers and eight multiple sclerosis (MS) patients was recruited to test the proposed sequence. Subtraction between two TEs was performed to extract uT2 signals. In addition, conventional longitudinal relaxation time (T1) weighted, T2-weighted, and PD-weighted MRI sequences were also acquired for comparison.

RESULTS

Typical PD-contrast was found in the second TE images, while uT2 signals were selectively captured in the first TE images. The subtraction images presented signals primarily originating from uT2 components, but only if the first TE is short enough. Lesions in the MS subjects showed hyperintense signals in the second TE images but were hypointense signals in the subtraction images. The lesions had significantly lower signal intensity in subtraction images than normal white matter (WM), which indicated a reduction of uT2 components likely associated with myelin.

CONCLUSION

3D isotropic sub-millimeter (0.94 mm) spatial resolution images were acquired with the novel bSSFP UTE sequence within 3 min. It provided easy extraction of uT2 signals and PD-contrast for reference within a single acquisition.

Fast measurement of the gradient system transfer function at 7 T

PURPOSE

In this work, a new method to determine the gradient system transfer function (GSTF) with high frequency resolution and high SNR is presented, using fast and simple phantom measurements. The GSTF is an effective instrument for hardware characterization and calibration, which can be used to correct for gradient distortions, or enhance gradient fidelity.

METHODS

The thin-slice approach for phantom-based measurements of the GSTF is expanded by adding excitations that are shifted after the application of the probing gradient, to capture long-lasting field fluctuations with high SNR. A physics-informed regularization procedure is implemented to derive high-quality transfer functions from a small number of measurements. The resulting GSTFs are evaluated by means of gradient time-course estimation and pre-emphasis of a trapezoidal test gradient on a 7T scanner.

RESULTS

The GSTFs determined with the proposed method capture sharp mechanical resonances with a high level of detail. The measured trapezoidal gradient progressions are authentically reproduced by the GSTF estimations on all three axes. The GSTF-based pre-emphasis considerably improves the gradient fidelity in the plateau phase of the test gradient and almost completely eliminates lingering field oscillations.

CONCLUSION

The presented approach allows fast and simple characterization of gradient field fluctuations caused by long-living eddy current and vibration effects, which become more pronounced at ultrahigh field strengths.

Diagnostic value of water-fat-separated images and CT-like susceptibility-weighted images extracted from a single ultrashort echo time sequence for the evaluation of vertebral fractures and degenerative changes of the spine

OBJECTIVES

To evaluate the performance of single-echo Dixon water-fat imaging and computed tomography (CT)-like imaging based on a single ultrashort echo time (sUTE) MR sequence for imaging of vertebral fractures as well as degenerative bone changes of the spine in comparison to conventional CT and MR sequences.

METHODS

Thirty patients with suspected acute vertebral fractures were examined using a 3-T MRI, including an sUTE sequence as well as short-tau inversion recovery (STIR) and T1-weighted sequences. During postprocessing, water-fat separation was performed by solving the smoothness-constrained inverse water-fat problem based on a single-complex UTE image. By removing the unwanted low-frequency phase terms, additional MR-based susceptibility-weighted-like (SW-like) images with CT-like contrast were created. Two radiologists evaluated semi-quantitative and quantitative features of fractures and degenerative changes independently and separately on CT and MR images.

RESULTS

In total, all 58 fractures were accurately detected of whom 24 were correctly classified as acute fractures with an edema detected on the water-fat-separated UTE images, using STIR and T1w sequences as standard of reference. For the morphological assessment of fractures and degenerative changes, the overall agreement between SW-like images and CT was substantial to excellent (e.g., Genant: κ 0.90 (95% confidence interval 0.54-1.00); AO/Magerl: κ 0.75 (95% confidence interval 0.43-1.00)). Overall inter-reader agreement for water-fat-separated UTE images and SW-like images was substantial to almost perfect.

CONCLUSION

Detection and assessment of vertebral fractures and degenerative bone changes of the spine were feasible and accurate using water-fat-separated images as well as SW-like images, both derived from the same sUTE-Dixon sequence.

KEY POINTS

• The detection of acute vertebral fractures was feasible using water-fat-separated images and CT-like images reconstructed from one sUTE sequence. • Assessment of the vertebral fractures using SW-like images with CT-like contrast was found to be comparable to conventional CT. • sUTE imaging of the spine can help reduce examination times and radiation exposure.

On quantification errors of R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction mapping in trabecularized bone marrow in the static dephasing regime

PURPOSE

To study the effect of field inhomogeneity distributions in trabecularized bone regions on the gradient echo (GRE) signal with short TEs and to characterize quantification errors on R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction (PDFF) maps when using a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model at short TEs.

METHODS

Field distortions were simulated based on a trabecular bone micro CT dataset. Simulations were performed for different bone volume fractions (BV/TV) and for different bone-fat composition values. A multi-TE UTE acquisition was developed to acquire multiple UTEs with random order to minimize eddy currents. The acquisition was validated in phantoms and applied in vivo in a volunteer's ankle and knee. Chemical shift encoded MRI (CSE-MRI) based on a Cartesian multi-TE GRE scan was acquired in the spine of patients with metastatic bone disease.

RESULTS

Simulations showed that signal deviations from the exponential signal decay at short TEs were more prominent for a higher BV/TV. UTE multi-TE measurements reproduced in vivo the simulation-based predicted behavior. In regions with high BV/TV, the presence of field inhomogeneities induced an R 2 * $$ {R}_2^{\ast } $$ underestimation in trabecularized bone marrow when using CSE-MRI at 3T with a short TE.

CONCLUSION

R 2 * $$ {R}_2^{\ast } $$ can be underestimated when using short TEs (<2 ms at 3 T) and a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model in multi-echo GRE acquisitions of trabecularized bone marrow.

Nonlinear droop compensation for current waveforms in MRI gradient systems

PURPOSE

Providing accurate gradient currents is challenging due to the gradient chain nonlinearities, arising from gradient power amplifiers and power supply stages. This work introduces a new characterization approach that takes the amplifier and power supply into account, resulting in a nonlinear model that compensates for the current droop.

METHODS

The gradient power amplifier and power supply stage were characterized by a modified state-space averaging technique. The resulting nonlinear model was inverted and used in feedforward to control the gradient coil current. A custom-built two-channel z-gradient coil was driven by high-switching (1 MHz), low-cost amplifiers (<$200) using linear and nonlinear controllers. High-resolution (<80 ps) pulse-width-modulation signals were used to drive the amplifiers. MRI experiments were performed to validate the nonlinear controller's effectiveness.

RESULTS

The simulation results validated the functionality of the state-space averaging method in characterizing the gradient system. The performance of linear and nonlinear controllers in generating a trapezoidal current waveform was compared in simulations and experiments. The integral errors between the desired waveform and waveforms generated by linear and nonlinear controllers were 1.9% and 0.13%, respectively, confirming the capability of the nonlinear controller to compensate for the current droop. Phantom images validated the nonlinear controller's ability to correct droop-induced distortions.

CONCLUSION

Benchtop measurements and MRI experiments demonstrated that the proposed nonlinear characterization and digitally implemented feedforward controller could drive gradient coils with droop-free current waveforms (without a feedback loop). In experiments, the nonlinear controller outperformed the linear controller by a 14-fold reduction in the integral error of a test waveform.

3D MR fingerprinting using Seiffert spirals

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Assessment of vertebral fractures and edema of the thoracolumbar spine based on water-fat and susceptibility-weighted images derived from a single ultra-short echo time scan

PURPOSE

To develop a methodology to simultaneously perform single echo Dixon water-fat imaging and susceptibility-weighted imaging (SWI) based on a single echo time (TE) ultra-short echo time (UTE) (sUTE) scan to assess vertebral fractures and degenerative bone changes in the thoracolumbar spine.

METHODS

A methodology was developed to solve the smoothness-constrained inverse water-fat problem to separate water and fat while removing unwanted low-frequency phase terms. Additionally, the corrected UTE phase was used for SWI. UTE imaging (TE: 0.14 ms, 3T MRI) was performed in the lumbar spine of nine patients with vertebral fractures and bone marrow edema (BME). All images were reviewed by two radiologists. Water- and fat-separated images were analyzed in comparison with short-tau inversion recovery (STIR) and with respect to BME visibility. The visibility of fracture lines and cortical outlining of the UTE magnitude images were analyzed in comparison with computed tomography.

RESULTS

Unwanted phase components, dominated by the B₁ phase, were removed from the UTE phase images. The rating of the diagnostic quality of BME visualization showed a high preference for the sUTE-Dixon water- and fat-separated images in comparison with STIR. The UTE magnitude images enabled better visualizing fracture lines compared with STIR and slightly better visibility of cortical outlining. With increasing SWI weighting osseous structures and fatty tissues were enhanced.

CONCLUSION

The proposed sUTE-Dixon-SWI methodology allows the removal of unwanted low-frequency phases and enables water-fat separation and SWI processing from a single complex UTE image. The methodology can be used for the simultaneous assessment of vertebral fractures and BME of the thoracolumbar spine.

Non-Cartesian k-space trajectory calculation based on concurrent reading of the gradient amplifiers' output currents

PURPOSE

Non-Cartesian imaging sequences involve sampling during rapid variation of the encoding field gradients. The quality of the reconstructed images often suffers from insufficient knowledge of the exact dynamics of the actual fields applied during sampling.

METHODS

We propose determination of the accurate field dynamics by measuring the currents at the gradient amplifier outputs using the amplifiers' internal sensors concurrently with imaging. The actual dynamic field evolution is then determined by convolution with the measured current-to-field impulse response function of the gradient coil. Integration of the gradient field evolution allows derivation of the k-space trajectory for reconstruction.

RESULTS

The current-based approach is investigated in spiral and ultrashort TE phantom imaging. In comparison with the model-based product reconstruction as well as a correction approach based on the conventional input waveform-to-field impulse response function, it provides slightly improved image quality. The improvement is ascribed to a better representation of eddy current and amplifier nonlinearity effects.

CONCLUSION

Trajectory calculation based on measured amplifier output currents offers a robust, purely measurement-based alternative to conventional model-based approaches. The implementation can mitigate gradient amplifier imperfections with no or little additional hardware effort.