Deep Learning Automated Background Phase Error Correction for Abdominopelvic 4D Flow MRI

Background Four-dimensional (4D) flow MRI has the potential to provide hemodynamic insights for a variety of abdominopelvic vascular diseases, but its clinical utility is currently impaired by background phase error, which can be challenging to correct. Purpose To assess the feasibility of using deep learning to automatically perform image-based background phase error correction in 4D flow MRI and to compare its effectiveness relative to manual image-based correction. Materials and Methods A convenience sample of 139 abdominopelvic 4D flow MRI acquisitions performed between January 2016 and July 2020 was retrospectively collected. Manual phase error correction was performed using dedicated imaging software and served as the reference standard. After reserving 40 examinations for testing, the remaining examinations were randomly divided into training (86% [85 of 99]) and validation (14% [14 of 99]) data sets to train a multichannel three-dimensional U-Net convolutional neural network. Flow measurements were obtained for the infrarenal aorta, common iliac arteries, common iliac veins, and inferior vena cava. Statistical analyses included Pearson correlation, Bland-Altman analysis, and F tests with Bonferroni correction. Results A total of 139 patients (mean age, 47 years ± 14 [standard deviation]; 108 women) were included. Inflow-outflow correlation improved after manual correction (ρ = 0.94, P < .001) compared with that before correction (ρ = 0.50, P < .001). Automated correction showed similar results (ρ = 0.91, P < .001) and demonstrated very strong correlation with manual correction (ρ = 0.98, P < .001). Both correction methods reduced inflow-outflow variance, improving mean difference from -0.14 L/min (95% limits of agreement: -1.61, 1.32) (uncorrected) to 0.05 L/min (95% limits of agreement: -0.32, 0.42) (manually corrected) and 0.05 L/min (95% limits of agreement: -0.38, 0.49) (automatically corrected). There was no significant difference in inflow-outflow variance between manual and automated correction methods (P = .10). Conclusion Deep learning automated phase error correction reduced inflow-outflow bias and variance of volumetric flow measurements in four-dimensional flow MRI, achieving results comparable with manual image-based phase error correction. © RSNA, 2021 See also the editorial by Roldán-Alzate and Grist in this issue.

Quantification of the Hemodynamic Changes of Cirrhosis with Free-Breathing Self-Navigated MRI

BACKGROUND

Non-invasive assessment of the hemodynamic changes of cirrhosis might help guide management of patients with liver disease but are currently limited.

PURPOSE

To determine whether free-breathing 4D flow MRI can be used to quantify the hemodynamic effects of cirrhosis and introduce hydraulic circuit indexes of severity.

STUDY TYPE

Retrospective.

POPULATION

Forty-seven patients including 26 with cirrhosis.

FIELD STRENGTH/SEQUENCE

3 T/free-breathing 4D flow MRI with soft gating and golden-angle view ordering.

ASSESSMENT

Measurements of the supra-celiac abdominal aorta, supra-renal abdominal aorta (SRA), celiac trunk (CeT), superior mesenteric artery (SMA), splenic artery (SpA), common hepatic artery (CHA), portal vein (PV), and supra-renal inferior vena cava (IVC) were made by two radiologists. Measures of hepatic vascular resistance (hepatic arterial relative resistance [HARR]; portal resistive index [PRI]) were proposed and calculated.

STATISTICAL ANALYSIS

Bland-Altman, Pearson's correlation, Tukey's multiple comparison, and Cohen's kappa. P < 0.05 was considered significant.

RESULTS

Forty-four of 47 studies yielded adequate image quality for flow quantification (94%). Arterial structures showed high inter-reader concordance (range; ρ = 0.948-0.987) and the IVC (ρ = 0.972), with moderate concordance in the PV (ρ = 0.866). Conservation of mass analysis showed concordance between large vessels (SRA vs. IVC; ρ = 0.806), small vessels (celiac vs. CHA + SpA; ρ = 0.939), and across capillary beds (CeT + SMA vs. PV; ρ = 0.862). Splanchnic flow was increased in patients with portosystemic shunting (PSS) relative to control patients and patients with cirrhosis without PSS (P < 0.05, difference range 0.11-0.68 liter/m). HARR was elevated and PRI was decreased in patients with PSS (3.55 and 1.49, respectively) compared to both the control (2.11/3.18) and non-PSS (2.11/2.35) cohorts.

DATA CONCLUSION

4D flow MRI with self-navigation was technically feasible, showing promise in quantifying the hemodynamic effects of cirrhosis. Proposed quantitative metrics of hepatic vascular resistance correlated with PSS.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 2.

Diffusion-weighted double-echo steady-state with a three-dimensional cones trajectory for non-contrast-enhanced breast MRI

The image quality limitations of echo-planar diffusion-weighted imaging (DWI) are an obstacle to its widespread adoption in the breast. Steady-state DWI is an alternative DWI method with more robust image quality but its contrast for imaging breast cancer is not well-understood. The aim of this study was to develop and evaluate diffusion-weighted double-echo steady-state imaging with a three-dimensional cones trajectory (DW-DESS-Cones) as an alternative to conventional DWI for non-contrast-enhanced MRI in the breast. This prospective study included 28 women undergoing clinically indicated breast MRI and six asymptomatic volunteers. In vivo studies were performed at 3 T and included DW-DESS-Cones, DW-DESS-Cartesian, DWI, and CE-MRI acquisitions. Phantom experiments (diffusion phantom, High Precision Devices) and simulations were performed to establish framework for contrast of DW-DESS-Cones in comparison to DWI in the breast. Motion artifacts of DW-DESS-Cones were measured with artifact-to-noise ratio in volunteers and patients. Lesion-to-fibroglandular tissue signal ratios were measured, lesions were categorized as hyperintense or hypointense, and an image quality observer study was performed in DW-DESS-Cones and DWI in patients. Effect of DW-DESS-Cones method on motion artifacts was tested by mixed-effects generalized linear model. Effect of DW-DESS-Cones on signal in phantom was tested by quadratic regression. Correlation was calculated between DW-DESS-Cones and DWI lesion-to-fibroglandular tissue signal ratios. Inter-observer agreement was assessed with Gwet's AC. Simulations predicted hyperintensity of lesions with DW-DESS-Cones but at a 3% to 67% lower degree than with DWI. Motion artifacts were reduced with DW-DESS-Cones versus DW-DESS-Cartesian (p < 0.05). Lesion-to-fibroglandular tissue signal ratios were not correlated between DW-DESS-Cones and DWI (r = 0.25, p = 0.38). Concordant hyperintensity/hypointensity was observed between DW-DESS-Cones and DWI in 11/14 lesions. DW-DESS-Cones improved sharpness, distortion, and overall image quality versus DWI. DW-DESS-Cones may be able to eliminate motion artifacts in the breast allowing for investigation of higher degrees of steady-state diffusion weighting. Malignant breast lesions in DW-DESS-Cones demonstrated hyperintensity with respect to surrounding tissue without an injection of contrast. LEVEL OF EVIDENCE: 2. TECHNICAL EFFICACY STAGE: 1.

Rosette Trajectories Enable Ungated, Motion-Robust, Simultaneous Cardiac and Liver T2 * Iron Assessment

BACKGROUND

Quantitative T₂ * MRI is the standard of care for the assessment of iron overload. However, patient motion corrupts T₂ * estimates.

PURPOSE

To develop and evaluate a motion-robust, simultaneous cardiac and liver T₂ * imaging approach using non-Cartesian, rosette sampling and a model-based reconstruction as compared to clinical-standard Cartesian MRI.

STUDY TYPE

Prospective.

PHANTOM/POPULATION

Six ferumoxytol-containing phantoms (26-288 μg/mL). Eight healthy subjects and 18 patients referred for clinically indicated iron overload assessment.

FIELD STRENGTH/SEQUENCE

1.5T, 2D Cartesian and rosette gradient echo (GRE) ASSESSMENT: GRE T₂ * values were validated in ferumoxytol phantoms. In healthy subjects, test-retest and spatial coefficient of variation (CoV) analysis was performed during three breathing conditions. Cartesian and rosette T₂ * were compared using correlation and Bland-Altman analysis. Images were rated by three experienced radiologists on a 5-point scale.

STATISTICAL TESTS

Linear regression, analysis of variance (ANOVA), and paired Student's t-testing were used to compare reproducibility and variability metrics in Cartesian and rosette scans. The Wilcoxon rank test was used to assess reader score comparisons and reader reliability was measured using intraclass correlation analysis.

RESULTS

