High Spatiotemporal Resolution 4D Flow MRI of Intracranial Aneurysms at 7T in 10 Minutes

L.M. Gottwald; J. Töger; K. Markenroth Bloch; E.S. Peper; B.F. Coolen; G.J. Strijkers; P. van Ooij; A.J. Nederveen

doi:10.3174/ajnr.A6603

Abstract

BACKGROUND AND PURPOSE: Patients with intracranial aneurysms may benefit from 4D flow MR imaging because the derived wall shear stress is considered a useful marker for risk assessment and growth of aneurysms. However, long scan times limit the clinical implementation of 4D flow MR imaging. Therefore, this study aimed to investigate whether highly accelerated, high resolution, 4D flow MR imaging at 7T provides reliable quantitative blood flow values in intracranial arteries and aneurysms.

MATERIALS AND METHODS: We used pseudospiral Cartesian undersampling with compressed sensing reconstruction to achieve high spatiotemporal resolution (0.5 mm isotropic, ∼30 ms) in a scan time of 10 minutes. We analyzed the repeatability of accelerated 4D flow scans and compared flow rates, stroke volume, and the pulsatility index with 2D flow and conventional 4D flow MR imaging in a flow phantom and 15 healthy subjects. Additionally, accelerated 4D flow MR imaging with high spatiotemporal resolution was acquired in 5 patients with aneurysms to derive wall shear stress.

RESULTS: Flow-rate bias compared with 2D flow was lower for accelerated than for conventional 4D flow MR imaging (0.31 ± 0.13, P = .22, versus 0.79 ± 0.17 mL/s, P < .01). Pulsatility index bias gave similar results. Stroke volume bias showed no difference for accelerated as well as for conventional 4D flow compared to 2D flow MR imaging. Repeatability for accelerated 4D flow was similar to that of 2D flow MR imaging. Increased temporal resolution for wall shear stress measurements in 5 intracranial aneurysms did not show a consistent effect for the wall shear stress but did show an effect for the oscillatory shear index.

CONCLUSIONS: Highly accelerated high spatiotemporal resolution 4D flow MR imaging at 7T in intracranial arteries and aneurysms provides repeatable and accurate quantitative flow values. Flow rate accuracy is significantly increased compared with conventional 4D flow scans.

ABBREVIATIONS:

CS: compressed sensing
N_card: number of reconstructed cardiac timeframes
OSI: oscillatory shear index
PROUD: PROspective Undersampling in multiple Dimensions
R: acceleration factor
SENSE: sensitivity encoding
UIA: unruptured intracranial aneurysm
WSS: wall shear stress
WSS_PS: peak systolic WSS
WSS_TA: time-averaged WSS

Intracranial aneurysms can threaten a patient’s life because of the risk of rupture followed by hemorrhagic stroke, which is associated with high morbidity and mortality.^1,2 Therefore, intracranial aneurysms need to be monitored regularly to assess whether they remain stable and asymptomatic or grow, with an increased risk of rupture. Treatment choices range from medical follow-up for stable aneurysms to coiling or clipping for growing aneurysms when the risk of rupture is considered too high.^1,2 For patients with unruptured intracranial aneurysms (UIAs) that need long-term follow-up and who have no contraindications for MR imaging, it is advised to consider MR imaging rather than CTA.² Typical MR imaging examinations use TOF or gadolinium contrast to visualize the aneurysms; however, advanced 4D flow MR imaging sequences can additionally visualize and quantify complex blood flow patterns in the aneurysm.³

In addition to measurements of basic blood flow volumes and velocities, 4D flow MR imaging allows the estimation of hemodynamic biomarkers, including wall shear stress (WSS).^4,5 Those have generated interest because abnormal WSS might be an indicator for aneurysm growth.^6,7 Moreover, it has been shown that the aneurysm wall is thin in regions of high WSS and vice versa.^8,9 This finding suggests that 4D flow MR imaging can provide added value in the clinical assessment of aneurysms. However, the implementation of 4D flow MR imaging in routine clinical practice has been hampered by long scan times required to achieve sufficient spatiotemporal resolution to image the complex blood flow patterns in the aneurysm.

A way to improve 4D flow MR imaging of intracranial aneurysms is the use of 7T MR imaging. The increased SNR compared with 3T can be exploited to obtain increased spatiotemporal resolution and/or reduced scan times.¹⁰ High spatiotemporal resolution is desirable because it has been shown that hemodynamic parameters such as mean flow, peak flow, flow waveform, WSS, and oscillatory shear index¹¹ (OSI) are influenced by spatial and temporal resolution of 4D flow MR imaging measurements.¹² A rule of thumb is that, to be clinically acceptable, the scan time should not exceed 10 minutes¹³ to fit into the clinical workflow and for patient comfort. Advanced accelerated 4D flow MR imaging techniques using pseudospiral Cartesian sampling with compressed sensing¹⁴ (CS) can be used to achieve such short scan times. These techniques have already proved effective to accelerate 4D flow MR imaging of the aorta and carotid arteries.^15,16 While we believe such a method could also be beneficial for intracranial 4D flow acquisitions at 7T, this benefit has not been thoroughly validated so far.

Therefore, the aim of this study was to investigate whether highly accelerated, high resolution 7T 4D flow MR imaging can provide reliable quantitative blood flow values in intracranial arteries. Specifically, we investigated flow accuracy in a flow phantom as well as healthy subjects. Furthermore, we aimed to explore the use of this accelerated sequence for deriving WSS in patients with unruptured intracranial aneurysms.

MATERIALS AND METHODS

Undersampled 4D flow MR Imaging Acquisition and CS Reconstruction

Accelerated 4D flow MR imaging was acquired on two similar 7T systems (Achieva 7T; Philips Healthcare), one located at the National 7T Facility, Lund, Sweden, and one at the Spinoza Center for Neuroimaging, Royal Netherlands Academy of Arts and Sciences, Amsterdam, the Netherlands. The scanners were equipped with a 32-channel receive head coil (2TX/32RX; NovaMedical) and an in-house-developed software modification called (PROspective Undersampling in multiple Dimensions [PROUD]), which enables a pseudospiral k_y/k_z-plane acquisition scheme designed for incoherent undersampling with a variable sampling density.^15,16 The acquisition was performed continuously without cardiac triggering. The signal of a peripheral pulse transducer was stored with the k-space data and used for retrospective binning in different cardiac timeframes. Due to physiologic heart rate variability, the resulting k-t space was randomly filled, which is favorable for CS reconstruction and enabled retrospective changes of the temporal resolution.¹⁶

