Free-breathing 3D cardiac MRI using iterative image-based respiratory motion correction

Inversion recovery and saturation recovery pulmonary vein MR angiography using an image based navigator fluoro trigger and variable-density 3D cartesian sampling with spiral-like order

Contrast enhanced pulmonary vein magnetic resonance angiography (PV CE-MRA) has value in atrial ablation pre-procedural planning. We aimed to provide high fidelity, ECG gated PV CE-MRA accelerated by variable density Cartesian sampling (VD-CASPR) with image navigator (iNAV) respiratory motion correction acquired in under 4 min. We describe its use in part during the global iodinated contrast shortage. VD-CASPR/iNAV framework was applied to ECG-gated inversion and saturation recovery gradient recalled echo PV CE-MRA in 65 patients (66 exams) using .15 mmol/kg Gadobutrol. Image quality was assessed by three physicians, and anatomical segmentation quality by two technologists. Left atrial SNR and left atrial/myocardial CNR were measured. 12 patients had CTA within 6 months of MRA. Two readers assessed PV ostial measurements versus CTA for intermodality/interobserver agreement. Inter-rater/intermodality reliability, reproducibility of ostial measurements, SNR/CNR, image, and anatomical segmentation quality was compared. The mean acquisition time was 3.58 ± 0.60 min. Of 35 PV pre-ablation datasets (34 patients), mean anatomical segmentation quality score was 3.66 ± 0.54 and 3.63 ± 0.55 as rated by technologists 1 and 2, respectively (p = 0.7113). Good/excellent anatomical segmentation quality (grade 3/4) was seen in 97% of exams. Each rated one exam as moderate quality (grade 2). 95% received a majority image quality score of good/excellent by three physicians. Ostial PV measurements correlated moderate to excellently with CTA (ICCs range 0.52-0.86). No difference in SNR was observed between IR and SR. High quality PV CE-MRA is possible in under 4 min using iNAV bolus timing/motion correction and VD-CASPR.

Retrospective motion correction through multi-average k-space data elimination (REMAKE) for free-breathing cardiac cine imaging

PURPOSE

To develop a motion-robust reconstruction technique for free-breathing cine imaging with multiple averages.

METHOD

Retrospective motion correction through multiple average k-space data elimination (REMAKE) was developed using iterative removal of k-space segments (from individual k-space samples) that contribute most to motion corruption while combining any remaining segments across multiple signal averages. A variant of REMAKE, termed REMAKE+, was developed to address any losses in SNR due to k-space information removal. With REMAKE+, multiple reconstructions using different initial conditions were performed, co-registered, and averaged. Both techniques were validated against clinical "standard" signal averaging reconstruction in a static phantom (with simulated motion) and 15 patients undergoing free-breathing cine imaging with multiple averages. Quantitative analysis of myocardial sharpness, blood/myocardial SNR, myocardial-blood contrast-to-noise ratio (CNR), as well as subjective assessment of image quality and rate of diagnostic quality images were performed.

RESULTS

In phantom, motion artifacts using "standard" (RMS error [RMSE]: 2.2 ± 0.5) were substantially reduced using REMAKE/REMAKE+ (RMSE: 1.5 ± 0.4/1.0 ± 0.4, p < 0.01). In patients, REMAKE/REMAKE+ led to higher myocardial sharpness (0.79 ± 0.09/0.79 ± 0.1 vs. 0.74 ± 0.12 for "standard", p = 0.004/0.04), higher image quality (1.8 ± 0.2/1.9 ± 0.2 vs. 1.6 ± 0.4 for "standard", p = 0.02/0.008), and a higher rate of diagnostic quality images (99%/100% vs. 94% for "standard"). Blood/myocardial SNR for "standard" (94 ± 30/33 ± 10) was higher vs. REMAKE (80 ± 25/28 ± 8, p = 0.002/0.005) and tended to be lower vs. REMAKE+ (105 ± 33/36 ± 12, p = 0.02/0.06). Myocardial-blood CNR for "standard" (61 ± 22) was higher vs. REMAKE (53 ± 19, p = 0.003) and lower vs. REMAKE+ (69 ± 24, p = 0.007).

CONCLUSIONS

Compared to "standard" signal averaging reconstruction, REMAKE and REMAKE+ provide improved myocardial sharpness, image quality, and rate of diagnostic quality images.

Accelerated 3D free-breathing high-resolution myocardial T1ρ mapping at 3 Tesla

