Analysis and correction of background velocity offsets in phase-contrast flow measurements using magnetic field monitoring

Beat phenomena of oscillating readouts

PURPOSE

To demonstrate slowly varying, erroneous magnetic field gradients for oscillating readouts due to the mechanically resonant behavior of gradient systems.

METHODS

Projections of a static phantom were acquired using a one-dimensional (1D) EPI sequence with varying EPI frequencies ranging from 1121 to 1580 Hz on clinical 3T systems (30 mT/m, 200 T/m/s). Phase due to static B0 inhomogeneities was eliminated by a complex division of two separate scans with different polarities of the EPI readout. The temporal evolution of phase was evaluated and related to the mechanical resonances of the gradient systems derived from the gradient modulation transfer function. Additionally, the impact of temporally varying mechanical resonance effects on EPI was evaluated using an echo-planar spectroscopic imaging sequence.

RESULTS

A beat phenomenon resulting in a slowly varying phase was observed. Its temporal frequency was given by the difference between the EPI frequency and the mechanical resonance frequency of the activated gradient axis. The maximum erroneous, oscillating phase during phase encoding was ±0.5 rad for an EPI frequency of 1281 Hz. Echo-planar spectroscopic imaging images showed the resulting time-dependent stretching/compression of the FOV.

CONCLUSION

Oscillating readouts such as those used in EPI can result in low-frequency, erroneous phase contributions, which are explained by the beat phenomenon. Therefore, EPI phase-correction approaches may need to include beat effects for accurate image reconstruction.

Flow evaluation software for four-dimensional flow MRI: a reliability and validation study

PURPOSE

Four-dimensional time-resolved phase-contrast cardiovascular magnetic resonance imaging (4D flow MRI) enables blood flow quantification in multiple vessels, which is crucial for patients with congenital heart disease (CHD). We investigated net flow volumes in the ascending aorta and pulmonary arteries by four different postprocessing software packages for 4D flow MRI in comparison with 2D cine phase-contrast measurements (2D PC).

MATERIAL AND METHODS

4D flow and 2D PC datasets of 47 patients with biventricular CHD (median age 16, range 0.6-52 years) were acquired at 1.5 T. Net flow volumes in the ascending aorta, the main, right, and left pulmonary arteries were measured using four different postprocessing software applications and compared to offset-corrected 2D PC data. Reliability of 4D flow postprocessing software was assessed by Bland-Altman analysis and intraclass correlation coefficient (ICC). Linear regression of internal flow controls was calculated. Interobserver reproducibility was evaluated in 25 patients.

RESULTS

Correlation and agreement of flow volumes were very good for all software compared to 2D PC (ICC ≥ 0.94; bias ≤ 5%). Internal controls were excellent for 2D PC (r ≥ 0.95, p < 0.001) and 4D flow (r ≥ 0.94, p < 0.001) without significant difference of correlation coefficients between methods. Interobserver reliability was good for all vendors (ICC ≥ 0.94, agreement bias < 8%).

CONCLUSION

Haemodynamic information from 4D flow in the large thoracic arteries assessed by four commercially available postprocessing applications matches routinely performed 2D PC values. Therefore, we consider 4D flow MRI-derived data ready for clinical use in patients with CHD.

Deep learning phase error correction for cerebrovascular 4D flow MRI

Background phase errors in 4D Flow MRI may negatively impact blood flow quantification. In this study, we assessed their impact on cerebrovascular flow volume measurements, evaluated the benefit of manual image-based correction, and assessed the potential of a convolutional neural network (CNN), a form of deep learning, to directly infer the correction vector field. With IRB waiver of informed consent, we retrospectively identified 96 MRI exams from 48 patients who underwent cerebrovascular 4D Flow MRI from October 2015 to 2020. Flow measurements of the anterior, posterior, and venous circulation were performed to assess inflow-outflow error and the benefit of manual image-based phase error correction. A CNN was then trained to directly infer the phase-error correction field, without segmentation, from 4D Flow volumes to automate correction, reserving from 23 exams for testing. Statistical analyses included Spearman correlation, Bland-Altman, Wilcoxon-signed rank (WSR) and F-tests. Prior to correction, there was strong correlation between inflow and outflow (ρ = 0.833-0.947) measurements with the largest discrepancy in the venous circulation. Manual phase error correction improved inflow-outflow correlation (ρ = 0.945-0.981) and decreased variance (p < 0.001, F-test). Fully automated CNN correction was non-inferior to manual correction with no significant differences in correlation (ρ = 0.971 vs ρ = 0.982) or bias (p = 0.82, Wilcoxon-Signed Rank test) of inflow and outflow measurements. Residual background phase error can impair inflow-outflow consistency of cerebrovascular flow volume measurements. A CNN can be used to directly infer the phase-error vector field to fully automate phase error correction.

Fully automated background phase correction using M-estimate SAmple consensus (MSAC)-Application to 2D and 4D flow

PURPOSE

Flow quantification by phase-contrast MRI is hampered by spatially varying background phase offsets. Correction performance by polynomial regression on stationary tissue may be affected by outliers such as wrap-around or constant flow. Therefore, we propose an alternative, M-estimate SAmple Consensus (MSAC) to reject outliers, and improve and fully automate background phase correction.

METHODS

The MSAC technique fits polynomials to randomly drawn small samples from the image. Over several trials, it aims to find the best consensus set of valid pixels by rejecting outliers to the fit and minimizing the residuals of the remaining pixels. The robustness of MSAC to its few parameters was investigated and verified using third-order polynomial correction fits on a total of 118 2D flow (97 with wrap-around) and 18 4D flow data sets (14 with wrap-around), acquired at 1.5 T and 3 T. Background phase was compared with standard stationary correction and phantom correction. Pulmonary/systemic flow ratios in 2D flow were derived, and exemplary 4D flow analysis was performed.

RESULTS

The MSAC technique is robust over a range of parameter choices, and a unique set of parameters is suitable for both 2D and 4D flow. In 2D flow, phase errors were significantly reduced by MSAC compared with stationary correction (p = 0.005), and stationary correction shows larger errors in pulmonary/systemic flow ratios compared with MSAC. In 4D flow, MSAC shows similar performance as stationary correction.

CONCLUSIONS

The MSAC method provides fully automated background phase correction to 2D and 4D flow data and shows improved robustness over stationary correction, especially with outliers present.

Integration of an RF coil and commercial field camera for ultrahigh-field MRI

PURPOSE

To develop an RF coil with an integrated commercial field camera for ultrahigh field (7T) neuroimaging. The RF coil would operate within a head-only gradient coil and be subject to the corresponding design constraints. The RF coil can thereafter be used for subject-specific correction of k-space trajectories-notably in gradient-sensitive sequences such as single-shot spiral imaging.

METHODS

The transmit and receive performance was evaluated before and after the integration of field probes, whereas field probes were evaluated when in an optimal configuration external to the coil and after their integration. Diffusion-weighted EPI and single-shot spiral acquisitions were employed to evaluate the efficacy of correcting higher order field perturbations and the consequent effect on image quality.

RESULTS

Field probes had a negligible effect on RF-coil performance, including the transmit efficiency, transmit uniformity, and mean SNR over the brain. Modest reductions in field-probe signal lifetimes were observed, caused primarily by nonidealities in the gradient and shim fields of the head-only gradient coil at the probe positions. The field-monitoring system could correct up to second-order field perturbations in single-shot spiral imaging.

CONCLUSION

The integrated RF coil and field camera was capable of concurrent-field monitoring within a 7T head-only scanner and facilitated the subsequent correction of k-space trajectories during spiral imaging.

