Convolution kernel design and efficient algorithm for sampling density correction

Impact of undersampling on preclinical lung T2* mapping with 3D radial UTE MRI at 7 T

Lung diseases are almost invariably heterogeneous and progressive, making it imperative to capture temporally and spatially explicit information to understand the disease initiation and progression. Imaging the lung with MRI-particularly in the preclinical setting-has historically been challenging because of relatively low lung tissue density, rapid cardiac and respiratory motion, and rapid transverse (T2*) relaxation. These limitations can largely be mitigated using ultrashort-echo-time (UTE) sequences, which are intrinsically robust to motion and avoid significant T2* decay. A significant disadvantage of common radial UTE sequences is that they require inefficient, center-out k-space sampling, resulting in long acquisition times relative to conventional Cartesian sequences. Therefore, pulmonary images acquired with radial UTE are often undersampled to reduce acquisition time. However, undersampling reduces image SNR, introduces image artifacts, and degrades true image resolution. The level of undersampling is further increased if offline gating techniques like retrospective gating are employed, because only a portion (∼40-50%) of the data is used in the final image reconstruction. Here, we explore the impact of undersampling on SNR and T2* mapping in mouse lung imaging using simulation and in-vivo data. Increased scatter in both metrics was noticeable at around 50% sampling. Parenchymal apparent SNR only decreased slightly (average decrease ∼ 1.4) with as little as 10% sampling. Apparent T2* remained similar across undersampling levels, but it became significantly increased (p < 0.05) below 80% sampling. These trends suggest that undersampling can generate quantifiable, but moderate changes in the apparent value of T2*. Moreover, these approaches to assess the impact of undersampling are straightforward to implement and can readily be expanded to assess the quantitative impact of other MR acquisition and reconstruction parameters.

Toward the use of MRI measurements of bound and pore water in fracture risk assessment

The current clinical assessment of fracture risk lacks information about the inherent quality of a person's bone tissue. Working toward an imaging-based approach to quantify both a bone tissue quality marker (tissue hydration as water bound to the matrix) and a bone microstructure marker (porosity as water in pores), we hypothesized that the concentrations of bound water (Cbw) are lower and concentrations of pore water (Cpw) are higher in patients with osteoporosis (OP) than in age- and sex-matched adults without the disease. Using recent developments in ultrashort echo time (UTE) magnetic resonance imaging (MRI), maps of Cbw and Cpw were acquired from the uninjured distal third radius (Study 1) of 20 patients who experienced a fragility fracture of the distal radius (Fx) and 20 healthy controls (Non-Fx) and from the tibia mid-diaphysis (Study 2) of 30 women with clinical OP (low T-scores) and 15 women without OP (normal T-scores). In Study 1, Cbw was significantly lower (p = 0.0018) and Cpw was higher (p = 0.0022) in the Fx than in the Non-Fx group. In forward stepwise, logistic regression models using Bayesian Information Criterion for selecting the best set of predictors (from imaging parameters, age, BMI, and DXA scanner type), the area-under-the-receiver operator characteristics-curve (AUC with 95 % confidence intervals) was 0.73 (0.56, 0.86) for hip aBMD (best predictors without MRI) and 0.86 (0.70, 0.95) for the combination of Cbw and Cpw (best predictors overall). In Study 2, Cbw was significantly lower (p = 0.0005) in women with OP (23.8 ± 4.3 1H mol/L) than in women without OP (29.9 ± 6.4 1H mol/L); Cpw was significantly higher by estimate of 2.9 1H mol/L (p = 0.0298) with clinical OP, but only when accounting for the type of UTE-MRI scan with 3D providing higher values than 2D (p < 0.0001). Lastly, Cbw, but not Cpw, was sensitive to bone forming osteoporosis medications over 12-months. UTE-MRI-derived measurements of bound and pore water concentrations are potential, aBMD-independent predictors of fracture risk.

Optimization in the space domain for density compensation with the nonuniform FFT

The non-uniform Discrete Fourier Transform algorithm has shown great utility for reconstructing images from non-uniformly spaced Fourier samples in several imaging modalities. Due to the non-uniform spacing, some correction for the variable density of the samples must be made. Common methods for generating density compensation values are either sub-optimal or only consider a finite set of points in the optimization. This manuscript presents an algorithm for generating density compensation values from a set of Fourier samples that takes into account the point spread function over an entire rectangular region in the image domain. We show that the reconstructed images using the density compensation values of this method are of superior quality when compared to other standard methods. Results are shown with a numerical phantom and with magnetic resonance images of the abdomen and the knee.

Finite element analysis of bone mechanical properties using MRI-derived bound and pore water concentration maps

Ultrashort echo time (UTE) MRI techniques can be used to image the concentration of water in bones. Particularly, quantitative MRI imaging of collagen-bound water concentration (Cbw) and pore water concentration (Cpw) in cortical bone have been shown as potential biomarkers for bone fracture risk. To investigate the effect of Cbw and Cpw on the evaluation of bone mechanical properties, MRI-based finite element models of cadaver radii were generated with tissue material properties derived from 3 D maps of Cbw and Cpw measurements. Three-point bending tests were simulated by means of the finite element method to predict bending properties of the bone and the results were compared with those from direct mechanical testing. The study results demonstrate that these MRI-derived measures of Cbw and Cpw improve the prediction of bone mechanical properties in cadaver radii and have the potential to be useful in assessing patient-specific bone fragility risk.

