Sparsity transform k-t principal component analysis for accelerating cine three-dimensional flow measurements

K-t PCA accelerated in-plane balanced steady-state free precession phase-contrast (PC-SSFP) for all-in-one diastolic function evaluation

PURPOSE

Diastolic function evaluation requires estimates of early and late diastolic mitral filling velocities (E and A) and of mitral annulus tissue velocity (e'). We aimed to develop an MRI method for simultaneous all-in-one diastolic function evaluation in a single scan by generating a 2D phase-contrast (PC) sequence with balanced steady-state free precession (bSSFP) contrast (PC-SSFP). E and A could then be measured with PC, and e' estimated by valve tracking on the magnitude images, using an established deep learning framework.

METHODS

Our PC-SSFP used in-plane flow-encoding, with zeroth and first moment nulling over each TR. For further acceleration, different k-t principal component analysis (PCA) methods were investigated with both retrospective and prospective undersampling. PC-SSFP was compared to separate balanced SSFP cine and PC-gradient echo acquisitions in phantoms and in 10 healthy subjects.

RESULTS

Phantom experiments showed that PC-SSFP measured accurate velocities compared to PC-gradient echo (r = 0.98 for a range of pixel-wise velocities -80 cm/s to 80 cm/s). In subjects, PC-SSFP generated high SNR and myocardium-blood contrast, and excellent agreement for E (limits of agreement [LOA] 0.8 ± 2.4 cm/s, r = 0.98), A (LOA 2.5 ± 4.1 cm/s, r = 0.97), and e' (LOA 0.3 ± 2.6 cm/s, r = 1.00), versus the standard methods. The best k-t PCA approach processed the complex difference data and substituted in raw k-space data. With prospective k-t PCA acceleration, higher frame rates were achieved (50 vs. 25 frames per second without k-t PCA), yielding a 13% higher e'.

CONCLUSION

The proposed PC-SSFP method achieved all-in-one diastolic function evaluation.

The role of 4-dimensional flow in the assessment of bicuspid aortic valve and its valvulo-aortopathies

Bicuspid aortic valve is the most common congenital cardiac malformation and the leading cause of aortopathy and aortic stenosis in younger patients. Aortic wall remodelling secondary to altered haemodynamic flow patterns, changes in peak velocity, and wall shear stress may be implicated in the development of aortopathy in the presence of bicuspid aortic valve and dysfunction. Assessment of these parameters as potential predictors of disease severity and progression is thus desirable. The anatomic and functional information acquired from 4D flow MRI can allow simultaneous visualisation and quantification of the pathological geometric and haemodynamic changes of the aorta. We review the current clinical utility of haemodynamic quantities including velocity, wall sheer stress and energy losses, as well as visual descriptors such as vorticity and helicity, and flow direction in assessing the aortic valve and associated aortopathies.

Altered Aortic Hemodynamics and Relative Pressure in Patients with Dilated Cardiomyopathy

Ventricular-vascular interaction is central in the adaptation to cardiovascular disease. However, cardiomyopathy patients are predominantly monitored using cardiac biomarkers. The aim of this study is therefore to explore aortic function in dilated cardiomyopathy (DCM). Fourteen idiopathic DCM patients and 16 controls underwent cardiac magnetic resonance imaging, with aortic relative pressure derived using physics-based image processing and a virtual cohort utilized to assess the impact of cardiovascular properties on aortic behaviour. Subjects with reduced left ventricular systolic function had significantly reduced aortic relative pressure, increased aortic stiffness, and significantly delayed time-to-pressure peak duration. From the virtual cohort, aortic stiffness and aortic volumetric size were identified as key determinants of aortic relative pressure. As such, this study shows how advanced flow imaging and aortic hemodynamic evaluation could provide novel insights into the manifestation of DCM, with signs of both altered aortic structure and function derived in DCM using our proposed imaging protocol.

Unlocking the Non-invasive Assessment of Conduit and Reservoir Function in the Aorta

Aortic surgeries in congenital conditions, such as hypoplastic left heart syndrome (HLHS), aim to restore and maintain the conduit and reservoir functions of the aorta. We proposed a method to assess these two functions based on 4D flow MRI, and we applied it to study the aorta in pre-Fontan HLHS. Ten pre-Fontan HLHS patients and six age-matched controls were studied to derive the advective pressure difference and viscous dissipation for conduit function, and pulse wave velocity and elastic modulus for reservoir function. The reconstructed neo-aorta in HLHS subjects achieved a good conduit function at a cost of an impaired reservoir function (69.7% increase of elastic modulus). The native descending HLHS aorta displayed enhanced reservoir (elastic modulus being 18.4% smaller) but impaired conduit function (three-fold increase in peak advection). A non-invasive and comprehensive assessment of aortic conduit and reservoir functions is feasible and has potentially clinical relevance in congenital vascular conditions.

Velocity quantification in 44 healthy volunteers using accelerated multi-VENC 4D flow CMR

BACKGROUND

To evaluate the feasibility of a k-t accelerated multi-VENC 4D phase contrast flow MRI acquisition of the main heart-surrounding vessels, its benefits over a traditional single-VENC acquisition and to present reference flow and velocity values in a large cohort of volunteers.

METHODS

44 healthy volunteers were examined on a 3 T MRI scanner (Ingenia, Philips, Best, The Netherlands). 4D flow measurements were obtained with a FOV including the aorta and the pulmonary arteries. VENC values were set to 40, 100 and 200 cm/s and unfolded based on an MRI signal model. Unfolded multi-VENC data was compared to the single-VENC with VENC 200 cm/s. Flow and velocity quantification was performed in several regions of interest (ROI) placed in the ascending aorta and in the main pulmonary artery. Conservation of mass analysis was performed for single- and multi-VENC datasets. Values for mean and maximal flow velocity and stroke volume were calculated and compared to the literature.

RESULTS

Mean scan time was 13.8 ± 4 min. Differences between stroke volumes between the ascending aorta and the main pulmonary artery were significantly lower in multi-VENC datasets compared to single-VENC datasets (9.6 ± 7.8 mL vs. 25.4 ± 26.4 mL, p < 0.001). This was also true for differences in stroke volume between up- and downstream ROIs in the ascending aorta and pulmonary artery. Values for mean and maximal velocities and stroke volume were in-line with previous studies. To highlight potential clinical applications two exemplary 4D flow measurements in patients with different pathologies are shown and compared to single-VENC datasets.

CONCLUSIONS

k-t accelerated multi-VENC 4D phase contrast flow MRI acquisition of the great vessels is feasible in a clinically acceptable scan duration. It offers improvements over traditional single-VENC 4D flow, expectedly being valuable when vessels with different flow velocities or complex flow phenomena are evaluated.

Validation of non-contrast multiple overlapping thin-slab 4D-flow cardiac magnetic resonance imaging

BACKGROUND

Cardiac magnetic resonance (CMR) flow quantification is typically performed using 2D phase-contrast (PC) imaging of a plane perpendicular to flow. 3D-PC imaging (4D-flow) allows offline quantification anywhere in a thick slab, but is often limited by suboptimal signal, potentially alleviated by contrast enhancement. We developed a non-contrast 4D-flow sequence, which acquires multiple overlapping thin slabs (MOTS) to minimize signal loss, and hypothesized that it could improve image quality, diagnostic accuracy, and aortic flow measurements compared to non-contrast single-slab approach.

METHODS

We prospectively studied 20 patients referred for transesophageal echocardiography (TEE), who underwent CMR (GE, 3 T). 2D-PC images of the aortic valve and three 4D-flow datasets covering the heart were acquired, including single-slab, pre- and post-contrast, and non-contrast MOTS. Each 4D-flow dataset was interpreted blindly for ≥moderate valve disease and compared to TEE. Flow visualization through each valve was scored (0 to 4), and aortic-valve flow measured on each 4D-flow dataset and compared to 2D-PC reference.

RESULTS

Diagnostic quality visualization was achieved with the pre- and post-contrast 4D-flow acquisitions in 25% and 100% valves, respectively (scores 0.9 ± 1.1 and 3.8 ± 0.5), and in 58% with the non-contrast MOTS (1.6 ± 1.1). Accuracy of detection of valve disease was 75%, 92% and 82%, respectively. Aortic flow measurements were possible in 53%, 95% and in 89% patients, respectively. The correlation between pre-contrast single-slab measurements and 2D-PC reference was weak (r = 0.21), but improved with both contrast enhancement (r = 0.71) and with MOTS (r = 0.67).

CONCLUSIONS

Although non-contrast MOTS 4D-flow improves valve function visualization and diagnostic accuracy, a significant proportion of valves cannot be accurately assessed. However, aortic flow measurements using non-contrast MOTS is feasible and reaches similar accuracy to that of contrast-enhanced 4D-flow.

Fully self-gated whole-heart 4D flow imaging from a 5-minute scan

PURPOSE

To develop and validate an acquisition and processing technique that enables fully self-gated 4D flow imaging with whole-heart coverage in a fixed 5-minute scan.

THEORY AND METHODS

The data are acquired continuously using Cartesian sampling and sorted into respiratory and cardiac bins using the self-gating signal. The reconstruction is performed using a recently proposed Bayesian method called ReVEAL4D. ReVEAL4D is validated using data from 8 healthy volunteers and 2 patients and compared with compressed sensing technique, L1-SENSE.

RESULTS

Healthy subjects-Compared with 2D phase-contrast MRI (2D-PC), flow quantification from ReVEAL4D shows no significant bias. In contrast, the peak velocity and peak flow rate for L1-SENSE are significantly underestimated. Compared with traditional parallel MRI-based 4D flow imaging, ReVEAL4D demonstrates small but significant biases in net flow and peak flow rate, with no significant bias in peak velocity. All 3 indices are significantly and more markedly underestimated by L1-SENSE. Patients-Flow quantification from ReVEAL4D agrees well with the 2D-PC reference. In contrast, L1-SENSE markedly underestimated peak velocity.

CONCLUSIONS

The combination of highly accelerated 5-minute Cartesian acquisition, self-gating, and ReVEAL4D enables whole-heart 4D flow imaging with accurate flow quantification.

Compressed-sensing accelerated 4D flow MRI of cerebrospinal fluid dynamics

BACKGROUND

4D flow magnetic resonance imaging (MRI) of CSF can make an important contribution to the understanding of hydrodynamic changes in various neurological diseases but remains limited in clinical application due to long acquisition times. The aim of this study was to evaluate the accuracy of compressed SENSE accelerated MRI measurements of the spinal CSF flow.

METHODS

In 20 healthy subjects 4D flow MRI of the CSF in the cervical spine was acquired using compressed sensitivity encoding [CSE, a combination of compressed sensing and parallel imaging (SENSE) provided by the manufacturer] with acceleration factors between 4 and 10. A conventional scan using SENSE was used as reference. Extracted parameters were peak velocity, absolute net flow, forward flow and backward flow. Bland-Altman analysis was performed to determine the scan-rescan reproducibility and the agreement between SENSE and compressed SENSE. Additionally, a time accumulated flow error was calculated. In one additional subject flow of the spinal canal at the level of the entire spinal cord was assessed.

RESULTS

Averaged acquisition times were 10:21 min (SENSE), 9:31 min (CSE4), 6:25 min (CSE6), 4:53 min (CSE8) and 3:51 min (CSE10). Acquisition of the CSF flow surrounding the entire spinal cord took 14:40 min. Bland-Altman analysis showed good agreement for peak velocity, but slight overestimations for absolute net flow, forward flow and backward flow (< 1 ml/min) in CSE4-8. Results of the accumulated flow error were similar for CSE4 to CSE8.

CONCLUSION

A quantitative analysis of acceleration factors CSE4-10 showed that CSE with an acceleration factor up to 6 is feasible. This allows a scan time reduction of 40% and enables the acquisition and analysis of the CSF flow dynamics surrounding the entire spinal cord within a clinically acceptable scan time.

4D Flow with MRI

Magnetic resonance imaging (MRI) has become an important tool for the clinical evaluation of patients with cardiac and vascular diseases. Since its introduction in the late 1980s, quantitative flow imaging with MRI has become a routine part of standard-of-care cardiothoracic and vascular MRI for the assessment of pathological changes in blood flow in patients with cardiovascular disease. More recently, time-resolved flow imaging with velocity encoding along all three flow directions and three-dimensional (3D) anatomic coverage (4D flow MRI) has been developed and applied to enable comprehensive 3D visualization and quantification of hemodynamics throughout the human circulatory system. This article provides an overview of the use of 4D flow applications in different cardiac and vascular regions in the human circulatory system, with a focus on using 4D flow MRI in cardiothoracic and cerebrovascular diseases.

Pseudo-spiral sampling and compressed sensing reconstruction provides flexibility of temporal resolution in accelerated aortic 4D flow MRI: A comparison with k-t principal component analysis

INTRODUCTION

Time-resolved three-dimensional phase contrast MRI (4D flow) of aortic blood flow requires acceleration to reduce scan time. Two established techniques for highly accelerated 4D flow MRI are k-t principal component analysis (k-t PCA) and compressed sensing (CS), which employ either regular or random k-space undersampling. The goal of this study was to gain insights into the quantitative differences between k-t PCA- and CS-derived aortic blood flow, especially for high temporal resolution CS 4D flow MRI.

METHODS

The scan protocol consisted of both k-t PCA and CS accelerated 4D flow MRI, as well as a 2D flow reference scan through the ascending aorta acquired in 15 subjects. 4D flow scans were accelerated with factor R = 8. For CS accelerated scans, we used a pseudo-spiral Cartesian sampling scheme, which could additionally be reconstructed at higher temporal resolution, resulting in R = 13. 4D flow data were compared with the 2D flow scan in terms of flow, peak flow and stroke volume. A 3D peak systolic voxel-wise velocity and wall shear stress (WSS) comparison between k-t PCA and CS 4D flow was also performed.

RESULTS

The mean difference in flow/peak flow/stroke volume between the 2D flow scan and the 4D flow CS with R = 8 and R = 13 was 4.2%/9.1%/3.0% and 5.3%/7.1%/1.9%, respectively, whereas for k-t PCA with R = 8 the difference was 9.7%/25.8%/10.4%. In the voxel-by-voxel 4D flow comparison we found 13.6% and 3.5% lower velocity and WSS values of k-t PCA compared with CS with R = 8, and 15.9% and 5.5% lower velocity and WSS values of k-t PCA compared with CS with R = 13.

CONCLUSION

Pseudo-spiral accelerated 4D flow acquisitions in combination with CS reconstruction provides a flexible choice of temporal resolution. We showed that our proposed strategy achieves better agreement in flow values with 2D reference scans compared with using k-t PCA accelerated acquisitions.

Improved acceleration of phase-contrast flow imaging with magnitude difference regularization

PURPOSE

To develop a regularized image reconstruction algorithm for improved scan acceleration of phase-contrast (PC) flow MRI.

METHODS

Based on the magnitude similarity between bipolar-encoded k-space data, magnitude-difference regularization was incorporated into the conventional compressed sensing (CS) reconstruction. The gradient of the magnitude regularization was derived so the reconstruction problem can be solved using non-linear conjugate gradient with backtracking line search. Phase contrast flow data obtained in the peripheral arteries of healthy and patient subjects were retrospectively undersampled for testing the proposed reconstruction method. Three-dimensional velocity-encoded PC flow MRI was performed with prospective 4-fold undersampling for measuring arotic flow velocity in a healthy volunteer.

RESULTS

In the femoral arteries of healthy volunteers, the root-mean-square (RMS) errors of mean velocities were 0.56 ± 0.09 cm/s with CS-only reconstruction and 0.46 ± 0.08 cm/s with addition of magnitude regularization for three-fold acceleration; 1.34 ± 0.17 cm/s (CS only) and 1.08 ± 0.15 cm/s (magnitude regularized) for four-fold acceleration. In the iliac arteries of the patient, the RMS errors of mean velocities were 0.72 ± 0.12 cm/s and 0.56 ± 0.10 for three-fold acceleration, and 1.75 ± 0.21 and 1.24 ± 0.19 cm/s for four-fold acceleration (in the order of CS-only and magnitude regularized reconstructions). In the popliteal arteries, the RMS errors were 0.61 ± 0.10 cm/s and 0.42 ± 0.11 for three-fold acceleration, and 1.41 ± 0.19 and 1.12 ± 0.17 cm/s for four-fold acceleration. The maximum through-plane mean flow velocities were measured as 63.2 cm/s and 84.5 cm/s in ascending and descending aortas, respectively.

CONCLUSION

The addition of magnitude-difference regularization into conventional CS reconstruction improves the accuracy of image reconstruction using highly undersampled phase-contrast flow MR data.

Highly accelerated 4D flow cardiovascular magnetic resonance using a pseudo-spiral Cartesian acquisition and compressed sensing reconstruction for carotid flow and wall shear stress

BACKGROUND

4D flow cardiovascular magnetic resonance (CMR) enables visualization of complex blood flow and quantification of biomarkers for vessel wall disease, such as wall shear stress (WSS). Because of the inherently long acquisition times, many efforts have been made to accelerate 4D flow acquisitions, however, no detailed analysis has been made on the effect of Cartesian compressed sensing accelerated 4D flow CMR at different undersampling rates on quantitative flow parameters and WSS.

METHODS

We implemented a retrospectively triggered 4D flow CMR acquisition with pseudo-spiral Cartesian k-space filling, which results in incoherent undersampling of k-t space. Additionally, this strategy leads to small jumps in k-space thereby minimizing eddy current related artifacts. The pseudo-spirals were rotated in a tiny golden-angle fashion, which provides optimal incoherence and a variable density sampling pattern with a fully sampled center. We evaluated this 4D flow protocol in a carotid flow phantom with accelerations of R = 2-20, as well as in carotids of 7 healthy subjects (27 ± 2 years, 4 male) for R = 10-30. Fully sampled 2D flow CMR served as a flow reference. Arteries were manually segmented and registered to enable voxel-wise comparisons of both velocity and WSS using a Bland-Altman analysis.

RESULTS

Magnitude images, velocity images, and pathline reconstructions from phantom and in vivo scans were similar for all accelerations. For the phantom data, mean differences at peak systole for the entire vessel volume in comparison to R = 2 ranged from - 2.3 to - 5.3% (WSS) and - 2.4 to - 2.2% (velocity) for acceleration factors R = 4-20. For the in vivo data, mean differences for the entire vessel volume at peak systole in comparison to R = 10 were - 9.9, - 13.4, and - 16.9% (WSS) and - 8.4, - 10.8, and - 14.0% (velocity), for R = 20, 25, and 30, respectively. Compared to single slice 2D flow CMR acquisitions, peak systolic flow rates of the phantom showed no differences, whereas peak systolic flow rates in the carotid artery in vivo became increasingly underestimated with increasing acceleration.

CONCLUSION

Acquisition of 4D flow CMR of the carotid arteries can be highly accelerated by pseudo-spiral k-space sampling and compressed sensing reconstruction, with consistent data quality facilitating velocity pathline reconstructions, as well as quantitative flow rate and WSS estimations. At an acceleration factor of R = 20 the underestimation of peak velocity and peak WSS was acceptable (< 10%) in comparison to an R = 10 accelerated 4D flow CMR reference scan. Peak flow rates were underestimated in comparison with 2D flow CMR and decreased systematically with higher acceleration factors.

Accelerated aortic 4D flow cardiovascular magnetic resonance using compressed sensing: applicability, validation and clinical integration

BACKGROUND

Three-dimensional time-resolved phase-contrast cardiovascular magnetic resonance (4D flow CMR) enables the quantification and visualisation of blood flow, but its clinical applicability remains hampered by its long scan time. The aim of this study was to evaluate the use of compressed sensing (CS) with on-line reconstruction to accelerate the acquisition and reconstruction of 4D flow CMR of the thoracic aorta.

METHODS

4D flow CMR of the thoracic aorta was acquired in 20 healthy subjects using CS with acceleration factors ranging from 4 to 10. As a reference, conventional parallel imaging (SENSE) with acceleration factor 2 was used. Flow curves, net flows, peak flows and peak velocities were extracted from six contours along the aorta. To measure internal data consistency, a quantitative particle trace analysis was performed. Additionally, scan-rescan, inter- and intraobserver reproducibility were assessed. Subsequently, 4D flow CMR with CS factor 6 was acquired in 3 patients with differing aortopathies. The flow patterns resulting from particle trace visualisation were qualitatively analysed.

RESULTS

All collected data were successfully acquired and reconstructed on-line. The average acquisition time including respiratory navigator efficiency with CS factor 6 was 5:02 ± 2:23 min while reconstruction took approximately 9 min. For CS factors of 8 or less, mean differences in net flow, peak flow and peak velocity as compared to SENSE were below 2.2 ± 7.8 ml/cycle, 4.6 ± 25.2 ml/s and - 7.9 ± 13.0 cm/s, respectively. For a CS factor of 10 differences reached 5.4 ± 8.0 ml/cycle, 14.4 ± 28.3 ml/s and - 4.0 ± 12.2 cm/s. Scan-rescan analysis yielded mean differences in net flow of - 0.7 ± 4.9 ml/cycle for SENSE and - 0.2 ± 8.5 ml/cycle for CS factor of 6.

CONCLUSIONS

A six- to eightfold acceleration of 4D flow CMR using CS is feasible. Up to a CS acceleration rate of 6, no statistically significant differences in measured flow parameters could be observed with respect to the reference technique. Acquisitions in patients with aortopathies confirm the potential to integrate the proposed method in a clinical routine setting, whereby its main benefits are scan-time savings and direct on-line reconstruction.

Respiration Dependency of Caval Blood Flow in Patients with Fontan Circulation: Quantification Using 5D Flow MRI

PURPOSE

To measure respiration-dependent blood flow in the total cavopulmonary connection (TCPC) of patients with Fontan circulation by using free-running, fully self-gated five-dimensional (5D) flow MRI.

MATERIALS AND METHODS

From July to November 2018, 10 volunteers (six female volunteers, mean age, 25.1 years ± 4.4 [standard deviation]) and six patients with Fontan circulation (two female patients, mean age, 19.7 years ± 7.5) with a TCPC were examined by using a cardiac- and respiration-resolved three-directional and three-dimensional phase-contrast MRI sequence (hereafter, 5D flow MRI). This prospective study was conducted with approval of the local ethics committee, and written informed consent was obtained from all participants and/or their representative. 5D flow data were acquired during free breathing. Data were reconstructed into 15-20 heart phases and four respiratory phases: end-expiration, inspiration, end-inspiration, and expiration. Respiration-dependent stroke volumes (SVs) and particle traces were analyzed from the caval circulation of volunteers and patients with Fontan circulation. Statistical analysis was performed by using parametric tests and scatterplots.

RESULTS

The respiration dependency of caval blood flow was evaluated in all participants and was significantly elevated in patients with Fontan circulation as compared with volunteers. In patients, SV in the inferior vena cava (IVC) showed variations of 120% between inspiration and expiration (P = .002). The flow distribution in the IVC and superior vena cava among the four respiratory phases was differentiated by 20% (range, 9%-30%) and 4% (range, 0%-13%), respectively.

CONCLUSION

Hemodynamic parameters (volume flow and blood flow distribution) throughout the cardiac and respiratory cycle can be measured using a single scan, potentially providing further insights into the Fontan circulation.© RSNA, 2019Supplemental material is available for this article.

On the limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy

PURPOSE

To investigate limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy (TKE).

METHODS

To assess the validity of partial Fourier reconstruction of TKE and phase images, computational fluid dynamics data of mean and turbulent velocities in a stenotic U-bend phantom was used. Partial Fourier acquisition with 75% k-space coverage was simulated and TKE data were reconstructed using zero-filling, homodyne reconstruction, and the method of projections onto convex sets (POCS). Results were compared to data from fully sampled k-space and 75% symmetric sampling. In addition, compressed sensing (CS) reconstruction was compared for a standard variable density sampling pattern and a variable density sampling pattern combined with 75% partial Fourier. For illustration purposes, in vivo examples of velocity magnitude and TKE maps of aortic flow reconstructed with the different methods are provided.

RESULTS

In accordance with theory, partial Fourier reconstruction of TKE maps from phase-contrast data results in artifacts relative to fully sampled data. It is demonstrated that neither homodyne reconstruction nor POCS can improve reconstruction of TKE data with respect to zero-filling reconstruction when compared to ground-truth (RMS error: 4.70%, 4.34%, and 2.45% for homodyne, POCS, and zero-filling reconstruction of in vivo data, respectively). CS reconstruction from data acquired with partial Fourier did not recover the resolution loss incurred by partial Fourier sampling.

CONCLUSION

Partial Fourier reconstruction of TKE maps from phase-contrast data does not yield a benefit over zero-filling reconstruction. In consequence, symmetric sampling is preferred over partial Fourier acquisition for a given number of phase-encodes in phase-contrast MRI.

A Bayesian approach for 4D flow imaging of aortic valve in a single breath-hold

PURPOSE

To develop and validate a data processing technique that allows phase-contrast MRI-based 4D flow imaging of the aortic valve in a single breath-hold.

THEORY AND METHODS

To regularize the ill-posed inverse problem, we extend a recently proposed 2D phase-contrast MRI method to 4D flow imaging. Adopting an empirical Bayes approach, spatial and temporal redundancies are exploited via sparsity in the wavelet domain, and the voxel-wise magnitude and phase structure across encodings is captured in a conditional mixture prior that applies regularizing constraints based on the presence of flow. We validate the proposed technique using data from a mechanical flow phantom and five healthy volunteers.

RESULTS

The flow parameters derived from the proposed technique are in good agreement with those derived from reference datasets for both in vivo and mechanical flow experiments at accelerations rates as high as R = 27. Additionally, the proposed technique outperforms kt SPARSE-SENSE and a method that exploits spatio-temporal sparsity but does not utilize signal structure across encodings.

CONCLUSIONS

Using the proposed technique, it is feasible to highly accelerate 4D flow acquisition and thus enable aortic valve imaging within a single breath-hold.

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms

OBJECTIVES

To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging.

MATERIALS AND METHODS

Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared.

RESULTS

Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: -3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, -1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6).