Rosette R2* (1/T₂ *) was linearly correlated with ferumoxytol concentration (r² = 1.00) and not significantly different than Cartesian values (P = 0.16). During breath-holding, ungated rosette liver and heart T₂ * had lower spatial CoV (liver: 18.4 ± 9.3% Cartesian, 8.8% ± 3.4% rosette, P = 0.02, heart: 37.7% ± 14.3% Cartesian, 13.4% ± 1.7% rosette, P = 0.001) and higher-quality scores (liver: 3.3 [3.0-3.6] Cartesian, 4.7 [4.1-4.9] rosette, P = 0.005, heart: 3.0 [2.3-3] Cartesian, 4.5 [3.8-5.0] rosette, P = 0.005) compared to Cartesian values. During free-breathing and failed breath-holding, Cartesian images had very poor to average image quality with significant artifacts, whereas rosette remained very good, with minimal artifacts (P = 0.001).

DATA CONCLUSION

Rosette k-sampling with a model-based reconstruction offers a clinically useful motion-robust T₂ * mapping approach for iron quantification. J. MAGN. RESON. IMAGING 2020;52:1688-1698.

Direct measurement of atrioventricular valve regurgitant jets using 4D flow cardiovascular magnetic resonance is accurate and reliable for children with congenital heart disease: a retrospective cohort study

BACKGROUND

3D-time resolved flow (4DF) cardiovascular magnetic resonance (CMR) with retrospective analysis of atrioventricular valve regurgitation (AVVR) allows for internal validation by multiple direct and indirect methods. Limited data exist on direct measurement of AVVR by 4DF CMR in pediatric congenital heart disease (CHD). We aimed to validate direct measurement of the AVVR jet as accurate and reliable compared to the volumetric method (clinical standard by 2D CMR) and as a superior method of internal validation than the annular inflow method.

METHODS

We identified 44 consecutive patients with diverse CHD referred for evaluation of AVVR by CMR. 1.5 T or 3 T scanners, intravenous contrast, and a combination of parallel imaging and compressed sensing were used. Four methods of measuring AVVR volume (RVol) were used: volumetric method (VOL; the clinical standard) = stroke volume by 2D balanced steady-state free precession - semilunar valve forward flow (SLFF); annular inflow method (AIM) = atrioventricular valve forward flow [AVFF] - semilunar valve net flow (SLNF); and direct measurement (JET). AVFF was measured using static and retrospective valve tracking planes. SLFF, SLNF, AVFF, and JET were measured by 4DF phase contrast. Regurgitant fraction was calculated as [RVol/(RVol+SLNF)]× 100. Statistical methods included Spearman, Wilcoxon rank sum test/Student paired t-test, Bland Altman analysis, and intra-class coefficient (ICC), where appropriate.

RESULTS

Regurgitant fraction by JET strongly correlated with the indirect methods (VOL and AIM) (ρ = 0.73-0.80, p < 0.001) and was similar to VOL with a median difference (interquartile range) of - 1.5% (- 8.3-7.2%; p = 0.624). VOL had weaker correlations with AIM and JET (ρ = 0.69-0.73, p < 0.001). AIM underestimated RF by 3.6-6.9% compared to VOL and JET, p < 0.03. Intra- and inter- observer reliability were excellent for all methods (ICC 0.94-0.99). The mean (±standard deviation) inter-observer difference for VOL was 2.4% (±5.1%), p < 0.05.

CONCLUSIONS

In a diverse cohort of pediatric CHD, measurement of AVVR using JET is accurate and reliable to VOL and is a superior method of internal validation compared to AIM. This study supports use of 4DF CMR for measurement of AVVR, obviating need for expert prospective prescription during image acquisition by 2D CMR.

Combined 5-minute double-echo in steady-state with separated echoes and 2-minute proton-density-weighted 2D FSE sequence for comprehensive whole-joint knee MRI assessment

BACKGROUND

Clinical knee MRI protocols require upwards of 15 minutes of scan time.

PURPOSE/HYPOTHESIS

To compare the imaging appearance of knee abnormalities depicted with a 5-minute 3D double-echo in steady-state (DESS) sequence with separate echo images, with that of a routine clinical knee MRI protocol. A secondary goal was to compare the imaging appearance of knee abnormalities depicted with 5-minute DESS paired with a 2-minute coronal proton-density fat-saturated (PDFS) sequence.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-six consecutive patients (19 male) referred for a routine knee MRI.

FIELD STRENGTH/SEQUENCES

DESS and PDFS at 3T.

ASSESSMENT

Five musculoskeletal radiologists evaluated all images for the presence of internal knee derangement using DESS, DESS+PDFS, and the conventional imaging protocol, and their associated diagnostic confidence of the reading.

STATISTICAL TESTS

Differences in positive and negative percent agreement (PPA and NPA, respectively) and 95% confidence intervals (CIs) for DESS and DESS+PDFS compared with the conventional protocol were calculated and tested using exact McNemar tests. The percentage of observations where DESS or DESS+PDFS had equivalent confidence ratings to DESS+Conv were tested with exact symmetry tests. Interreader agreement was calculated using Krippendorff's alpha.

RESULTS

DESS had a PPA of 90% (88-92% CI) and NPA of 99% (99-99% CI). DESS+PDFS had increased PPA of 99% (95-99% CI) and NPA of 100% (99-100% CI) compared with DESS (both P < 0.001). DESS had equivalent diagnostic confidence to DESS+Conv in 94% of findings, whereas DESS+PDFS had equivalent diagnostic confidence in 99% of findings (both P < 0.001). All readers had moderate concordance for all three protocols (Krippendorff's alpha 47-48%).

DATA CONCLUSION

Both 1) 5-minute 3D-DESS with separated echoes and 2) 5-minute 3D-DESS paired with a 2-minute coronal PDFS sequence depicted knee abnormalities similarly to a routine clinical knee MRI protocol, which may be a promising technique for abbreviated knee MRI.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018.

Targeted rapid knee MRI exam using T2 shuffling

BACKGROUND

MRI is commonly used to evaluate pediatric musculoskeletal pathologies, but same-day/near-term scheduling and short exams remain challenges.

PURPOSE

To investigate the feasibility of a targeted rapid pediatric knee MRI exam, with the goal of reducing cost and enabling same-day MRI access.

STUDY TYPE

A cost effectiveness study done prospectively.

SUBJECTS

Forty-seven pediatric patients.

FIELD STRENGTH/SEQUENCE

3T. The 10-minute protocol was based on T2 Shuffling, a four-dimensional acquisition and reconstruction of images with variable T2 contrast, and a T1 2D fast spin-echo (FSE) sequence. A distributed, compressed sensing-based reconstruction was implemented on a four-node high-performance compute cluster and integrated into the clinical workflow.

ASSESSMENT

In an Institutional Review Board-approved study with informed consent/assent, we implemented a targeted pediatric knee MRI exam for assessing pediatric knee pain. Pediatric patients were subselected for the exam based on insurance plan and clinical indication. Over a 2-year period, 47 subjects were recruited for the study and 49 MRIs were ordered. Date and time information was recorded for MRI referral, registration, and completion. Image quality was assessed from 0 (nondiagnostic) to 5 (outstanding) by two readers, and consensus was subsequently reached.

STATISTICAL TESTS

A Wilcoxon rank-sum test assessed the null hypothesis that the targeted exam times compared with conventional knee exam times were unchanged.

RESULTS

Of the 49 cases, 20 were completed on the same day as exam referral. Median time from registration to exam completion was 18.7 minutes. Median reconstruction time for T2 Shuffling was reduced from 18.9 minutes to 95 seconds using the distributed implementation. Technical fees charged for the targeted exam were one-third that of the routine clinical knee exam. No subject had to return for additional imaging.

DATA CONCLUSION

The targeted knee MRI exam is feasible and reduces the imaging time, cost, and barrier to same-day MRI access for pediatric patients.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 6 J. Magn. Reson. Imaging 2019.

Flow dynamics of type II endoleaks can determine sac expansion after endovascular aneurysm repair using four-dimensional flow-sensitive magnetic resonance imaging analysis

OBJECTIVE

The objective of this study was to investigate the hemodynamic parameters of type II endoleaks (T2ELs) to predict sac expansion using four-dimensional flow-sensitive magnetic resonance imaging (4D-flow MRI) analysis.

METHODS

Patients who underwent endovascular aneurysm repair (EVAR) and were diagnosed with a T2EL were included in the study. Using 4D-flow MRI at 7 days, the peak flow velocity and amplitude of dynamics of blood flow per minute were measured in each T2EL vessel. The peak flow velocity was defined as the maximum of the absolute value of the blood flow velocity. The amplitude of dynamics of blood flow in the tributary arteries was defined as the sum of the absolute values of the inflow and outflow volume in each vessel. The amplitude of dynamics of blood flow in the tributary arteries per sac was calculated in each sac. The aneurysm sac diameter was measured by computed tomography (CT) at 1 year. The patients were divided into two groups according to the presence or absence of sac expansion.