View-sharing for 4D magnetic resonance imaging with randomized projection-encoding enables improvements of respiratory motion imaging for treatment planning in abdominothoracic radiotherapy

Background and Purpose

The accuracy and precision of radiation therapy are dependent on the characterization of organ-at-risk and target motion. This work aims to demonstrate a 4D magnetic resonance imaging (MRI) method for improving spatial and temporal resolution in respiratory motion imaging for treatment planning in abdominothoracic radiotherapy.

Materials and Methods

The spatial and temporal resolution of phase-resolved respiratory imaging is improved by considering a novel sampling function based on quasi-random projection-encoding and peripheral k-space view-sharing. The respiratory signal is determined directly from k-space, obviating the need for an external surrogate marker. The average breathing curve is used to optimize spatial resolution and temporal blurring by limiting the extent of data sharing in the Fourier domain. Improvements in image quality are characterized by evaluating changes in signal-to-noise ratio (SNR), resolution, target detection, and level of artifact. The method is validated in simulations, in a dynamic phantom, and in-vivo imaging.

Results

Sharing of high-frequency k-space data, driven by the average breathing curve, improves spatial resolution and reduces artifacts. Although equal sharing of k-space data improves resolution and SNR in stationary features, phases with large temporal changes accumulate significant artifacts due to averaging of high frequency features. In the absence of view-sharing, no averaging and detection artifacts are observed while spatial resolution is degraded.

Conclusions

The use of a quasi-random sampling function, with view-sharing driven by the average breathing curve, provides a feasible method for self-navigated 4D-MRI at improved spatial resolution.

Preclinical MRI to Quantify Pulmonary Disease Severity and Trajectories in Poorly Characterized Mouse Models: A Pedagogical Example Using Data from Novel Transgenic Models of Lung Fibrosis

Structural remodeling in lung disease is progressive and heterogeneous, making temporally and spatially explicit information necessary to understand disease initiation and progression. While mouse models are essential to elucidate mechanistic pathways underlying disease, the experimental tools commonly available to quantify lung disease burden are typically invasive (e.g., histology). This necessitates large cross-sectional studies with terminal endpoints, which increases experimental complexity and expense. Alternatively, magnetic resonance imaging (MRI) provides information noninvasively, thus permitting robust, repeated-measures statistics. Although lung MRI is challenging due to low tissue density and rapid apparent transverse relaxation (T₂* <1 ms), various imaging methods have been proposed to quantify disease burden. However, there are no widely accepted strategies for preclinical lung MRI. As such, it can be difficult for researchers who lack lung imaging expertise to design experimental protocols-particularly for novel mouse models. Here, we build upon prior work from several research groups to describe a widely applicable acquisition and analysis pipeline that can be implemented without prior preclinical pulmonary MRI experience. Our approach utilizes 3D radial ultrashort echo time (UTE) MRI with retrospective gating and lung segmentation is facilitated with a deep-learning algorithm. This pipeline was deployed to assess disease dynamics over 255 days in novel, transgenic mouse models of lung fibrosis based on disease-associated, loss-of-function mutations in Surfactant Protein-C. Previously identified imaging biomarkers (tidal volume, signal coefficient of variation, etc.) were calculated semi-automatically from these data, with an objectively-defined high signal volume identified as the most robust metric. Beyond quantifying disease dynamics, we discuss common pitfalls encountered in preclinical lung MRI and present systematic approaches to identify and mitigate these challenges. While the experimental results and specific pedagogical examples are confined to lung fibrosis, the tools and approaches presented should be broadly useful to quantify structural lung disease in a wide range of mouse models.

Rapid single scan ramped hybrid-encoding for bicomponent T2* mapping in a human knee joint: A feasibility study

The purpose of this study is to determine the feasibility of using a single scan ramped hybrid-encoding (RHE) method for rapid bicomponent T2* analysis of the human knee joint. The proposed method utilizes RHE to acquire ultrashort echo time (UTE) and subsequent gradient echo images at 16 different echo times ranging between 40 μs and 30 ms in a single scan. In the proposed RHE technique, UTE imaging was followed by acquisition of 14 gradient recalled echo images, where an additional UTE image was obtained within the first readout by oversampling single point imaging (SPI) encoding. The single scan RHE method with a 9-minute scan time was performed on human cadaveric knee joints from six donors and in vivo knee joints from four healthy volunteers at 3 T. A bicomponent signal model was used to characterize the short T2* and long T2* water components. Mean bicomponent T2* parameters for patellar tendon, anterior cruciate ligament (ACL), posterior cruciate ligament (PCL) and meniscus were calculated. In the experimental results, the RHE technique provided bicomponent T2* parameter estimations of tendon, ACL, PCL and meniscus, which were similar to previously reported values in the literature. In conclusion, the proposed single scan RHE technique provides rapid bicomponent T2* analysis of the human knee joint with a total scan time of less than 9 minutes.

Short-T2 MRI: Principles and recent advances