CONCLUSION

The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

Exploring kinetic energy as a new marker of cardiac function in the single ventricle circulation

Ventricular volumetric ejection fraction (VV EF) is often normal in patients with single ventricle circulations despite them experiencing symptoms related to circulatory failure. We sought to determine if kinetic energy (KE) could be a better marker of ventricular performance. KE was prospectively quantified using four-dimensional flow MRI in 41 patients with a single ventricle circulation (aged 0.5-28 yr) and compared with 43 healthy volunteers (aged 1.5-62 yr) and 14 patients with left ventricular (LV) dysfunction (aged 28-79 yr). Intraventricular end-diastolic blood was tracked through systole and divided into ejected and residual blood components. Two ejection fraction (EF) metrics were devised based on the KE of the ejected component over the total of both the ejected and residual components using 1) instantaneous peak KE to assess KE EF or 2) summating individual peak particle energy (PE) to assess PE EF. KE EF and PE EF had a smaller range than VV EF in healthy subjects (97.9 ± 0.8 vs. 97.3 ± 0.8 vs. 60.1 ± 5.2%). LV dysfunction caused a fall in KE EF ( P = 0.01) and PE EF ( P = 0.0001). VV EF in healthy LVs and single ventricle hearts was equivalent; however, KE EF and PE EF were lower ( P < 0.001) with a wider range indicating a spectrum of severity. Those reporting the greatest symptomatic impairment (New York Heart Association II) had lower PE EF than asymptomatic subjects ( P = 0.0067). KE metrics are markers of healthy cardiac function. PE EF may be useful in grading dysfunction. NEW & NOTEWORTHY Kinetic energy (KE) represents the useful work of the heart in ejecting blood. This article details the utilization of KE indexes to assess cardiac function in health and a variety of pathophysiological conditions. KE ejection fraction and particle energy ejection fraction (PE EF) showed a narrow range in health and a lower wider range in disease representing a spectrum of severity. PE EF was altered by functional status potentially offering the opportunity to grade dysfunction.