A subject-specific assessment of measurement errors and their correction in cerebrospinal fluid velocity maps using 4D flow MRI

PURPOSE

Phase-contrast MRI (PC-MRI) of cerebrospinal fluid (CSF) velocity is used to evaluate the characteristics of intracranial diseases, such as normal-pressure hydrocephalus (NPH). Nevertheless, PC-MRI has several potential error sources, with eddy-current-based phase offset error being non-negligible in CSF measurement. In this study, we assess the measurement error of CSF velocity maps obtained using 4D flow MRI and evaluate correction methods.

METHODS

CSF velocity maps of 10 patients with NPH were acquired using 4D flow MRI (velocity-encoding = 5 cm/s). Distributed phase offset error was estimated for a whole 3D background field by polynomial fitting using robust regression analysis. This estimated phase offset error was then used to correct the CSF velocity maps. The estimated error profiles were compared with those obtained using an existing 2D correction approach involving local background information near the region of interest.

RESULTS

The residual standard error of the polynomial fitting against the phase offset error extracted from the measured velocities was within 0.2 cm/s. The spatial dependencies of the phase offset errors showed similar tendencies in all cases, but sufficient differences in these values were found to indicate requirement of velocity correction. Differences of the estimated errors among other correction approaches were in the order of 10^-2 cm/s, and the estimated errors were in good agreement with those obtained using existing approaches.

CONCLUSION

Our method is capable of estimating the measurement error of CSF velocity maps obtained from 4D flow MRI and provides quantitatively reasonable characteristics for the main CSF profile in the cerebral aqueduct in patients with NPH.

The clinical impact of phase offset errors and different correction methods in cardiovascular magnetic resonance phase contrast imaging: a multi-scanner study

BACKGROUND

Cardiovascular magnetic resonance (CMR) phase contrast (PC) flow measurements suffer from phase offset errors. Background subtraction based on stationary phantom measurements can most reliably be used to overcome this inaccuracy. Stationary tissue correction is an alternative and does not require additional phantom scanning. The aim of this study was 1) to compare measurements with and without stationary tissue correction to phantom corrected measurements on different GE Healthcare CMR scanners using different software packages and 2) to evaluate the clinical implications of these methods.

METHODS

CMR PC imaging of both the aortic and pulmonary artery flow was performed in patients on three different 1.5 T CMR scanners (GE Healthcare) using identical scan parameters. Uncorrected, first, second and third order stationary tissue corrected flow measurement were compared to phantom corrected flow measurements, our reference method, using Medis QFlow, Circle cvi42 and MASS software. The optimal (optimized) stationary tissue order was determined per scanner and software program. Velocity offsets, net flow, clinically significant difference (deviation > 10% net flow), and regurgitation severity were assessed.

RESULTS

Data from 175 patients (28 (17-38) years) were included, of which 84% had congenital heart disease. First, second and third order and optimized stationary tissue correction did not improve the velocity offsets and net flow measurements. Uncorrected measurements resulted in the least clinically significant differences in net flow compared to phantom corrected data. Optimized stationary tissue correction per scanner and software program resulted in net flow differences (> 10%) in 19% (MASS) and 30% (Circle cvi42) of all measurements compared to 18% (MASS) and 23% (Circle cvi42) with no correction. Compared to phantom correction, regurgitation reclassification was the least common using uncorrected data. One CMR scanner performed worse and significant net flow differences of > 10% were present both with and without stationary tissue correction in more than 30% of all measurements.

CONCLUSION

Phase offset errors had a significant impact on net flow quantification, regurgitation assessment and varied greatly between CMR scanners. Background phase correction using stationary tissue correction worsened accuracy compared to no correction on three GE Healthcare CMR scanners. Therefore, careful assessment of phase offset errors at each individual scanner is essential to determine whether routine use of phantom correction is necessary.

TRIAL REGISTRATION

Observational Study.

Fold-Over Oversampling Effects in the Measurements of Cerebral Cerebrospinal Fluid and Blood Flows with 2D Cine Phase-Contrast MRI

This prospective study investigated the effects of fold-over oversampling on phase-offset background errors with 2D-Cine phase contrast (Cine-PC) magnetic resonance imaging (MRI). It was performed on brain MRI and compared to conventional Full-field of view FOV coverage and it was tested with two different velocity encoding (Venc) values. We chose Venc = 100 mm/s to encode cerebrospinal fluid (CSF) flows in the aqueduct and 600 mm/s to encode blood flow in the carotid artery. Cine-PC was carried out on 10 healthy adult volunteers followed simultaneously by an acquisition on static agar-gel phantom to measure the phase-offset background errors. Pixel-wise correction of both the CSF and the blood flows was calculated through 32 points of the cardiac-cycle. We compared the velocity-to-noise ratio, the section area, the absolute and the corrected velocity (peak; mean and minimum), the net flow, and the stroke volume before and after correction. We performed the statistical T-test to compare Full-FOV and fold-over and Bland–Altman plots to analyze their differences. Our results showed that following phase-offset error correction, the blood stroke-volume was significantly higher with Full-FOV compared to fold-over. We observed a significantly higher CSF mean velocity and net flow values in the fold-over option. Compared to Full-FOV, fold-over provides a significantly larger section area and significantly lower peak velocity-offset in the aqueduct. No significant difference between the two coverages was reported before and after phase-offset in blood flow measurements. In conclusion, fold-over oversampling can be chosen as an alternative to increase spatial resolution and accurate cerebral flow quantification in Cine-PC.

Limits to flow detection in phase contrast MRI

Pulsed gradient spin echo (PGSE) complex signal behavior becomes dominated by attenuation rather than oscillation when displacements due to flow are similar or less than diffusive displacements. In this "slow-flow" regime, the optimal displacement encoding parameter q for phase contrast velocimetry depends on the diffusive length scale q s l o w = 1 / l D = 1 / 2 D Δ rather than the velocity encoding parameter v _enc = π/(qΔ). The minimum detectable mean velocity using the difference between the phase at +q _slow and -q _slow is 〈 v m i n 〉 = 1 / SNR D / Δ . These theories are then validated and applied to MRI by performing PGSE echo planar imaging experiments on water flowing through a column with a bulk region and a beadpack region at controlled flow rates. Velocities as slow as 6 μm/s are detected with velocimetry. Theories, MRI experimental protocols, and validation on a controlled phantom help to bridge the gap between porous media NMR and pre-clinical phase contrast and diffusion MRI.

Evaluation of self-calibrated non-linear phase-contrast correction in pediatric and congenital cardiovascular magnetic resonance imaging

BACKGROUND

The need for background error correction in phase-contrast flow analysis has historically posed a challenge in cardiac magnetic resonance (MR) imaging. While previous studies have shown that phantom correction improves flow measurements, it impedes scanner workflow.

OBJECTIVE

To evaluate the efficacy of self-calibrated non-linear phase-contrast correction on flows in pediatric and congenital cardiac MR compared to phantom correction as the standard.

MATERIALS AND METHODS

We retrospectively identified children who had great-vessel phase-contrast and static phantom sequences acquired between January 2015 and June 2015. We applied a novel correction method to each phase-contrast sequence post hoc. Uncorrected, non-linear, and phantom-corrected flows were compared using intraclass correlation. We used paired t-tests to compare how closely non-linear and uncorrected flows approximated phantom-corrected flows. In children without intra- or extracardiac shunts or significant semilunar valvular regurgitation, we used paired t-tests to compare how closely the uncorrected pulmonary-to-systemic flow ratio (Qp:Qs) and non-linear Qp:Qs approximated phantom-corrected Qp:Qs.