RESULTS

Of 155 patients who underwent EVAR, both CT angiography and 4D-flow MRI were performed in 107 patients at 7 days after EVAR. Among them, 39 (36.4%) were found to have a T2EL, of whom 28 were re-evaluated with CT angiography and 4D-flow at 1 year; 7 patients had expanding sacs (expanding group), whereas 21 had nonexpanding sacs (not-expanding group). At 7 days, 28 patients had 80 T2EL vessels detected by 4D-flow MRI, of which 39 vessels (48.8%) had stopped flowing at 1 year (transient vessels); 41 vessels (51.3%) had sustained flow (persistent vessels). The persistent vessels had significantly larger peak flow velocity and amplitude of dynamics of blood flow. The comprehensive analysis of T2EL vessels per sac identified that the amplitude of dynamics of blood flow in the tributary arteries per sac was significantly higher in the expanding group than in the not-expanding group. A receiver operating characteristic curve analysis revealed that the sensitivity and specificity of sac enlargement at a cutoff value of 3750 mm³/min were 85.7% and 76.2%, respectively.

CONCLUSIONS

The fate of aneurysm sacs with T2ELs after EVAR has remained difficult to predict. A comprehensive analysis of concurrent multiple T2EL vessels using 4D-flow MRI analysis may enable prediction of the sac expansion after EVAR.

4D flow MRI quantification of mitral and tricuspid regurgitation: Reproducibility and consistency relative to conventional MRI

BACKGROUND

In patients with mitral or tricuspid valve regurgitation, evaluation of regurgitant severity is essential for determining the need for surgery. While transthoracic echocardiography is widely accessible, it has limited reproducibility for grading inlet valve regurgitation. Multiplanar cardiac MRI is the quantitative standard but requires specialized local expertise, and is thus not widely available. Volumetric 4D flow MRI has potential for quantitatively grading the severity of inlet valve regurgitation in adult patients.

PURPOSE

To evaluate the accuracy and reproducibility of volumetric 4D flow MRI for quantification of inlet valvular regurgitation compared to conventional multiplanar MRI, which may simplify and improve accessibility of cardiac MRI.

STUDY TYPE

This retrospective, HIPAA-compliant imaging-based comparison study was conducted at a single institution.

SUBJECTS

Twenty-one patients who underwent concurrent multiplanar and 4D flow cardiac MRI between April 2015 and January 2017.

FIELD STRENGTH/SEQUENCES

3T; steady-state free-precession (SSFP), 2D phase contrast (2D-PC), and postcontrast 4D flow.

ASSESSMENT

We evaluated the intertechnique (4D flow vs. 2D-PC), intermethod (direct vs. indirect measurement), interobserver and intraobserver reproducibility of measurements of regurgitant flow volume (RFV), fraction (RF), and volume (RVol).

STATISTICAL TESTS

Statistical analysis included Pearson correlation, Bland-Altman statistics, and intraclass correlation coefficients.

RESULTS

There was high concordance between 4D flow and multiplanar MRI, whether using direct or indirect methods of quantifying regurgitation (r = 0.813-0.985). Direct interrogation of the regurgitant jet with 4D flow showed high intraobserver consistency (r = 0.976-0.999) and interobserver consistency (r = 0.861-0.992), and correlated well with traditional indirect measurements obtained as the difference between stroke volume and forward outlet valve flow.

DATA CONCLUSION

4D flow MRI provides highly reproducible measurements of mitral and tricuspid regurgitant volume, and may be used in place of conventional multiplanar MRI.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;48:1147-1158.

High-resolution 3D volumetric contrast-enhanced MR angiography with a blood pool agent (ferumoxytol) for diagnostic evaluation of pediatric brain arteriovenous malformations

OBJECTIVE Patients with brain arteriovenous malformations (AVMs) often require repeat imaging with MRI or MR angiography (MRA), CT angiography (CTA), and digital subtraction angiography (DSA). The ideal imaging modality provides excellent vascular visualization without incurring added risks, such as radiation exposure. The purpose of this study is to evaluate the performance of ferumoxytol-enhanced MRA using a high-resolution 3D volumetric sequence (fe-SPGR) for visualizing and grading pediatric brain AVMs in comparison with CTA and DSA, which is the current imaging gold standard. METHODS In this retrospective cohort study, 21 patients with AVMs evaluated by fe-SPGR, CTA, and DSA between April 2014 and August 2017 were included. Two experienced raters graded AVMs using Spetzler-Martin criteria on all imaging studies. Lesion conspicuity (LC) and diagnostic confidence (DC) were assessed using a 5-point Likert scale, and interrater agreement was determined. The Kruskal-Wallis test was performed to assess the raters' grades and scores of LC and DC, with subsequent post hoc pairwise comparisons to assess for statistically significant differences between pairs of groups at p < 0.05. RESULTS Assigned Spetzler-Martin grades for AVMs on DSA, fe-SPGR, and CTA were not significantly different (p = 0.991). LC and DC scores were higher with fe-SPGR than with CTA (p < 0.05). A significant difference in LC scores was found between CTA and fe-SPGR (p < 0.001) and CTA and DSA (p < 0.001) but not between fe-SPGR and DSA (p = 0.146). A significant difference in DC scores was found among DSA, fe-SPGR, and CTA (p < 0.001) and between all pairs of the groups (p < 0.05). Interrater agreement was good to very good for all image groups (κ = 0.77-1.0, p < 0.001). CONCLUSIONS Fe-SPGR performed robustly in the diagnostic evaluation of brain AVMs, with improved visual depiction of AVMs compared with CTA and comparable Spetzler-Martin grading relative to CTA and DSA.

Volumetric segmentation-free method for rapid visualization of vascular wall shear stress using 4D flow MRI

PURPOSE

To develop a rapid segmentation-free method to visualize and compute wall shear stress (WSS) throughout the aorta using 4D Flow MRI data. WSS is the drag force-per-area the vessel endothelium exerts on luminal blood; abnormal levels of WSS are associated with cardiovascular pathologies. Previous methods for computing WSS are bottlenecked by labor-intensive manual segmentation of vessel boundaries. A rapid automated segmentation-free method for computing WSS is presented.

THEORY AND METHODS

Shear stress is the dot-product of the viscous stress tensor and the inward normal vector. The inward normal vectors are approximated as the gradient of fluid speed at every voxel. Subsequently, a 4D map of shear stress is computed as the partial derivatives of velocity with respect to the inward normal vectors. We highlight the shear stress near the wall by fusing visualization with edge-emphasized anatomical data.

RESULTS

As a proof-of-concept, four cases with aortic pathologies are presented. Visualization allows for rapid localization of pathologic WSS. Subsequent analysis of these pathological regions enables quantification of WSS. Average WSS during peak systole measures approximately 50-60 cPa in nonpathological regions of the aorta and is elevated in regions of stenosis, coarctation, and dissection. WSS is reduced in regions of aneurysm.

CONCLUSION

A volumetric technique for calculation and visualization of WSS from 4D Flow MRI data is presented. Traditional labor-intensive methods for WSS rely on explicit manual segmentation of vessel boundaries before visualization. This automated volumetric strategy for visualization and quantification of WSS may facilitate its clinical translation.

T2 shuffling: Sharp, multicontrast, volumetric fast spin-echo imaging

PURPOSE

A new acquisition and reconstruction method called T2 Shuffling is presented for volumetric fast spin-echo (three-dimensional [3D] FSE) imaging. T2 Shuffling reduces blurring and recovers many images at multiple T2 contrasts from a single acquisition at clinically feasible scan times (6-7 min).

THEORY AND METHODS

The parallel imaging forward model is modified to account for temporal signal relaxation during the echo train. Scan efficiency is improved by acquiring data during the transient signal decay and by increasing echo train lengths without loss in signal-to-noise ratio (SNR). By (1) randomly shuffling the phase encode view ordering, (2) constraining the temporal signal evolution to a low-dimensional subspace, and (3) promoting spatio-temporal correlations through locally low rank regularization, a time series of virtual echo time images is recovered from a single scan. A convex formulation is presented that is robust to partial voluming and radiofrequency field inhomogeneity.

RESULTS

Retrospective undersampling and in vivo scans confirm the increase in sharpness afforded by T2 Shuffling. Multiple image contrasts are recovered and used to highlight pathology in pediatric patients. A proof-of-principle method is integrated into a clinical musculoskeletal imaging workflow.