Reconstructions were performed off-line in Matlab R2016a (MathWorks) using MRecon (GyroTools) in combination with the Berkeley Advanced Reconstruction Toolbox (https://mrirecon.github.io/bart/).¹⁷ A nonlinear parallel imaging and CS reconstruction were performed using a total variation sparsifying operator in time as previously described.¹⁶

Flow Phantom

A flow phantom was constructed, consisting of two merging plastic tubes surrounded by water (On-line Fig 1). Table salt (∼20 g/L) was added to the water to improve B₁ field homogeneity. The inlets were connected to a custom-made MR imaging–compatible mechanical pulsatile pump placed on the scanner table. A linear actuator, in the form of a servo motor (CM2-X-56B20A; MUSCLE CORPORATION) coupled to a ball screw actuator (VLACT55-12-0150; THK CO., LTD.), was used to produce a pulsatile flow (peak flow rate, 2 mL/s; pump output, 0.175 L/min; heart rate, 57 beats per minute). The pump software provided a trigger signal to the MR imaging system for retrospective binning. The flow phantom was scanned at the National 7T Facility, Lund, Sweden.

At both inlets as well as the outlet, 2D flow MR imaging was acquired and considered as a reference (spatial resolution, 0.3 × 0.3 × 3.0 mm³; temporal resolution, 26.0 ms). The 2D flow scans were performed at the beginning and end of the experiment to verify flow stability. Two 4D flow scans were obtained with sensitivity encoding (SENSE) with an acceleration factor R = 3 (spatial resolution, 0.5 mm isotropic; temporal resolution, 65.8–87.7 ms). Moreover, six PROUD 4D flow scans were obtained with nominal acceleration factors of R = 4, 8, 12, 16, 32 (spatial resolution, 0.5 mm isotropic). These scans were reconstructed in 15 different combinations of effective acceleration factors and number of cardiac frames (N_card) (temporal resolution, 26.3–87.7 ms). Detailed scan parameters are listed in the On-line Table. All measurements were performed in the same session without interrupting the flow in the phantom.

Healthy Subjects

Fifteen healthy subjects (9 men/6 women; mean age, 25.4 ± 1.2 years) were included in the study. Ethical approval was obtained at the Regional Ethical Review Board in Lund, Sweden, or waived by the local ethics committee of the Amsterdam University Medical Centers, University of Amsterdam, Amsterdam, the Netherlands. Written informed consent was provided by all subjects.

After scout scans to locate the intracranial arteries in the circle of Willis, five blocks of flow MR imaging scans were performed in a randomized order per subject: 1) SENSE 4D flow, 2) PROUD 4D flow, 3) PROUD 4D flow rescan, 4) 2D flow, and 5) 2D flow rescan. The 4D flow scans (spatial resolution, 0.5 mm isotropic) were placed in a transversal plane to include the greatest possible part of the vasculature. The 2D flow scans were placed in a sagittal plane orthogonal to both left and right middle cerebral arteries and were considered as the reference for flow measurements (spatial resolution, 0.5 × 0.5 × 3.0 mm³; mean temporal resolution, 31.3 ± 0.6 ms; N_card, 30). The PROUD 4D flow scans were reconstructed with three different numbers of cardiac frames (N_card = 7, 15, 30), resulting in different mean temporal resolutions (134.0 ± 2.7 ms, 62.6 ± 1.2 ms, 31.3 ± 0.6 ms) and mean acceleration factors R = 5.7 ± 0.1, 10.6 ± 0.2, 19.9 ± 0.3, respectively. The SENSE 4D flow scans were reconstructed into 7 cardiac frames. Detailed scan parameters are shown in the On-line Table.

Patients with Unruptured Intracranial Aneurysms

Five patients with UIAs (2 men/3 women; mean age, 63.1 ± 3.6 years) were scanned. Subjects provided written informed consent before the start of the study, which was approved by the local ethics committee of the Amsterdam University Medical Centers, University of Amsterdam, Amsterdam, the Netherlands.

A PROUD 4D flow scan was added to the examination protocol with a scan time of around 10 minutes. These scans were reconstructed with N_card = 10, 30 and R = 9.9 ± 0.2, 26.1 ± 0.6, respectively (spatial resolution, 0.5 mm isotropic; mean temporal resolution, 94.7 ± 3.6 and 31.6 ± 1.2 ms). Detailed scan parameters are shown in the On-line Table.

Flow Analysis