PURPOSE

To develop a fast free-breathing whole-heart high-resolution myocardial T_1ρ mapping technique with robust spin-lock preparation that can be performed at 3 Tesla.

METHODS

An adiabatically excited continuous-wave spin-lock module, insensitive to field inhomogeneities, was implemented with an electrocardiogram-triggered low-flip angle spoiled gradient echo sequence with variable-density 3D Cartesian undersampling at a 3 Tesla whole-body scanner. A saturation pulse was performed at the beginning of each cardiac cycle to null the magnetization before T_1ρ preparation. Multiple T_1ρ -weighted images were acquired with T_1ρ preparations with different spin-lock times in an interleaved fashion. Respiratory self-gating approach was adopted along with localized autofocus to enable 3D translational motion correction of the data acquired in each heartbeat. After motion correction, multi-contrast locally low-rank reconstruction was performed to reduce undersampling artifacts. The accuracy and feasibility of the 3D T_1ρ mapping technique was investigated in phantoms and in vivo in 10 healthy subjects compared with the 2D T_1ρ mapping.

RESULTS

The 3D T_1ρ mapping technique provided similar phantom T_1ρ measurements in the range of 25-120 ms to the 2D T_1ρ mapping reference over a wide range of simulated heart rates. With the robust adiabatically excited continuous-wave spin-lock preparation, good quality 2D and 3D in vivo T_1ρ -weighted images and T_1ρ maps were obtained. Myocardial T_1ρ values with the 3D T_1ρ mapping were slightly longer than 2D breath-hold measurements (septal T_1ρ : 52.7 ± 1.4 ms vs. 50.2 ± 1.8 ms, P < 0.01).

CONCLUSION

A fast 3D free-breathing whole-heart T_1ρ mapping technique was proposed for T_1ρ quantification at 3 T with isotropic spatial resolution (2 mm³ ) and short scan time of ∼4.5 min.

Whole-heart non-rigid motion corrected coronary MRA with autofocus virtual 3D iNAV

PURPOSE

Respiratory motion-corrected coronary MR angiography (CMRA) has shown promise for assessing coronary disease. By incorporating coronal 2D image navigators (iNAVs), respiratory motion can be corrected for in a beat-to-beat basis using translational correction in the foot-head (FH) and right-left (RL) directions and in a bin-to-bin basis using non-rigid motion correction addressing the remaining FH, RL and anterior-posterior (AP) motion. However, with this approach beat-to-beat AP motion is not corrected for. In this work we investigate the effect of remaining beat-to-beat AP motion and propose a virtual 3D iNAV that exploits autofocus motion correction to enable beat-to-beat AP and improved RL intra-bin motion correction.

METHODS

Free-breathing 3D whole-heart CMRA was acquired using a 3-fold undersampled variable-density Cartesian trajectory. Beat-to-beat 3D translational respiratory motion was estimated from the 2D iNAVs in FH and RL directions, and in AP direction with autofocus assuming a linear relationship between FH and AP movement of the heart. Furthermore, motion in RL was also refined using autofocus. This virtual 3D (v3D) iNAV was incorporated in a non-rigid motion correction (NRMC) framework. The proposed approach was tested in 12 cardiac patients, and visible vessel length and vessel sharpness for the right (RCA) and left (LAD) coronary arteries were compared against 2D iNAV-based NRMC.

RESULTS

Average vessel sharpness and length in v3D iNAV NRMC was improved compared to 2D iNAV NRMC (vessel sharpness: RCA: 56 ± 1% vs 52 ± 11%, LAD: 49 ± 8% vs 49 ± 7%; visible vessel length: RCA: 5.98 ± 1.37 cm vs 5.81 ± 1.62 cm, LAD: 5.95 ± 1.85 cm vs 4.83 ± 1.56 cm), however these improvements were not statistically significant.

CONCLUSION

The proposed virtual 3D iNAV NRMC reconstruction further improved NRMC CMRA image quality by reducing artefacts arising from residual AP motion, however the level of improvement was subject-dependent.

Faster 3D saturation-recovery based myocardial T1 mapping using a reduced number of saturation points and denoising

PURPOSE

To accelerate the acquisition of free-breathing 3D saturation-recovery-based (SASHA) myocardial T1 mapping by acquiring fewer saturation points in combination with a post-processing 3D denoising technique to maintain high accuracy and precision.

METHODS