CONCLUSION

The proposed T2 Shuffling method improves the diagnostic utility of 3D FSE by reducing blurring and producing multiple image contrasts from a single scan. Magn Reson Med 77:180-195, 2017. © 2016 Wiley Periodicals, Inc.

Feasibility of ferumoxytol-enhanced neonatal and young infant cardiac MRI without general anesthesia

PURPOSE

To assess the feasibility of ferumoxytol-enhanced anesthesia-free cardiac MRI in neonates and young infants for complex congenital heart disease (CHD).

MATERIALS AND METHODS

With Institutional Review Board approval, 21 consecutive neonates and young infants (1 day to 11 weeks old; median age of 3 days) who underwent a rapid two-sequence (MR angiography [MRA] and four-dimensional [4D] flow) MRI protocol with intravenous ferumoxytol without sedation (n = 17) or light sedation (n = 4) at 3 Tesla (T) (except one case at 1.5T) between June 2014 and February 2016 were retrospectively identified. Medical records were reviewed for indication, any complications, if further diagnostic imaging was performed after MRI, and surgical findings. Two radiologists scored the images in two sessions on a 5-point scale for overall image quality and delineation of various anatomical structures. Confidence interval of proportions for likelihood of requiring additional diagnostic imaging after MRI was determined. For the possibility of reducing the protocol to a single rapid sequence, Wilcoxon-rank sum test was used to assess whether 4D flow and MRA significantly differed in anatomical delineation.

RESULTS

One of 21 patients (4.8%, 80% confidence interval 0-11%) required additional imaging, a computed tomography angiography to assess lung parenchyma and peripheral pulmonary arteries. Only 1 of 13 patients (7.7%) with operative confirmation had a minor discrepancy between radiology and operative reports (80% confidence interval 0-17%). 4D flow was significantly superior to MRA (P < 0.05) for the evaluation of systemic arteries, valves, ventricular trabeculae, and overall quality. Using Cohen's kappa coefficient, there was good interobserver agreement for the evaluation of systemic arteries by 4D flow (κ = 0.782), and systemic veins and pulmonary arteries by MRA (κ > 0.6). Overall 4D flow measurements (mean κ = 0.64-0.74) had better internal agreement compared with MRA (mean κ = 0.30-0.64).

CONCLUSION

Ferumoxytol-enhanced cardiac MRI, without anesthesia, is feasible for the evaluation of complex CHD in neonates and young infants, with a low likelihood of need for additional diagnostic studies. The decreased risk by avoiding anesthesia must be balanced against the potential for adverse reactions with ferumoxytol.

LEVEL OF EVIDENCE

2 J. MAGN. RESON. IMAGING 2017;45:1407-1418.

Non-invasive evaluation of pulmonary arterial blood flow and wall shear stress in pulmonary arterial hypertension with 3D phase contrast magnetic resonance imaging

Background

Recently, time-resolved 3D phase contrast magnetic resonance imaging (4D-flow) allows flow dynamics in patients with pulmonary arterial hypertension to be measured. Abnormal flow dynamics, such as vortex blood flow pattern in the pulmonary artery (PA), may reflect progression of pulmonary arterial hypertension (PAH). Some reports suggested that abnormal blood flow parameters including wall shear stress (WSS) could be markers of PAH. However, it was not fully assessed clinical usefulness of these variables. We aimed to assess whether these flow dynamic parameters, such as vortex formation time (VFT) and WSS, were associated with right ventricular (RV) function.

Results

Fifteen subjects, nine with PAH and six healthy volunteers, underwent 4D-flow. Differences of Blood flow patterns, blood flow velocities and WSS between PAH patients and healthy volunteers were evaluated. We also assessed the association between VFT, WSS and RV function in PAH patients. Both vortex blood flow patterns and early systolic retrograde flow in the main PA were observed in all patients with PAH. The PA flow velocities and WSS in patients with PAH were lower than those in healthy volunteers, but that blood flow volumes in the MPA, RPA and LPA and SV in the MPA were broadly comparable between the groups. The mean VFT was 35.0 ± 16.6 % of the cardiac cycle. The VFT significantly correlated with RV ejection fraction, RV end systolic volume, and RV end systolic volume index (RVEF = 75.1 + (−85.7)·VFT, p = 0.003, RVESV = 12.4 + 181.8·VFT, p = 0.037 and RVESVI = 10.6 + 114.8·VFT, p = 0.038, respectively) in PAH patients, whereas WSS did not correlate with RV function.

Conclusions

We confirmed that abnormal blood flow dynamics, including the vortex formation and the early onset of retrograde flow, low WSS in the PA were characteristics of PAH. The VFT may be associated with right ventricular dysfunction, whereas WSS was not. Our results suggest that 4D-flow is an effective means of detecting right heart failure as well as diagnosing PAH.

Clinical trial registration URL: https://upload.umin.ac.jp/cgi-open-bin/ctr/ctr.cgi. Unique identifier: UMIN000011128

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-2755-7) contains supplementary material, which is available to authorized users.

Resolving phase ambiguity in dual-echo dixon imaging using a projected power method

PURPOSE

To develop a fast and robust method to resolve phase ambiguity in dual-echo Dixon imaging.

METHODS

A major challenge in dual-echo Dixon imaging is to estimate the phase error resulting from field inhomogeneity. In this work, a binary quadratic optimization program was formulated to resolve the phase ambiguity. A projected power method was developed to efficiently solve the optimization problem. Both the 1-peak fat model and 6-peak fat model were applied to three-dimensional (3D) datasets. Additionally, the proposed method was extended to dynamic magnetic resonance imaging (MRI) applications using the 6-peak fat model. With institutional review board (IRB) approval and patient consent/assent, the proposed method was evaluated and compared with region growing on 29 consecutive 3D high-resolution patient datasets.

RESULTS

Fast and robust water/fat separation was achieved by the proposed method in different representative 3D datasets and dynamic 3D datasets. Superior water/fat separation was achieved using the 6-peak fat model compared with the 1-peak fat model. Compared to region growing, the proposed method reduced water/fat swaps from 76 to 7% of the patient cohort.

CONCLUSION

The proposed method can achieve fast and robust phase error estimation in dual-echo Dixon imaging. Magn Reson Med 77:2066-2076, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Assessment of the precision and reproducibility of ventricular volume, function, and mass measurements with ferumoxytol-enhanced 4D flow MRI

PURPOSE

To compare the precision and interobserver agreement of ventricular volume, function, and mass quantification by 3D time-resolved (4D) flow MRI relative to cine steady-state free precession (SSFP).

MATERIALS AND METHODS

With Institutional Research Board approval, informed consent, and HIPAA compliance, 22 consecutive patients with congenital heart disease (CHD) (10 males, 6.4 ± 4.8 years) referred for 3T ferumoxytol-enhanced cardiac MRI were prospectively recruited. Complete ventricular coverage with standard 2D short-axis cine SSFP and whole chest coverage with axial 4D flow were obtained. Two blinded radiologists independently segmented images for left ventricular (LV) and right ventricular (RV) myocardium at end systole (ES) and end diastole (ED). Statistical analysis included linear regression, analysis of variance (ANOVA), Bland-Altman (BA) analysis, and intraclass correlation (ICC).

RESULTS

Significant positive correlations were found between 4D flow and SSFP for ventricular volumes (r = 0.808-0.972, P < 0.001), ejection fraction (EF) (r = 0.900-928, P < 0.001), and mass (r = 0.884-0.934, P < 0.001). BA relative limits of agreement for both ventricles were between -52% to 34% for volumes, -29% to 27% for EF, and -41% to 48% for mass, with wider limits of agreement for the RV compared to the LV. There was no significant difference between techniques with respect to mean square difference of ED-ES mass for either LV (F = 2.05, P = 0.159) or RV (F = 0.625, P = 0.434). Interobserver agreement was moderate to good with both 4D flow (ICC 0.523-0.993) and SSFP (ICC 0.619-0.982), with overlapping confidence intervals.

CONCLUSION

Quantification of ventricular volume, function, and mass can be accomplished with 4D flow MRI with precision and interobserver agreement comparable to that of cine SSFP. J. Magn. Reson. Imaging 2016;44:383-392.

High-resolution diffusion-weighted imaging of the breast with multiband 2D radiofrequency pulses and a generalized parallel imaging reconstruction

PURPOSE

To develop a technique for high-resolution diffusion-weighted imaging (DWI) and to compare it with standard DWI methods.

METHODS