Among current modalities of biomedical and diagnostic imaging, MRI stands out by virtue of its versatile contrast obtained without ionizing radiation. However, in various cases, e.g., water protons in tissues such as bone, tendon, and lung, MRI performance is limited by the rapid decay of resonance signals associated with short transverse relaxation times T₂ or T₂*. Efforts to address this shortcoming have led to a variety of specialized short-T₂ techniques. Recent progress in this field expands the choice of methods and prompts fresh considerations with regard to instrumentation, data acquisition, and signal processing. In this review, the current status of short-T₂ MRI is surveyed. In an attempt to structure the growing range of techniques, the presentation highlights overarching concepts and basic methodological options. The most frequently used approaches are described in detail, including acquisition strategies, image reconstruction, hardware requirements, means of introducing contrast, sources of artifacts, limitations, and applications.

Diffusion MRI using two-dimensional single-shot radial imaging (2D ss-rDWI) with variable flip angle and random view ordering

PURPOSE

The main objective of this study is to develop a 2D single-shot radial-DWI (2D ss-rDWI) technique to reduce motion artifacts and geometric distortion in DW images.

METHOD

A diffusion-preparation module is developed and applied prior to the data acquisition. Because the diffusion-prepared longitudinal magnetization is measured over multiple RF excitations in each shot, 2D ss-rDWI is subject to low signal-to-noise ratio (SNR). We used variable-flip angle (VFA), random view ordering (RVO), and sliding spokes, and compared the performances to constant flip angle (CFA), smooth view ordering (SVO), and identical spoke averaging, respectively. For each technique, we performed numerical simulation and MRI experiments on a fluid phantom as well as in-vivo human brain studies with a 3 T MRI system.

RESULTS

Using VFA, optimal SNR was acquired for 2D ss-rDWI. Using SVO, the high signal is clustered at specific quadrant in 2D k-space: the first quadrant using high initial flip angle or the last quadrant using the low flip angle. This clustered signal in k-space led to geometric distortion in image space. 2D ss-rDWI using RVO spreads the high signaled spokes over all angular directions and removes the view-order-related distortion. The in-vivo images using 2D ss-rDWI with VFA and RVO show no geometric distortion at the skull base brain, but greatly reduced SNR compared with those using 2D ss-DWEPI.

CONCLUSION

2D ss-rDWI is optimized by using VFA with RVO. The resultant DWI using 2D ss-rDWI is insensitive to motion-induced artifacts and geometric distortion. Even with low SNR, it may be useful for DWI of organs limited by severe susceptibility-induced geometric distortion.

Longitudinal free-breathing MRI measurement of murine lung physiology in a progressive model of lung fibrosis

To longitudinally monitor progressive fibrosis in the transforming growth factor-α (TGF-α) transgenic mouse model of lung fibrosis, we used retrospective self-gating ultrashort echo time (UTE) magnetic resonance imaging (MRI) to image mouse lung at baseline and after 4 and 8 wk of fibrosis initiation via doxycycline administration. Only bitransgenic mice were used in this study and divided into two cohorts: six mice were fed doxycycline food to induce lung fibrosis (referred to as Dox cohort), and five other mice were fed normal food (referred to as control cohort). Lung mechanics, histology, and hydroxyproline were assessed after the final MRI. A linear mixed-effects model was used to analyze MRI-derived longitudinal lung-function parameters. Tidal volume decreased at a rate of -0.016 ± 0.002 ml/week [χ²(1) = 16.48, P < 0.001] for Dox cohort and increased at a rate of 0.010 ± 0.003 ml/week [χ²(1) = 6.37, P = 0.01] for control cohort. Minute ventilation decreased at a rate of -1.71 ± 0.26 ml·min^-1·wk^-1 [χ²(1) = 14.04, P < 0.001] for Dox cohort but did not change significantly over time for control cohort. High-density lung volume percentage increased at a rate of 3.9 ± 0.7%/wk for Dox cohort [χ²(1) = 11.47, P < 0.001] but did not change significantly over time for control cohort. MRI-derived lung structure and function parameters were strongly correlated with pleural thickness, hydroxyproline content, lung compliance, airway resistance, and airway elastance. We conclude that self-gating UTE MRI could be used to longitudinally monitor lung fibrosis in the TGF-α transgenic mouse model. NEW & NOTEWORTHY Self-gating UTE MRI was used to monitor morphology and physiology in lung fibrosis in a transforming growth factor-α transgenic mouse model. Tidal volume was shown for the first time to correlate strongly with conventional metrics of fibrosis such as hydroxyproline and pleural thickness.

Whole-heart coronary MR angiography using a 3D cones phyllotaxis trajectory

PURPOSE

To develop a 3D cones steady-state free precession sequence with improved robustness to respiratory motion while mitigating eddy current artifacts for free-breathing whole-heart coronary magnetic resonance angiography.

METHOD

The proposed sequence collects cone interleaves using a phyllotaxis pattern, which allows for more distributed k-space sampling for each heartbeat compared to a typical sequential collection pattern. A Fibonacci number of segments is chosen to minimize eddy current effects with the trade-off of an increased number of acquisition heartbeats. For verification, phyllotaxis-cones is compared to sequential-cones through simulations, phantom studies, and in vivo coronary scans with 8 subjects using 2D image-based navigators for retrospective motion correction.