3D SASHA T1 mapping acquires nine T1-weighted images along the saturation recovery curve, resulting in long acquisition times. In this work, we propose to accelerate conventional cardiac T1 mapping by reducing the number of saturation points. High T1 accuracy and low standard deviation (as a surrogate for precision) is maintained by applying a 3D denoising technique to the T1-weighted images prior to pixel-wise T1 fitting. The proposed approach was evaluated on a T1 phantom and 20 healthy subjects, by varying the number of T1-weighted images acquired between three and nine, both prospectively and retrospectively. Following the results from the healthy subjects, three patients with suspected cardiovascular disease were acquired using five T1-weighted images. T1 accuracy and precision was determined for all the acquisitions before and after denoising.

RESULTS

In the T1 phantom, no statistical difference was found in terms of accuracy and precision for the different number of T1-weighted images before or after denoising (P = 0.99 and P = 0.99 for accuracy, P = 0.64 and P = 0.42 for precision, respectively). In vivo, both prospectively and retrospectively, the precision improved considerably with the number of T1-weighted images employed before denoising (P<0.05) but was independent on the number of T1-weighted images after denoising.

CONCLUSION

We demonstrate the feasibility of accelerating 3D SASHA T1 mapping by reducing the number of acquired T1-weighted images in combination with an efficient 3D denoising, without affecting accuracy and precision of T1 values.

Optimized respiratory-resolved motion-compensated 3D Cartesian coronary MR angiography

Purpose

To develop a robust and efficient reconstruction framework that provides high‐quality motion‐compensated respiratory‐resolved images from free‐breathing 3D whole‐heart Cartesian coronary magnetic resonance angiography (CMRA) acquisitions.

Methods

Recently, XD‐GRASP (eXtra‐Dimensional Golden‐angle RAdial Sparse Parallel MRI) was proposed to achieve 100% scan efficiency and provide respiratory‐resolved 3D radial CMRA images by exploiting sparsity in the respiratory dimension. Here, a reconstruction framework for Cartesian CMRA imaging is proposed, which provides respiratory‐resolved motion‐compensated images by incorporating 2D beat‐to‐beat translational motion information to increase sparsity in the respiratory dimension. The motion information is extracted from interleaved image navigators and is also used to compensate for 2D translational motion within each respiratory phase. The proposed Optimized Respiratory‐resolved Cartesian Coronary MR Angiography (XD‐ORCCA) method was tested on 10 healthy subjects and 2 patients with cardiovascular disease, and compared against XD‐GRASP.

Results

The proposed XD‐ORCCA provides high‐quality respiratory‐resolved images, allowing clear visualization of the right and left coronary arteries, even for irregular breathing patterns. Compared with XD‐GRASP, the proposed method improves the visibility and sharpness of both coronaries. Significant differences (p < .05) in visible vessel length and proximal vessel sharpness were found between the 2 methods. The XD‐GRASP method provides good‐quality images in the absence of intraphase motion. However, motion blurring is observed in XD‐GRASP images for respiratory phases with larger motion amplitudes and subjects with irregular breathing patterns.

Conclusion

A robust respiratory‐resolved motion‐compensated framework for Cartesian CMRA has been proposed and tested in healthy subjects and patients. The proposed XD‐ORCCA provides high‐quality images for all respiratory phases, independently of the regularity of the breathing pattern.

Improved respiratory self-navigation for 3D radial acquisitions through the use of a pencil-beam 2D-T2 -prep for free-breathing, whole-heart coronary MRA

PURPOSE

In respiratory self-navigation (SN), signal from static structures, such as the chest wall, may complicate motion detection or introduce post-correction artefacts. Suppressing signal from superfluous tissues may therefore improve image quality. We thus test the hypothesis that SN whole-heart coronary magnetic resonance angiography (MRA) will benefit from an outer-volume suppressing 2D-T2 -Prep and present both phantom and in vivo results.

METHODS

A 2D-T2 -Prep and a conventional T2 -Prep were used prior to a free-breathing 3D-radial SN sequence. Both techniques were compared by imaging a home-built moving cardiac phantom and by performing coronary MRA in nine healthy volunteers. Reconstructions were performed using both a reference-based and a reference-independent approach to motion tracking, along with several coil combinations. Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were compared, along with vessel sharpness (VS).

RESULTS