Multiple in-plane bands of magnetization were simultaneously excited by identically phase modulating each subpulse of a two-dimensional (2D) RF pulse. Several excitations with the same multiband pattern progressively shifted in the phase-encode direction were used to cover the prescribed field of view (FOV). The phase-encoded FOV was limited to the width of a single band to reduce off-resonance-induced distortion and blurring. Parallel imaging (PI) techniques were used to resolve aliasing from the other bands and to combine the different excitations. Following validation in phantoms and healthy volunteers, a preliminary study in breast cancer patients (N=14) was performed to compare the proposed method to conventional DWI with PI and to reduced-FOV DWI.

RESULTS

The proposed method gave high-resolution diffusion-weighted images with minimal artifacts at the band intersections. Compared to PI alone, higher phase-encoded FOV-reduction factors and reduced noise amplification were obtained, which translated to higher resolution images than conventional (non-multiband) DWI. The same resolution and image quality achievable over targeted regions using existing reduced-FOV methods was obtained, but the proposed method also enables complete bilateral coverage.

CONCLUSION

We developed an in-plane multiband technique for high-resolution DWI and compared its performance with other standard DWI methods. Magn Reson Med 77:209-220, 2017. © 2016 Wiley Periodicals, Inc.

Comprehensive motion-compensated highly accelerated 4D flow MRI with ferumoxytol enhancement for pediatric congenital heart disease

PURPOSE

To develop and evaluate motion-compensation and compressed-sensing techniques in 4D flow MRI for anatomical assessment in a comprehensive ferumoxytol-enhanced congenital heart disease (CHD) exam.

MATERIALS AND METHODS

A Cartesian 4D flow sequence was developed to enable intrinsic navigation and two variable-density sampling schemes: VDPoisson and VDRad. Four compressed-sensing methods were developed: A) VDPoisson scan reconstructed using spatial wavelets; B) added temporal total variation to A; C) VDRad scan using the same reconstruction as in B; and D) added motion compensation to C. With Institutional Review Board (IRB) approval and Health Insurance Portability and Accountability Act (HIPAA) compliance, 23 consecutive patients (eight females, mean 6.3 years) referred for ferumoxytol-enhanced CHD 3T MRI were recruited. Images were acquired and reconstructed using methods A-D. Two cardiovascular radiologists independently scored the images on a 5-point scale. These readers performed a paired wall motion and functional assessment between method D and 2D balanced steady-state free precession (bSSFP) CINE for 16 cases.

RESULTS

Method D had higher diagnostic image quality for most anatomical features (mean 3.8-4.8) compared to A (2.0-3.6), B (2.2-3.7), and C (2.9-3.9) with P < 0.05 with good interobserver agreement (κ ≥ 0.49). Method D had similar or better assessment of myocardial borders and cardiac motion compared to 2D bSSFP (P < 0.05, κ ≥ 0.77). All methods had good internal agreement in comparing aortic with pulmonic flow (BA mean < 0.02%, r > 0.85) and compared to method A (BA mean < 0.13%, r > 0.84) with P < 0.01.

CONCLUSION

Flow, functional, and anatomical assessment in CHD with ferumoxytol-enhanced 4D flow is feasible and can be significantly improved using motion compensation and compressed sensing. J. Magn. Reson. Imaging 2016;43:1355-1368.

Improved quantification and mapping of anomalous pulmonary venous flow with four-dimensional phase-contrast MRI and interactive streamline rendering

BACKGROUND

Cardiac MRI is routinely performed for quantification of shunt flow in patients with anomalous pulmonary veins, but can be technically-challenging to perform. Four-dimensional phase-contrast (4D-PC) MRI has potential to simplify this exam. We sought to determine whether 4D-PC may be a viable clinical alternative to conventional 2D phase-contrast MR imaging.

METHODS

With institutional review board approval and HIPAA-compliance, we retrospectively identified all patients with anomalous pulmonary veins who underwent cardiac MRI at either 1.5 Tesla (T) or 3T with parallel-imaging compressed-sensing (PI-CS) 4D-PC between April, 2011 and October, 2013. A total of 15 exams were included (10 male, 5 female). Algorithms for interactive streamline visualization were developed and integrated into in-house software. Blood flow was measured at the valves, pulmonary arteries and veins, cavae, and any associated shunts. Pulmonary veins were mapped to their receiving atrial chamber with streamlines. The intraobserver, interobserver, internal consistency of flow measurements, and consistency with conventional MRI were then evaluated with Pearson correlation and Bland-Altman analysis.

RESULTS

Triplicate measurements of blood flow from 4D-PC were highly consistent, particularly at the aortic and pulmonary valves (cv 2-3%). Flow measurements were reproducible by a second observer (ρ = 0.986-0.999). Direct measurements of shunt volume from anomalous veins and intracardiac shunts matched indirect estimates from the outflow valves (ρ = 0.966). Measurements of shunt fraction using 4D-PC using any approach were more consistent with ventricular volumetric displacements than conventional 2D-PC (ρ = 0.972-0.991 versus 0.929).

CONCLUSION

Shunt flow may be reliably quantified with 4D-PC MRI, either indirectly or with detailed delineation of flow from multiple shunts. The 4D-PC may be a more accurate alternative to conventional MRI.

Congenital heart disease assessment with 4D flow MRI

With improvements in surgical and medical management, patients with congenital heart disease (CHD) are often living well into adulthood. MRI provides critical data for diagnosis and monitoring of these patients, yielding information on cardiac anatomy, blood flow, and cardiac function. Though historically these exams have been complex and lengthy, four-dimensional (4D) flow is emerging as a single fast technique for comprehensive assessment of CHD. The 4D flow consists of a volumetric time-resolved acquisition that is gated to the cardiac cycle, providing a time-varying vector field of blood flow as well as registered anatomic images. In this article, we provide an overview of MRI evaluation of congenital heart disease by means of example of three relatively common representative conditions: tetralogy of Fallot, aortic coarctation, and anomalous pulmonary venous drainage. Then 4D flow data acquisition, data correction, and postprocessing techniques are reviewed. We conclude with several examples that highlight the comprehensive nature of the evaluation of congenital heart disease with 4D flow.

Fast pediatric 3D free-breathing abdominal dynamic contrast enhanced MRI with high spatiotemporal resolution

PURPOSE

To develop a method for fast pediatric 3D free-breathing abdominal dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI) and investigate its clinical feasibility.

MATERIALS AND METHODS

A combined locally low rank parallel imaging method with soft gating is proposed for free-breathing DCE MRI acquisition. With Institutional Review Board (IRB) approval and informed consent/assent, 23 consecutive pediatric patients were recruited for this study. Free-breathing DCE MRI with ∼1 mm(3) spatial resolution and a 6.5-sec frame rate was acquired on a 3T scanner. Undersampled data were reconstructed with a compressed sensing method without motion correction (FB-CS) and the proposed method (FB-LR). A follow-up respiratory-triggered acquisition (RT-CS) was performed as a reference standard. The reconstructed images were evaluated independently by two radiologists. Wilcoxon tests were performed to test the hypothesis that there was no significant difference between different reconstructions. Quantitative evaluation of contrast dynamics was also performed.

RESULTS

The mean score of overall image quality of FB-LR was 4.0 on a 5-point scale, significantly better (P < 0.05) than FB-CS reconstruction (mean score 2.9), and similar to RT-CS (mean score 4.1). FB-LR also matched the temporal fidelity of contrast dynamics with a root mean square error less than 5%.

CONCLUSION

Fast 3D free-breathing DCE MRI with high scan efficiency and image quality similar to respiratory-triggered acquisition is feasible in a pediatric clinical setting.

Free-breathing pediatric MRI with nonrigid motion correction and acceleration

PURPOSE

To develop and assess motion correction techniques for high-resolution pediatric abdominal volumetric magnetic resonance images acquired free-breathing with high scan efficiency.

MATERIALS AND METHODS

First, variable-density sampling and radial-like phase-encode ordering were incorporated into the 3D Cartesian acquisition. Second, intrinsic multichannel butterfly navigators were used to measure respiratory motion. Lastly, these estimates are applied for both motion-weighted data-consistency in a compressed sensing and parallel imaging reconstruction, and for nonrigid motion correction using a localized autofocusing framework. With Institutional Review Board approval and informed consent/assent, studies were performed on 22 consecutive pediatric patients. Two radiologists independently scored the images for overall image quality, degree of motion artifacts, and sharpness of hepatic vessels and the diaphragm. The results were assessed using paired Wilcoxon test and weighted kappa coefficient for interobserver agreements.

RESULTS