RESULTS

We included 211 diagnostic-quality phase-contrast sequences (93 aorta, 74 main pulmonary artery [MPA], 21 left pulmonary artery [LPA], 23 right pulmonary artery [RPA]) from 108 children (median age 15 years, interquartile range 11-18 years). Intraclass correlation showed strong agreement between non-linear and phantom-corrected flow measurements but also between uncorrected and phantom-corrected flow measurements. Non-linear flow measurements did not more closely approximate phantom-corrected measurements than did uncorrected measurements for any vessel. In 39 children without significant shunting or regurgitation, mean non-linear Qp:Qs (1.07; 95% confidence interval [CI] = 1.01, 1.13) was no closer than mean uncorrected Qp:Qs (1.06; 95% CI = 1.00, 1.13) to mean phantom-corrected Qp:Qs (1.02; 95% CI = 0.98, 1.06).

CONCLUSION

Despite strong agreement between self-calibrated non-linear and phantom correction, cardiac flows and shunt calculations with non-linear correction were no closer to phantom-corrected measurements than those without background correction. However, phantom-corrected flows also demonstrated minimal differences from uncorrected flows. These findings suggest that in the current era, more accurate phase-contrast flow measurements might limit the need for background correction. Further investigation of the clinical impact and optimal methods of background correction in the pediatric and congenital cardiac population is needed.

An MRI-based switched gradient impulse response characterization method with uniform eigenmode excitation

Switching gradients generate eddy currents and mechanical vibrations of the gradient assembly causing errors in the gradient time integrals. This results in image distortions in k-space and inaccuracies in q-space imaging. The purpose of this work is to develop an MRI based unbiased measurement of the switched gradient impulse response function (sGIRF). A new gradient pattern, called the Tukey windowed Shifted Sine-Integral (Tw-SSI) pulse, is introduced to excite the gradient eigenmodes uniformly over a user-defined bandwidth. A 3D MRI-based method with Hadamard encoding was developed to map the spatiotemporal magnetic field generated after the excitation pulse to obtain the sGIRF for all the three gradient axes simultaneously. Compared to an energy-equivalent traditional trapezoidal pulse, the Tw-SSI pulse is able to excite the weak bandlimited cross-terms of the sGIRF by uniformly distributing the energy across eigenmodes. The developed field mapping method is sensitive enough to capture both the direct and cross-terms in the sGIRF. The various mechanical resonant modes of the gradient coils are also revealed, which were found to last longer than eddy currents in the shielded gradient coil studied. Tunable Tw-SSI pulse offers the flexibility to perform unbiased sGIRF measurements over a bandwidth of interest. Rapid MRI field mapping can be easily implemented in any MRI system. The method may be used to perform gradient pre-emphasis, to evaluate new gradient coil designs, and to characterize higher order shims.

Stationary tissue background correction increases the precision of clinical evaluation of intra-cardiac shunts by cardiovascular magnetic resonance

We aimed to evaluate the clinical utility of stationary tissue background phase correction for affecting precision in the measurement of Qp/Qs by cardiovascular magnetic resonance (CMR). We enrolled consecutive patients (n = 91) referred for CMR at 1.5T without suspicion of cardiac shunt, and patients (n = 10) with verified cardiac shunts in this retrospective study. All patients underwent phase contrast flow quantification in the ascending aorta and pulmonary trunk. Flow was quantified using two semi-automatic software platforms (SyngoVia VA30, Vendor 1; Segment 2.0R4534, Vendor 2). Measurements were performed both uncorrected and corrected for linear (Vendor 1 and Vendor 2) or quadratic (Vendor 2) background phase. The proportion of patients outside the normal range of Qp/Qs was compared using the McNemar’s test. Compared to uncorrected measurements, there were fewer patients with a Qp/Qs outside the normal range following linear correction using Vendor 1 (10% vs 18%, p < 0.001), and Vendor 2 (10% vs 18%, p < 0.001), and following quadratic correction using Vendor 2 (7% vs 18%, p < 0.001). No patient with known shunt was reclassified as normal following stationary background correction. Therefore, we conclude that stationary tissue background correction reduces the number of patients with a Qp/Qs ratio outside the normal range in a consecutive clinical population, while simultaneously not reclassifying any patient with known cardiac shunts as having a normal Qp/Qs. Stationary tissue background correction may be used in clinical patients to increase diagnostic precision.

Gradient Response Harvesting for Continuous System Characterization During MR Sequences

MRI gradient systems are required to generate magnetic field gradient waveforms with very high fidelity. This is commonly implemented by gradient system calibration and pre-emphasis. However, a number of mechanisms, particularly thermal changes, cause variation in the gradient response over time, which cannot be addressed by calibration approaches. To overcome this limitation, we present a novel method termed gradient response harvesting, where the gradient response is continuously characterized during the course of a normal MR sequence. Snippets of field measurements are repeatedly acquired during an MR sequence, and from these multiple field measurements and the known nominal MR sequence gradients, the gradient response and gradient/field offsets are calculated. The calculation is implemented in a model-based and a model-free variant. The method is demonstrated for EPI with high gradient duty-cycle, where the continuous gradient characterization is used to obtain k-space trajectory estimates that are employed in the subsequent image reconstruction. During the course of the MR sequence, changes in both the envelope and the phase of the gradient response functions were observed, including shifts of mechanical resonances. The gradient response changes were also reflected in the calculated uninterrupted gradient waveforms and thus in the k-space trajectories. Using the updated encoding information in the image reconstruction removed ghosting artifacts, that otherwise impaired the image quality. We introduced the concept of gradient response harvesting and demonstrated its feasibility. The obtained gradient response functions may be used for quality assurance/preventive maintenance, real-time adaptation of gradient pre-emphasis or to calculate uninterrupted gradient field evolutions.

Comparison of 4D Flow MRI to 2D Flow MRI in the pulmonary arteries in healthy volunteers and patients with pulmonary hypertension

Purpose

4D and 2D phase-contrast MRI (2D Flow MRI, 4D Flow MRI, respectively) are increasingly being used to noninvasively assess pulmonary hypertension (PH). The goals of this study were i) to evaluate whether established quantitative parameters in 2D Flow MRI associated with pulmonary hypertension can be assessed using 4D Flow MRI; ii) to compare results from 4D Flow MRI on a digital broadband 3T MR system with data from clinically established MRI-techniques as well as conservation of mass analysis and phantom correction and iii) to elaborate on the added value of secondary flow patterns in detecting PH.

Methods

11 patients with PH (4f, 63 ± 16y), 15 age-matched healthy volunteers (9f, 56 ± 11y), and 20 young healthy volunteers (13f, 23 ± 2y) were scanned on a 3T MR scanner (Philips Ingenia). Subjects were examined with a 4D Flow, a 2D Flow and a bSSFP sequence. For extrinsic comparison, quantitative parameters measured with 4D Flow MRI were compared to i) a static phantom, ii) 2D Flow acquisitions and iii) stroke volume derived from a bSSFP sequence. For intrinsic comparison conservation of mass-analysis was employed. Dedicated software was used to extract various flow, velocity, and anatomical parameters. Visualization of blood flow was performed to detect secondary flow patterns.

Results

Overall, there was good agreement between all techniques, 4D Flow results revealed a considerable spread. Data improved after phantom correction. Both 4D and 2D Flow MRI revealed concordant results to differentiate patients from healthy individuals, especially based on values derived from anatomical parameters. The visualization of a vortex, indicating the presence of PH was achieved in 9 /11 patients and 2/35 volunteers.

