Motion compensated whole-heart coronary cardiovascular magnetic resonance angiography using focused navigation (fNAV)

Respiratory-resolved 5D Flow MRI: in-vivo validation and respiratory dependent flow changes in healthy volunteers and patients with congenital heart disease

BACKGROUND

This study aimed to validate respiratory-resolved 5D flow MRI against real-time 2D phase contrast MRI, assess the impact of number of respiratory states, and measure the impact of respiration on hemodynamics in congenital heart disease (CHD) patients.

METHODS

Respiratory-resolved 5D flow MRI derived net and peak flow measurements were compared to real-time 2D phase contrast MRI derived measurements in 10 healthy volunteers. Pulmonary to systemic flow ratios (Qp:Qs) were measured in 19 CHD patients and aortopulmonary collateral burden was measured in 5 Fontan patients. Additionally, the impact of number of respiratory states on measured respiratory-driven net flow changes was investigated in 10 healthy volunteers and 19 CHD patients (shunt physiology, n=11, single ventricle disease (SVD), n=8).

RESULTS

There was good agreement between 5D flow MRI and real-time 2D phase contrast derived net and peak flow. Respiratory driven changes had good correlation (rho=0.64, p<0.001). In healthy volunteers, fewer than four respiratory states reduced measured respiratory driven flow changes in veins (5.2mL/cycle, p<0.001) and arteries (1.7mL/cycle, p=0.05). Respiration drove substantial venous net flow changes in SVD (64% change) and shunt patients (57% change). Respiration had significantly greater impact in SVD patients compared to shunt patients in the right and left pulmonary arteries (46% vs 15%, p=0.003 & 59% vs 20%, p=0.002). Qp:Qs varied by 37±24% over respiration in SVD patients and 12±20% in shunt patients. Aortopulmonary collateral burden varied by 118±84% over respiration in Fontan patients. The smallest collateral burden was measured during active inspiration in all patients and the greatest burden was during active expiration in 4 of 5 patients. Reduced respiratory resolution blunted measured flow changes in the caval veins of shunt and SVD patients (p<0.005).

CONCLUSIONS

Respiratory-resolved 5D flow MRI measurements agree with real-time 2D phase contrast. Venous measurements are sensitive to number of respiratory states, whereas arterial measurements are more robust. Respiration has substantial impact on caval vein flow, Qp:Qs, and collateral burden in CHD patients.

Stop moving: MR motion correction as an opportunity for artificial intelligence

Subject motion is a long-standing problem of magnetic resonance imaging (MRI), which can seriously deteriorate the image quality. Various prospective and retrospective methods have been proposed for MRI motion correction, among which deep learning approaches have achieved state-of-the-art motion correction performance. This survey paper aims to provide a comprehensive review of deep learning-based MRI motion correction methods. Neural networks used for motion artifacts reduction and motion estimation in the image domain or frequency domain are detailed. Furthermore, besides motion-corrected MRI reconstruction, how estimated motion is applied in other downstream tasks is briefly introduced, aiming to strengthen the interaction between different research areas. Finally, we identify current limitations and point out future directions of deep learning-based MRI motion correction.

Similarity-driven motion-resolved reconstruction for ferumoxytol-enhanced whole-heart MRI in congenital heart disease

A similarity-driven multi-dimensional binning algorithm (SIMBA) reconstruction of free-running cardiac magnetic resonance imaging data was previously proposed. While very efficient and fast, the original SIMBA focused only on the reconstruction of a single motion-consistent cluster, discarding the remaining data acquired. However, the redundant data clustered by similarity may be exploited to further improve image quality. In this work, we propose a novel compressed sensing (CS) reconstruction that performs an effective regularization over the clustering dimension, thanks to the integration of inter-cluster motion compensation (XD-MC-SIMBA). This reconstruction was applied to free-running ferumoxytol-enhanced datasets from 24 patients with congenital heart disease, and compared to the original SIMBA, the same XD-MC-SIMBA reconstruction but without motion compensation (XD-SIMBA), and a 5D motion-resolved CS reconstruction using the free-running framework (FRF). The resulting images were compared in terms of lung-liver and blood-myocardium sharpness, blood-myocardium contrast ratio, and visible length and sharpness of the coronary arteries. Moreover, an automated image quality score (IQS) was assigned using a pretrained deep neural network. The lung-liver sharpness and blood-myocardium sharpness were significantly higher in XD-MC-SIMBA and FRF. Consistent with these findings, the IQS analysis revealed that image quality for XD-MC-SIMBA was improved in 18 of 24 cases, compared to SIMBA. We successfully tested the hypothesis that multiple motion-consistent SIMBA clusters can be exploited to improve the quality of ferumoxytol-enhanced cardiac MRI when inter-cluster motion-compensation is integrated as part of a CS reconstruction.