The complete procedure yielded significantly better overall image quality (mean score of 4.7 out of 5) when compared to using no correction (mean score of 3.4, P < 0.05) and to using motion-weighted accelerated imaging (mean score of 3.9, P < 0.05). With an average scan time of 28 seconds, the proposed method resulted in comparable image quality to conventional prospective respiratory-triggered acquisitions with an average scan time of 91 seconds (mean score of 4.5).

CONCLUSION

With the proposed methods, diagnosable high-resolution abdominal volumetric scans can be obtained from free-breathing data acquisitions.

Robust 4D flow denoising using divergence-free wavelet transform

PURPOSE

To investigate four-dimensional flow denoising using the divergence-free wavelet (DFW) transform and compare its performance with existing techniques.

THEORY AND METHODS

DFW is a vector-wavelet that provides a sparse representation of flow in a generally divergence-free field and can be used to enforce "soft" divergence-free conditions when discretization and partial voluming result in numerical nondivergence-free components. Efficient denoising is achieved by appropriate shrinkage of divergence-free wavelet and nondivergence-free coefficients. SureShrink and cycle spinning are investigated to further improve denoising performance.

RESULTS

DFW denoising was compared with existing methods on simulated and phantom data and was shown to yield better noise reduction overall while being robust to segmentation errors. The processing was applied to in vivo data and was demonstrated to improve visualization while preserving quantifications of flow data.

CONCLUSION

DFW denoising of four-dimensional flow data was shown to reduce noise levels in flow data both quantitatively and visually.

Inlet and outlet valve flow and regurgitant volume may be directly and reliably quantified with accelerated, volumetric phase-contrast MRI

PURPOSE

To determine whether it is feasible to use solely an accelerated 4D phase-contrast magnetic resonance imaging (4D-PC MRI) acquisition to quantify net and regurgitant flow volume through each of the cardiac valves.

MATERIALS AND METHODS

Accelerated, 4D-PC MRI examinations performed between March 2010 through June 2011 as part of routine MRI examinations for congenital, structural heart disease were retrospectively reviewed and analyzed using valve-tracking visualization and quantification algorithms developed in Java and OpenGL. Excluding patients with transposition or single ventricle physiology, a total of 34 consecutive pediatric patients (19 male, 15 female; mean age 6.9 years; age range 10 months to 15 years) were identified. 4D-PC flow measurements were compared at each valve and against routine measurements from conventional cardiac MRI using Bland-Altman and Pearson correlation analysis.

RESULTS

Inlet and outlet valve net flow were highly correlated between all valves (P = 0.940-0.985). The sum of forward flow at the outlet valve and regurgitant flow at the inlet valve were consistent with volumetric displacements in each ventricle (P = 0.939-0.948). These were also highly consistent with conventional planar MRI measurements with net flow (P = 0.923-0.935) and regurgitant fractions (P = 0.917-0.972) at the outlet valve and ventricular volumes (P = 0.925-0.965).

CONCLUSION

It is possible to obtain consistent measurements of net and regurgitant blood flow across the inlet and outlet valves relying solely on accelerated 4D-PC. This may facilitate more efficient clinical quantification of valvular regurgitation.

High-resolution, three-dimensional diffusion-weighted breast imaging using DESS

PURPOSE

To evaluate the use of the double-echo steady-state (DESS) sequence for acquiring high-resolution breast images with diffusion and T2 weighting.

MATERIALS AND METHODS

Phantom scans were used to verify the T2 and diffusion weighting of the DESS sequence. Image distortion was evaluated in volunteers by comparing DESS images and conventional diffusion-weighted images (DWI) to spoiled gradient-echo images. The DESS sequence was added to a standard clinical protocol, and the resulting patient images were used to evaluate overall image quality and image contrast in lesions.

RESULTS

The diffusion weighting of the DESS sequence can be easily modulated by changing the spoiler gradient area and flip angle. Radiologists rated DESS images as having higher resolution and less distortion than conventional DWI. Lesion-to-tissue contrast ratios are strongly correlated between DWI and DESS images (R=0.83) and between T2-weighted fast spin-echo and DESS images (R=0.80).

CONCLUSION

The DESS sequence is able to acquire high-resolution 3D diffusion- and T2-weighted images in short scan times, with image quality that facilitates morphological assessment of lesions.

Clinical performance of contrast enhanced abdominal pediatric MRI with fast combined parallel imaging compressed sensing reconstruction

PURPOSE

To deploy clinically, a combined parallel imaging compressed sensing method with coil compression that achieves a rapid image reconstruction, and assess its clinical performance in contrast-enhanced abdominal pediatric MRI.

MATERIALS AND METHODS

With Institutional Review Board approval and informed patient consent/assent, 29 consecutive pediatric patients were recruited. Dynamic contrast-enhanced MRI was acquired on a 3 Tesla scanner using a dedicated 32-channel pediatric coil and a three-dimensional SPGR sequence, with pseudo-random undersampling at a high acceleration (R = 7.2). Undersampled data were reconstructed with three methods: a traditional parallel imaging method and a combined parallel imaging compressed sensing method with and without coil compression. The three sets of images were evaluated independently and blindly by two radiologists at one siting, for overall image quality and delineation of anatomical structures. Wilcoxon tests were performed to test the hypothesis that there was no significant difference in the evaluations, and interobserver agreement was analyzed.

RESULTS

Fast reconstruction with coil compression did not deteriorate image quality. The mean score of structural delineation of the fast reconstruction was 4.1 on a 5-point scale, significantly better (P < 0.05) than traditional parallel imaging (mean score 3.1). Fair to substantial interobserver agreement was reached in structural delineation assessment.

CONCLUSION

A fast combined parallel imaging compressed sensing method is feasible in a pediatric clinical setting. Preliminary results suggest it may improve structural delineation over parallel imaging.

Evaluation of valvular insufficiency and shunts with parallel-imaging compressed-sensing 4D phase-contrast MR imaging with stereoscopic 3D velocity-fusion volume-rendered visualization

PURPOSE

To assess the potential of compressed-sensing parallel-imaging four-dimensional (4D) phase-contrast magnetic resonance (MR) imaging and specialized imaging software in the evaluation of valvular insufficiency and intracardiac shunts in patients with congenital heart disease.

MATERIALS AND METHODS

Institutional review board approval was obtained for this HIPAA-compliant study. Thirty-four consecutive retrospectively identified patients in whom a compressed-sensing parallel-imaging 4D phase-contrast sequence was performed as part of routine clinical cardiac MR imaging between March 2010 and August 2011 and who had undergone echocardiography were included. Multiplanar, volume-rendered, and stereoscopic three-dimensional velocity-fusion visualization algorithms were developed and implemented in Java and OpenGL. Two radiologists independently reviewed 4D phase-contrast studies for each of 34 patients (mean age, 6 years; age range, 10 months to 21 years) and tabulated visible shunts and valvular regurgitation. These results were compared with color Doppler echocardiographic and cardiac MR imaging reports, which were generated without 4D phase-contrast visualization. Cohen κ statistics were computed to assess interobserver agreement and agreement with echocardiographic results.

RESULTS

The 4D phase-contrast acquisitions were performed, on average, in less than 10 minutes. Among 123 valves seen in 34 4D phase-contrast studies, 29 regurgitant valves were identified, with good agreement between observers (k=0.85). There was also good agreement with the presence of at least mild regurgitation at echocardiography (observer 1, κ=0.76; observer 2, κ=0.77) with high sensitivity (observer 1, 75%; observer 2, 82%) and specificity (observer 1, 97%; observer 2, 95%) relative to the reference standard. Eight intracardiac shunts were identified, four of which were not visible with conventional cardiac MR imaging but were detected with echocardiography. No intracardiac shunts were found with echocardiography alone.

CONCLUSION

With velocity-fusion visualization, the compressed-sensing parallel-imaging 4D phase-contrast sequence can augment conventional cardiac MR imaging by improving sensitivity for and depiction of hemodynamically significant shunts and valvular regurgitation.

Nonrigid motion correction in 3D using autofocusing with localized linear translations

MR scans are sensitive to motion effects due to the scan duration. To properly suppress artifacts from nonrigid body motion, complex models with elements such as translation, rotation, shear, and scaling have been incorporated into the reconstruction pipeline. However, these techniques are computationally intensive and difficult to implement for online reconstruction. On a sufficiently small spatial scale, the different types of motion can be well approximated as simple linear translations. This formulation allows for a practical autofocusing algorithm that locally minimizes a given motion metric--more specifically, the proposed localized gradient-entropy metric. To reduce the vast search space for an optimal solution, possible motion paths are limited to the motion measured from multichannel navigator data. The novel navigation strategy is based on the so-called "Butterfly" navigators, which are modifications of the spin-warp sequence that provides intrinsic translational motion information with negligible overhead. With a 32-channel abdominal coil, sufficient number of motion measurements were found to approximate possible linear motion paths for every image voxel. The correction scheme was applied to free-breathing abdominal patient studies. In these scans, a reduction in artifacts from complex, nonrigid motion was observed.