Influence of principal component analysis acceleration factor on velocity measurement in 2D and 4D PC-MRI

OBJECTIVE

The objective of the study was to determine how to optimize 2D and 4D phase-contrast magnetic resonance imaging (PC-MRI) acquisitions to acquire flow velocities in millimetric vessels. In particular, we search for the best compromise between acquisition time and accuracy and assess the influence of the principal component analysis (PCA).

MATERIALS AND METHODS

2D and 4D PC-MRI measurements are conducted within two in vitro vessel phantoms: a Y-bifurcation phantom, the branches of which range from 2 to 5 mm in diameter, and a physiological subject-specific phantom of the carotid bifurcation. The same sequences are applied in vivo in carotid vasculature.

RESULTS

For a vessel oriented in the axial direction, both 2D and axial 4D PC-MRI provided accuracy measurements regardless of the k-t PCA factor, while the acquisition time is reduced by a factor 6 for k-t PCA maximum value. The in vivo measurements show that the proposed sequences are adequate to acquire 2D and 4D velocity fields in millimetric vessels and with clinically realistic time durations.

CONCLUSION

The study shows the feasibility of conducting fast, high-resolution PC-MRI flow measurements in millimetric vessels and that it is worth maximizing the k-t PCA factor to reduce the acquisition time in the case of 2D and 4D axial acquisitions.