Quantitative flow analysis was performed using Segment,¹⁸ Version 2.2 (Medviso; http://medviso.com/).

Flow Phantom.

All 4D flow MR imaging scans were resliced to the positions of the three 2D flow MR imaging scans (outlet, inlet right, and inlet left). A ROI was drawn on the cross-section of each tube, and these ROIs were copied to all datasets. Background phase correction was performed by defining a donut-shaped region as static tissue around the actual ROI. 4D flow scans were compared with the 2D flow scans, and mean differences over the three ROIs were measured of the peak flow rate, systolic stroke volume, and diastolic flow rate standard error. The systolic part (half of a cardiac cycle) was used for flow quantification, and the diastolic part for diastolic flow rate standard error calculation.

The 2D flow MR imaging scans were analyzed for repeatability by comparing the scan and rescan values.

Healthy Subjects.

All 4D flow MR imaging scans were resliced to the positions of the two 2D flow MR imaging scans. ROIs were drawn on the cross-section of both MCAs in every dataset. Background phase correction was performed as in the flow phantom experiment. The first 4D flow MR imaging scans were compared with the first 2D flow MR imaging scan in peak flow rate, stroke volume, and the pulsatility index. The pulsatility index was calculated as the difference between the cross-sectional mean velocities at systole/maximum and diastole/minimum velocity divided by the mean velocity over the cardiac cycle.¹⁹ Repeatability of the 2D flow MR imaging scan and the PROUD 4D flow MR imaging scan was analyzed by comparing the scan/rescan variability.

Patients with UIAs.

The PROUD 4D flow MR imaging scans were postprocessed into a phase-contrast MRA,²⁰ which was used with the magnitude images to segment the intracranial vasculature and the aneurysms. WSS was calculated using in-house-developed software²¹ in Matlab R2016a by multiplying the wall shear rate by the dynamic viscosity of blood (3.2 × 10⁻³Pa·s). We performed two types of WSS calculations: WSS calculated at peak systole (WSS_PS) and time-averaged WSS (WSS_TA), in which the WSS is expressed as the average over the cardiac cycle. The WSS_TA, WSS_PS, and OSI were calculated for five patients with UIAs and expressed as the spatial mean value over the two ROIs: UIA only or surrounding vessel. Streamline flow visualizations of the aneurysms were made in GTFlow, Version 3.2.4 (GyroTools).

Statistical Analysis.

Differences between two flow scans (A and B) were tested by orthogonal regression analysis and Spearman rank ρ as well as Bland-Altman plots (X = A, Y = A-B). Bland-Altman results were presented as the bias or mean difference ± standard error of the mean and limits of agreement. A Wilcoxon signed rank test was used for statistical evaluations with a P < .05 level of significance.

RESULTS

Flow Phantom

Flow phantom measurements showed good accuracy even at high acceleration factors (peak flow rate difference, <10%). The best trade-off between acceleration and flow rate error was obtained for R = 20–30 and N_card ≥ 30. For a small peak flow rate error, at least 30 cardiac frames were needed.

Flow rates measured at the outlet ROI for low, moderate, and high temporal resolution are depicted in Fig 1. The differences between 2D flow and 4D flow MR imaging (2D flow – 4D flow) over all ROIs and scan sequences are listed in Table 1, and more flow rates are shown in On-line Fig 2.

FIG 1.

A, Flow rates of the outlet ROI of exemplary low (12 frames), moderate (16–24 frames), and high temporal resolution (30–40 frames). B, Zoomed view into the peak flow rate of the outlet ROI.

View this table:

Table 1:

Results of the flow phantom analysis in terms of peak flow rate difference, stroke volume difference, and the SD of the diastolic flow

Compared with 2D flow MR imaging, PROUD 4D flow MR imaging with N_card ≥ 30 showed the smallest peak flow rate underestimation. Furthermore, PROUD 4D flow MR imaging with N_card = 24 underestimated peak flow rate, comparable with conventional SENSE-accelerated 4D flow MR imaging. All PROUD 4D flow scans with N_card = 12 underestimated peak flow rate to a greater extent than 4D flow scans with N_card = 24. Stroke volume differences had no clear trend over the temporal resolutions, and the error was never higher than 15%. Diastolic flow rate standard error showed higher signal variation for the 4D flow than the 2D flow scans, but no clear trend with acceleration could be observed.

The 2D flow MR imaging repeatability analysis showed steady pump conditions (Bland-Altman: bias, 0.002 ± 0.007 mL/s; limits of agreement, 0.14 mL/s; regression, y = 0.98x + 0.00, ρ = 0.99).

Healthy Subjects

Flow measurements in healthy subjects showed good agreement between highly accelerated (R = 20) high temporal resolution 4D flow and 2D flow scans.

Examples of magnitude and phase images of the three 4D flow scans are shown in On-line Fig 3. Image quality was comparable, but more blurring was visible for high acceleration, particularly in the magnitude images. Peak flow rate differences between 2D flow and 4D flow MR imaging are shown in Table 2 and On-line Fig 4. It can be seen that decreased for higher temporal resolutions. For N_card = 7, 15, values were significantly different, but not for N_card = 30. Orthogonal regression analysis revealed an improved correlation between 2D and 4D flow MR imaging when using a higher number of cardiac frames. Similar peak flow rate repeatability (bias and limits of agreement) was observed for PROUD 4D flow compared with the 2D flow scans (Table 2 and On-line Fig 5). No significant differences in stroke volumes were observed for the four 4D flow scans compared with the 2D flow scans (Table 2 and On-line Fig 6). The stroke volume in PROUD 4D flow and 2D flow scans also had similar repeatability (Table 2 and On-line Fig 7). Pulsatility index differences were significant for 4D flow MR imaging with N_card = 7, 15, compared with the 2D flow MR imaging, but not for 4D flow MR imaging with N_card = 30 (Table 2 and On-line Fig 8). However, the Spearman ρ did not indicate a high correlation for the pulsatility index between 2D flow and 4D flow MR imaging.

View this table:

Table 2:

Results of the healthy subject analysis in terms of peak flow rate difference, peak flow rate repeatability, stroke volume difference, stroke volume repeatability, and pulsatility index difference

Patients with UIAs

The 4D flow MR imaging scans in the patients with UIAs showed a trend toward higher OSI for increased temporal resolution, but a similar trend could not be found in the WSS of all aneurysms.

Vortices were seen in the streamline visualizations for patients 1, 4, and 5 (Fig 2), but not for patients 2 and 3. Streamlines over the cardiac cycle of all patients are visualized in the On-line Video. Visualizations of WSS_TA, WSS_PS, and OSI of patient 5 for N_card = 10, 30 are shown in Fig 3. With increased temporal resolution, the WSS_PS and the OSI increased as well in some areas of the aneurysm wall. In Fig 4, the mean values of WSS_TA, WSS_PS, and OSI per UIA are presented as well as the mean values of WSS_TA and WSS_PS of their in- and outflow vessels. In the UIAs of patients 2 and 5 as well as in the vessel in all patients, WSS_PS was higher than WSS_TA, regardless of the temporal resolution. In the UIAs of patients 3 and 5, WSS_TA and WSS_PS were higher for N_card = 30 than for N_card = 10 frames. In the UIAs of patients 1, 2, and 4, no clear difference or trend between the two temporal resolutions was observed. In all UIAs, the OSI was higher for N_card = 30 than for N_card = 10 frames.

FIG 2.

Streamlines of patient 5 with an UIA to visualize blood flow patterns in the aneurysm (red arrow) of 10-mm diameter located at the left MCA. A video of streamlines of all five patients with UIAs over the cardiac cycle can be seen in the On-line Video. MCA indicates middle cerebral artery.

FIG 3.

Example WSS and OSI analysis of patient 5. Time-averaged WSS on the left, peak-systolic WSS in the middle, and OSI on the right for datasets with either 10 (upper row) or 30 (lower row) cardiac frames. MCA indicates middle cerebral artery; a.u., arbitrary units.

FIG 4.

The changes in time-averaged WSS, peak systolic WSS, and OSI in the aneurysm and in the surrounding vessel when increasing the temporal resolution from 10 to 30 cardiac frames, or from 95 to 32 ms, respectively. ACOM indicates anterior communicating artery; ICA, internal carotid artery; MCA, middle cerebral artery; a.u., arbitrary units.

DISCUSSION

This study investigated whether highly accelerated high resolution 4D flow MR imaging in intracranial arteries and aneurysms at 7T can provide reliable quantitative flow values in a clinically acceptable scan time of about 10 minutes using a pseudospiral Cartesian sampling scheme with CS reconstruction. Results showed that stroke volumes and peak flow rate values from 4D flow MR imaging, with R = 20–30 acceleration and a high temporal resolution of 30 cardiac frames, were similar to those obtained with 2D flow MR imaging.

Pulsatile flow phantom experiments were used to investigate the trade-off between scan acceleration and accuracy of flow measurements. The error in the peak flow rate decreased with an increasing number of cardiac frames, and stroke volume was underestimated in an acceptable level even for high temporal resolution (<3× difference of 2D rescan). The diastolic flow rate standard error in 4D flow scans showed 2–4 times larger variation than the 2D reference scans, and there was no clear trend observable for changes in acceleration or temporal resolution. Thus, we interpret this finding as acceptable random reconstruction- or undersampling-based artifacts. Taken together, we used R = 30 as maximum acceleration factor for the in vivo scans when reconstructing the data into 30 cardiac frames. In vivo measurements then confirmed that highly accelerated 4D flow MR imaging with a high temporal resolution resulted in small errors compared with 2D flow MR imaging. Furthermore, PROUD 4D flow MR imaging was superior to conventional 4D flow MR imaging, especially for those reconstructed in 30 cardiac frames, which resulted in the lowest peak flow rate difference. When we compared PROUD with conventional 4D flow MR imaging, the flow rate accuracy (peak flow rate difference to 2D flow MR imaging) was increased by 60%, while both had the same scan duration.

Reduced scan time is important for patients, which is why high acceleration is needed. However, there are only a few reports on acceleration factors above 10, especially at 7T field strength. Comparable high acceleration is only reported for acquisition and reconstruction similar to methods used in this study, which all deal with large multidimensional datasets. Cheng et al²² reported R ≈ 19 for accelerated 4D flow MR imaging in the heart at 3T and also used temporal total variation in the reconstruction. Walheim et al²³ reported R = 19 for accelerated 5D flow MR imaging in the aorta at 3T. Rich et al²⁴ demonstrated 4D flow imaging of the aortic valve at 1.5T with acceleration up to a factor 27. Gu et al²⁵ also acquired 4D flow data in the intracranial arteries, but at 1.5T, and reported a usable acceleration up to a factor of 30.

Considering the aneurysm and surrounding intracranial vasculature, two main observations could be made from the derived WSS. First, WSS_PS and WSS_TA as well as OSI in the vessels were higher with increased temporal resolution. Second, WSS_PS was higher than WSS_TA regardless of the temporal resolution. These two observations can be explained by the increased peak flow rate accuracy observed in the phantom experiment and healthy subjects. Measured WSS in the aneurysms was lower than in the surrounding intracranial vasculature, which is in line with previous results.²⁶ Increased temporal resolution resulted in increased OSI in all aneurysms, which has also been reported previously.¹² Most interesting, increased temporal resolution resulted in increased WSS in only two aneurysms but rather constant WSS in three aneurysms. Moreover, WSS_PS was higher than WSS_TA in only two aneurysms. Thus, we did not observe a clear trend for the WSS in the aneurysms. A likely explanation can be found in the differences in size and geometry of the aneurysms; for a UIA with low and constant flow, the improved temporal resolution will not change WSS estimations.²⁶ For cases of complex flow patterns, this finding suggests that flow and WSS estimations in aneurysms will benefit from the improved temporal resolution provided by our accelerated protocols. We will pursue studies in larger cohorts of patients to further investigate the added value of high temporal resolution WSS in aneurysm rupture risk predictions.

Previous studies on WSS do not agree on which WSS calculation method should be used,^{4,13,26⇓-28} but temporal fluctuations in WSS were associated with aneurysm risk of growth or formation.^11,27 Thus, indicators like the OSI¹¹ or the aneurysm formation indicator²⁷ were developed to express these temporal fluctuations in WSS. Most of these models use a reference vector, commonly the time-averaged WSS vector, based on the assumption that the endothelial cells are preferentially aligned in the direction of this vector.²⁷ When one uses these metrics, it is likely that a higher temporal resolution will give a more accurate outcome of temporal fluctuations in WSS.

A possible key application for high temporal resolution 4D flow MR imaging data is the combination of 4D flow MR imaging with computational fluid dynamics simulations, which have been shown to be a useful clinical tool for the prediction of the initiation, growth, and rupture of aneurysms.^13,26 Prior research combining fluid dynamics and 4D flow showed better agreement than calculations alone.^4,29 MR imaging–guided computational fluid dynamics can be used to correct the MR imaging–measured flow field, forcing it to satisfy the fluid mechanics equations. It can also be used as a means of reducing imperfections in the 4D flow MR imaging measurements and may improve the ability to accurately derive clinically relevant secondary parameters such as WSS and pressure gradients at a much higher level of detail and confidence than was previously possible.^29,30

Because highly accelerated 4D flow MR imaging with high temporal resolution provided reliable quantitative flow values, we encourage using 4D flow MR imaging more often for UIA monitoring and growth prediction in future clinically focused studies.^6,26 The techniques must be investigated in a larger study of patients with UIAs to determine whether the derived WSS is affected by either the calculation method or the temporal resolution or both. If it turns out that high temporal resolution does not add any value for clinical outcome, then the acceleration of the 4D flow MR imaging scan can be used to improve coverage and resolution or to further shorten the scan duration. Eventually, a consensus on 4D flow MR imaging followed by WSS calculation and analysis for intracranial aneurysm is desirable.

Limitations

A more complex phantom could mimic the intracranial arteries and their flow distribution better. Moreover, a pump with combined systolic pulsatile and diastolic steady flow could reduce the signal variation during the diastolic period. Furthermore, temporal fluctuations in WSS or OSI are sensitive to deviations induced by high acceleration. This limitation might be reduced when combining 4D flow MR imaging with computational fluid dynamics, which performs well in de-noising and preserving details in the velocity profiles. Finally, the study population to investigate the effect of increased temporal resolution for WSS calculations in patients with UIAs was too small to draw a significant conclusion.

CONCLUSIONS

Highly accelerated high spatiotemporal resolution 4D flow MR imaging in intracranial arteries and aneurysms at 7T provides repeatable quantitative flow values in a clinically acceptable scan time of ∼10 minutes, using pseudospiral Cartesian sampling with compressed sensing reconstruction. In vivo measurements showed that the accuracy of the flow rate was significantly increased compared with conventional low temporal resolution 4D flow scans. Data from patients with unruptured intracranial aneurysms show that the improved temporal resolution influences the oscillatory shear index.

Acknowledgment

We would like to thank Eva L. Leemans and Bart M.W. Cornelissen for the collection of the patient data, which was funded by the Stichting Toegepast Wetenschappelijk Instituut voor Neuromodulatie.

Footnotes

Funding was obtained by the Netherlands Organization for Scientific Research grant No. 13928 (HTSM2014).
Disclosures: Lukas M. Gottwald—RELATED: Grant: Stichting voor de Technische Wetenschappen, Comments: Stichting voor de Technische Wetenschappen, grant 13928. The Foundation for the Technical Sciences (Stichting voor de Technische Wetenschappen) is a Dutch foundation that was founded in 1981 and aims to realize the transfer of knowledge between technical sciences and users. To this end, Stichting voor de Technische Wetenschappen brings researchers and potential users together and finances technical-scientific research.* *Money paid to the institution.

Indicates open access to non-subscribers at www.ajnr.org

References

1.↵
1. Steiner T,
2. Juvela S,
3. Unterberg A, et al
. European stroke organization guidelines for the management of intracranial aneurysms and subarachnoid haemorrhage. Cerebrovasc Dis 2013;35:93–112 doi:10.1159/000346087 pmid:23406828
CrossRef PubMed
2.↵
1. Thompson BG,
2. Brown RD,
3. Amin-Hanjani S, et al
; American Heart Association Stroke Council, Council on Cardiovascular and Stroke Nursing, and Council on Epidemiology and Prevention; American Heart Association; American Stroke Association. Guidelines for the management of patients with unruptured intracranial aneurysms. Stroke 2015;46:2368–2400 doi:10.1161/STR.0000000000000070 pmid:26089327
Abstract/FREE Full Text
3.↵
1. Markl M,
2. Frydrychowicz A,
3. Kozerke S, et al
. 4D flow MRI. J Magn Reson Imaging 2012;36:1015–36 doi:10.1002/jmri.23632 pmid:23090914
CrossRef PubMed
4.↵
1. van Ooij P,
2. Potters WV,
3. Guédon A, et al
. Wall shear stress estimated with phase contrast MRI in an in vitro and in vivo intracranial aneurysm. J Magn Reson Imaging 2013;38:876–84 doi:10.1002/jmri.2405]1 pmid:23417769
CrossRef PubMed
5.↵
1. Coolen BF,
2. Calcagno C,
3. van Ooij P, et al
. Vessel wall characterization using quantitative MRI: what’s in a number? MAGMA 2018;31:201–22 doi:10.1007/s10334-017-0644-x pmid:28808823
CrossRef PubMed
6.↵
1. Boussel L,
2. Rayz V,
3. McCulloch C, et al
. Aneurysm growth occurs at region of low wall shear stress: patient-specific correlation of hemodynamics and growth in a longitudinal study. Stroke 2008;39:2997–3002 doi:10.1161/STROKEAHA.108.521617 pmid:18688012
Abstract/FREE Full Text
7.↵
1. Chung BJ,
2. Mut F,
3. Putman CM, et al
. Identification of hostile hemodynamics and geometries of cerebral aneurysms: a case-control study. AJNR Am J Neuroradiol 2018;39:1860–66 doi:10.3174/ajnr.A5764 pmid:30166431
Abstract/FREE Full Text
8.↵
1. Blankena R,
2. Kleinloog R,
3. Verweij BH, et al
. Thinner regions of intracranial aneurysm wall correlate with regions of higher wall shear stress: a 7T MRI study. AJNR Am J Neuroradiol 2016;37:1310–17 doi:10.3174/ajnr.A4734] pmid:26892986
Abstract/FREE Full Text
9.↵
1. Cebral JR,
2. Detmer F,
3. Chung BJ, et al
. Local hemodynamic conditions associated with focal changes in the intracranial aneurysm wall. AJNR Am J Neuroradiol 2019;40:510–16 doi:10.3174/ajnr.A5970 30733253 pmid:30733253
Abstract/FREE Full Text
10.↵
1. van Ooij P,
2. Zwanenburg JJM,
3. Visser F, et al
. Quantification and visualization of flow in the circle of Willis: time-resolved three-dimensional phase contrast MRI at 7 T compared with 3 T. Magn Reson Med 2013;69:868–76 doi:10.1002/mrm.24317 pmid:22618854
CrossRef PubMed
11.↵
1. Ku DN,
2. Giddens DP,
3. Phillips DJ, et al
. Hemodynamics of the normal human carotid bifurcation: in vitro and in vivo studies. Ultrasound Med Biol 1985;11:13–26 pmid:3160152
CrossRef PubMed
12.↵
1. Cibis M,
2. Potters WV,
3. Gijsen FJ, et al
. The effect of spatial and temporal resolution of cine phase contrast MRI on wall shear stress and oscillatory shear index assessment. PLoS One 2016;11:e0163316–15 doi:10.1371/journal.pone.0163316 pmid:27669568
CrossRef PubMed
13.↵
1. Boussel L,
2. Rayz V,
3. Martin A, et al
. Phase-contrast magnetic resonance imaging measurements in intracranial aneurysms in vivo of flow patterns, velocity fields, and wall shear stress: comparison with computational fluid dynamics. Magn Reson Med 2009;6:409–17 doi:10.1002/mrm.21861 pmid:19161132
CrossRef PubMed
14.↵
1. Lustig M,
2. Donoho D,
3. Pauly JM
. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med 2007;58:1182–95 doi:10.1002/mrm.21391 pmid:17969013
CrossRef PubMed
15.↵
1. Peper ES,
2. Gottwald LM,
3. Zhang Q, et al
. Highly accelerated 4D flow cardiovascular magnetic resonance using a pseudo-spiral Cartesian acquisition and compressed sensing reconstruction for carotid flow and wall shear stress. J Cardiovasc Magn Reson 2020;22:7 doi:10.1186/s12968-019-0582-z pmid:31959203
CrossRef PubMed
16.↵
1. Gottwald LM,
2. Peper ES,
3. Zhang Q, et al
. 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. NMR Biomed 2020;33:e4255 doi:10.1002/nbm.4255 pmid:31957927
CrossRef PubMed
17.↵
1. Uecker M,
2. Ong F,
3. Tamir JI, et al
. Berkeley advanced reconstruction toolbox. In: Proceedings of the International Society for Magnetic Resonance in Medicine, May 30 to June 5, 2015; Toronto, Ontario, Canada; 2486
18.↵
1. Heiberg E,
2. Sjögren J,
3. Ugander M, et al
. Design and validation of Segment: freely available software for cardiovascular image analysis. BMC Med Imaging 2010;10:1 doi:10.1186/1471-2342-10-1 pmid:20064248
CrossRef PubMed
19.↵
1. Bouvy WH,
2. Geurts LJ,
3. Kuijf HJ, et al
. Assessment of blood flow velocity and pulsatility in cerebral perforating arteries with 7-T quantitative flow MRI. NMR Biomed 2016;29:1295–1304 doi:10.1002/nbm.3306 pmid:25916399
CrossRef PubMed
20.↵
1. Bock J,
2. Frydrychowicz A,
3. Stalder AF, et al
. 4D phase contrast MRI at 3 T: effect of standard and blood-pool contrast agents on SNR, PC-MRA, and blood flow visualization. Magn Reson Med 2010;63:330–38 doi:10.1002/mrm.22199 pmid:20024953
CrossRef PubMed
21.↵
1. Potters WV,
2. van Ooij P,
3. Marquering H, et al
. Volumetric arterial wall shear stress calculation based on cine phase contrast MRI. J Magn Reson Imaging 2015;41:505–16 doi:10.1002/jmri.24560 pmid:24436246
CrossRef PubMed
22.↵
1. Cheng JY,
2. Hanneman K,
3. Zhang T, et al
. Comprehensive motion-compensated highly accelerated 4D flow MRI with ferumoxytol enhancement for pediatric congenital heart disease. J Magn Reson Imaging 2016;43:1355–68 doi:10.1002/jmri.25106 pmid:26646061
CrossRef PubMed
23.↵
1. Walheim J,
2. Dillinger H,
3. Kozerke S
. Multipoint 5D flow cardiovascular magnetic resonance - accelerated cardiac- and respiratory-motion resolved mapping of mean and turbulent velocities. J Cardiovasc Magn Reson 2019;21:42 doi:10.1186/s12968-019-0549-0 pmid:31331353
CrossRef PubMed
24.↵
1. Rich A,
2. Potter LC,
3. Jin N, et al
. A Bayesian model for highly accelerated phase-contrast MRI. Magn Reson Med 2016;76:689–701 doi:10.1002/mrm.25904 pmid:26444911
CrossRef PubMed
25.↵
1. Gu T,
2. Korosec FR,
3. Block WF, et al
. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol 2005;26:743–49 pmid:15814915
Abstract/FREE Full Text
26.↵
1. Shojima M,
2. Oshima M,
3. Takagi K, et al
. Magnitude and role of wall shear stress on cerebral aneurysm: computational fluid dynamic study of 20 middle cerebral artery aneurysms. Stroke 2004;35:2500–05 doi:10.1161/01.STR.0000144648.89172.0f pmid:15514200
Abstract/FREE Full Text
27.↵
1. Mantha A,
2. Karmonik C,
3. Benndorf G, et al
. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. AJNR Am J Neuroradiol 2006;27:1113–18 pmid:16687554
Abstract/FREE Full Text
28.↵
1. Schnell S,
2. Wu C,
3. Ansari SA
. Four-dimensional MRI flow examinations in cerebral and extracerebral vessels: ready for clinical routine? Curr Opin Neurol 2016;29:19–28 doi:10.1097/WCO.0000000000000341 pmid:27262148
CrossRef PubMed
29.↵
1. Töger J,
2. Zahr MJ,
3. Aristokleous N, et al
. Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data. Magn Reson Med 2020 Apr 8. [Epub ahead of print] doi:10.1002/mrm.28269 pmid:32270549
CrossRef PubMed
30.↵
1. Bakhshinejad A,
2. Baghaie A,
3. Vali A, et al
. Merging computational fluid dynamics and 4D-Flow MRI using proper orthogonal decomposition and ridge regression. J Biomech 2017;58:162–73 doi:10.1016/j.jbiomech.2017.05.004 pmid:2857790
CrossRef PubMed