In phantoms, using the 2D-T2 -Prep increased SNR by 16% to 53% and mean VS by 8%; improved motion tracking precision was also achieved. In volunteers, SNR increased by an average of 29% to 33% in the blood pool and by 15% to 25% in the myocardium, depending on the choice of reconstruction coils and algorithm, and VS increased by 34%.

CONCLUSION

A 2D-T2 -Prep significantly improves image quality in both phantoms and volunteers when performing SN coronary MRA. Magn Reson Med 79:1293-1303, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Motion robust high resolution 3D free-breathing pulmonary MRI using dynamic 3D image self-navigator

PURPOSE

To achieve motion robust high resolution 3D free-breathing pulmonary MRI utilizing a novel dynamic 3D image navigator derived directly from imaging data.

METHODS

Five-minute free-breathing scans were acquired with a 3D ultrashort echo time (UTE) sequence with 1.25 mm isotropic resolution. From this data, dynamic 3D self-navigating images were reconstructed under locally low rank (LLR) constraints and used for motion compensation with one of two methods: a soft-gating technique to penalize the respiratory motion induced data inconsistency, and a respiratory motion-resolved technique to provide images of all respiratory motion states.

RESULTS

Respiratory motion estimation derived from the proposed dynamic 3D self-navigator of 7.5 mm isotropic reconstruction resolution and a temporal resolution of 300 ms was successful for estimating complex respiratory motion patterns. This estimation improved image quality compared to respiratory belt and DC-based navigators. Respiratory motion compensation with soft-gating and respiratory motion-resolved techniques provided good image quality from highly undersampled data in volunteers and clinical patients.

CONCLUSION

An optimized 3D UTE sequence combined with the proposed reconstruction methods can provide high-resolution motion robust pulmonary MRI. Feasibility was shown in patients who had irregular breathing patterns in which our approach could depict clinically relevant pulmonary pathologies. Magn Reson Med 79:2954-2967, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Highly-accelerated self-gated free-breathing 3D cardiac cine MRI: validation in assessment of left ventricular function

OBJECTIVE

This work presents a highly-accelerated, self-gated, free-breathing 3D cardiac cine MRI method for cardiac function assessment.

MATERIALS AND METHODS

A golden-ratio profile based variable-density, pseudo-random, Cartesian undersampling scheme was implemented for continuous 3D data acquisition. Respiratory self-gating was achieved by deriving motion signal from the acquired MRI data. A multi-coil compressed sensing technique was employed to reconstruct 4D images (3D+time). 3D cardiac cine imaging with self-gating was compared to bellows gating and the clinical standard breath-held 2D cine imaging for evaluation of self-gating accuracy, image quality, and cardiac function in eight volunteers. Reproducibility of 3D imaging was assessed.

RESULTS

Self-gated 3D imaging provided an image quality score of 3.4 ± 0.7 vs 4.0 ± 0 with the 2D method (p = 0.06). It determined left ventricular end-systolic volume as 42.4 ± 11.5 mL, end-diastolic volume as 111.1 ± 24.7 mL, and ejection fraction as 62.0 ± 3.1%, which were comparable to the 2D method, with bias ± 1.96 × SD of -0.8 ± 7.5 mL (p = 0.90), 2.6 ± 3.3 mL (p = 0.84) and 1.4 ± 6.4% (p = 0.45), respectively.

CONCLUSION

The proposed 3D cardiac cine imaging method enables reliable respiratory self-gating performance with good reproducibility, and provides comparable image quality and functional measurements to 2D imaging, suggesting that self-gated, free-breathing 3D cardiac cine MRI framework is promising for improved patient comfort and cardiac MRI scan efficiency.

4D respiratory motion-compensated image reconstruction of free-breathing radial MR data with very high undersampling

PURPOSE

To develop four-dimensional (4D) respiratory time-resolved MRI based on free-breathing acquisition of radial MR data with very high undersampling.

METHODS

We propose the 4D joint motion-compensated high-dimensional total variation (4D joint MoCo-HDTV) algorithm, which alternates between motion-compensated image reconstruction and artifact-robust motion estimation at multiple resolution levels. The algorithm is applied to radial MR data of the thorax and upper abdomen of 12 free-breathing subjects with acquisition times between 37 and 41 s and undersampling factors of 16.8. Resulting images are compared with compressed sensing-based 4D motion-adaptive spatio-temporal regularization (MASTeR) and 4D high-dimensional total variation (HDTV) reconstructions.

RESULTS