RESULTS

Simulated point spread functions and moving phantom results show less coherent motion artifacts for phyllotaxis-cones compared to sequential-cones. Assessment of the right and left coronary arteries using reader scores and the image edge profile acutance vessel sharpness metric indicate superior image quality and sharpness for phyllotaxis-cones.

CONCLUSION

Phyllotaxis 3D cones results in improved qualitative image scores and coronary vessel sharpness for free-breathing whole-heart coronary magnetic resonance angiography compared to standard sequential ordering when using a steady-state free precession sequence.

Murine pulmonary imaging at 7T: T2* and T1 with anisotropic UTE

PURPOSE

To measure the T2* and T₁ of mouse lung at 7T via anisotropic-resolution radial ultrashort echo-time imaging with ellipsoidal k-space coverage.

METHODS

Ellipsoidal field-of-view was created by expanding uniform spherical k-space coverage. The effects of T2* and ellipsoidal sampling on image resolution were investigated by using point-spread-function analysis and resolution phantoms. Finally, this ellipsoidal sampling approach was used to measure the lung T2* and T₁ of healthy C57BL/6 mice at the increasingly common preclinical field strength of 7T.

RESULTS

Lung parenchyma T2* of 17- to 23-week-old mice at 7T was 0.395 ± 0.033 ms. T₁ of lung and left- and right-heart ventricles was 1452.5 ± 87.0 ms, 1810.5 ± 54.6 ms, and 1602.6 ± 120.9 ms, respectively. Ellipsoidal k-space sampling provides enhanced resolution for a fixed scanning time or provides equivalent (although anisotropic) spatial resolution with reduced scanning times, while simultaneously avoiding fold-in artifacts.

CONCLUSION

Using these techniques, the first T2* and T₁ measures of mouse lung at 7T are reported. Ultrashort echo-time imaging with ellipsoidal k-space coverage significantly increases in-plane resolution without increasing scanning time, or equivalently, decreases scanning time while maintaining equivalent in-plane resolution. Magn Reson Med 79:2254-2264, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A rapid and robust gradient measurement technique using dynamic single-point imaging

PURPOSE

We propose a new gradient measurement technique based on dynamic single-point imaging (SPI), which allows simple, rapid, and robust measurement of k-space trajectory.

METHODS

To enable gradient measurement, we utilize the variable field-of-view (FOV) property of dynamic SPI, which is dependent on gradient shape. First, one-dimensional (1D) dynamic SPI data are acquired from a targeted gradient axis, and then relative FOV scaling factors between 1D images or k-spaces at varying encoding times are found. These relative scaling factors are the relative k-space position that can be used for image reconstruction. The gradient measurement technique also can be used to estimate the gradient impulse response function for reproducible gradient estimation as a linear time invariant system.

RESULTS

The proposed measurement technique was used to improve reconstructed image quality in 3D ultrashort echo, 2D spiral, and multi-echo bipolar gradient-echo imaging. In multi-echo bipolar gradient-echo imaging, measurement of the k-space trajectory allowed the use of a ramp-sampled trajectory for improved acquisition speed (approximately 30%) and more accurate quantitative fat and water separation in a phantom.

CONCLUSION

The proposed dynamic SPI-based method allows fast k-space trajectory measurement with a simple implementation and no additional hardware for improved image quality. Magn Reson Med 78:950-962, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Ramped hybrid encoding for improved ultrashort echo time imaging

PURPOSE

We propose a new acquisition to minimize the per-excitation encoding duration and improve the imaging capability for short T2 * species.

METHODS

In the proposed ramped hybrid encoding (RHE) technique, gradients are applied before the radiofrequency (RF) pulse as in pointwise encoding time reduction with radial acquisition (PETRA) and zero echo time (ZTE) imaging. However, in RHE, gradients are rapidly ramped after RF excitation to the maximum amplitude to minimize encoding duration. To acquire central k-space data not measured during RF deadtime, RHE uses a hybrid encoding scheme similar to PETRA. A new gradient calibration method based on single-point imaging was developed to estimate the k-space trajectory and enable robust and high quality reconstruction.

RESULTS

RHE enables a shorter per-excitation encoding time and provides the highest spatial resolution among ultrashort T2 * imaging methods. In phantom and in vivo experiments, RHE exhibited robust imaging with negligible chemical shift or blurriness caused by T2 * decay and unwanted slice selection.

CONCLUSION

RHE allows the shortest per-excitation encoding time for ultrashort T2 * imaging, which alleviates the impact of fast T2 * decay occurring during encoding, and enables improved spatial resolution. Magn Reson Med 76:814-825, 2016. © 2015 Wiley Periodicals, Inc.

High-resolution variable-density 3D cones coronary MRA

PURPOSE

To improve the spatial/temporal resolution of whole-heart coronary MR angiography by developing a variable-density (VD) 3D cones acquisition suitable for image reconstruction with parallel imaging and compressed sensing techniques.

METHODS

A VD 3D cones trajectory design incorporates both radial and spiral trajectory undersampling techniques to achieve higher resolution. This design is used to generate a VD 3D cones trajectory with 0.8 mm/66 ms isotropic spatial/temporal resolution, using a similar number of readouts as our previous fully sampled cones trajectory (1.2 mm/100 ms). Scans of volunteers and patients are performed to evaluate the performance of the VD trajectory, using non-Cartesian L1 -ESPIRiT for high-resolution image reconstruction.