Comprehensive Multi-Dimensional MRI for the Simultaneous Assessment of Cardiopulmonary Anatomy and Physiology

Diagnostic testing often assesses the cardiovascular or respiratory systems in isolation, ignoring the major pathophysiologic interactions between the systems in many diseases. When both systems are assessed currently, multiple modalities are utilized in costly fashion with burdensome logistics and decreased accessibility. Thus, we have developed a new acquisition and reconstruction paradigm using the flexibility of MRI to enable a comprehensive exam from a single 5-15 min scan. We constructed a compressive-sensing approach to pseudo-randomly acquire highly subsampled, multi-dimensionally-encoded and time-stamped data from which we reconstruct volumetric cardiac and respiratory motion phases, contrast-agent dynamics, and blood flow velocity fields. The proposed method, named XD flow, is demonstrated for (a) evaluating congenital heart disease, where the impact of bulk motion is reduced in a non-sedated neonatal patient and (b) where the observation of the impact of respiration on flow is necessary for diagnostics; (c) cardiopulmonary imaging, where cardiovascular flow, function, and anatomy information is needed along with pulmonary perfusion quantification; and in (d) renal function imaging, where blood velocities and glomerular filtration rates are simultaneously measured, which highlights the generality of the technique. XD flow has the ability to improve quantification and to provide additional data for patient diagnosis for comprehensive evaluations.