Received February 12, 2020.
Accepted after revision April 21, 2020.

View Abstract

In this issue

Download PDF

Email Article

Citation Tools

Purchase

Cited By...

More in this TOC Section

Show more ADULT BRAIN

[1] 1.↵
Steiner T,
Juvela S,
Unterberg A, et al
. European stroke organization guidelines for the management of intracranial aneurysms and subarachnoid haemorrhage. Cerebrovasc Dis 2013;35:93–112 doi:10.1159/000346087 pmid:23406828
CrossRef PubMed

[2] Steiner T,

[3] Juvela S,

[4] Unterberg A, et al

[5] 2.↵
Thompson BG,
Brown RD,
Amin-Hanjani S, et al
; American Heart Association Stroke Council, Council on Cardiovascular and Stroke Nursing, and Council on Epidemiology and Prevention; American Heart Association; American Stroke Association. Guidelines for the management of patients with unruptured intracranial aneurysms. Stroke 2015;46:2368–2400 doi:10.1161/STR.0000000000000070 pmid:26089327
Abstract/FREE Full Text

[6] Thompson BG,

[7] Brown RD,

[8] Amin-Hanjani S, et al

[9] 3.↵
Markl M,
Frydrychowicz A,
Kozerke S, et al
. 4D flow MRI. J Magn Reson Imaging 2012;36:1015–36 doi:10.1002/jmri.23632 pmid:23090914
CrossRef PubMed