RESULTS

With gridding reconstruction, the high-resolution scans experience an expected drop in signal-to-noise and contrast-to-noise ratios, but with L1 -ESPIRiT, the apparent noise is substantially reduced. Compared with 1.2 mm images, in each volunteer, the L1 -ESPIRiT 0.8 mm images exhibit higher vessel sharpness values in the right and left anterior descending arteries.

CONCLUSION

Coronary MR angiography with isotropic submillimeter spatial resolution and high temporal resolution can be performed with VD 3D cones to improve the depiction of coronary arteries.

Graphical programming interface: A development environment for MRI methods

PURPOSE

To introduce a multiplatform, Python language-based, development environment called graphical programming interface for prototyping MRI techniques.

METHODS

The interface allows developers to interact with their scientific algorithm prototypes visually in an event-driven environment making tasks such as parameterization, algorithm testing, data manipulation, and visualization an integrated part of the work-flow. Algorithm developers extend the built-in functionality through simple code interfaces designed to facilitate rapid implementation.

RESULTS

This article shows several examples of algorithms developed in graphical programming interface including the non-Cartesian MR reconstruction algorithms for PROPELLER and spiral as well as spin simulation and trajectory visualization of a FLORET example.

CONCLUSION

The graphical programming interface framework is shown to be a versatile prototyping environment for developing numeric algorithms used in the latest MR techniques.

An augmented Lagrangian based compressed sensing reconstruction for non-Cartesian magnetic resonance imaging without gridding and regridding at every iteration

Background

Non-Cartesian trajectories are used in a variety of fast imaging applications, due to the incoherent image domain artifacts they create when undersampled. While the gridding technique is commonly utilized for reconstruction, the incoherent artifacts may be further removed using compressed sensing (CS). CS reconstruction is typically done using conjugate-gradient (CG) type algorithms, which require gridding and regridding to be performed at every iteration. This leads to a large computational overhead that hinders its applicability.

Methods

We sought to develop an alternative method for CS reconstruction that only requires two gridding and one regridding operation in total, irrespective of the number of iterations. This proposed technique is evaluated on phantom images and whole-heart coronary MRI acquired using 3D radial trajectories, and compared to conventional CS reconstruction using CG algorithms in terms of quantitative vessel sharpness, vessel length, computation time, and convergence rate.

Results

Both CS reconstructions result in similar vessel length (P = 0.30) and vessel sharpness (P = 0.62). The per-iteration complexity of the proposed technique is approximately 3-fold lower than the conventional CS reconstruction (17.55 vs. 52.48 seconds in C++). Furthermore, for in-vivo datasets, the convergence rate of the proposed technique is faster (60±13 vs. 455±320 iterations) leading to a ∼23-fold reduction in reconstruction time.

Conclusions

The proposed reconstruction provides images of similar quality to the conventional CS technique in terms of removing artifacts, but at a much lower computational complexity.

Iterative method for predistortion of MRI gradient waveforms

The purpose of this work is to correct for transient gradient waveform errors in magnetic resonance imaging (MRI), whether from eddy currents, group delay, or gradient amplifier nonlinearities, which are known to affect image quality. An iterative method is proposed to minimize error between desired and measured gradient waveforms, whose success does not depend on accurate knowledge of the gradient system impulse response. The method was applied to half-pulse excitation for 2-D ultra-short echo time (UTE) imaging on a small animal MRI system and to spiral 2-D excitation on a human 7T MRI system. Predistorted gradient waveforms reduced temporal signal variation caused by excitation gradient trajectory errors in 2-D UTE, and improved the quality of excitation patterns produced by spiral excitation pulses. Iterative gradient predistortion is useful for minimizing transient gradient errors without requiring accurate characterization of the gradient system impulse response.

Simple and robust saturation-based slice selection for ultrashort echo time MRI

PURPOSE

To present a new method for localizing signal within a two-dimensional (2D) slice suitable for ultrashort echo time (UTE) imaging, called saturation-based UTE (sat-UTE). The new method digitally subtracts two acquisitions that are nonselectively excited with and without selective saturation of the slice of interest.

METHODS

Sat-UTE was compared with half-pulse and double-half pulse excited UTE within phantoms, as well as 3D-UTE within ex vivo femur and in vivo tibia. Numerical simulations were also used to quantify the effects of slice profile broadening and signal component amplitudes for quantitative UTE.

RESULTS

Sat-UTE is robust to suppress out-of-slice signal, and produces short T2 signal decay curves comparable to 3D-UTE, but has a lower signal to noise ratio efficiency compared with other slice-selective methods.

CONCLUSION

The proposed method is useful for fast, quantitative evaluation of short T2 signals, and is insensitive to gradient performance.

Real-time correction of rigid body motion-induced phase errors for diffusion-weighted steady-state free precession imaging

PURPOSE

