Robust 4D flow denoising using divergence-free wavelet transform

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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

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

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

BACKGROUND

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

OBJECTIVE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

BACKGROUND

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

OBJECTIVE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Accelerated slice encoding for metal artifact correction

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

Breathheld autocalibrated phase-contrast imaging

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Clinical evaluation of aortic coarctation with 4D flow MR imaging

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

Improved pediatric MR imaging with compressed sensing

PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSION

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

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

INTRODUCTION

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

Evaluation of intracranial stenoses and aneurysms with accelerated 4D flow

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

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

SUMMARY

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