Intra-bin correction and inter-bin compensation of respiratory motion in free-running five-dimensional whole-heart magnetic resonance imaging

BACKGROUND

Free-running cardiac and respiratory motion-resolved whole-heart five-dimensional (5D) cardiovascular magnetic resonance (CMR) can reduce scan planning and provide a means of evaluating respiratory-driven changes in clinical parameters of interest. However, respiratory-resolved imaging can be limited by user-defined parameters which create trade-offs between residual artifact and motion blur. In this work, we develop and validate strategies for both correction of intra-bin and compensation of inter-bin respiratory motion to improve the quality of 5D CMR.

METHODS

Each component of the reconstruction framework was systematically validated and compared to the previously established 5D approach using simulated free-running data (N = 50) and a cohort of 32 patients with congenital heart disease. The impact of intra-bin respiratory motion correction was evaluated in terms of image sharpness while inter-bin respiratory motion compensation was evaluated in terms of reconstruction error, compression of respiratory motion, and image sharpness. The full reconstruction framework (intra-acquisition correction and inter-acquisition compensation of respiratory motion [IIMC] 5D) was evaluated in terms of image sharpness and scoring of image quality by expert reviewers.

RESULTS

Intra-bin motion correction provides significantly (p < 0.001) sharper images for both simulated and patient data. Inter-bin motion compensation results in significant (p < 0.001) lower reconstruction error, lower motion compression, and higher sharpness in both simulated (10/11) and patient (9/11) data. The combined framework resulted in significantly (p < 0.001) sharper IIMC 5D reconstructions (End-expiration (End-Exp): 0.45 ± 0.09, End-inspiration (End-Ins): 0.46 ± 0.10) relative to the previously established 5D implementation (End-Exp: 0.43 ± 0.08, End-Ins: 0.39 ± 0.09). Similarly, image scoring by three expert reviewers was significantly (p < 0.001) higher using IIMC 5D (End-Exp: 3.39 ± 0.44, End-Ins: 3.32 ± 0.45) relative to 5D images (End-Exp: 3.02 ± 0.54, End-Ins: 2.45 ± 0.52).

CONCLUSION

The proposed IIMC reconstruction significantly improves the quality of 5D whole-heart MRI. This may be exploited for higher resolution or abbreviated scanning. Further investigation of the diagnostic impact of this framework and comparison to gold standards is needed to understand its full clinical utility, including exploration of respiratory-driven changes in physiological measurements of interest.

Combined free-running four-dimensional anatomical and flow magnetic resonance imaging with native contrast using Synchronization of Neighboring Acquisitions by Physiological Signals

BACKGROUND

Four-dimensional (4D) flow magnetic resonance imaging (MRI) often relies on the injection of gadolinium- or iron-oxide-based contrast agents to improve vessel delineation. In this work, a novel technique is developed to acquire and reconstruct 4D flow data with excellent dynamic visualization of blood vessels but without the need for contrast injection. Synchronization of Neighboring Acquisitions by Physiological Signals (SyNAPS) uses pilot tone (PT) navigation to retrospectively synchronize the reconstruction of two free-running three-dimensional radial acquisitions, to create co-registered anatomy and flow images.

METHODS

Thirteen volunteers and two Marfan syndrome patients were scanned without contrast agent using one free-running fast interrupted steady-state (FISS) sequence and one free-running phase-contrast MRI (PC-MRI) sequence. PT signals spanning the two sequences were recorded for retrospective respiratory motion correction and cardiac binning. The magnitude and phase images reconstructed, respectively, from FISS and PC-MRI, were synchronized to create SyNAPS 4D flow datasets. Conventional two-dimensional (2D) flow data were acquired for reference in ascending (AAo) and descending aorta (DAo). The blood-to-myocardium contrast ratio, dynamic vessel area, net volume, and peak flow were used to compare SyNAPS 4D flow with Native 4D flow (without FISS information) and 2D flow. A score of 0-4 was given to each dataset by two blinded experts regarding the feasibility of performing vessel delineation.