Diffusion contrast in diffusion-weighted steady-state free precession magnetic resonance imaging (MRI) is generated through the constructive addition of signal from many coherence pathways. Motion-induced phase causes destructive interference which results in loss of signal magnitude and diffusion contrast. In this work, a three-dimensional (3D) navigator-based real-time correction of the rigid body motion-induced phase errors is developed for diffusion-weighted steady-state free precession MRI.

METHODS

The efficacy of the real-time prospective correction method in preserving phase coherence of the steady state is tested in 3D phantom experiments and 3D scans of healthy human subjects.

RESULTS

In nearly all experiments, the signal magnitude in images obtained with proposed prospective correction was higher than the signal magnitude in images obtained with no correction. In the human subjects, the mean magnitude signal in the data was up to 30% higher with prospective motion correction than without. Prospective correction never resulted in a decrease in mean signal magnitude in either the data or in the images.

CONCLUSIONS

The proposed prospective motion correction method is shown to preserve the phase coherence of the steady state in diffusion-weighted steady-state free precession MRI, thus mitigating signal magnitude losses that would confound the desired diffusion contrast.

Time-resolved angiography using inflow subtraction (TRAILS)

PURPOSE

A novel pseudo-continuous arterial spin labeling based angiographic method called Time-Resolved Angiography using InfLow Subtraction is introduced and used to acquire time-resolved whole-head angiographic data sets in healthy volunteers in a clinical feasible scan time of less than 5 min.

METHODS

Using this new method, in conjunction with a sliding window reconstruction, a temporal resolution of 7.2 ms with a low temporal footprint of 432 ms can be achieved.

RESULTS

Excellent vessel delineation compared to a time-of-flight MRA was demonstrated. Normal variations of the vascular system including the Circle of Willis (CoW) were identified using Time-Resolved Angiography Using Inflow Subtraction. Signal intensities were measured in various vascular segments to quantify the blood transit time.

CONCLUSION

In this feasibility study, we showed that Time-Resolved Angiography using InfLow Subtraction can be used to acquire hemodynamic information of the whole head in healthy volunteers with a high temporal and spatial resolution. Further studies in patients that suffer from vascular diseases to explore various flow patterns including longer transit time are needed.

Hybrid-SPRITE MRI

In a FID based frequency encoding MRI experiment the central part of k-space is not generally accessible due to the probe dead time. This portion of k-space is however crucial for image reconstruction. SPRITE (Single Point Ramped Imaging with T1 Enhancement), SPI with a linearly ramped phase encode gradient, has been employed to image short relaxation time systems for many years with great success. It is a robust imaging method in significant measure because it provides acquisition of high quality k-space origin data. We propose a new sampling scheme, termed hybrid-SPRITE, combining phase and frequency encoding to ensure high quality images with reduced acquisition times, reduced gradient duty cycle and increased sensitivity. In hybrid-SPRITE, numerous time domain points are collected to assist image reconstruction. An Inverse Non-uniform Discrete Fourier Transform (INDFT) is employed in 1D applications. A pseudo-polar grid is exploited in 2D hybrid-SPRITE for rapid and accurate image reconstruction.

A comparison of radial keyhole strategies for high spatial and temporal resolution 4D contrast-enhanced MRI in small animal tumor models

PURPOSE

Dynamic contrast-enhanced (DCE) MRI has been widely used as a quantitative imaging method for monitoring tumor response to therapy. The simultaneous challenges of increasing temporal and spatial resolution in a setting where the signal from the much smaller voxel is weaker have made this MR technique difficult to implement in small-animal imaging. Existing protocols employed in preclinical DCE-MRI acquire a limited number of slices resulting in potentially lost information in the third dimension. This study describes and compares a family of four-dimensional (3D spatial + time), projection acquisition, radial keyhole-sampling strategies that support high spatial and temporal resolution.

METHODS

The 4D method is based on a RF-spoiled, steady-state, gradient-recalled sequence with minimal echo time. An interleaved 3D radial trajectory with a quasi-uniform distribution of points in k-space was used for sampling temporally resolved datasets. These volumes were reconstructed with three different k-space filters encompassing a range of possible radial keyhole strategies. The effect of k-space filtering on spatial and temporal resolution was studied in a 5 mM CuSO(4) phantom consisting of a meshgrid with 350-μm spacing and in 12 tumors from three cell lines (HT-29, LoVo, MX-1) and a primary mouse sarcoma model (three tumors∕group). The time-to-peak signal intensity was used to assess the effect of the reconstruction filters on temporal resolution. As a measure of heterogeneity in the third dimension, the authors analyzed the spatial distribution of the rate of transport (K(trans)) of the contrast agent across the endothelium barrier for several different types of tumors.

RESULTS

Four-dimensional radial keyhole imaging does not degrade the system spatial resolution. Phantom studies indicate there is a maximum 40% decrease in signal-to-noise ratio as compared to a fully sampled dataset. T1 measurements obtained with the interleaved radial technique do not differ significantly from those made with a conventional Cartesian spin-echo sequence. A bin-by-bin comparison of the distribution of the time-to-peak parameter shows that 4D radial keyhole reconstruction does not cause significant temporal blurring when a temporal resolution of 9.9 s is used for the subsamples of the keyhole data. In vivo studies reveal substantial tumor heterogeneity in the third spatial dimension that may be missed with lower resolution imaging protocols.

CONCLUSIONS