[10] Markl M,

[11] Frydrychowicz A,

[12] Kozerke S, et al

[13] 4.↵
van Ooij P,
Potters WV,
Guédon A, et al
. Wall shear stress estimated with phase contrast MRI in an in vitro and in vivo intracranial aneurysm. J Magn Reson Imaging 2013;38:876–84 doi:10.1002/jmri.2405]1 pmid:23417769
CrossRef PubMed

[14] van Ooij P,

[15] Potters WV,

[16] Guédon A, et al

[17] 5.↵
Coolen BF,
Calcagno C,
van Ooij P, et al
. Vessel wall characterization using quantitative MRI: what’s in a number? MAGMA 2018;31:201–22 doi:10.1007/s10334-017-0644-x pmid:28808823
CrossRef PubMed

[18] Coolen BF,

[19] Calcagno C,

[20] van Ooij P, et al

[21] 6.↵
Boussel L,
Rayz V,
McCulloch C, et al
. Aneurysm growth occurs at region of low wall shear stress: patient-specific correlation of hemodynamics and growth in a longitudinal study. Stroke 2008;39:2997–3002 doi:10.1161/STROKEAHA.108.521617 pmid:18688012
Abstract/FREE Full Text

[22] Boussel L,

[23] Rayz V,

[24] McCulloch C, et al