Accelerated phase contrast flow imaging with direct complex difference reconstruction

PURPOSE

To propose and evaluate a new model-based reconstruction method for highly accelerated phase-contrast magnetic resonance imaging (PC-MRI) with sparse sampling.

THEORY AND METHODS

This work presents a new constrained reconstruction method based on low-rank and sparsity constraints to accelerate PC-MRI. More specifically, we formulate the image reconstruction problem into separate reconstructions of flow-reference image sequence and complex differences. We then utilize the joint partial separability and sparsity constraints to enable high quality reconstruction from highly undersampled (k,t)-space data. We further integrate the proposed method with ESPIRiT based parallel imaging model to effectively handle multichannel acquisition.

RESULTS

The proposed method was evaluated with in vivo data acquired from both 2D and 3D PC flow imaging experiments, and compared with several state-of-the-art methods. Experimental results demonstrate that the proposed method leads to more accurate velocity reconstruction from highly undersampled (k,t)-space data, and particularly superior capability of capturing the peak velocity of blood flow. In terms of flow visualization, blood flow patterns obtained from the proposed reconstruction also exhibit better agreement with those obtained from the fully sampled reference.

CONCLUSION

The proposed method achieves improved accuracy over several state-of-the-art methods for velocity reconstruction with highly accelerated (k,t)-space data. Magn Reson Med 77:1036-1048, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