Volumetric keyhole imaging with projection acquisition provides a means to increase spatiotemporal resolution and coverage over that provided by existing 2D Cartesian protocols. Furthermore, there is no difference in temporal resolution between the higher spatial resolution keyhole reconstruction and the undersampled projection data. The technique allows one to measure complex heterogeneity of kinetic parameters with isotropic, microscopic spatial resolution.

Blind retrospective motion correction of MR images

PURPOSE

Subject motion can severely degrade MR images. A retrospective motion correction algorithm, Gradient-based motion correction, which significantly reduces ghosting and blurring artifacts due to subject motion was proposed. The technique uses the raw data of standard imaging sequences; no sequence modifications or additional equipment such as tracking devices are required. Rigid motion is assumed.

METHODS

The approach iteratively searches for the motion trajectory yielding the sharpest image as measured by the entropy of spatial gradients. The vast space of motion parameters is efficiently explored by gradient-based optimization with a convergence guarantee.

RESULTS

The method has been evaluated on both synthetic and real data in two and three dimensions using standard imaging techniques. MR images are consistently improved over different kinds of motion trajectories. Using a graphics processing unit implementation, computation times are in the order of a few minutes for a full three-dimensional volume.

CONCLUSION

The presented technique can be an alternative or a complement to prospective motion correction methods and is able to improve images with strong motion artifacts from standard imaging sequences without requiring additional data.

Compressed sensing reconstruction for whole-heart imaging with 3D radial trajectories: a graphics processing unit implementation

A disadvantage of three-dimensional (3D) isotropic acquisition in whole-heart coronary MRI is the prolonged data acquisition time. Isotropic 3D radial trajectories allow undersampling of k-space data in all three spatial dimensions, enabling accelerated acquisition of the volumetric data. Compressed sensing (CS) reconstruction can provide further acceleration in the acquisition by removing the incoherent artifacts due to undersampling and improving the image quality. However, the heavy computational overhead of the CS reconstruction has been a limiting factor for its application. In this article, a parallelized implementation of an iterative CS reconstruction method for 3D radial acquisitions using a commercial graphics processing unit is presented. The execution time of the graphics processing unit-implemented CS reconstruction was compared with that of the C++ implementation, and the efficacy of the undersampled 3D radial acquisition with CS reconstruction was investigated in both phantom and whole-heart coronary data sets. Subsequently, the efficacy of CS in suppressing streaking artifacts in 3D whole-heart coronary MRI with 3D radial imaging and its convergence properties were studied. The CS reconstruction provides improved image quality (in terms of vessel sharpness and suppression of noise-like artifacts) compared with the conventional 3D gridding algorithm, and the graphics processing unit implementation greatly reduces the execution time of CS reconstruction yielding 34-54 times speed-up compared with C++ implementation.

Fast direct fourier reconstruction of radial and PROPELLER MRI data using the chirp transform algorithm on graphics hardware

PURPOSE

To develop and test a new algorithm for fast direct Fourier transform (DrFT) reconstruction of MR data on non-Cartesian trajectories composed of lines with equally spaced points.

THEORY AND METHODS

The DrFT, which is normally used as a reference in evaluating the accuracy of other reconstruction methods, can reconstruct images directly from non-Cartesian MR data without interpolation. However, DrFT reconstruction involves substantially intensive computation, which makes the DrFT impractical for clinical routine applications. In this article, the Chirp transform algorithm was introduced to accelerate the DrFT reconstruction of radial and Periodically Rotated Overlapping ParallEL Lines with Enhanced Reconstruction (PROPELLER) MRI data located on the trajectories that are composed of lines with equally spaced points. The performance of the proposed Chirp transform algorithm-DrFT algorithm was evaluated by using simulation and in vivo MRI data.

RESULTS

After implementing the algorithm on a graphics processing unit, the proposed Chirp transform algorithm-DrFT algorithm achieved an acceleration of approximately one order of magnitude, and the speed-up factor was further increased to approximately three orders of magnitude compared with the traditional single-thread DrFT reconstruction.

CONCLUSION

Implementation the Chirp transform algorithm-DrFT algorithm on the graphics processing unit can efficiently calculate the DrFT reconstruction of the radial and PROPELLER MRI data.

Distributed spirals: a new class of three-dimensional k-space trajectories

This work presents a new class of three-dimensional spiral based-trajectories for sampling magnetic resonance data. The distributed spirals trajectory efficiently traverses a cylinder or sphere or intermediate shape in k-space. The trajectory is shown to be nearly as efficient as a conventional stack of spirals trajectory in terms of scan time and signal-to-noise ratio, while reducing coherent aliasing in all three spatial directions and reducing Gibbs ringing due to the nature of collecting data from a sphere in k-space. The trajectory uses a single two-dimensional spiral waveform with the addition of a single orthogonal waveform which is scaled with each repetition, making it relatively easy to implement. Blurring from off-resonance only occurs in two dimensions due to the temporal nature of the sampling.

Simple method for MR gradient system characterization and k-space trajectory estimation