Discussion

This study confirms that quantitative parameters used for characterizing pulmonary hypertension can be gathered using 4D Flow MRI within clinically reasonable limits of agreement. Despite its unfavorable spatial and lesser temporal resolution and a non-neglible spread of results, the identification of diseased study participants was possible.

In-vivo validation of interpolation-based phase offset correction in cardiovascular magnetic resonance flow quantification: a multi-vendor, multi-center study

BACKGROUND

A velocity offset error in phase contrast cardiovascular magnetic resonance (CMR) imaging is a known problem in clinical assessment of flow volumes in vessels around the heart. Earlier studies have shown that this offset error is clinically relevant over different systems, and cannot be removed by protocol optimization. Correction methods using phantom measurements are time consuming, and assume reproducibility of the offsets which is not the case for all systems. An alternative previously published solution is to correct the in-vivo data in post-processing, interpolating the velocity offset from stationary tissue within the field-of-view. This study aims to validate this interpolation-based offset correction in-vivo in a multi-vendor, multi-center setup.

METHODS

Data from six 1.5 T CMR systems were evaluated, with two systems from each of the three main vendors. At each system aortic and main pulmonary artery 2D flow studies were acquired during routine clinical or research examinations, with an additional phantom measurement using identical acquisition parameters. To verify the phantom acquisition, a region-of-interest (ROI) at stationary tissue in the thorax wall was placed and compared between in-vivo and phantom measurements. Interpolation-based offset correction was performed on the in-vivo data, after manually excluding regions of spatial wraparound. Correction performance of different spatial orders of interpolation planes was evaluated.

RESULTS

A total of 126 flow measurements in 82 subjects were included. At the thorax wall the agreement between in-vivo and phantom was - 0.2 ± 0.6 cm/s. Twenty-eight studies were excluded because of a difference at the thorax wall exceeding 0.6 cm/s from the phantom scan, leaving 98. Before correction, the offset at the vessel as assessed in the phantom was - 0.4 ± 1.5 cm/s, which resulted in a - 5 ± 16% error in cardiac output. The optimal order of the interpolation correction plane was 1st order, except for one system at which a 2nd order plane was required. Application of the interpolation-based correction revealed a remaining offset velocity of 0.1 ± 0.5 cm/s and 0 ± 5% error in cardiac output.

CONCLUSIONS

This study shows that interpolation-based offset correction reduces the offset with comparable efficacy as phantom measurement phase offset correction, without the time penalty imposed by phantom scans.

TRIAL REGISTRATION

The study was registered in The Netherlands National Trial Register (NTR) under TC 4865 . Registered 19 September 2014. Retrospectively registered.

Segmental differences of cervical spinal cord motion: advancing from confounders to a diagnostic tool

Increased cranio-caudal spinal cord motion is associated with clinical impairment in degenerative cervical myelopathy. However, whether spinal cord motion holds potential as a neuroimaging biomarker requires further validation. Different confounders (i.e. subject characteristics, methodological problems such as phase drift, etc.) on spinal cord motion readouts have to be considered. Twenty-two healthy subjects underwent phase contrast MRI, a subset of subjects (N = 9) had repeated scans. Parameters of interest included amplitude of velocity signal, maximum cranial respectively maximum caudal velocity, displacement (=area under curve of the velocity signal). The cervical spinal cord showed pulse synchronic oscillatory motions with significant differences in all readouts across cervical segments, with a maximum at C5. The Inter-rater reliability was excellent for all readouts. The test-retest reliability was excellent for all parameters at C2 to C6, but not for maximum cranial velocity at C6 and all readouts at C7. Spinal cord motion was correlated with spinal canal size, heart rate and body size. This is the first study to propose a standardized MRI measurement of spinal cord motion for further clinical implementation based on satisfactory phase drift correction and excellent reliability. Understanding the influence of confounders (e.g. structural conditions of the spine) is essential for introducing cord motion into the diagnostic work up.

Influence of respiration-induced B0 variations in real-time phase-contrast echo planar imaging of the cervical cerebrospinal fluid

PURPOSE

Respiration induces temporal variations of the main magnetic field B₀ along the spinal cord. These variations are typically not compensated for in velocity quantifications using phase-contrast MRI. The goal of this study was to analyze errors caused by respiration-induced B₀ variations in real-time phase-contrast echo planar imaging (PCEPI) of cervical cerebrospinal fluid (CSF) velocity measurements and to evaluate this effect for various sequence parameters using numerical simulations.

METHODS

Real-time B₀ measurements with double gradient echo sequence and PCEPI measurements were acquired in the cervical CSF of 10 healthy subjects. Dynamic phase offsets attributed to respiration-induced B₀ variations were analyzed by quantifying amplitudes and comparing the temporal behavior with respiratory signals. In experiments and simulations, the influence of the echo time (TE) and the delay between PCEPI images (Δt) with respect to respiration on the dynamic phase offsets were investigated.

RESULTS

A good agreement was found between phase offsets extracted from both acquisition types. Furthermore, respiratory signals qualitatively matched the temporal behavior of the measured phase offsets showing a dependency on subject-dependent local B₀ distribution and respiration physiology. Simulations revealed residual background phases in PCEPI velocity quantification varying with TE and Δt.

CONCLUSION

Respiration-induced B₀ variations result in dynamic background phases in real-time PCEPI velocity quantifications of the CSF in the cervical spine. The current work underlines that these background phases need to be corrected to avoid confounding effects.

Phase contrast MRI of creeping flows using stimulated echo

Creeping flows govern many important physiological phenomena such as elevated interstitial fluid flows in tumors, glymphatic flows in the brain, among other applications. However, few methods exist to measure such slow flows non-invasively in optically opaque biological tissues in vivo. Phase-contrast MRI is a velocimetry technique routinely used in the clinic to measure fast flows in biological tissues, such as blood and cerebrospinal fluid (CSF), in the order of cm/s. Use of this technique to encode slower flows is hampered by diffusion weighting and phase error introduced by gradient hardware imperfections. In this study, a new PC-MRI technique is developed using stimulated echo preparation to overcome these challenges. Flows as slow as 1 μm/s are measured and validated using controlled water flow through a pipe at 4.7 T. The error in measured flow rate obtained by integrating the measured velocity over the cross-sectional area of the pipe is less than 10%. The developed method was also able to capture slow natural convection flows appearing in liquids placed inside a horizontal bore magnet. Monitoring the 4D velocity vector field revealed that the natural convection flows decay exponentially with time. This method could be applied in future to study creeping flows, e.g. in tissue.

Effect of upper airway on tracheobronchial fluid dynamics

The upper airways play a significant role in the tracheal flow dynamics. Despite many previous studies, however, the effect of the upper airways on the ventilation distribution in distal airways has remained a challenge. The aim of this study is to experimentally and computationally investigate the dynamic behaviour in the intratracheal flow induced by the upper respiratory tract and to assess its influence on the subsequent tributaries. Patient-specific images from 2 different modalities (magnetic resonance imaging of the upper airways and computed tomography of the lower airways) were segmented and combined. An experimental phantom of patient-specific airways (including the oral cavity, larynx, trachea, down to generations 6-8) was generated using 3D printing. The flow velocities in this phantom model were measured by the flow-sensitised phase contrast magnetic resonance imaging technique and compared with the computational fluid dynamics simulations. Both experimental and computational results show a good agreement in the time-averaged velocity fields as well as fluctuating velocity. The flows in the proximal trachea were complex and unsteady under both lower- and higher-flow rate conditions. Computational fluid dynamics simulations were also performed with an airways model without the upper airways. Although the flow near the carina remained unstable only when the inflow rate was high, the influence of the upper airways caused notable changes in distal flow distributions when the 2 airways models were compared with and without the upper airways. The results suggest that the influence of the upper airways should be included in the respiratory flow assessment as the upper airways extensively affect the flows in distal airways and consequent ventilation distribution in the lungs.