For all subjects, 4D joint MoCo-HDTV achieves higher similarity in terms of normalized mutual information and cross-correlation than 4D MASTeR and 4D HDTV when compared with reference 4D gated gridding reconstructions with 8.4 ± 1.1 times longer acquisition times. In a qualitative assessment of artifact level and image sharpness by two radiologists, 4D joint MoCo-HDTV reveals higher scores (P < 0.05) than 4D HDTV and 4D MASTeR at the same undersampling factor and the reference 4D gated gridding reconstructions, respectively.

CONCLUSIONS

4D joint MoCo-HDTV enables time-resolved image reconstruction of free-breathing radial MR data with undersampling factors of 16.8 while achieving low-streak artifact levels and high image sharpness. Magn Reson Med 77:1170-1183, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

An iterative approach to respiratory self-navigated whole-heart coronary MRA significantly improves image quality in a preliminary patient study

PURPOSE

In respiratory self-navigated coronary MRA, the selection of a reference position may have a direct effect on image quality. While end-expiration is commonly used as reference, it may be ill defined in cases of irregular breathing. Here, an iterative self-navigation approach that operates without a reference position was implemented and tested in healthy volunteers and patients.

METHODS

Data were acquired in 15 healthy volunteers and in 23 patients. Images obtained with end-expiratory self-navigation were compared with those obtained with the iterative approach that incorporates cross-correlation to iteratively minimize a global measure of respiratory displacement. Vessel sharpness, length, signal-to-noise ratio (SNR), and contrast-to-noise ratio (CNR) were evaluated while differences in breathing patterns between the two sub-groups were assessed, too.

RESULTS

Vessel sharpness and length were similar for both methods in healthy volunteers. In patients, a significant improvement in vessel sharpness and length was obtained using the iterative approach. SNR and CNR remained constant. While end-expiration was the most frequent respiratory phase in healthy volunteers (57.6 ± 16.2%), intermediate respiratory phases (43.4 ± 30.1%) were predominantly found in patients.

CONCLUSION

An iterative approach to respiratory motion correction in self-navigation may lead to significant improvements in coronary MRA image quality in patients with a less consistent end-expiratory respiratory phase.

Nonrigid autofocus motion correction for coronary MR angiography with a 3D cones trajectory

PURPOSE

To implement a nonrigid autofocus motion correction technique to improve respiratory motion correction of free-breathing whole-heart coronary magnetic resonance angiography acquisitions using an image-navigated 3D cones sequence.

METHODS

2D image navigators acquired every heartbeat are used to measure superior-inferior, anterior-posterior, and right-left translation of the heart during a free-breathing coronary magnetic resonance angiography scan using a 3D cones readout trajectory. Various tidal respiratory motion patterns are modeled by independently scaling the three measured displacement trajectories. These scaled motion trajectories are used for 3D translational compensation of the acquired data, and a bank of motion-compensated images is reconstructed. From this bank, a gradient entropy focusing metric is used to generate a nonrigid motion-corrected image on a pixel-by-pixel basis. The performance of the autofocus motion correction technique is compared with rigid-body translational correction and no correction in phantom, volunteer, and patient studies.

RESULTS

Nonrigid autofocus motion correction yields improved image quality compared to rigid-body-corrected images and uncorrected images. Quantitative vessel sharpness measurements indicate superiority of the proposed technique in 14 out of 15 coronary segments from three patient and two volunteer studies.

CONCLUSION

The proposed technique corrects nonrigid motion artifacts in free-breathing 3D cones acquisitions, improving image quality compared to rigid-body motion correction.

Advanced respiratory motion compensation for coronary MR angiography

Despite technical advances, respiratory motion remains a major impediment in a substantial amount of patients undergoing coronary magnetic resonance angiography (CMRA). Traditionally, respiratory motion compensation has been performed with a one-dimensional respiratory navigator positioned on the right hemi-diaphragm, using a motion model to estimate and correct for the bulk respiratory motion of the heart. Recent technical advancements has allowed for direct respiratory motion estimation of the heart, with improved motion compensation performance. Some of these new methods, particularly using image-based navigators or respiratory binning, allow for more advanced motion correction which enables CMRA data acquisition throughout most or all of the respiratory cycle, thereby significantly reducing scan time. This review describes the three components typically involved in most motion compensation strategies for CMRA, including respiratory motion estimation, gating and correction, and how these processes can be utilized to perform advanced respiratory motion compensation.