Fast imaging trajectories are used in MRI to speed up the acquisition process, but imperfections in the gradient system create artifacts in the reconstructed images. Artifacts result from the deviation between k-space trajectories achieved on the scanner and their original prescription. Measuring or approximating actual k-space trajectories with predetermined gradient timing delays reduces the artifacts, but are generally based on a specific trajectory and scan orientation. A single linear time-invariant characterization of the gradient system provides a method to predict k-space trajectories scanned in arbitrary orientations through convolution. This is done efficiently, by comparing the Fourier transforms of the input and measured waveforms of a single high-bandwidth test gradient waveform. This new method is tested for spiral, interleaved echo-planar, and three-dimensional cones imaging, demonstrating its ability to reduce reconstructed image artifacts for various k-space trajectories.

Description of chewing and food intake over the course of a meal

While the average frequency of chewing and food intake have been reported before, a detailed description of the pattern of chewing and the cumulative intake of food over the course of a meal have not. In order to achieve this goal, video recording of the maxillary-mandibular region of women eating food from a plate was synchronized with video recording of the plate and computer recording of the weight-loss of the plate. Video recording of chewing correlated strongly with chewing identified by magnetic tracking of jaw displacement in a test with chewing gum at three different frequencies, thus ensuring the validity of video recording of chewing. Weight-loss data were corrected by convolution algorithms, validated against human correction, using sliding window filtering to correct errors with video events as reference points. By use of this method, women ate on average 264 g of food over 114 min, they took an average of 51 mouthfuls during the meal and displayed on average 794 chews with 15 chews per chewing sequence. The number of mouthfuls decreased and the duration of the pauses after each mouthful increased in the middle of the meal and these measures were then restored. The ratio between chewing sequences and subsequent pauses remained stable although the weight of each mouthful decreased by the end of the meal, a measure that is hypothesized to be reflected in a decelerated speed of eating. The method allows this hypothesis to be tested and its implication for clinical intervention to be examined.

Efficient sample density estimation by combining gridding and an optimized kernel

The reconstruction of non-Cartesian k-space trajectories often requires the estimation of nonuniform sampling density. Particularly for 3D, this calculation can be computationally expensive. The method proposed in this work combines an iterative algorithm previously proposed by Pipe and Menon (Magn Reson Med 1999;41:179-186) with the optimal kernel design previously proposed by Johnson and Pipe (Magn Reson Med 2009;61:439-447). The proposed method shows substantial time reductions in estimating the densities of center-out trajectories, when compared with that of Johnson. It is demonstrated that, depending on the trajectory, the proposed method can provide reductions in execution time by factors of 12 to 85. The method is also shown to be robust in areas of high trajectory overlap, when compared with two analytical density estimation methods, producing a 10-fold increase in accuracy in one case. Initial conditions allow the proposed method to converge in fewer iterations and are shown to be flexible in terms of the accuracy of information supplied. The proposed method is not only one of the fastest and most accurate algorithms, it is also completely generic, allowing any arbitrary trajectory to be density compensated extemporaneously. The proposed method is also simple and can be implemented on parallel computing platforms in a straightforward manner.

Rigid body motion compensation for spiral projection imaging

Spiral projection imaging (SPI) is a 3D, spiral based magnetic resonance imaging (MRI) acquisition scheme that allows for self-navigated motion estimation of all six degrees-of-freedom. The trajectory, a set of spiral planes, is enhanced to accommodate motion tracking by adding orthogonal planes. Rigid-body motion tracking is accomplished by comparing the overlapping data and deducing the motion that is consistent with the comparisons. The accuracy of the proposed method is quantified for simulated data and for data collected using both a phantom and a volunteer. These tests were repeated to measure the effect of off-resonance blurring, coil sensitivity, gradient warping, undersampling, and nonrigid motion (e.g., neck). The artifacts of off-resonance, coils sensitivity, and gradient warping impose an unnotable effect on the accuracy of motion estimation. The worst mean accuracy is 0.15° and 0.20 mm for the phantom while the worst mean accuracy is 0.48° and 0.34 mm when imaging a brain, indicating that the nonrigid component in human subjects slightly degrades accuracy. When applied to in vivo motion, the proposed technique considerably reduces motion artifact.

Combined prospective and retrospective motion correction to relax navigator requirements

Prospective motion correction can prevent motion artifacts in magnetic resonance imaging of the brain. However, for high-resolution imaging, the technique relies on precise tracking of head motion. This precision is often limited by tracking noise, which leads to residual errors in the prospectively-corrected k-space data and artifacts in the image. This work shows that it is possible to estimate these tracking errors, and hence the true k-space sample locations, by applying a two-sided filter to the tracking data after imaging. A conjugate gradient reconstruction is compared to gridding as a means of using this information to retrospectively correct for the effects of the residual errors.

Fast, simple gradient delay estimation for spiral MRI

Timing delays between data acquisition and gradient transmission result in image degradation. This is especially true in spiral MRI, where delays can alter data in a nonuniform manner, generating significant artifact in the reconstructed data. The many methods that exist to mitigate these delays or measure the k-space coordinates require long measurement times, complicated analysis, specialized phantoms or hardware, or significant changes to the sequence of interest. A fast and simple method is proposed to measure delays on each gradient channel. It requires only minimal modification to an existing spiral sequence and can be used to measure independent delays on three gradient channels and any scan subject within six sequence repetition times. The effectiveness and accuracy of this method are analyzed.