MRI-based biomechanical parameters for carotid artery plaque vulnerability assessment

Carotid atherosclerotic plaques are a major cause of ischaemic stroke. The biomechanical environment to which the arterial wall and plaque is subjected to plays an important role in the initiation, progression and rupture of carotid plaques. MRI is frequently used to characterize the morphology of a carotid plaque, but new developments in MRI enable more functional assessment of carotid plaques. In this review, MRI based biomechanical parameters are evaluated on their current status, clinical applicability, and future developments. Blood flow related biomechanical parameters, including endothelial wall shear stress and oscillatory shear index, have been shown to be related to plaque formation. Deriving these parameters directly from MRI flow measurements is feasible and has great potential for future carotid plaque development prediction. Blood pressure induced stresses in a plaque may exceed the tissue strength, potentially leading to plaque rupture. Multi-contrast MRI based stress calculations in combination with tissue strength assessment based on MRI inflammation imaging may provide a plaque stress-strength balance that can be used to assess the plaque rupture risk potential. Direct plaque strain analysis based on dynamic MRI is already able to identify local plaque displacement during the cardiac cycle. However, clinical evidence linking MRI strain to plaque vulnerability is still lacking. MRI based biomechanical parameters may lead to improved assessment of carotid plaque development and rupture risk. However, better MRI systems and faster sequences are required to improve the spatial and temporal resolution, as well as increase the image contrast and signal-to-noise ratio.

Age-related changes in intraventricular kinetic energy: a physiological or pathological adaptation?

Measuring intracardiac kinetic energy using four-dimensionl flow cardiac magnetic resonance provides important information on the decline in the early diastolic kinetic energy of blood with aging. The decline is comparable with that seen in those with heart failure and may be a marker of cardiac function.

Aging has important deleterious effects on the cardiovascular system. We sought to compare intraventricular kinetic energy (KE) in healthy subjects of varying ages with subjects with ventricular dysfunction to understand if changes in energetic momentum may predispose individuals to heart failure. Four-dimensional flow MRI was acquired in 35 healthy subjects (age: 1–67 yr) and 10 patients with left ventricular (LV) dysfunction (age: 28–79 yr). Healthy subjects were divided into age quartiles (1st quartile: <16 yr, 2nd quartile: 17–32 yr, 3rd quartile: 33–48 yr, and 4th quartile: 49–64 yr). KE was measured in the LV throughout the cardiac cycle and indexed to ventricular volume. In healthy subjects, two large peaks corresponding to systole and early diastole occurred during the cardiac cycle. A third smaller peak was seen during late diastole in eight adults. Systolic KE (P = 0.182) and ejection fraction (P = 0.921) were preserved through all age groups. Older adults showed a lower early peak diastolic KE compared with children (P < 0.0001) and young adults (P = 0.025). Subjects with LV dysfunction had reduced ejection fraction (P < 0.001) and compared with older healthy adults exhibited a similar early peak diastolic KE (P = 0.142) but with the addition of an elevated KE in diastasis (P = 0.029). In healthy individuals, peak diastolic KE progressively decreases with age, whereas systolic peaks remain constant. Peak diastolic KE in the oldest subjects is comparable to those with LV dysfunction. Unique age-related changes in ventricular diastolic energetics might be physiological or herald subclinical pathology.

A Bayesian model for highly accelerated phase-contrast MRI

PURPOSE

Phase-contrast magnetic resonance imaging is a noninvasive tool to assess cardiovascular disease by quantifying blood flow; however, low data acquisition efficiency limits the spatial and temporal resolutions, real-time application, and extensions to four-dimensional flow imaging in clinical settings. We propose a new data processing approach called Reconstructing Velocity Encoded MRI with Approximate message passing aLgorithms (ReVEAL) that accelerates the acquisition by exploiting data structure unique to phase-contrast magnetic resonance imaging.

THEORY AND METHODS

The proposed approach models physical correlations across space, time, and velocity encodings. The proposed Bayesian approach exploits the relationships in both magnitude and phase among velocity encodings. A fast iterative recovery algorithm is introduced based on message passing. For validation, prospectively undersampled data are processed from a pulsatile flow phantom and five healthy volunteers.

RESULTS

The proposed approach is in good agreement, quantified by peak velocity and stroke volume (SV), with reference data for acceleration rates R≤10. For SV, Pearson r≥0.99 for phantom imaging (n = 24) and r≥0.96 for prospectively accelerated in vivo imaging (n = 10) for R≤10.

CONCLUSION