Magnetic Resonance Imaging technology-bridging the gap between noninvasive human imaging and optical microscopy

Technological advances in Magnetic Resonance Imaging (MRI) have provided substantial gains in the sensitivity and specificity of functional neuroimaging. Mounting evidence demonstrates that the hemodynamic changes utilized in functional MRI can be far more spatially and thus neuronally specific than previously believed. This has motivated a push toward novel, high-resolution MR imaging strategies that can match this biological resolution limit while recording from the entire human brain. Although sensitivity increases are a necessary component, new MR encoding technologies are required to convert improved sensitivity into higher resolution. These new sampling strategies improve image acquisition efficiency and enable increased image encoding in the time-frame needed to follow hemodynamic changes associated with brain activation.

Radiofrequency saturation induced bias in aqueductal cerebrospinal fluid flow quantification obtained using two-dimensional cine phase contrast magnetic resonance imaging

PURPOSE

To explore the extent of bias in cerebrospinal fluid flow estimates due to radiofrequency saturation, and its possible impact on the use of two-dimensional cine phase contrast magnetic resonance imaging in the diagnosis and characterization of normal pressure hydrocephalus in patients.

THEORY AND METHODS

Theoretical signal equations were generated to describe saturation dependence on velocity. An experimental set of phase contrast magnetic resonance imaging scans with two different flip angles was used to show bias in flow estimates in a flow phantom, and in six different healthy volunteers. The cerebral aqueduct was targeted as the flow region of interest.

RESULTS

Data from a constant flow phantom showed a spatial distribution of voxels with significant bias in flow at the periphery of the flow region. The velocity difference (bias) maps of the cerebral aqueduct correlated with the spatial velocity gradients around peak systole and peak diastole, and high correlation with temporal velocity gradients during transition between systole and diastole. The aqueductal stroke volume for θ = 30° were found to be significantly higher than for θ = 10° using a Wilcoxon signed rank test.

CONCLUSION

This work shows the extent of bias in cerebrospinal fluid flow quantification due to radiofrequency saturation effects. This clinical relevance of this error was presented with respect to shunt responsiveness among normal pressure hydrocephalus patients. Magn Reson Med 79:2067-2076, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Comparison of radial 4D Flow-MRI with perivascular ultrasound to quantify blood flow in the abdomen and introduction of a porcine model of pre-hepatic portal hypertension

OBJECTIVES

Objectives of this study were to compare radial time-resolved phase contrast magnetic resonance imaging (4D Flow-MRI) with perivascular ultrasound (pvUS) and to explore a porcine model of acute pre-hepatic portal hypertension (PHTN).

METHODS

Abdominal 4D Flow-MRI and pvUS in portal and splenic vein, hepatic and both renal arteries were performed in 13 pigs of approximately 60 kg. In six pigs, measurements were repeated after partial portal vein (PV) ligature. Inter- and intra-reader comparisons and statistical analysis including Bland-Altman (BA) comparison, paired Student's t tests and linear regression were performed.

RESULTS

PvUS and 4D Flow-MRI measurements agreed well; flow before partial PV ligature was 322 ± 30 ml/min in pvUS and 297 ± 27 ml/min in MRI (p = 0.294), and average BA difference was 25 ml/min [-322; 372]. Inter- and intra-reader results differed very little, revealed excellent correlation (R ² = 0.98 and 0.99, respectively) and resulted in BA differences of -5 ml/min [-161; 150] and -2 ml/min [-28; 25], respectively. After PV ligature, PV flow decreased from 356 ± 50 to 298 ± 61 ml/min (p = 0.02), and hepatic arterial flow increased from 277 ± 36 to 331 ± 65 ml/min (p = n.s.).

CONCLUSION

The successful in vivo comparison of radial 4D Flow-MRI to perivascular ultrasound revealed good agreement of abdominal blood flow although with considerable spread of results. A model of pre-hepatic PHTN was successfully introduced and acute responses monitored.

KEY POINTS

• Radial 4D Flow-MRI in the abdomen was successfully compared to perivascular ultrasound. • Inter- and intra-reader testing demonstrated excellent reproducibility of upper abdominal 4D Flow-MRI. • A porcine model of acute pre-hepatic portal hypertension was successfully introduced. • 4D Flow-MRI successfully monitored acute changes in a model of portal hypertension.

Acquisition of electrocardiogram signals during magnetic resonance imaging

The recording of the electrocardiogram (ECG) during magnetic resonance imaging (MRI) acquisition is of great interest and importance. Firstly, MRI acquisition is a relatively slow process, which therefore complicates the imaging of moving organs. Cardiac MRI requires the development of strategies for acquiring high quality images, which is mainly achieved by synchronising the image acquisition with a specific time during the cardiac cycle. The ECG is used to monitor the heart's activity, and the detection of the largest and steepest peak in the cardiac cycle (the QRS complex) triggers the acquisition of slices of the k-space. Secondly, patients undergoing an MRI examination need to be monitored for safety during the procedure, and therefore ECG signals are used to track their cardiovascular state in real time. However, there are significant barriers to the accurate observation and processing of the ECG during MRI acquisition. In particular, the flow of charged blood particles through the large applied magnetic field leads to an extra current source, known as the magnetohdrodymanic (MHD) effect. This review article discusses these barriers and state-of-the-art solutions. An overview of the relevant technology including hardware and applications are described. The development of new software tools for the processing of the ECG signals acquired during MRI is also detailed. These developments include the design of specific QRS detection algorithms, which are able to distinguish QRS complexes from the MHD effect but also the gradient artefacts. Different techniques for the suppression of the gradient artefacts are also presented as well as the most challenging problem to-date-the problem of separating the MHD effect from the ECG. The article concludes by summarising the advantages of using ECG signals during MRI, but also presents the current limitations of modern analysis techniques in this domain. The most promising avenues of research are also discussed and suggestions for new methodological analyses for the development of this field are given.

Accuracy of 4D Flow Measurement of Cerebrospinal Fluid Dynamics in the Cervical Spine: An In Vitro Verification Against Numerical Simulation

Abnormal alterations in cerebrospinal fluid (CSF) flow are thought to play an important role in pathophysiology of various craniospinal disorders such as hydrocephalus and Chiari malformation. Three directional phase contrast MRI (4D Flow) has been proposed as one method for quantification of the CSF dynamics in healthy and disease states, but prior to further implementation of this technique, its accuracy in measuring CSF velocity magnitude and distribution must be evaluated. In this study, an MR-compatible experimental platform was developed based on an anatomically detailed 3D printed model of the cervical subarachnoid space and subject specific flow boundary conditions. Accuracy of 4D Flow measurements was assessed by comparison of CSF velocities obtained within the in vitro model with the numerically predicted velocities calculated from a spatially averaged computational fluid dynamics (CFD) model based on the same geometry and flow boundary conditions. Good agreement was observed between CFD and 4D Flow in terms of spatial distribution and peak magnitude of through-plane velocities with an average difference of 7.5 and 10.6% for peak systolic and diastolic velocities, respectively. Regression analysis showed lower accuracy of 4D Flow measurement at the timeframes corresponding to low CSF flow rate and poor correlation between CFD and 4D Flow in-plane velocities.