RESULTS

Blood-to-myocardium contrast ratio for SyNAPS 4D flow magnitude images (1.5 ± 0.3) was significantly higher than for Native 4D flow (0.7 ± 0.1, p < 0.01) and was comparable to 2D flow (2.3 ± 0.9, p = 0.02). Image quality scores of SyNAPS 4D flow from the experts (M.P.: 1.9 ± 0.3, E.T.: 2.5 ± 0.5) were overall significantly higher than the scores from Native 4D flow (M.P.: 1.6 ± 0.6, p = 0.03, E.T.: 0.8 ± 0.4, p < 0.01) but still significantly lower than the scores from the reference 2D flow datasets (M.P.: 2.8 ± 0.4, p < 0.01, E.T.: 3.5 ± 0.7, p < 0.01). The Pearson correlation coefficient between the dynamic vessel area measured on SyNAPS 4D flow and that from 2D flow was 0.69 ± 0.24 for the AAo and 0.83 ± 0.10 for the DAo, whereas the Pearson correlation between Native 4D flow and 2D flow measurements was 0.12 ± 0.48 for the AAo and 0.08 ± 0.39 for the DAo. Linear correlations between SyNAPS 4D flow and 2D flow measurements of net volume (r2 = 0.83) and peak flow (r2 = 0.87) were larger than the correlations between Native 4D flow and 2D flow measurements of net volume (r2 = 0.79) and peak flow (r2 = 0.76).

CONCLUSION

The feasibility and utility of SyNAPS were demonstrated for joint whole-heart anatomical and flow MRI without requiring electrocardiography gating, respiratory navigators, or contrast agents. Using SyNAPS, a high-contrast anatomical imaging sequence can be used to improve 4D flow measurements that often suffer from poor delineation of vessel boundaries in the absence of contrast agents.

Reproducibility of 3D thoracic aortic displacement from 3D cine balanced SSFP at 3 T without contrast enhancement

PURPOSE

Aortic motion has direct impact on the mechanical stresses acting on the aorta. In aortic disease, increased stiffness of the aorta may lead to decreased aortic motion over time, which could be a predictor for aortic dissection or rupture. This study investigates the reproducibility of obtaining 3D displacement and diameter maps quantified using accelerated 3D cine MRI at 3 T.

METHODS

A noncontrast-enhanced, free-breathing 3D cine sequence based on balanced SSFP and pseudo-spiral undersampling with high spatial isotropic resolution was developed (spatial/temporal resolution [1.6 mm]3 /67 ms). The thoracic aorta of 14 healthy volunteers was prospectively scanned three times at 3 T: twice on the same day and a third time 2 weeks later. Aortic displacement was calculated using iterative closest point nonrigid registration of manual segmentations of the 3D aorta at end-systole and mid-diastole. Interexamination and interobserver regional analysis of mean displacement for five regions of interest was performed using Bland-Altman analysis. Additionally, a complementary voxel-by-voxel analysis was done, allowing a more local inspection of the method.

RESULTS

No significant differences were found in mean and maximum displacement for any of the regions of interest for the interexamination and interobserver analysis. The maximum displacement measured in the lower half of the ascending aorta was 11.0 ± 3.4 mm (range: 3.0-17.5 mm) for the first scan. The smallest detectable change in mean displacement in the lower half of the ascending aorta was 3 mm.

CONCLUSION

Detailed 3D cine balanced SSFP at 3 T allows for reproducible quantification of systolic-diastolic mean aortic displacement within acceptable limits.

Rapid three-dimensional quantification of high-intensity plaques from coronary atherosclerosis T1-weighted characterization to predict periprocedural myocardial injury

BACKGROUND

High-intensity plaque (HIP) on magnetic resonance imaging (MRI) has been documented as a powerful predictor of periprocedural myocardial injury (PMI) following percutaneous coronary intervention (PCI). Despite the recent proposal of three-dimensional HIP quantification to enhance the predictive capability, the conventional pulse sequence, which necessitates the separate acquisition of anatomical reference images, hinders accurate three-dimensional segmentation along the coronary vasculature. Coronary atherosclerosis T1-weighted characterization (CATCH) enables the simultaneous acquisition of inherently coregistered dark-blood plaque and bright-blood coronary artery images. We aimed to develop a novel HIP quantification approach using CATCH and to ascertain its superior predictive performance compared to the conventional two-dimensional assessment based on plaque-to-myocardium signal intensity ratio (PMR).