[25] 7.↵
Chung BJ,
Mut F,
Putman CM, et al
. Identification of hostile hemodynamics and geometries of cerebral aneurysms: a case-control study. AJNR Am J Neuroradiol 2018;39:1860–66 doi:10.3174/ajnr.A5764 pmid:30166431
Abstract/FREE Full Text

[26] Chung BJ,

[27] Mut F,

[28] Putman CM, et al

[29] 8.↵
Blankena R,
Kleinloog R,
Verweij BH, et al
. Thinner regions of intracranial aneurysm wall correlate with regions of higher wall shear stress: a 7T MRI study. AJNR Am J Neuroradiol 2016;37:1310–17 doi:10.3174/ajnr.A4734] pmid:26892986
Abstract/FREE Full Text

[30] Blankena R,

[31] Kleinloog R,

[32] Verweij BH, et al

[33] 9.↵
Cebral JR,
Detmer F,
Chung BJ, et al
. Local hemodynamic conditions associated with focal changes in the intracranial aneurysm wall. AJNR Am J Neuroradiol 2019;40:510–16 doi:10.3174/ajnr.A5970 30733253 pmid:30733253
Abstract/FREE Full Text

[34] Cebral JR,

[35] Detmer F,

[36] Chung BJ, et al

[37] 10.↵
van Ooij P,
Zwanenburg JJM,
Visser F, et al
. Quantification and visualization of flow in the circle of Willis: time-resolved three-dimensional phase contrast MRI at 7 T compared with 3 T. Magn Reson Med 2013;69:868–76 doi:10.1002/mrm.24317 pmid:22618854
CrossRef PubMed