Advanced flow MRI: emerging techniques and applications

Magnetic resonance imaging (MRI) techniques provide non-invasive and non-ionising methods for the highly accurate anatomical depiction of the heart and vessels throughout the cardiac cycle. In addition, the intrinsic sensitivity of MRI to motion offers the unique ability to acquire spatially registered blood flow simultaneously with the morphological data, within a single measurement. In clinical routine, flow MRI is typically accomplished using methods that resolve two spatial dimensions in individual planes and encode the time-resolved velocity in one principal direction, typically oriented perpendicular to the two-dimensional (2D) section. This review describes recently developed advanced MRI flow techniques, which allow for more comprehensive evaluation of blood flow characteristics, such as real-time flow imaging, 2D multiple-venc phase contrast MRI, four-dimensional (4D) flow MRI, quantification of complex haemodynamic properties, and highly accelerated flow imaging. Emerging techniques and novel applications are explored. In addition, applications of these new techniques for the improved evaluation of cardiovascular (aorta, pulmonary arteries, congenital heart disease, atrial fibrillation, coronary arteries) as well as cerebrovascular disease (intra-cranial arteries and veins) are presented.

Phase Error Correction in Time-Averaged 3D Phase Contrast Magnetic Resonance Imaging of the Cerebral Vasculature

Purpose

Volume flow rate (VFR) measurements based on phase contrast (PC)-magnetic resonance (MR) imaging datasets have spatially varying bias due to eddy current induced phase errors. The purpose of this study was to assess the impact of phase errors in time averaged PC-MR imaging of the cerebral vasculature and explore the effects of three common correction schemes (local bias correction (LBC), local polynomial correction (LPC), and whole brain polynomial correction (WBPC)).

Methods

Measurements of the eddy current induced phase error from a static phantom were first obtained. In thirty healthy human subjects, the methods were then assessed in background tissue to determine if local phase offsets could be removed. Finally, the techniques were used to correct VFR measurements in cerebral vessels and compared statistically.

Results

In the phantom, phase error was measured to be <2.1 ml/s per pixel and the bias was reduced with the correction schemes. In background tissue, the bias was significantly reduced, by 65.6% (LBC), 58.4% (LPC) and 47.7% (WBPC) (p < 0.001 across all schemes). Correction did not lead to significantly different VFR measurements in the vessels (p = 0.997). In the vessel measurements, the three correction schemes led to flow measurement differences of -0.04 ± 0.05 ml/s, 0.09 ± 0.16 ml/s, and -0.02 ± 0.06 ml/s. Although there was an improvement in background measurements with correction, there was no statistical difference between the three correction schemes (p = 0.242 in background and p = 0.738 in vessels).

Conclusions

While eddy current induced phase errors can vary between hardware and sequence configurations, our results showed that the impact is small in a typical brain PC-MR protocol and does not have a significant effect on VFR measurements in cerebral vessels.

Cardiac index after acute ST-segment elevation myocardial infarction measured with phase-contrast cardiac magnetic resonance imaging

OBJECTIVES

Phase-contrast CMR (PC-CMR) might provide a fast and robust non-invasive determination of left ventricular function in patients after ST-segment elevation myocardial infarction (STEMI).

METHODS