METHODS

In this prospective study, CATCH MRI was conducted before elective stent implantation in 137 lesions from 125 patients. On CATCH images, dedicated software automatically generated tubular three-dimensional volumes of interest on the dark-blood plaque images along the coronary vasculature, based on the precisely matched bright-blood coronary artery images, and subsequently computed PMR and HIP volume (HIPvol). Specifically, HIPvol was calculated as the volume of voxels with signal intensity exceeding that of the myocardium, weighted by their respective signal intensities. PMI was defined as post-PCI cardiac troponin-T > 5 × the upper reference limit.

RESULTS

The entire analysis process was completed within 3 min per lesion. PMI occurred in 44 lesions. Based on the receiver operating characteristic curve analysis, HIPvol outperformed PMR for predicting PMI (C-statistics, 0.870 [95% CI, 0.805-0.936] vs. 0.787 [95% CI, 0.706-0.868]; p = 0.001). This result was primarily driven by the higher sensitivity HIPvol offered: 0.886 (95% CI, 0.754-0.962) vs. 0.750 for PMR (95% CI, 0.597-0.868; p = 0.034). Multivariable analysis identified HIPvol as an independent predictor of PMI (odds ratio, 1.15 per 10-μL increase; 95% CI, 1.01-1.30, p = 0.035).

CONCLUSIONS

Our semi-automated method of analyzing coronary plaque using CATCH MRI provided rapid HIP quantification. Three-dimensional assessment using this approach had a better ability to predict PMI than conventional two-dimensional assessment.

Comparison between diaphragmatic-navigated and self-navigated coronary magnetic resonance angiography at 3T in pediatric patients with congenital coronary artery anomalies

Background

Coronary magnetic resonance angiography (CMRA) is being increasingly used in pediatric patients with congenital coronary artery anomalies (CAAs). However, the data on the free-breathing self-navigation technique, which has the potential to simplify the acquisition plan with a high success rate at 3T, remain scarce. This study investigated the clinical application value of self-navigated (sNAV) CMRA at 3T in pediatric patients with suspected CAAs and compared it to conventional diaphragmatic-navigated (dNAV) CMRA.

Methods

From April 2019 to March 2022, we enrolled 65 pediatric patients (38 males and 27 females; mean age 8.5±4.4 years) with suspected CAAs in this prospective study. All patients underwent both dNAV and sNAV sequences in random order with gradient recalled echo (GRE) sequence during free breathing, with 39 (20 males and 19 females; mean age 10.2±3.6 years) of them additionally undergoing coronary computed tomography angiography (CCTA) or invasive coronary angiography (ICA). We measured and compared the success rate, scan time, visual score of the 9 main coronary artery segments, vessel sharpness, and vessel length between the two sequences. The diagnostic accuracy was compared using CCTA or ICA as a reference.

Results

The success rate of sNAV-CMRA (65/65, 100%) was higher than that of dNAV-CMRA (61/65, 93.8%) (P<0.001), and the scan time of sNAV-CMRA (7.3±2.5 min) was significantly shorter than that of dNAV-CMRA (9.1±3.6 min) (P=0.002). The acquisition efficiency of dNAV-CMRA was 40.5%±12.9%, while for sNAV-CMRA, 100% acquisition efficiency was achieved. There was no significant difference in vessel length of any of the coronary arteries, visual score, or vessel sharpness of the left circumflex coronary artery (LCX) between the two sequences (all P values >0.050). The visual score and vessel sharpness of the right coronary artery and left anterior descending coronary artery (LAD) were significantly improved in dNAV-CMRA compared with sNAV-CMRA (all P values <0.050). The sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) for the detection of CAAs were not significantly different between the two sequences (all P values >0.050).

Conclusions

Our findings demonstrated that both sNAV and dNAV in CMRA provide clinical application value in pediatric patients with CAAs and have similar diagnostic performance. Although the image quality of sNAV-CMRA is slightly inferior compared to that of dNAV-CMRA, sNAV-CMRA allows for a simpler scanning procedure.

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.