[38] van Ooij P,

[39] Zwanenburg JJM,

[40] Visser F, et al

[41] 11.↵
Ku DN,
Giddens DP,
Phillips DJ, et al
. Hemodynamics of the normal human carotid bifurcation: in vitro and in vivo studies. Ultrasound Med Biol 1985;11:13–26 pmid:3160152
CrossRef PubMed

[42] Ku DN,

[43] Giddens DP,

[44] Phillips DJ, et al

[45] 12.↵
Cibis M,
Potters WV,
Gijsen FJ, et al
. The effect of spatial and temporal resolution of cine phase contrast MRI on wall shear stress and oscillatory shear index assessment. PLoS One 2016;11:e0163316–15 doi:10.1371/journal.pone.0163316 pmid:27669568
CrossRef PubMed

[46] Cibis M,

[47] Potters WV,

[48] Gijsen FJ, et al

[49] 13.↵
Boussel L,
Rayz V,
Martin A, et al
. Phase-contrast magnetic resonance imaging measurements in intracranial aneurysms in vivo of flow patterns, velocity fields, and wall shear stress: comparison with computational fluid dynamics. Magn Reson Med 2009;6:409–17 doi:10.1002/mrm.21861 pmid:19161132
CrossRef PubMed

[50] Boussel L,

[51] Rayz V,

[52] Martin A, et al

[53] 14.↵
Lustig M,
Donoho D,
Pauly JM
. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med 2007;58:1182–95 doi:10.1002/mrm.21391 pmid:17969013
CrossRef PubMed

[54] Lustig M,

[55] Donoho D,

[56] Pauly JM

[57] 15.↵
Peper ES,
Gottwald LM,
Zhang Q, et al
. Highly accelerated 4D flow cardiovascular magnetic resonance using a pseudo-spiral Cartesian acquisition and compressed sensing reconstruction for carotid flow and wall shear stress. J Cardiovasc Magn Reson 2020;22:7 doi:10.1186/s12968-019-0582-z pmid:31959203
CrossRef PubMed

[58] Peper ES,

[59] Gottwald LM,

[60] Zhang Q, et al

[61] 16.↵
Gottwald LM,
Peper ES,
Zhang Q, et al
. 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. NMR Biomed 2020;33:e4255 doi:10.1002/nbm.4255 pmid:31957927
CrossRef PubMed

[62] Gottwald LM,

[63] Peper ES,

[64] Zhang Q, et al

[65] 17.↵
Uecker M,
Ong F,
Tamir JI, et al
. Berkeley advanced reconstruction toolbox. In: Proceedings of the International Society for Magnetic Resonance in Medicine, May 30 to June 5, 2015; Toronto, Ontario, Canada; 2486

[66] Uecker M,

[67] Ong F,

[68] Tamir JI, et al

[69] 18.↵
Heiberg E,
Sjögren J,
Ugander M, et al
. Design and validation of Segment: freely available software for cardiovascular image analysis. BMC Med Imaging 2010;10:1 doi:10.1186/1471-2342-10-1 pmid:20064248
CrossRef PubMed

[70] Heiberg E,

[71] Sjögren J,

[72] Ugander M, et al