Simultaneous estimation of T(2) and apparent diffusion coefficient in human articular cartilage in vivo with a modified three-dimensional double echo steady state (DESS) sequence at 3 T

T(2) mapping and diffusion-weighted imaging complement morphological imaging for assessing cartilage disease and injury. The double echo steady state sequence has been used for morphological imaging and generates two echoes with markedly different T(2) and diffusion weighting. Modifying the spoiler gradient area and flip angle of the double echo steady state sequence allows greater control of the diffusion weighting of both echoes. Data from two acquisitions with different spoiler gradient areas and flip angles are used to simultaneously estimate the T(2) and apparent diffusion coefficient of each voxel. This method is verified in phantoms and validated in vivo in the knee; estimates from different regions of interest in the phantoms and cartilage are compared to those obtained using standard spin-echo methods. The Pearson correlations were 0.984 for T(2) (∼2% relative difference between spin-echo and double echo steady state estimates) and 0.997 for apparent diffusion coefficient (˜1% relative difference between spin-echo and double echo steady state estimates) for the phantom study and 0.989 for T(2) and 0.987 for apparent diffusion coefficient in regions of interest in the human knee in vivo. High accuracy for simultaneous three-dimensional T(2) and apparent diffusion coefficient measurements are demonstrated, while also providing morphologic three-dimensional images without blurring or distortion in reasonable scan times.

Combined respiratory and cardiac triggering improves blood pool contrast-enhanced pediatric cardiovascular MRI

BACKGROUND

Contrast-enhanced cardiac MRA suffers from cardiac motion artifacts and often requires a breath-hold.

OBJECTIVE

This work develops and evaluates a blood pool contrast-enhanced combined respiratory- and ECG-triggered MRA method.

MATERIALS AND METHODS

An SPGR sequence was modified to enable combined cardiac and respiratory triggering on a 1.5-T scanner. Twenty-three consecutive children referred for pediatric heart disease receiving gadofosveset were recruited in HIPAA-compliant fashion with IRB approval and informed consent. Children underwent standard non-triggered contrast-enhanced MRA with or without suspended respiration. Additionally, a free-breathing-triggered MRA was acquired. Triggered and non-triggered studies were presented in blinded random order independently to two radiologists twice. Anatomical structure delineation was graded for each triggered and non-triggered acquisition and the visual quality on triggered MRA was compared directly to that on non-triggered MRA.

RESULTS

Triggered images received higher scores from each radiologist for all anatomical structures on each of the two reading sessions (Wilcoxon rank sum test, P < 0.05). In direct comparison, triggered images were preferred over non-triggered images for delineating cardiac structures, with most comparisons reaching statistical significance (binomial test, P < 0.05).

CONCLUSION

Combined cardiac and respiratory triggering, enabled by a blood pool contrast agent, improves delineation of most anatomical structures in pediatric cardiovascular MRA.

Improved cardiovascular flow quantification with time-resolved volumetric phase-contrast MRI

BACKGROUND

Cardiovascular flow is commonly assessed with two-dimensional, phase-contrast MRI (2-D PC-MRI). However, scan prescription and acquisition over multiple planes is lengthy, often requires direct physician oversight and has inconsistent results. Time-resolved volumetric PC-MRI (4-D flow) may address these limitations.

OBJECTIVE

We assess the degree of agreement and internal consistency between 2-D and 4-D flow quantification in our clinical population.

MATERIALS AND METHODS

Software enabling flow calculation from 4-D flow was developed in Java. With IRB approval and HIPAA compliance, 18 consecutive patients without shunts were identified who underwent both (1) conventional 2-D PC-MRI of the aorta and main pulmonary artery and (2) 4-D flow imaging. Aortic and pulmonary flow rates were assessed with both techniques.

RESULTS

Both methods showed general agreement in flow rates (ρ: 0.87-0.90). Systemic and pulmonary arterial flow rates were well-correlated (ρ: 4-D 0.98-0.99, 2-D 0.93), but more closely matched with 4-D (P < 0.05, Brown-Forsythe). Pulmonary flow rates were lower than systemic rates for 2-D (P < 0.05, two-sample t-test). In a sub-analysis of patients without pulmonary or aortic regurgitation, 2-D showed improved correlation of flow rates while 4-D phase-contrast remained tightly correlated (ρ: 4-D 0.99-1.00, 2-D 0.99).

CONCLUSION

4-D PC-MRI demonstrates greater consistency than conventional 2-D PC-MRI for flow quantification.

Custom-fitted 16-channel bilateral breast coil for bidirectional parallel imaging

A 16-channel receive-only, closely fitted array coil is described and tested in vivo for bilateral breast imaging at 3 T. The primary purpose of this coil is to provide high signal-to-noise ratio and parallel imaging acceleration in two directions for breast MRI. Circular coil elements (7.5-cm diameter) were placed on a closed "cup-shaped" platform, and nearest neighbor coils were decoupled through geometric overlap. Comparisons were made between the 16-channel custom coil and a commercially available 8-channel coil. SENSitivity Encoding (SENSE) parallel imaging noise amplification (g-factor) was evaluated in phantom scans. In healthy volunteers, we compared signal-to-noise ratio, parallel imaging in one and two directions, Autocalibrating Reconstruction for Cartesian sampling (ARC) g-factor, and high spatial resolution imaging. When compared with a commercially available 8-channel coil, the 16-channel custom coil shows 3.6× higher mean signal-to-noise ratio in the breast and higher quality accelerated images. In patients, the 16-channel custom coil has facilitated high-quality, high-resolution images with bidirectional acceleration of R = 6.3.

Accelerated slice encoding for metal artifact correction

PURPOSE

To demonstrate accelerated imaging with both artifact reduction and different contrast mechanisms near metallic implants.

MATERIALS AND METHODS

Slice-encoding for metal artifact correction (SEMAC) is a modified spin echo sequence that uses view-angle tilting and slice-direction phase encoding to correct both in-plane and through-plane artifacts. Standard spin echo trains and short-TI inversion recovery (STIR) allow efficient PD-weighted imaging with optional fat suppression. A completely linear reconstruction allows incorporation of parallel imaging and partial Fourier imaging. The signal-to-noise ratio (SNR) effects of all reconstructions were quantified in one subject. Ten subjects with different metallic implants were scanned using SEMAC protocols, all with scan times below 11 minutes, as well as with standard spin echo methods.

RESULTS

The SNR using standard acceleration techniques is unaffected by the linear SEMAC reconstruction. In all cases with implants, accelerated SEMAC significantly reduced artifacts compared with standard imaging techniques, with no additional artifacts from acceleration techniques. The use of different contrast mechanisms allowed differentiation of fluid from other structures in several subjects.

CONCLUSION

SEMAC imaging can be combined with standard echo-train imaging, parallel imaging, partial-Fourier imaging, and inversion recovery techniques to offer flexible image contrast with a dramatic reduction of metal-induced artifacts in scan times under 11 minutes.

Breathheld autocalibrated phase-contrast imaging

PURPOSE

To compare generalized autocalibrating partially parallel acquisitions (GRAPPA), modified sensitivity encoding (mSENSE), and SENSE in phase-contrast magnetic resonance imaging (PC-MRI) applications.

MATERIALS AND METHODS

Aliasing of the torso can occur in PC-MRI applications. If the data are further undersampled for parallel imaging, SENSE can be problematic in correctly unaliasing signals due to coil sensitivity maps that do not match that of the aliased volume. Here, a method for estimating coil sensitivities in flow applications is described. Normal volunteers (n = 5) were scanned on a 1.5 T MRI scanner and underwent PC-MRI scans using GRAPPA, mSENSE, SENSE, and conventional PC-MRI acquisitions. Peak velocity and flow through the aorta and pulmonary artery were evaluated.

RESULTS

Bland-Altman statistics for flow in the aorta and pulmonary artery acquired with mSENSE and GRAPPA methods (R = 2 and R = 3 cases) have comparable mean differences to flow acquired with conventional PC-MRI. GRAPPA and mSENSE PC-MRI have more robust measurements than SENSE when there is aliasing artifact caused by insufficient coil sensitivity maps. For peak velocity, there are no considerable differences among the mSENSE, GRAPPA, and SENSE reconstructions and are comparable to conventional PC-MRI.