Cine sequences in the left-ventricular (LV) short-axis and free-breathing, retrospectively gated PC-CMR were performed in 90 patients with first acute STEMI and 15 healthy volunteers. Inter- and intra-observer agreement was determined. The correlations of clinical variables (age, gender, ejection fraction, NT pro-brain natriuretic peptide [NT-proBNP] with cardiac index (CI) were calculated.

RESULTS

For CI, there was a strong agreement of cine CMR with PC-CMR in healthy volunteers (r: 0.82, mean difference: -0.14 l/min/m(2), error ± 23 %). Agreement was lower in STEMI patients (r: 0.61, mean difference: -0.17 l/min/m(2), error ± 32 %). In STEMI patients, CI measured with PC-CMR showed lower intra-observer (1 % vs. 9 %) and similar inter-observer variability (9 % vs. 12 %) compared to cine CMR. CI was significantly correlated with age, ejection fraction and NT-proBNP values in STEMI patients.

DISCUSSION

The agreement of PC-CMR and cine CMR for the determination of CI is lower in STEMI patients than in healthy volunteers. After acute STEMI, CI measured with PC-CMR decreases with age, LV ejection fraction and higher NT-proBNP.

KEY POINTS

• Cine CMR and PC-CMR correlate well in healthy volunteers. • Agreement is lower in STEMI patients. • Cardiac Output should be measured with one method longitudinally. • Cardiac output decreases with age after myocardial infarction.

Effects of contour propagation and background corrections in different MRI flow software packages

Background

Velocity-encoded magnetic resonance imaging (VENC-MRI) is a commonly used technique in cardiac examinations. This technique utilizes the phase shift properties of protons moving along a magnetic field gradient. VENC-MRI offers a unique way of measuring the severity of valve regurgitation by directly quantifying the regurgitation flow volume.

Purpose

To compare flow analysis results of different software programs and to assess the effect of background correction in sample patient cases.

Material and Methods

A phantom was built out of Polymethyl methacrylate (PMMA) which provides tubes of different diameters. These tubes can be connected to an external water circuit to generate a water flow inside the tubes. Expected absolute flow quantities inside the tubes were determined from preset tube- and flow-parameters. Different flow conditions were measured with a VENC-MRI sequence and the images evaluated using different software packages. In a second step six randomly selected patients showing different degrees of aortic insufficiency were evaluated in clinical terms.

Results

The contour propagation algorithms used in the software packages performed differently even on static phantom geometry. In terms of clinical evaluation the software packages performed similarly. Enabling background correction or leaving out manual correction of propagated contours changed results for severity of aortic insufficiency.

Conclusion

Turning on background correction and manual correction of propagated contours in MRI flow volume measurements is strongly recommended.

4D flow cardiovascular magnetic resonance consensus statement

Pulsatile blood flow through the cavities of the heart and great vessels is time-varying and multidirectional. Access to all regions, phases and directions of cardiovascular flows has formerly been limited. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has enabled more comprehensive access to such flows, with typical spatial resolution of 1.5×1.5×1.5 - 3×3×3 mm(3), typical temporal resolution of 30-40 ms, and acquisition times in the order of 5 to 25 min. This consensus paper is the work of physicists, physicians and biomedical engineers, active in the development and implementation of 4D Flow CMR, who have repeatedly met to share experience and ideas. The paper aims to assist understanding of acquisition and analysis methods, and their potential clinical applications with a focus on the heart and greater vessels. We describe that 4D Flow CMR can be clinically advantageous because placement of a single acquisition volume is straightforward and enables flow through any plane across it to be calculated retrospectively and with good accuracy. We also specify research and development goals that have yet to be satisfactorily achieved. Derived flow parameters, generally needing further development or validation for clinical use, include measurements of wall shear stress, pressure difference, turbulent kinetic energy, and intracardiac flow components. The dependence of measurement accuracy on acquisition parameters is considered, as are the uses of different visualization strategies for appropriate representation of time-varying multidirectional flow fields. Finally, we offer suggestions for more consistent, user-friendly implementation of 4D Flow CMR acquisition and data handling with a view to multicenter studies and more widespread adoption of the approach in routine clinical investigations.

Cardiovascular magnetic resonance phase contrast imaging

Cardiovascular magnetic resonance (CMR) phase contrast imaging has undergone a wide range of changes with the development and availability of improved calibration procedures, visualization tools, and analysis methods. This article provides a comprehensive review of the current state-of-the-art in CMR phase contrast imaging methodology, clinical applications including summaries of past clinical performance, and emerging research and clinical applications that utilize today's latest technology.

Accelerating 4D flow MRI by exploiting vector field divergence regularization

PURPOSE

To improve velocity vector field reconstruction from undersampled four-dimensional (4D) flow MRI by penalizing divergence of the measured flow field.

THEORY AND METHODS

Iterative image reconstruction in which magnitude and phase are regularized separately in alternating iterations was implemented. The approach allows incorporating prior knowledge of the flow field being imaged. In the present work, velocity data were regularized to reduce divergence, using either divergence-free wavelets (DFW) or a finite difference (FD) method using the ℓ1-norm of divergence and curl. The reconstruction methods were tested on a numerical phantom and in vivo data. Results of the DFW and FD approaches were compared with data obtained with standard compressed sensing (CS) reconstruction.

RESULTS

Relative to standard CS, directional errors of vector fields and divergence were reduced by 55-60% and 38-48% for three- and six-fold undersampled data with the DFW and FD methods. Velocity vector displays of the numerical phantom and in vivo data were found to be improved upon DFW or FD reconstruction.

CONCLUSION

Regularization of vector field divergence in image reconstruction from undersampled 4D flow data is a valuable approach to improve reconstruction accuracy of velocity vector fields.

Analysis of temperature dependence of background phase errors in phase-contrast cardiovascular magnetic resonance

BACKGROUND

The accuracy of phase-contrast cardiovascular magnetic resonance (PC-CMR) can be compromised by background phase errors. It is the objective of the present work to provide an analysis of the temperature dependence of background phase errors in PC-CMR by means of gradient mount temperature sensing and magnetic field monitoring.

METHODS

Background phase errors were measured for various temperatures of the gradient mount using magnetic field monitoring and validated in a static phantom. The effect of thermal changes during k-space acquisition was simulated and confirmed with measurements in a stationary phantom.

RESULTS

The temperature of the gradient mount was found to increase by 20-30 K during PC-CMR measurements of 6-12 min duration. Associated changes in background phase errors of up to 11% or 0.35 radian were measured at 10 cm from the magnet's iso-center as a result of first order offsets. Zeroth order phase errors exhibited little thermal dependence.

CONCLUSIONS

It is concluded that changes in gradient mount temperature significantly modify background phase errors during PC-CMR with high gradient duty cycle. Since temperature increases significantly during the first minutes of scanning the results presented are also of relevance for single-slice or multi-slice PC-CMR scans. The findings prompt for further studies to investigate advanced correction methods taking into account gradient temperature and/or the use of concurrent field-monitoring to map gradient-induced fields throughout the scan.

Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring

PURPOSE

Ultrashort echo time (UTE) MRI has been proposed as a way to produce segmented attenuation maps for PET, as it provides contrast between bone, air, and soft tissue. However, UTE sequences require samples to be acquired during rapidly changing gradient fields, which makes the resulting images prone to eddy current artifacts. In this work it is demonstrated that this can lead to misclassification of tissues in segmented attenuation maps (AC maps) and that these effects can be corrected for by measuring the true k-space trajectories using a magnetic field camera.

METHODS

The k-space trajectories during a dual echo UTE sequence were measured using a dynamic magnetic field camera. UTE images were reconstructed using nominal trajectories and again using the measured trajectories. A numerical phantom was used to demonstrate the effect of reconstructing with incorrect trajectories. Images of an ovine leg phantom were reconstructed and segmented and the resulting attenuation maps were compared to a segmented map derived from a CT scan of the same phantom, using the Dice similarity measure. The feasibility of the proposed method was demonstrated in in vivo cranial imaging in five healthy volunteers. Simulated PET data were generated for one volunteer to show the impact of misclassifications on the PET reconstruction.

RESULTS

Images of the numerical phantom exhibited blurring and edge artifacts on the bone-tissue and air-tissue interfaces when nominal k-space trajectories were used, leading to misclassification of soft tissue as bone and misclassification of bone as air. Images of the tissue phantom and the in vivo cranial images exhibited the same artifacts. The artifacts were greatly reduced when the measured trajectories were used. For the tissue phantom, the Dice coefficient for bone in MR relative to CT was 0.616 using the nominal trajectories and 0.814 using the measured trajectories. The Dice coefficients for soft tissue were 0.933 and 0.934 for the nominal and measured cases, respectively. For air the corresponding figures were 0.991 and 0.993. Compared to an unattenuated reference image, the mean error in simulated PET uptake in the brain was 9.16% when AC maps derived from nominal trajectories was used, with errors in the SUV max for simulated lesions in the range of 7.17%-12.19%. Corresponding figures when AC maps derived from measured trajectories were used were 0.34% (mean error) and -0.21% to +1.81% (lesions).

CONCLUSIONS

Eddy current artifacts in UTE imaging can be corrected for by measuring the true k-space trajectories during a calibration scan and using them in subsequent image reconstructions. This improves the accuracy of segmented PET attenuation maps derived from UTE sequences and subsequent PET reconstruction.

Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring

Diffusion tensor imaging (DTI) of moving organs is gaining increasing attention but robust performance requires sequence modifications and dedicated correction methods to account for system imperfections. In this study, eddy currents in the "unipolar" Stejskal-Tanner and the velocity-compensated "bipolar" spin-echo diffusion sequences were investigated and corrected for using a magnetic field monitoring approach in combination with higher-order image reconstruction. From the field-camera measurements, increased levels of second-order eddy currents were quantified in the unipolar sequence relative to the bipolar diffusion sequence while zeroth and linear orders were found to be similar between both sequences. Second-order image reconstruction based on field-monitoring data resulted in reduced spatial misalignment artifacts and residual displacements of less than 0.43 mm and 0.29 mm (in the unipolar and bipolar sequences, respectively) after second-order eddy-current correction. Results demonstrate the need for second-order correction in unipolar encoding schemes but also show that bipolar sequences benefit from second-order reconstruction to correct for incomplete intrinsic cancellation of eddy-currents.

Towards highly accelerated Cartesian time-resolved 3D flow cardiovascular magnetic resonance in the clinical setting

BACKGROUND

The clinical applicability of time-resolved 3D flow cardiovascular magnetic resonance (CMR) remains compromised by the long scan times associated with phase-contrast imaging. The present work demonstrates the applicability of 8-fold acceleration of Cartesian time-resolved 3D flow CMR in 10 volunteers and in 9 patients with different congenital heart diseases (CHD). It is demonstrated that accelerated 3D flow CMR data acquisition and image reconstruction using k-t PCA (principal component analysis) can be implemented into clinical workflow and results are sufficiently accurate relative to conventional 2D flow CMR to permit for comprehensive flow quantification in CHD patients.

METHODS

The fidelity of k-t PCA was first investigated on retrospectively undersampled data for different acceleration factors and compared to k-t SENSE and fully sampled reference data. Subsequently, k-t PCA with 8-fold nominal undersampling was applied on 10 healthy volunteers and 9 CHD patients on a clinical 1.5 T MR scanner. Quantitative flow validation was performed in vessels of interest on the 3D flow datasets and compared to 2D through-plane flow acquisitions. Particle trace analysis was used to qualitatively visualise flow patterns in patients.

RESULTS

Accelerated time-resolved 3D flow data were successfully acquired in all subjects with 8-fold nominal scan acceleration. Nominal scan times excluding navigator efficiency were on the order of 6 min and 7 min in patients and volunteers. Mean differences in stroke volume in selected vessels of interest were 2.5 ± 8.4 ml and 1.63 ± 4.8 ml in volunteers and patients, respectively. Qualitative flow pattern analysis in the time-resolved 3D dataset revealed valuable insights into hemodynamics including circular and helical patterns as well as flow distributions and origin in the Fontan circulation.

CONCLUSION

Highly accelerated time-resolved 3D flow using k-t PCA is readily applicable in clinical routine protocols of CHD patients. Nominal scan times of 6 min are well tolerated and allow for quantitative and qualitative flow assessment in all great vessels.

Volumetric velocity measurements in restricted geometries using spiral sampling: a phantom study

OBJECT

The aim of this study was to evaluate the accuracy of maximum velocity measurements using volumetric phase-contrast imaging with spiral readouts in a stenotic flow phantom.

MATERIALS AND METHODS

In a phantom model, maximum velocity, flow, pressure gradient, and streamline visualizations were evaluated using volumetric phase-contrast magnetic resonance imaging (MRI) with velocity encoding in one (extending on current clinical practice) and three directions (for characterization of the flow field) using spiral readouts. Results of maximum velocity and pressure drop were compared to computational fluid dynamics (CFD) simulations, as well as corresponding low-echo-time (TE) Cartesian data. Flow was compared to 2D through-plane phase contrast (PC) upstream from the restriction.

RESULTS

Results obtained with 3D through-plane PC as well as 4D PC at shortest TE using a spiral readout showed excellent agreements with the maximum velocity values obtained with CFD (<1 % for both methods), while larger deviations were seen using Cartesian readouts (-2.3 and 13 %, respectively). Peak pressure drop calculations from 3D through-plane PC and 4D PC spiral sequences were respectively 14 and 13 % overestimated compared to CFD.

CONCLUSION

Identification of the maximum velocity location, as well as the accurate velocity quantification can be obtained in stenotic regions using short-TE spiral volumetric PC imaging.

Postsurgical hemodynamics of the aortic valve bypass operation evaluated with phase contrast magnetic resonance

PURPOSE

To characterize the postsurgical hemodynamics in aortic valve bypass (AVB) patients, and to determine the relationship between presurgical native aortic valve pressure gradient and postsurgical hemodynamics.

MATERIALS AND METHODS

Twenty patients scheduled for AVB surgery underwent presurgical transthoracic Doppler echocardiography to assess the degree of aortic stenosis and postsurgical cardiac magnetic resonance imaging (MRI) to acquire phase contrast magnetic resonance (PCMR) flow values along the ascending and descending aorta, and in the conduit. Net flow values were calculated from the PCMR images and compared to presurgical aortic valve pressure gradient measurements.

RESULTS

PCMR showed that: 1) The blood flow split between the aorta and the conduit was 35%:65% of cardiac output and 2) 60% of patients had net retrograde blood flow in the superior thoracic aorta over the cardiac cycle. Patients with presurgical pressure gradient (ΔP) > 45 mmHg had significantly less blood flow out of the native aorta than patients with ΔP < 45 mmHg, and had significantly more retrograde flow in the superior thoracic aorta postsurgery.