The proposed approach enables accurate quantification of blood flow from highly undersampled data. The technique is extensible to four-dimensional flow imaging, where higher acceleration may be possible due to additional redundancy. Magn Reson Med 76:689-701, 2016. © 2015 Wiley Periodicals, Inc.

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.

Retrospectively gated intracardiac 4D flow MRI using spiral trajectories

PURPOSE

To develop and evaluate retrospectively gated spiral readout four-dimensional (4D) flow MRI for intracardiac flow analysis.

METHODS

Retrospectively gated spiral 4D flow MRI was implemented on a 1.5-tesla scanner. The spiral sequence was compared against conventional Cartesian 4D flow (SENSE [sensitivity encoding] 2) in seven healthy volunteers and three patients (only spiral). In addition to comparing flow values, linear regression was used to assess internal consistency of aortic versus pulmonary net volume flows and left ventricular inflow versus outflow using quantitative pathlines analysis.

RESULTS

Total scan time with spiral 4D flow was 44% ± 6% of the Cartesian counterpart (13 ± 3 vs. 31 ± 7 min). Aortic versus pulmonary flow correlated strongly for the spiral sequence (P < 0.05, slope = 1.03, R(2) = 0.88, N = 10), whereas the linear relationship for the Cartesian sequence was not significant (P = 0.06, N = 7). Pathlines analysis indicated good data quality for the spiral (P < 0.05, slope = 1.02, R(2) = 0.90, N = 10) and Cartesian sequence (P < 0.05, slope = 1.10, R(2) = 0.93, N = 7). Spiral and Cartesian peak flow rate (P < 0.05, slope = 0.96, R(2) = 0.72, N = 14), peak velocity (P < 0.05, slope = 1.00, R(2) = 0.81, N = 14), and pathlines flow components (P < 0.05, slope = 1.04, R(2) = 0.87, N = 28) correlated well.

CONCLUSION

Retrospectively gated spiral 4D flow MRI permits more than two-fold reduction in scan time compared to conventional Cartesian 4D flow MRI, while maintaining similar data quality.

Accelerated three-dimensional cine phase contrast imaging using randomly undersampled echo planar imaging with compressed sensing reconstruction

The aim of this study was to implement and evaluate an accelerated three-dimensional (3D) cine phase contrast MRI sequence by combining a randomly sampled 3D k-space acquisition sequence with an echo planar imaging (EPI) readout. An accelerated 3D cine phase contrast MRI sequence was implemented by combining EPI readout with randomly undersampled 3D k-space data suitable for compressed sensing (CS) reconstruction. The undersampled data were then reconstructed using low-dimensional structural self-learning and thresholding (LOST). 3D phase contrast MRI was acquired in 11 healthy adults using an overall acceleration of 7 (EPI factor of 3 and CS rate of 3). For comparison, a single two-dimensional (2D) cine phase contrast scan was also performed with sensitivity encoding (SENSE) rate 2 and approximately at the level of the pulmonary artery bifurcation. The stroke volume and mean velocity in both the ascending and descending aorta were measured and compared between two sequences using Bland-Altman plots. An average scan time of 3 min and 30 s, corresponding to an acceleration rate of 7, was achieved for 3D cine phase contrast scan with one direction flow encoding, voxel size of 2 × 2 × 3 mm(3) , foot-head coverage of 6 cm and temporal resolution of 30 ms. The mean velocity and stroke volume in both the ascending and descending aorta were statistically equivalent between the proposed 3D sequence and the standard 2D cine phase contrast sequence. The combination of EPI with a randomly undersampled 3D k-space sampling sequence using LOST reconstruction allows a seven-fold reduction in scan time of 3D cine phase contrast MRI without compromising blood flow quantification.

4D flow imaging: current status to future clinical applications

4D flow MRI permits a comprehensive in-vivo assessment of three-directional blood flow within 3-dimensional vascular structures throughout the cardiac cycle. Given the large coverage permitted from a 4D flow acquisition, the distribution of vessel wall and flow parameters along an entire vessel of interest can thus be derived from a single measurement without being dependent on multiple predefined 2D acquisitions. In addition to qualitative 3D visualizations of complex cardiac and vascular flow patterns, quantitative flow analysis can be performed and is complemented by the ability to compute sophisticated hemodynamic parameters, such as wall shear stress or 3D pressure difference maps. These metrics can provide information previously unavailable with conventional modalities regarding the impact of cardiovascular disease or therapy on global and regional changes in hemodynamics. This review provides an introduction to the methodological aspects of 4D flow MRI to assess vascular hemodynamics and describes its potential for the assessment and understanding of altered hemodynamics in the presence of cardiovascular disease.

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.

Highly accelerated aortic 4D flow MR imaging with variable-density random undersampling

PURPOSE

To investigate an effective time-resolved variable-density random undersampling scheme combined with an efficient parallel image reconstruction method for highly accelerated aortic 4D flow MR imaging with high reconstruction accuracy.

MATERIALS AND METHODS

Variable-density Poisson-disk sampling (vPDS) was applied in both the phase-slice encoding plane and the temporal domain to accelerate the time-resolved 3D Cartesian acquisition of flow imaging. In order to generate an improved initial solution for the iterative self-consistent parallel imaging method (SPIRiT), a sample-selective view sharing reconstruction for time-resolved random undersampling (STIRRUP) was introduced. The performance of different undersampling and image reconstruction schemes were evaluated by retrospectively applying those to fully sampled data sets obtained from three healthy subjects and a flow phantom.

RESULTS

Undersampling pattern based on the combination of time-resolved vPDS, the temporal sharing scheme STIRRUP, and parallel imaging SPIRiT, were able to achieve 6-fold accelerated 4D flow MRI with high accuracy using a small number of coils (N=5). The normalized root mean square error between aorta flow waveforms obtained with the acceleration method and the fully sampled data in three healthy subjects was 0.04±0.02, and the difference in peak-systolic mean velocity was -0.29±2.56cm/s.

CONCLUSION

Qualitative and quantitative evaluation of our preliminary results demonstrate that time-resolved variable-density random sampling is efficient for highly accelerating 4D flow imaging while maintaining image reconstruction accuracy.

Estimation of the CSA-ODF using Bayesian compressed sensing of multi-shell HARDI

PURPOSE