Free-running cardiac and respiratory motion-resolved 5D whole-heart coronary cardiovascular magnetic resonance angiography in pediatric cardiac patients using ferumoxytol

BACKGROUND

Coronary cardiovascular magnetic resonance angiography (CCMRA) of congenital heart disease (CHD) in pediatric patients requires accurate planning, adequate sequence parameter adjustments, lengthy scanning sessions, and significant involvement from highly trained personnel. Anesthesia and intubation are commonplace to minimize movements and control respiration in younger subjects. To address the above concerns and provide a single-click imaging solution, we applied our free-running framework for fully self-gated (SG) free-breathing 5D whole-heart CCMRA to CHD patients after ferumoxytol injection. We tested the hypothesis that spatial and motion resolution suffice to visualize coronary artery ostia in a cohort of CHD subjects, both for intubated and free-breathing acquisitions.

METHODS

In 18 pediatric CHD patients, non-electrocardiogram (ECG) triggered 5D free-running gradient echo CCMRA with whole-heart 1 mm3 isotropic spatial resolution was performed in seven minutes on a 1.5T CMR scanner. Eleven patients were anesthetized and intubated, while seven were breathing freely without anesthesia. All patients were slowly injected with ferumoxytol (4 mg/kg) over 15 minutes. Cardiac and respiratory motion-resolved 5D images were reconstructed with a fully SG approach. To evaluate the performance of motion resolution, visibility of coronary artery origins was assessed. Intubated and free-breathing patient sub-groups were compared for image quality using coronary artery length and conspicuity as well as lung-liver interface sharpness.

RESULTS

Data collection using the free-running framework was successful in all patients in less than 8 min; scan planning was very simple without the need for parameter adjustments, while no ECG lead placement and triggering was required. From the resulting SG 5D motion-resolved reconstructed images, coronary artery origins could be retrospectively extracted in 90% of the cases. These general findings applied to both intubated and free-breathing pediatric patients (no difference in terms of lung-liver interface sharpness), while image quality and coronary conspicuity between both cohorts was very similar.

CONCLUSIONS

A simple-to-use push-button framework for 5D whole-heart CCMRA was successfully employed in pediatric CHD patients with ferumoxytol injection. This approach, working without any external gating and for a wide range of heart rates and body sizes provided excellent definition of cardiac anatomy for both intubated and free-breathing patients.

Continuous cardiac thermometry via simultaneous catheter tracking and undersampled radial golden angle acquisition for radiofrequency ablation monitoring

The complexity of the MRI protocol is one of the factors limiting the clinical adoption of MR temperature mapping for real-time monitoring of cardiac ablation procedures and a push-button solution would ease its use. Continuous gradient echo golden angle radial acquisition combined with intra-scan motion correction and undersampled temperature determination could be a robust and more user-friendly alternative than the ultrafast GRE-EPI sequence which suffers from sensitivity to magnetic field susceptibility artifacts and requires ECG-gating. The goal of this proof-of-concept work is to establish the temperature uncertainty as well as the spatial and temporal resolutions achievable in an Agar-gel phantom and in vivo using this method. GRE radial golden angle acquisitions were used to monitor RF ablations in a phantom and in vivo in two sheep hearts with different slice orientations. In each case, 2D rigid motion correction based on catheter micro-coil signal, tracking its motion, was performed and its impact on the temperature imaging was assessed. The temperature uncertainty was determined for three spatial resolutions (1 × 1 × 3 mm³, 2 × 2 × 3 mm³, and 3 × 3 × 3 mm³) and three temporal resolutions (0.48, 0.72, and 0.97 s) with undersampling acceleration factors ranging from 2 to 17. The combination of radial golden angle GRE acquisition, simultaneous catheter tracking, intra-scan 2D motion correction, and undersampled thermometry enabled temperature monitoring in the myocardium in vivo during RF ablations with high temporal (< 1 s) and high spatial resolution. The temperature uncertainty ranged from 0.2 ± 0.1 to 1.8 ± 0.2 °C for the various temporal and spatial resolutions and, on average, remained superior to the uncertainty of an EPI acquisition while still allowing clinical monitoring of the RF ablation process. The proposed method is a robust and promising alternative to EPI acquisition to monitor in vivo RF cardiac ablations. Further studies remain required to improve the temperature uncertainty and establish its clinical applicability.

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.