[73] 19.↵
Bouvy WH,
Geurts LJ,
Kuijf HJ, et al
. Assessment of blood flow velocity and pulsatility in cerebral perforating arteries with 7-T quantitative flow MRI. NMR Biomed 2016;29:1295–1304 doi:10.1002/nbm.3306 pmid:25916399
CrossRef PubMed

[74] Bouvy WH,

[75] Geurts LJ,

[76] Kuijf HJ, et al

[77] 20.↵
Bock J,
Frydrychowicz A,
Stalder AF, et al
. 4D phase contrast MRI at 3 T: effect of standard and blood-pool contrast agents on SNR, PC-MRA, and blood flow visualization. Magn Reson Med 2010;63:330–38 doi:10.1002/mrm.22199 pmid:20024953
CrossRef PubMed

[78] Bock J,

[79] Frydrychowicz A,

[80] Stalder AF, et al

[81] 21.↵
Potters WV,
van Ooij P,
Marquering H, et al
. Volumetric arterial wall shear stress calculation based on cine phase contrast MRI. J Magn Reson Imaging 2015;41:505–16 doi:10.1002/jmri.24560 pmid:24436246
CrossRef PubMed

[82] Potters WV,

[83] van Ooij P,

[84] Marquering H, et al

[85] 22.↵
Cheng JY,
Hanneman K,
Zhang T, et al
. Comprehensive motion-compensated highly accelerated 4D flow MRI with ferumoxytol enhancement for pediatric congenital heart disease. J Magn Reson Imaging 2016;43:1355–68 doi:10.1002/jmri.25106 pmid:26646061
CrossRef PubMed

[86] Cheng JY,

[87] Hanneman K,

[88] Zhang T, et al

[89] 23.↵
Walheim J,
Dillinger H,
Kozerke S
. Multipoint 5D flow cardiovascular magnetic resonance - accelerated cardiac- and respiratory-motion resolved mapping of mean and turbulent velocities. J Cardiovasc Magn Reson 2019;21:42 doi:10.1186/s12968-019-0549-0 pmid:31331353
CrossRef PubMed

[90] Walheim J,

[91] Dillinger H,

[92] Kozerke S

[93] 24.↵
Rich A,
Potter LC,
Jin N, et al
. A Bayesian model for highly accelerated phase-contrast MRI. Magn Reson Med 2016;76:689–701 doi:10.1002/mrm.25904 pmid:26444911
CrossRef PubMed

[94] Rich A,

[95] Potter LC,

[96] Jin N, et al

[97] 25.↵
Gu T,
Korosec FR,
Block WF, et al
. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol 2005;26:743–49 pmid:15814915
Abstract/FREE Full Text

[98] Gu T,

[99] Korosec FR,

[100] Block WF, et al

[101] 26.↵
Shojima M,
Oshima M,
Takagi K, et al
. Magnitude and role of wall shear stress on cerebral aneurysm: computational fluid dynamic study of 20 middle cerebral artery aneurysms. Stroke 2004;35:2500–05 doi:10.1161/01.STR.0000144648.89172.0f pmid:15514200
Abstract/FREE Full Text

[102] Shojima M,

[103] Oshima M,

[104] Takagi K, et al

[105] 27.↵
Mantha A,
Karmonik C,
Benndorf G, et al
. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. AJNR Am J Neuroradiol 2006;27:1113–18 pmid:16687554
Abstract/FREE Full Text

[106] Mantha A,

[107] Karmonik C,

[108] Benndorf G, et al

[109] 28.↵
Schnell S,
Wu C,
Ansari SA
. Four-dimensional MRI flow examinations in cerebral and extracerebral vessels: ready for clinical routine? Curr Opin Neurol 2016;29:19–28 doi:10.1097/WCO.0000000000000341 pmid:27262148
CrossRef PubMed

[110] Schnell S,

[111] Wu C,

[112] Ansari SA

[113] 29.↵
Töger J,
Zahr MJ,
Aristokleous N, et al
. Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data. Magn Reson Med 2020 Apr 8. [Epub ahead of print] doi:10.1002/mrm.28269 pmid:32270549
CrossRef PubMed

[114] Töger J,

[115] Zahr MJ,

[116] Aristokleous N, et al

[117] 30.↵
Bakhshinejad A,
Baghaie A,
Vali A, et al
. Merging computational fluid dynamics and 4D-Flow MRI using proper orthogonal decomposition and ridge regression. J Biomech 2017;58:162–73 doi:10.1016/j.jbiomech.2017.05.004 pmid:2857790
CrossRef PubMed

[118] Bakhshinejad A,

[119] Baghaie A,

[120] Vali A, et al

Main menu

User menu

Search

American Journal of Neuroradiology

High Spatiotemporal Resolution 4D Flow MRI of Intracranial Aneurysms at 7T in 10 Minutes

Abstract

ABBREVIATIONS:

MATERIALS AND METHODS

Undersampled 4D flow MR Imaging Acquisition and CS Reconstruction

Flow Phantom

Healthy Subjects

Patients with Unruptured Intracranial Aneurysms

Flow Analysis

Flow Phantom.

Healthy Subjects.

Patients with UIAs.

Statistical Analysis.

RESULTS

Flow Phantom

Healthy Subjects

Patients with UIAs

DISCUSSION

Limitations

CONCLUSIONS

Acknowledgment

Footnotes

References

In this issue

Citation Manager Formats

Related Articles

Cited By...

More in this TOC Section

Similar Articles

News and Updates

Resources

Opportunities

American Society of Neuroradiology

Main menu

User menu

Search

High Spatiotemporal Resolution 4D Flow MRI of Intracranial Aneurysms at 7T in 10 Minutes

Abstract

ABBREVIATIONS:

MATERIALS AND METHODS

Undersampled 4D flow MR Imaging Acquisition and CS Reconstruction

Flow Phantom

Healthy Subjects

Patients with Unruptured Intracranial Aneurysms

Flow Analysis

Flow Phantom.

Healthy Subjects.

Patients with UIAs.

Statistical Analysis.

RESULTS

Flow Phantom

Healthy Subjects

Patients with UIAs

DISCUSSION

Limitations

CONCLUSIONS

Acknowledgment

Footnotes

References

In this issue

Citation Manager Formats

Jump to section

Related Articles

Cited By...

More in this TOC Section

Similar Articles

News and Updates

Resources

Opportunities

American Society of Neuroradiology