Diffusion MRI provides important information about the brain white matter structures and has opened new avenues for neuroscience and translational research. However, acquisition time needed for advanced applications can still be a challenge in clinical settings. There is consequently a need to accelerate diffusion MRI acquisitions.

METHODS

A multi-task Bayesian compressive sensing (MT-BCS) framework is proposed to directly estimate the constant solid angle orientation distribution function (CSA-ODF) from under-sampled (i.e., accelerated image acquisition) multi-shell high angular resolution diffusion imaging (HARDI) datasets, and accurately recover HARDI data at higher resolution in q-space. The proposed MT-BCS approach exploits the spatial redundancy of the data by modeling the statistical relationships within groups (clusters) of diffusion signal. This framework also provides uncertainty estimates of the computed CSA-ODF and diffusion signal, directly computed from the compressive measurements. Experiments validating the proposed framework are performed using realistic multi-shell synthetic images and in vivo multi-shell high angular resolution HARDI datasets.

RESULTS

Results indicate a practical reduction in the number of required diffusion volumes (q-space samples) by at least a factor of four to estimate the CSA-ODF from multi-shell data.

CONCLUSION

This work presents, for the first time, a multi-task Bayesian compressive sensing approach to simultaneously estimate the full posterior of the CSA-ODF and diffusion-weighted volumes from multi-shell HARDI acquisitions. It demonstrates improvement of the quality of acquired datasets by means of CS de-noising, and accurate estimation of the CSA-ODF, as well as enables a reduction in the acquisition time by a factor of two to four, especially when "staggered" q-space sampling schemes are used. The proposed MT-BCS framework can naturally be combined with parallel MR imaging to further accelerate HARDI acquisitions.

Radial k-t SPIRiT: autocalibrated parallel imaging for generalized phase-contrast MRI

PURPOSE

To extend SPIRiT to additionally exploit temporal correlations for highly accelerated generalized phase-contrast MRI and to compare the performance of the proposed radial k-t SPIRiT method relative to frame-by-frame SPIRiT and radial k-t GRAPPA reconstruction for velocity and turbulence mapping in the aortic arch.

THEORY AND METHODS

Free-breathing navigator-gated two-dimensional radial cine imaging with three-directional multi-point velocity encoding was implemented and fully sampled data were obtained in the aortic arch of healthy volunteers. Velocities were encoded with three different first gradient moments per axis to permit quantification of mean velocity and turbulent kinetic energy. Velocity and turbulent kinetic energy maps from up to 14-fold undersampled data were compared for k-t SPIRiT, frame-by-frame SPIRiT, and k-t GRAPPA relative to the fully sampled reference.

RESULTS

Using k-t SPIRiT, improvements in magnitude and velocity reconstruction accuracy were found. Temporally resolved magnitude profiles revealed a reduction in spatial blurring with k-t SPIRiT compared with frame-by-frame SPIRiT and k-t GRAPPA for all velocity encodings, leading to improved estimates of turbulent kinetic energy.

CONCLUSION

k-t SPIRiT offers improved reconstruction accuracy at high radial undersampling factors and hence facilitates the use of generalized phase-contrast MRI for routine use.

Free-breathing phase contrast MRI with near 100% respiratory navigator efficiency using k-space-dependent respiratory gating

PURPOSE

To investigate the efficacy of a novel respiratory motion scheme, where only the center of k-space is gated using respiratory navigators, versus a fully respiratory-gated acquisition for three-dimensional flow imaging.

METHODS

Three-dimensional flow images were acquired axially using a gradient echo sequence in a volume, covering the ascending and descending aorta, and the pulmonary artery bifurcation in 12 healthy subjects (33.2 ± 15.8 years; five men). For respiratory motion compensation, two gating and tracking strategies were used with a 7-mm gating window: (1) All of k-space acquired within the gating window (fully gated) and (2) central k-space acquired within the gating window, and the remainder of k-space acquired without any gating (center gated). Each scan was repeated twice. Stroke volume, mean flow, peak velocity, and signal-to-noise-ratio measurements were performed both on the ascending and on the descending aorta for all acquisitions, which were compared using a linear mixed-effects model and Bland-Altman analysis.

RESULTS

There were no statistical differences between the fully gated and the center-gated strategies for the quantification of stroke volume, peak velocity, and mean flow, as well as the signal-to-noise-ratio measurements. Furthermore, the proposed center-gated strategy had significantly shorter acquisition time compared to the fully gated strategy (13:19 ± 3:02 vs. 19:35 ± 5:02, P < 0.001).

CONCLUSIONS

The proposed novel center-gated strategy for three-dimensional flow MRI allows for markedly shorter acquisition time without any systematic variation in quantitative flow measurements in this small group of healthy volunteers.

Mapping mean and fluctuating velocities by Bayesian multipoint MR velocity encoding-validation against 3D particle tracking velocimetry

PURPOSE

To validate Bayesian multipoint MR velocity encoding against particle tracking velocimetry for measuring velocity vector fields and fluctuating velocities in a realistic aortic model.

METHODS

An elastic cast of a human aortic arch equipped with an 80 or 64% stenotic section was driven by a pulsatile pump. Peak velocities and peak turbulent kinetic energies of more than 3 m/s and 1000 J/m(3) could be generated. Velocity vector fields and fluctuating velocities were assessed using Bayesian multipoint MR velocity encoding with varying numbers of velocity encoding points and particle tracking velocimetry in the ascending aorta.

RESULTS

Velocities and turbulent kinetic energies measured with 5-fold k-t undersampled 10-point MR velocity encoding and particle tracking velocimetry were found to reveal good correlation with mean differences of -4.8 ± 13.3 cm/s and r(2) = 0.98 for velocities and -21.8 ± 53.9 J/m(3) and r(2) = 0.98 for turbulent kinetic energies, respectively. Three-dimensional velocity patterns of fast flow downstream of the stenoses and regions of elevated velocity fluctuations were found to agree well.

CONCLUSION

Accelerated Bayesian multipoint MR velocity encoding has been demonstrated to be accurate for assessing mean and fluctuating velocities against the reference standard particle tracking velocimetry. The MR method holds considerable potential to map velocity vector fields and turbulent kinetic energies in clinically feasible exam times of <15 min.