Coronary Flow Assessment Using Accelerated 4D Flow MRI With Respiratory Motion Correction

Magnetic resonance imaging (MRI) can potentially be used for non-invasive screening of patients with stable angina pectoris to identify probable obstructive coronary artery disease. MRI-based coronary blood flow quantification has to date only been performed in a 2D fashion, limiting its clinical applicability. In this study, we propose a framework for coronary blood flow quantification using accelerated 4D flow MRI with respiratory motion correction and compressed sensing image reconstruction. We investigate its feasibility and repeatability in healthy subjects at rest. Fourteen healthy subjects received 8 times-accelerated 4D flow MRI covering the left coronary artery (LCA) with an isotropic spatial resolution of 1.0 mm³. Respiratory motion correction was performed based on 1) lung-liver navigator signal, 2) real-time monitoring of foot-head motion of the liver and LCA by a separate acquisition, and 3) rigid image registration to correct for anterior-posterior motion. Time-averaged diastolic LCA flow was determined, as well as time-averaged diastolic maximal velocity (V_MAX) and diastolic peak velocity (V_PEAK). 2D flow MRI scans of the LCA were acquired for reference. Scan-rescan repeatability and agreement between 4D flow MRI and 2D flow MRI were assessed in terms of concordance correlation coefficient (CCC) and coefficient of variation (CV). The protocol resulted in good visibility of the LCA in 11 out of 14 subjects (six female, five male, aged 28 ± 4 years). The other 3 subjects were excluded from analysis. Time-averaged diastolic LCA flow measured by 4D flow MRI was 1.30 ± 0.39 ml/s and demonstrated good scan-rescan repeatability (CCC/CV = 0.79/20.4%). Time-averaged diastolic V_MAX (17.2 ± 3.0 cm/s) and diastolic V_PEAK (24.4 ± 6.5 cm/s) demonstrated moderate repeatability (CCC/CV = 0.52/19.0% and 0.68/23.0%, respectively). 4D flow- and 2D flow-based diastolic LCA flow agreed well (CCC/CV = 0.75/20.1%). Agreement between 4D flow MRI and 2D flow MRI was moderate for both diastolic V_MAX and V_PEAK (CCC/CV = 0.68/20.3% and 0.53/27.0%, respectively). In conclusion, the proposed framework of accelerated 4D flow MRI equipped with respiratory motion correction and compressed sensing image reconstruction enables repeatable diastolic LCA flow quantification that agrees well with 2D flow MRI.

Pilot tone navigation for respiratory and cardiac motion-resolved free-running 5D flow MRI

Purpose

In this work, we integrated the pilot tone (PT) navigation system into a reconstruction framework for respiratory and cardiac motion‐resolved 5D flow. We tested the hypotheses that PT would provide equivalent respiratory curves, cardiac triggers, and corresponding flow measurements to a previously established self‐gating (SG) technique while being independent from changes to the acquisition parameters.

Methods

Fifteen volunteers and 9 patients were scanned with a free‐running 5D flow sequence, with PT integrated. Respiratory curves and cardiac triggers from PT and SG were compared across all subjects. Flow measurements from 5D flow reconstructions using both PT and SG were compared to each other and to a reference electrocardiogram‐gated and respiratory triggered 4D flow acquisition. Radial trajectories with variable readouts per interleave were also tested in 1 subject to compare cardiac trigger quality between PT and SG.

Results

The correlation between PT and SG respiratory curves were 0.95 ± 0.06 for volunteers and 0.95 ± 0.04 for patients. Heartbeat duration measurements in volunteers and patients showed a bias to electrocardiogram measurements of, respectively, 0.16 ± 64.94 ms and 0.01 ± 39.29 ms for PT versus electrocardiogram and of 0.24 ± 63.68 ms and 0.09 ± 32.79 ms for SG versus electrocardiogram. No significant differences were reported for the flow measurements between 5D flow PT and from 5D flow SG. A decrease in the cardiac triggering quality of SG was observed for increasing readouts per interleave, whereas PT quality remained constant.

Conclusion

PT has been successfully integrated in 5D flow MRI and has shown equivalent results to the previously described 5D flow SG technique, while being completely acquisition‐independent.