CONCLUSION

Flow measurements of images reconstructed with autocalibration techniques have comparable agreement with conventional PC-MRI and provide robust measurements in the presence of wraparound.

Three-dimensional T(1), T(2) and proton density mapping with inversion recovery balanced SSFP

By combining a balanced steady-state free precession (bSSFP) readout with an initial inversion pulse, all three contrast parameters, T(1), T(2) and proton density (M(0)), may be rapidly calculated from the signal progression in time. However, here it is shown that this technique is quite sensitive to variation in the applied transmit RF (B(1)) field, leading to pronounced errors in calculated values. Two-dimensional (2D) acquisitions are taxed to accurately quantify the relaxation, as the short RF pulses required by SSFP's rapid TR contain a broad spectrum of excitation angles. A 3D excitation using a large diameter excitation coil was able to correctly quantify the parameters. While the extreme B(1) sensitivity was previously problematic and has precluded use of IR-bSSFP for relaxometry, in this work these obstacles were significantly reduced, allowing the rapid quantification of T(1), T(2) and M(0). The results may further be used to simulate image contrast from common sequences, such as a T(1)-weighted or fluid-attenuated inversion recovery (FLAIR) examination.

Clinical evaluation of aortic coarctation with 4D flow MR imaging

PURPOSE

To show that 4D Flow is a clinically viable tool for evaluation of collateral blood flow and demonstration of distorted blood flow patterns in patients with treated and untreated aortic coarctation.

MATERIALS AND METHODS

Time-resolved, 3D phase contrast magnetic resonance imaging (MRI) (4D Flow) was used to assess blood flow in the thoracic aorta of 34 individuals: 26 patients with coarctation (22 after surgery or stent placement) and eight healthy volunteers.

RESULTS

Direct comparison of blood flow calculated with 2D and 4D phase contrast data at standard levels for analysis in coarctation patients showed good correlation and agreement (correlation coefficient r = 0.99, limits of agreement = -20% to 20% for collateral blood flow calculations). Abnormal blood flow patterns were demonstrated at peak systole with 4D Flow visualization techniques in the descending thoracic aorta of patients but not volunteers. Marked helical flow was seen in 9 of 13 patients with angulated aortic arch geometries after coarctation repair. Vortical flow was seen in regions of poststenotic dilation.

CONCLUSION

4D Flow is a fast and reliable means of evaluating collateral blood flow in patients with aortic coarctation in order to establish hemodynamic significance. It also can detect distorted blood flow patterns in the descending aorta after coarctation repair.

Improved pediatric MR imaging with compressed sensing

PURPOSE

To develop a method that combines parallel imaging and compressed sensing to enable faster and/or higher spatial resolution magnetic resonance (MR) imaging and show its feasibility in a pediatric clinical setting.

MATERIALS AND METHODS

Institutional review board approval was obtained for this HIPAA-compliant study, and informed consent or assent was given by subjects. A pseudorandom k-space undersampling pattern was incorporated into a three-dimensional (3D) gradient-echo sequence; aliasing then has an incoherent noiselike pattern rather than the usual coherent fold-over wrapping pattern. This k-space-sampling pattern was combined with a compressed sensing nonlinear reconstruction method that exploits the assumption of sparsity of medical images to permit reconstruction from undersampled k-space data and remove the noiselike aliasing. Thirty-four patients (15 female and 19 male patients; mean age, 8.1 years; range, 0-17 years) referred for cardiovascular, abdominal, and knee MR imaging were scanned with this 3D gradient-echo sequence at high acceleration factors. Obtained k-space data were reconstructed with both a traditional parallel imaging algorithm and the nonlinear method. Both sets of images were rated for image quality, radiologist preference, and delineation of specific structures by two radiologists. Wilcoxon and symmetry tests were performed to test the hypothesis that there was no significant difference in ratings for image quality, preference, and delineation of specific structures.

RESULTS

Compressed sensing images were preferred more often, had significantly higher image quality ratings, and greater delineation of anatomic structures (P < .001) than did images obtained with the traditional parallel reconstruction method.

CONCLUSION

A combination of parallel imaging and compressed sensing is feasible in a clinical setting and may provide higher resolution and/or faster imaging, addressing the challenge of delineating anatomic structures in pediatric MR imaging.

In vivo hemodynamic analysis of intracranial aneurysms obtained by magnetic resonance fluid dynamics (MRFD) based on time-resolved three-dimensional phase-contrast MRI

INTRODUCTION

Hemodynamics is thought to play a very important role in the initiation, growth, and rupture of intracranial aneurysms. The purpose of our study was to perform in vivo hemodynamic analysis of unruptured intracranial aneurysms of magnetic resonance fluid dynamics using time-resolved three-dimensional phase-contrast MRI (4D-Flow) at 1.5 T and to analyze relationships between hemodynamics and wall shear stress (WSS) and oscillatory shear index (OSI).

METHODS

This study included nine subjects with 14 unruptured aneurysms. 4D-Flow was performed by a 1.5-T magnetic resonance scanner with a head coil. We calculated in vivo streamlines, WSS, and OSI of intracranial aneurysms based on 4D-Flow with our software. We evaluated the number of spiral flows in the aneurysms and compared the differences in WSS or OSI between the vessel and aneurysm and between whole aneurysm and the apex of the spiral flow.

RESULTS

3D streamlines, WSS, and OSI distribution maps in arbitrary direction during the cardiac phase were obtained for all intracranial aneurysms. Twelve aneurysms had one spiral flow each, and two aneurysms had two spiral flows each. The WSS was lower and the OSI was higher in the aneurysm compared to the vessel. The apex of the spiral flow had a lower WSS and higher OSI relative to the whole aneurysm.

CONCLUSION

Each intracranial aneurysm in this study had at least one spiral flow. The WSS was lower and OSI was higher at the apex of the spiral flow than the whole aneurysmal wall.

Hemodynamics in Normal Cerebral Arteries: Qualitative Comparison of 4D Phase-Contrast Magnetic Resonance and Image-Based Computational Fluid Dynamics

Detailed knowledge of the hemodynamic conditions in normal cerebral arteries is important for a better understanding of the underlying mechanisms leading to the initiation and progression of cerebrovascular diseases. Information about the baseline values of hemodynamic variables such as wall shear stresses is necessary for comparison to pathological conditions such as in cerebral aneurysms or arterial stenoses. The purpose of this study was to compare the blood flow patterns in cerebral arteries of normal subjects determined by 4D phase-contrast magnetic resonance and image-based computational fluid dynamics techniques in order to assess their consistency and to highlight their differences. The goal was not to validate (or disprove) any of the two methodologies but rather to identify regions where disagreements are to be expected and to provide guidance when interpreting the data produced by each technique.

Evaluation of intracranial stenoses and aneurysms with accelerated 4D flow

The aim of this study was to evaluate intracranial arterial stenoses and aneurysms with accelerated time-resolved three-dimensional (3D) phase-contrast MRI or 4D flow. The 4D flow technique was utilized to image four normal volunteers, two patients with intracranial stenoses and two patients with intracranial aneurysms. In order to reduce scan time, parallel imaging was combined with an acquisition strategy that eliminates the corners of k-space. In the two patients with intracranial stenoses, 4D flow velocity measurements showed that one patient had normal velocity profiles in agreement with a previous magnetic resonance angiogram (MRA), while the second showed increased velocities that indicated a less significant narrowing than suspected on a previous MRA, as confirmed by catheter angiography. This result may have prevented an invasive angiogram. In the two patients with 4-mm intracranial aneurysm, one had a stable helical flow pattern with a large jet, while the other had a temporally unstable flow pattern with a more focal jet possibly indicating that the second aneurysm may have a higher likelihood of rupture. Accelerated 4D flow provides time-resolved 3D velocity data in an 8- to 10-min scan. In the stenosis patients, the addition of 4D flow to a traditional MRA adds the velocity data provided from transcranial Doppler ultrasound (TCD) possibly allowing for more accurate grading of stenoses. In the aneurysm patients, visualization of flow patterns may help to provide prognostic information about future risk of rupture.

Complete intracranial arterial and venous blood flow evaluation with 4D flow MR imaging

SUMMARY

Time-resolved, 3D velocity-encoded MR imaging (4D Flow) allows for the acquisition of dynamic, multidirectional data on blood flow and has recently been used for the evaluation of intracranial arterial flow. Using a 3T system with optimization of both temporal resolution and k-space subsampling with a combination of parallel imaging and cut-corner acquisition, we present the clinical assessment of a patient with an arteriovenous malformation by providing complete intracranial arterial and venous coverage in a reasonable scan time.