CONCLUSION

In patients undergoing AVB, presurgical aortic valve pressure gradient is associated with the volume of blood flow out the aorta and the direction of blood flow in the superior thoracic aorta after conduit addition as measured by PCMR.

Field camera measurements of gradient and shim impulse responses using frequency sweeps

PURPOSE

Applications of dynamic shimming require high field fidelity, and characterizing the shim field dynamics is therefore necessary. Modeling the system as linear and time-invariant, the purpose of this work was to measure the impulse response function with optimal sensitivity.

THEORY AND METHODS

Frequency-swept pulses as inputs are analyzed theoretically, showing that the sweep speed is a key factor for the measurement sensitivity. By adjusting the sweep speed it is possible to achieve any prescribed noise profile in the measured system response. Impulse response functions were obtained for the third-order shim system of a 7 Tesla whole-body MR scanner. Measurements of the shim fields were done with a dynamic field camera, yielding also cross-term responses.

RESULTS

The measured shim impulse response functions revealed system characteristics such as response bandwidth, eddy currents and specific resonances, possibly of mechanical origin. Field predictions based on the shim characterization were shown to agree well with directly measured fields, also in the cross-terms.

CONCLUSION

Frequency sweeps provide a flexible tool for shim or gradient system characterization. This may prove useful for applications involving dynamic shimming by yielding accurate estimates of the shim fields and a basis for setting shim pre-emphasis.

A multi-center inter-manufacturer study of the temporal stability of phase-contrast velocity mapping background offset errors

BACKGROUND

Phase-contrast velocity images often contain a background or baseline offset error, which adds an unknown offset to the measured velocities. For accurate flow measurements, this offset must be shown negligible or corrected. Some correction techniques depend on replicating the clinical flow acquisition using a uniform stationary phantom, in order to measure the baseline offset at the region of interest and subtract it from the clinical study. Such techniques assume that the background offset is stable over the time of a patient scan, or even longer if the phantom scans are acquired later, or derived from pre-stored background correction images. There is no published evidence regarding temporal stability of the background offset.

METHODS

This study assessed the temporal stability of the background offset on 3 different manufacturers' scanners over 8 weeks, using a retrospectively-gated phase-contrast cine acquisition with fixed parameters and at a fixed location, repeated 5 times in rapid succession each week. A significant offset was defined as 0.6 cm/s within 50 mm of isocenter, based upon an accuracy of 10% in a typical cardiac shunt measurement.

RESULTS

Over the 5 repeated cine acquisitions, temporal drift in the baseline offset was insignificant on two machines (0.3 cm/s, 0.2 cm/s), and marginally insignificant on the third machine (0.5 cm/s) due to an apparent heating effect. Over a longer timescale of 8 weeks, insignificant drift (0.4 cm/s) occurred on one, with larger drifts (0.9 cm/s, 0.6 cm/s) on the other machines.

CONCLUSIONS

During a typical patient study, background drift was insignificant. Extended high gradient power scanning with work requires care to avoid drift on some machines. Over the longer term of 8 weeks, significant drift is likely, preventing accurate correction by delayed phantom corrections or derivation from pre-stored background offset data.

Gradient system characterization by impulse response measurements with a dynamic field camera

This work demonstrates a fast, sensitive method of characterizing the dynamic performance of MR gradient systems. The accuracy of gradient time-courses is often compromised by field imperfections of various causes, including eddy currents and mechanical oscillations. Characterizing these perturbations is instrumental for corrections by pre-emphasis or post hoc signal processing. Herein, a gradient chain is treated as a linear time-invariant system, whose impulse response function is determined by measuring field responses to known gradient inputs. Triangular inputs are used to probe the system and response measurements are performed with a dynamic field camera consisting of NMR probes. In experiments on a whole-body MR system, it is shown that the proposed method yields impulse response functions of high temporal and spectral resolution. Besides basic properties such as bandwidth and delay, it also captures subtle features such as mechanically induced field oscillations. For validation, measured response functions were used to predict gradient field evolutions, which was achieved with an error below 0.2%. The field camera used records responses of various spatial orders simultaneously, rendering the method suitable also for studying cross-responses and dynamic shim systems. It thus holds promise for a range of applications, including pre-emphasis optimization, quality assurance, and image reconstruction.