Influence of Spatial Resolution and Compressed SENSE Acceleration Factor on Flow Quantification with 4D Flow MRI at 3 Tesla

Four-dimensional (4D) flow MRI allows quantifying flow in blood vessels–non invasively and in vivo. The clinical use of 4D flow MRI in small vessels, however, is hampered by long examination times and limited spatial resolution. Compressed SENSE (CS-SENSE) is a technique that can accelerate 4D flow dramatically. Here, we investigated the effect of spatial resolution and CS acceleration on flow measurements by using 4D flow MRI in small vessels in vitro at 3 T. We compared the flow in silicon tubes (inner diameters of 2, 3, 4, and 5 mm) measured with 4D flow MRI, accelerated with four CS factors (CS = 2.5, 4.5, 6.5, and 13) and three voxel sizes (0.5, 1, and 1.5 mm3) to 2D flow MRI and a flow sensor. Additionally, the velocity field in an aneurysm model acquired with 4D flow MRI was compared to the one simulated with computational fluid dynamics (CFD). A strong correlation was observed between flow sensor, 2D flow MRI, and 4D flow MRI (rho > 0.94). The use of fewer than seven voxels per vessel diameter (nROI) resulted in an overestimation of flow in more than 5% of flow measured with 2D flow MRI. A negative correlation (rho = −0.81) between flow error and nROI were found for CS = 2.5 and 4.5. No statistically significant impact of CS factor on differences in flow rates was observed. However, a trend of increased flow error with increased CS factor was observed. In an aneurysm model, the peak velocity and stagnation zone were detected by CFD and all 4D flow MRI variants. The velocity difference error in the aneurysm sac did not exceed 11% for CS = 4.5 in comparison to CS = 2.5 for all spatial resolutions. Therefore, CS factors from 2.5–4.5 can appear suitable to improve spatial or temporal resolution for accurate quantification of flow rate and velocity. We encourage reporting the number of voxels per vessel diameter to standardize 4D flow MRI protocols.


Introduction
Blood flow is altered in a variety of intracranial diseases, such as arterial stenosis [1,2], intracranial aneurysms [3,4], and vascular malformations [5]. The study of hemodynamics in general can potentially improve diagnostic capabilities, therapeutic planning, and patient follow-up.
Four-dimensional (4D) phase-contrast magnetic resonance imaging (4D flow MRI) is a noninvasive technique and can be used to quantitatively and qualitatively analyze the hemodynamics of blood vessels [6]. It can provide time-resolved 3D visualization of complex flow patterns and quantify flow rates, velocities, and more advanced parameters, such as wall shear stress (WSS) and oscillatory shear index (OSI) [7][8][9][10]. Four-dimensional flow MRI is widely used for cardiac examinations, where the diameter of the arteries (e.g., thoracic or abdominal aorta) is relatively large. In neurovascular applications, 4D flow MRI has been used primarily for research purposes to study intracranial aneurysms, arteriovenous malformation, and Alzheimer's disease, among others [11]. Furthermore, neurovascular 4D flow MRI can potentially be used to monitor the effect of aneurysm therapy [12][13][14].
Although the number of neuroradiological applications for 4D flow MRI is increasing [15], the technique is rarely used as part of standard clinical care, as compared to the cardiovascular applications. Factors impeding widespread clinical use of 4D flow MRI include insufficient spatial resolution to assess the flow in small vessels and velocities near the vessel wall, long examination times, and high sensitivity to velocity-encoding parameters (Venc) [16]. Similar to reconstruction, quantification, and visualization of velocity, these datasets have needed to be evaluated offline, too, often by physicists. Nowadays, the vendors offer online reconstruction of the datasets, as well as basic data post-processing, for example, flow peak velocity, but not kinetic energy or vorticity [17].
Nevertheless, some initial studies have demonstrated that certain types of 4D flow MRI examinations can be obtained in about 2-8 min [18][19][20][21]. Therefore, we believe that implementation in clinical routine is inevitable. First, however, 4D flow MRI protocols need to be standardized and awareness of confounders is essential to achieve a broad clinical application.
Thus, we aimed to study the effect of spatial resolution in combination with acceleration on flow rates and velocity measured using 4D flow MRI. Three voxel sizes (0.5, 1.0, and 1.5 mm 3 ) and four compressed SENSE (CS) acceleration factors (2.5, 4.5, 6.5, and 13) were considered. We investigated the matter in vitro using silicone tubes of different diameters (2,3,4, and 5 mm) and a patient-derived 3D-printed aneurysm model. The results were compared to flow sensor, 2D flow MRI, and CFD simulations.

Flow Models and Circulation Setup
A set of silicone tubes with inner diameters (ID) in the range of 2-5 mm was used to create a model of tortuous neurovascular vessels (Figure 1a). To do so, tubes were wound around a 50-mL falcon tube. Chosen diameters cover the dimensions reported for human intracranial arteries, such as internal carotid artery (ICA, ID = 3.6 ± 0.9 mm), middle carotid artery (MCA, ID = 2.5 ± 0.6 mm), anterior cerebral artery (ACA, ID = 1.8 ± 0.6 mm), basilar artery (BA, ID = 3.3 ± 1.0 mm), and posterior cerebral arteries (PCA, ID = 2 ± 0.5 mm) [29]. Tubes were connected with a pulsatile pump (PD-1100, BDC Laboratories, Wheat Ridge, CO, USA). The glycerol-water mixture (40/60% by volume) circulated in the flow setup. The testing fluid had a viscosity of 3.72 cP according to Segur et al. [25], which corresponds to the viscosity of blood (3.5 cP [26]). Pulsatile flow at mean flow rates of 0.6, 1.2, 2.3, and 3.1 mL/s was applied to tubes with an ID of 2-5 mm, respectively, which correspond to the flow observed in vivo in the ICA (4.3 ± 0.8 mL/s), MCA (2.4 ± 0.5 mL/s), BA (2.4 ± 0.7 mL/s), and PCA (0.9 ± 0.2 mL/s) [30]. The flow peak velocity was about 60 cm/s in all four tubes. Therefore, the same flow velocity-encoding parameter (60 cm/s) was set during all MRI measurements (details below). The flow in the tubes was measured before MRI experiments one by one with a transonic flow sensor (US sensor, Transonic System Inc., Ithaca, NY, USA) located 10 cm downstream from the model outlets. The flow setup proved to be stable over at least 3 h (details in Figure S1). cm/s in all four tubes. Therefore, the same flow velocity-encoding parameter (60 cm/s) was set during all MRI measurements (details below). The flow in the tubes was measured before MRI experiments one by one with a transonic flow sensor (US sensor, Transonic System Inc., Ithaca, NY, USA) located 10 cm downstream from the model outlets. The flow setup proved to be stable over at least 3 h (details in Figure S1). A patient-specific aneurysm model was designed and 3D printed (Form 3, Formlabs Inc., Somerville, MA, USA) in-house [31] (Figure 1b). Aneurysm geometry was segmented from 3D rotational angiographic data of a patient with a 17.5-mm infraophthalmic extradural internal carotid artery (ICA) aneurysm. A mixture of methacrylic acid esters and a photoinitiator (Clear Photoreactive Resin, Formlabs Inc., Somerville, MA, USA) was used for 3D printing of an aneurysm model with a rigid wall. The inlet and outlet flow and A patient-specific aneurysm model was designed and 3D printed (Form 3, Formlabs Inc., Somerville, MA, USA) in-house [31] (Figure 1b). Aneurysm geometry was segmented from 3D rotational angiographic data of a patient with a 17.5-mm infraophthalmic extradural internal carotid artery (ICA) aneurysm. A mixture of methacrylic acid esters and a photoinitiator (Clear Photoreactive Resin, Formlabs Inc., Somerville, MA, USA) was used for 3D printing of an aneurysm model with a rigid wall. The inlet and outlet flow and pressure were measured approx. 10-15 cm apart from the model inlet and the outlets with the US flow and pressure (PRESS-N-000, PendoTech, Princeton, NJ, USA) sensors, respectively. For simplicity, tap water was circulated at the flow setup at a flow rate of 3.8 mL/s (Ismaltec MCP Standart, Cole Parmer, Vernon Hills, IL, USA). The chosen flow rate corresponds to a flow observed in vivo [30]. The velocity field was analyzed at the 2D plane (Figure 1c, left) and over the 3D volumetric aneurysm ROI (ROI, Figure 1c, right).
In addition, 3% (weight/weight) agarose gel was poured around the silicone tubes and aneurysm model before the MRI experiment to increase the overall MR signal of the in vitro models and to eliminate the air-tissue interface artifacts. Reinforced, polyvinyl chloride, braided tubing with an inner diameter of 6 mm and 2 m in length was used to connect the models to the pump placed outside of the MRI room.
All measurements were acquired at room temperature (21 • C). The silicone tubes and aneurysm model were supplied with different pumps due to the unavailability of the same circulation setup during each experiment. However, this fact should not impact the general applicability of these results.

Magnetic Resonance Imaging
A whole-body 3 T MRI clinical system equipped with a 32-channel volume head coil (Ingenia CX, R5 V6.1, Philips Healthcare, Best, the Netherlands) was used for all MRI experiments. The imaging protocol comprised: time-of-flight (TOF), 2D flow, and 4D flow MRI. All sequences were accelerated by using a compressed SENSE technique based on sensitive encoding (CS-SENSE) implemented by the vendor (Philips) [32,33]. CS-SENSE uses incoherent undersampling of variable density together with iterative (non-linear) reconstruction. CS-SENSE consists of a mix between compress sensing theory [34] and SENSE technology (parallel imaging technique) [35].
Both 2D and 4D flow MRI were acquired with the spoiled gradient-echo sequence with Cartesian sampling [15]. The maximum velocity was estimated by performing a series of 2D flow MRI measurements with varied velocity-encoding parameters (Venc) in the range from 30-150 cm/s. The final Venc values were chosen to be 10% above the maximum detected velocity and were equal to 60 and 80 cm/s for the experiment with silicone tubes and an aneurysm model, respectively, for both 2D and 4D flow MRI. One planar velocity encoding was used for 2D flow MRI and three velocity encoding directions were used for 4D flow MRI. The cardiac cycle was equal to 800 and 816 ms for the silicone tubes and aneurysm model, respectively, and 24 measurement points over a cardiac cycle were obtained. The resulting temporal resolution was 33-34 ms (acquired heart phases were set at 100%). Four-dimensional flow MR images were acquired using a balanced symmetric four-point phase-contrast velocity-encoding scheme (Hadamard) [36]. Acquisition planes of the 2D flow MRI were positioned perpendicular to the flow direction, as shown in Figure 1a (ROI A-C).
Isotropic voxel size (0.5, 1.0, and 1.5 mm 3 ) and CS acceleration factor were varied by the operator; all other parameters were kept the same or were set automatically by the system. TR and TE were automatically set to the shortest values allowed by the MRI system. MRI protocol parameters are summarized in Table 1.

Time-Resolved Hemodynamic Simulations
To enable a comparison between the measured aneurysm velocity field and the highly resolved in silico velocity results, image-based blood flow simulations were carried out. The numerical geometry was segmented from the TOF MRI data in vitro. Before performing the corresponding simulation, the model was spatially discretized using STAR-CCM+ 2020.01 (Siemens Product Lifecycle Management Software Inc., Plano, TX, USA) with a cell base size of ∆x = 0.1 mm, resulting in a total number of 5.5 million cells (polyhedral and prism layers). The flow rate measured using a flow sensor was applied to the inlet cross-section, and the outlet boundary condition was set to the pressure values provided by the pressure sensor placed in the aneurysm outlet. Similar to the in vitro model, all vessels were assumed to be rigid, and no-slip boundary conditions were applied. Finally, identical fluid properties (water) were defined for density and viscosity. The maximum Reynolds number (Re) was calculated with the mean velocity (u = 20 cm/s) and diameter (ID = 0.5 cm) of the ICA as Re = (u × ID × ρ)/µ, [37], with the water viscosity and density of µ = 0.890 cP and ρ = 997 kg/m 3 , respectively. This resulted in Re ≈ 1120; therefore, laminar flow conditions were considered.
To obtain a periodic solution, the time-dependent blood flow simulation constitutes three cardiac cycles (time step size ∆t = 0.001 s). However, the first two cycles were discarded, and only the last one was included in the analysis. The final post-processing of the simulation data was carried out in EnSight 10.2 (ANSYS Inc., Canonsburg, PA, USA).

Data Processing
The 2D and 4D flow MR images were reconstructed on the MRI console. Subsequently, the 2D and 4D flow MRI data were pre-processed as follows (GTflow, Version 3.1.12, Gyrotools, Switzerland): 1.
The linear offset phase correction was conducted on each slice individually to correct for the presence of eddy currents. The fit was calculated at the reference heart phase (at the peak flow time) and then applied to all heart phases. A phase correction to compensate for concomitant gradients (Maxwell terms) and geometry correction to compensate for inhomogeneities of the main magnetic field and non-linearity of the gradient fields was performed on MR systems as part of the standard phase-contrast MR image reconstruction.

2.
Velocities in voxels outside of the flow lumen were nulled based on a magnitude intensity threshold. 3.
The data were inspected against phase-aliasing and manually corrected if necessary.
Next, the following processing steps were performed: 1.
ROIs were created manually on MRI magnitude data. First, contours around the tube's lumen were drawn using a b-spline curve (feature in GTflow) on 2D flow MRI. Note that the 2D flow MRI acquisition planes are already perpendicular to the flow direction. Next, the resulting 2D flow ROIs were translated to the 4D flow MRI data, ensuring identical placement of ROIs on the 2D and 4D flow datasets.

2.
In a given ROI, the flow of all voxels was summed up for each time point where i indicates the voxel and t the temporal point. The flow was spatially averaged over ROI A-C, as follows: The number of voxels per ROI diameter (nROI) was calculated to obtain a measurement not depending on the voxel size and vessel diameter as The time-dependent difference between flow values obtained with 4D and 2D flow . Similarly, the difference between flow values obtained with 4D flow MRI and US sensor was calculated.

5.
The normalized root-mean-square (RMS) error was used to assess the accuracy of flow quantification. RMS was calculated as the sum of squared differences between 2D and 4D flow MRI data over the time steps and normalized by the time-averaged flow acquired with 2D flow MRI: , where t indicates the temporal measurement point and Num is the number of temporal points. Similarly, RMS between flow values obtained with 4D flow MRI and US sensor wasc alculated. 6.
The velocity magnitude for each voxel in the aneurysm 3D ROI (Figure 1d, right) was calculated and averaged over time v , where i indicates voxel, t temporal point, and Num the number of temporal points. Aneurysm velocity distributions were presented using histograms similar to Schnell et al. [38]. 7.
Time-averaged velocity magnitude in the evaluation plane across the aneurysm was visualized pixel-wise on a color-coded representation (MATLAB R2019a, MathWork, Natick, MA, USA). 8.
Repeatability of flow measurements with 2D flow MRI was assessed with repeatability coefficient (RC) as follows: (1)  Here, we mostly aimed to assess flow in silicone tubes measured with 4D flow MRI against 2D flow MRI. A detailed comparison to the US sensor is presented in the paper supplement. In the main manuscript, we focus on the comparison with 2D flow MRI, because the location of 4D and 2D flow MRI were identical while the US sensor was located at a considerable distance between MR imaging volume. The distance between flow measurement positions has no impact on total flow at the tubes in the closed circuit. However, it might impact individual time-depended flow values, for example making the flow curve flatter, but keeping the total flow volume the same. As we want to assess time-depended flow curves, we limited the comparison of 4D flow MRI at the main text to the method that can be performed at the same location, namely 2D flow MRI here.

Statistical Analysis
Normality (Gaussian distribution) was tested for all variables (flow rate, flow difference) using the Shapiro-Wilk test (Tables S1 and S2) [40]. The median and range of variables were reported as summary statistics for all variables due to the severe skewness of the respective distributions. Spearman's rank correlation coefficient (rho) and linear fit coefficient were calculated to assess the similarity of flow rate values obtained by US sensor and 2D flow, US sensor and 4D flow MRI, 2D and 4D flow MRI measured in silicone tubes, and to assess the relationship between flow differences and RMS on nROI. A two-sided paired Wilcoxon signed-rank test was conducted for the flow rate values measured with a US sensor, 2D flow MRI, and 4D flow MRI in silicone tubes to assess the statistical difference between datasets. The same test was conducted to assess the significance of flow differences from 5 and 10% of median flow rate values measured with US sensor and 2D flow MRI. To assess the impact of the CS factor on RMS and median flow differences, the Kruskal-Wallis test was performed where CS was as a grouping factor. The impact of spatial resolution on velocity distribution in an aneurysm model was assessed with the Kruskal-Wallis test too, where voxel size was as a grouping factor. A two-sided Wilcoxon rank-sum test was conducted to compare velocity distribution obtained using CS = 2.5 with CS = 4.5 and 6.5 for each spatial resolution independently. The level of significance for all tests was set as p-value = 0.05. MATLAB R2020a (MathWork, Natick, MA, USA) was used for statistical analysis. The graphic representation was generated using Microsoft Excel 2016 (Microsoft Corporation, Redmond, WA, USA).

Flow in Silicone Tubes
The flow patterns obtained using the US sensor and 2D flow MRI were qualitatively similar to those obtained with 4D flow MRI for all voxel sizes and CS acceleration factors ( Figure 2).
First, as a sanity check, we compared flow measured with US sensor and the 2D flow MRI. Data were strongly correlated (rho = 0.96, Figure S2). Median flow difference was 12.94, −2.97, −2.90, and −6.76% for tube ID = 2-5 mm, respectively. Only for tube ID = 5 mm, median flow difference was significantly different from 0% (p = 0.03), but not different from 5% (p = 0.61). Therefore, the flow measured with US sensor and 2D flow MRI was in good agreement. Furthermore, we report the results of the comparison between 4D and 2D flow, while the quantitative comparison between flow obtained with a US sensor and 4D flow MRI is presented in a supplement to this paper (Figures S2 and S3,  Tables S3 and S4).
Quantitatively, a strong correlation (rho > 0.97) was observed between 2D and 4D flow MRI for all voxel sizes and CS factors ( Figure 3). A strong correlation (rho > 0.94) was also observed between US sensor and 4D flow MRI ( Figure S3). In addition, a linear fit was performed on 2D and 4D flow MRI flow values (R 2 > 0.96, Table 2). The linear slope was close to 1 for 0.5 and 1 mm 3 voxel size, regardless of the CS factor. For a voxel size of 1.5 mm 3 , the linear slope was close to 1 only for CS = 2.5. At higher acceleration factors (CS = 4.5-13), the flow values with 4D flow MRI were overestimated in comparison to 2D flow MRI (1.21 < slope < 1.29).     (Figure 4 left). For example, the flow median was statistically overestimated by voxel size equal to 1.5 mm 3 and CS = 2.5 and 4.5 for all tubes (p < 0.05). Summary statistics of flow rate values for each tube ID, voxel size, and CS factor are summarized in the supporting information (Table S3).
Flow difference between 4D and 2D flow ranged from -23.75% (tube ID = 4 mm, voxel size = 0.5 mm 3 , CS = 6.5) to +18.31% (tube ID = 2 mm, voxel size = 1.5 mm 3 , CS = 2.5). However, the median of flow difference calculated between 4D and 2D flow MRI was statistically higher than 10% only for voxel size of 1.5 mm 3 and CS = 2.5 for tube ID = 2, 4,  Overall, the flow was overestimated in all experiments by an average of 1.33 ± 8.31% and 0.24 ± 10.73% in comparison to 2D flow MRI and US sensor, respectively, for all voxel sizes and CS acceleration factors used for 4D flow MRI. A two-sided paired Wilcoxon signed-rank test revealed the statistical difference between flow values measured with 4D and 2D flow MRI for almost half of the used combinations of spatial resolution and CS factors (Figure 4 left). For example, the flow median was statistically overestimated by voxel size equal to 1.5 mm 3 and CS = 2.5 and 4.5 for all tubes (p < 0.05). Summary statistics of flow rate values for each tube ID, voxel size, and CS factor are summarized in the supporting information (Table S3). and 5 mm (p < 0.001), and for a voxel size of 0.5 mm 3 and CS = 6.5 for tube ID = 4 mm (p < 0.001, Figure 4 right). Summary statistics of flow differences for each tube ID, voxel size, and CS factor is summarized in supporting information (Tables S4 and S5). Flow values and flow differences were obtained by 4D flow MRI in comparison to US sensor presented in Figure S4. in statistically different flow median values (p < 0.001, star sign, left, a two-sided paired Wilcoxon signed-rank test). However, the difference was significantly higher than 10% only for four datasets (p < 0.001, red circle, right, a two-sided paired Wilcoxon signed-rank test).  ). However, the difference was significantly higher than 10% only for four datasets (p < 0.001, red circle, right, a two-sided paired Wilcoxon signed-rank test).
Flow difference between 4D and 2D flow ranged from -23.75% (tube ID = 4 mm, voxel size = 0.5 mm 3 , CS = 6.5) to +18.31% (tube ID = 2 mm, voxel size = 1.5 mm 3 , CS = 2.5). However, the median of flow difference calculated between 4D and 2D flow MRI was statistically higher than 10% only for voxel size of 1.5 mm 3 and CS = 2.5 for tube ID = 2, 4, and 5 mm (p < 0.001), and for a voxel size of 0.5 mm 3 and CS = 6.5 for tube ID = 4 mm (p < 0.001, Figure 4 right). Summary statistics of flow differences for each tube ID, voxel size, and CS factor is summarized in supporting information (Tables S4 and S5). Flow values and flow differences were obtained by 4D flow MRI in comparison to US sensor presented in Figure S4.
There might be an interplay between the size of the tube diameter and voxel size. Therefore, we constructed a dimensionless surrogate representing both voxel size and tube, termed here voxel number per vessel diameter (nROI). Increased nROI leads to a decreased median flow difference (Figure 5a). A weak correlation was observed between median flow differences and nROI (rho = −0.55. p < 0.01) when all CS factors were considered. A strong negative correlation was found for individual acceleration factors CS = 2.5 and 4.5 between median difference and nROI, rho = −0.81 (p = 0.02) and rho = −0.85 (p < 0.01), respectively, while no correlation was present for CS = 6.5 and 13 (rho < −0.44). Therefore, nROI plays the dominant role in the flow accuracy when CS ≤ 4.5, while when CS ≥ 6.5, some other factors might play a greater role. The trend of decreased median flow differences was observed with increased CS factor ( Figure 5b); however, no significant impact of CS factor was found according to the Kruskal-Wallis test (p = 0.18). ered. A strong negative correlation was found for individual acceleration factors CS = 2.5 and 4.5 between median difference and nROI, rho = −0.81 (p = 0.02) and rho = −0.85 (p < 0.01), respectively, while no correlation was present for CS = 6.5 and 13 (rho < −0.44). Therefore, nROI plays the dominant role in the flow accuracy when CS ≤ 4.5, while when CS ≥ 6.5, some other factors might play a greater role. The trend of decreased median flow differences was observed with increased CS factor (Figure 5b); however, no significant impact of CS factor was found according to the Kruskal-Wallis test (p = 0.18).
Next, we assessed RMS value to evaluate an absolute flow error. Overall, RMS was equal by an average of 0.18 ± 0.08 for all voxel sizes and CS acceleration factors used. RMS ranged from 0.07 (nROI = 4.8, CS = 4.5) to 0.38 (nROI = 1.66, CS = 4.5). As a result, we observed that increased nROI leads to decreased RMS values (Figure 5c). If more than 7 voxels are used, the RMS value tends to plateau. No correlation between RMS and nROI was found when the RMS values from all CS factors were considered (rho = −0.38, p = 0.007). However, for individual CS factors = 2.5 and 4.5, a strong negative correlation was found between RMS and nROI, rho = −0.63 (p = 0.03) and rho = −0.87 (p < 0.01), respectively. While for CS factors = 6.5 and 13, no correlation was found (rho < −0.15). Here, again nROI plays the dominant role in the flow accuracy when CS ≤ 4.5, while when CS ≥ 6.5, some other factors might play a greater role. The trend of increased RMS was observed with increased CS factor from 2.5 to 6.5 (Figure 5d). Statistical difference between CS = 2.5 and CS = 6.5 was found according to Kruskal-Wallis test (p = 0.01).
RMS and median flow difference depending on nROI calculated for 4D flow MRI in comparison to the US sensor is presented in Figure S5.  (d) RMS tended to increase with increased CS factor from 2.5 to 6.5. RMS obtained with CS = 2.5 was significantly different from CS = 6.5. For creating box plots in (b,d), RMS and flow difference values, obtained by 4D flow MRI with various nROIs but with the same CS factor, were joined together to make one dataset. Therefore, the dependency of data on nROI was excluded.
Next, we assessed RMS value to evaluate an absolute flow error. Overall, RMS was equal by an average of 0.18 ± 0.08 for all voxel sizes and CS acceleration factors used. RMS ranged from 0.07 (nROI = 4.8, CS = 4.5) to 0.38 (nROI = 1.66, CS = 4.5). As a result, we observed that increased nROI leads to decreased RMS values (Figure 5c). If more than 7 voxels are used, the RMS value tends to plateau. No correlation between RMS and nROI was found when the RMS values from all CS factors were considered (rho = −0.38, p = 0.007). However, for individual CS factors = 2.5 and 4.5, a strong negative correlation was found between RMS and nROI, rho = −0.63 (p = 0.03) and rho = −0.87 (p < 0.01), respectively. While for CS factors = 6.5 and 13, no correlation was found (rho < −0.15). Here, again nROI plays the dominant role in the flow accuracy when CS ≤ 4.5, while when CS ≥ 6.5, some other factors might play a greater role. The trend of increased RMS was observed with increased CS factor from 2.5 to 6.5 (Figure 5d). Statistical difference between CS = 2.5 and CS = 6.5 was found according to Kruskal-Wallis test (p = 0.01).
RMS and median flow difference depending on nROI calculated for 4D flow MRI in comparison to the US sensor is presented in Figure S5.
The effect of high CS factors = 6.5 and 13 is not clear based on our results, although the trend of changing flow differences and RMS values was observed with an increasing CS factor. The flow curves still were very similar to each other. However, severe distortions were apparent in the MRI magnitude, and phase images obtained with 4D flow MRI with CS = 13 ( Figure 6). The robustness of flow curves to the acceleration factor is likely due to a great difference in the voxel intensity originating from the silicone tube wall and lumen. Thus, the definition of the flow region of interest was not affected by signal distortion.
Tomography 2022, 8, FOR PEER REVIEW 12 increasing voxel number per diameter. (b) Median flow difference tended to decrease with an increased CS factor. (c) RMS tended to decrease with increasing voxel number per diameter. (d) RMS tended to increase with increased CS factor from 2.5 to 6.5. RMS obtained with CS = 2.5 was significantly different from CS = 6.5. For creating box plots in (b,d), RMS and flow difference values, obtained by 4D flow MRI with various nROIs but with the same CS factor, were joined together to make one dataset. Therefore, the dependency of data on nROI was excluded.
The effect of high CS factors = 6.5 and 13 is not clear based on our results, although the trend of changing flow differences and RMS values was observed with an increasing CS factor. The flow curves still were very similar to each other. However, severe distortions were apparent in the MRI magnitude, and phase images obtained with 4D flow MRI with CS = 13 ( Figure 6). The robustness of flow curves to the acceleration factor is likely due to a great difference in the voxel intensity originating from the silicone tube wall and lumen. Thus, the definition of the flow region of interest was not affected by signal distortion.

Figure 6.
A representative cross-section of the silicone tubes with a diameter of 4 and 5 mm that were imaged with 4D flow MRI. The sequence was accelerated factors of 2.5-, 4.5-, 6.5-, and 13 (from left to right), and an isotropic voxel size of 1 mm 3 was used. Note the severe image distortion when an acceleration factor of 13 is applied. Phase image with velocity encoding in foot to head direction, which is collinear with a flow direction, is shown.
It is important to note that on some 4D flow datasets, phase wrapping was observed, which was manually corrected within GTflow. The datasets affected with phase wrapping were the ones acquired with a spatial resolution of 0.5 mm 3 and CS = 2.5-13, 1.0 mm 3 and CS = 4.5 and 13, and 1.5 mm 3 and CS = 6.5. Flow calculated in ROI C at the tube with ID = 2 mm was most affected, while ROI A and B were almost intact.

Velocity in an Aneurysm Model
To give an overview of the flow pattern inside the aneurysm model, velocity streamlines were generated from 4D flow MRI data (Figure 7a). Blood enters in the distal part of the ostium, impinging on the aneurysm wall, while the velocity is continuously reduced until the blood reenters the parent vessel. The aneurysm has a regular shape, e.g., with no bulbs. Consequently, a single vortex forms, and a stagnation zone with very low-velocity values develops.
Flow in parental vessel and aneurysm was pulsatile at a time-averaged flow rate of 4.03 ± 1.34 and 2.99 ± 0.74 mL/s, respectively, as measured with 2D flow MRI ( Figure S6).

Figure 6.
A representative cross-section of the silicone tubes with a diameter of 4 and 5 mm that were imaged with 4D flow MRI. The sequence was accelerated factors of 2.5-, 4.5-, 6.5-, and 13 (from left to right), and an isotropic voxel size of 1 mm 3 was used. Note the severe image distortion when an acceleration factor of 13 is applied. Phase image with velocity encoding in foot to head direction, which is collinear with a flow direction, is shown.
It is important to note that on some 4D flow datasets, phase wrapping was observed, which was manually corrected within GTflow. The datasets affected with phase wrapping were the ones acquired with a spatial resolution of 0.5 mm 3 and CS = 2.5-13, 1.0 mm 3 and CS = 4.5 and 13, and 1.5 mm 3 and CS = 6.5. Flow calculated in ROI C at the tube with ID = 2 mm was most affected, while ROI A and B were almost intact.

Velocity in an Aneurysm Model
To give an overview of the flow pattern inside the aneurysm model, velocity streamlines were generated from 4D flow MRI data (Figure 7a). Blood enters in the distal part of the ostium, impinging on the aneurysm wall, while the velocity is continuously reduced until the blood reenters the parent vessel. The aneurysm has a regular shape, e.g., with no bulbs. Consequently, a single vortex forms, and a stagnation zone with very low-velocity values develops.
The flow was measured with 2D flow MRI experiments five times at parental vessel and aneurysm to test repeatability ( Figure S7). The time-averaged repeatability coefficient was equal to 8.6 ± 4.1% and 6.6 ± 2.6% for parental vessel and aneurysm, respectively.
The peak velocity and stagnation zone were detected by CFD, too, and were reproduced qualitatively by all combinations of voxel sizes and CS acceleration factors of 4D flow MRI (Figure 7b). The medians of velocity distributions measured in 3D aneurysm ROI (Figure 1d left) with 4D flow MRI among three voxel sizes, regardless of CS factor, were found to be close to each other-10.4, 10.1, and 9.7 cm/s for a voxel size of 0.5, 1.0, and 1.5 mm 3 , respectively (Figure 8a, left). However, the velocity distributions were statistically different (Kruskal-Wallis test, p-value << 0.01, Figure 8a, right). The decreasing median velocity value was observed for an increasing acceleration factor while keeping the voxel size the same (Figure 8b). In particular, (1) a voxel size of 0.5 mm 3 yielded a median velocity of 11.2, 10.5, and 9.4 cm/s for acceleration factors 2.5, 4.5, and 6.5, respectively; (2) a voxel size of 1.0 mm 3 resulted in a median velocity of 10.9, 9.8, and 9.8 cm/s; and (3) a voxel size of 1.5 mm 3 yielded a median velocity equal to 10.7, 9.7, and 8.9 cm/s. Note that the CS factor affected median values stronger than voxel size. For each spatial resolution, the velocity obtained  Figure S7). The time-averaged repeatability coefficient was equal to 8.6 ± 4.1% and 6.6 ± 2.6% for parental vessel and aneurysm, respectively.
The peak velocity and stagnation zone were detected by CFD, too, and were reproduced qualitatively by all combinations of voxel sizes and CS acceleration factors of 4D flow MRI (Figure 7b).
The medians of velocity distributions measured in 3D aneurysm ROI (Figure 1d left) with 4D flow MRI among three voxel sizes, regardless of CS factor, were found to be close to each other-10.4, 10.1, and 9.7 cm/s for a voxel size of 0.5, 1.0, and 1.5 mm 3 , respectively (Figure 8a, left). However, the velocity distributions were statistically different (Kruskal-Wallis test, p-value << 0.01, Figure 8a, right). The decreasing median velocity value was observed for an increasing acceleration factor while keeping the voxel size the same (Figure 8b). In particular, (1) a voxel size of 0.5 mm 3 yielded a median velocity of 11.2, 10.5, and 9.4 cm/s for acceleration factors 2.5, 4.5, and 6.5, respectively; (2) a voxel size of 1.0 mm 3 resulted in a median velocity of 10.9, 9.8, and 9.8 cm/s; and (3) a voxel size of 1.5 mm 3 yielded a median velocity equal to 10.7, 9.7, and 8.9 cm/s. Note that the CS factor affected median values stronger than voxel size. For each spatial resolution, the velocity obtained with CS = 4.5 and 6.5 were compared against the one obtained with CS = 2.5. Statistical difference was observed for all cases (Wilcoxon rank-sum test, p << 0.01, Table S7). However, the changes did not exceed 11% for CS = 4.5 and 17% for CS = 6.5 in comparison to CS = 2.5 (Table S7). In addition, the median of velocity distribution obtained using 4D flow MRI with voxel size 0.5 mm 3 and CS = 2.5 was compared to all other combinations of CS factors and voxel sizes. In this case, the difference between medians ranged from −20.82% (voxel size = 1.5 mm 3 , CS = 6.5) to −3.09 % (voxel size = 1 mm 3 , CS = 2.5). Overall, the median velocity obtained with voxel size 0.5 and CS = 2.5 was underestimated by average on −11.42 ± 6.1% in comparison to all other combinations of voxel size and CS factors used. with CS = 4.5 and 6.5 were compared against the one obtained with CS = 2.5. Statistical difference was observed for all cases (Wilcoxon rank-sum test, p << 0.01, Table S7). However, the changes did not exceed 11% for CS = 4.5 and 17% for CS = 6.5 in comparison to CS = 2.5 (Table S7). In addition, the median of velocity distribution obtained using 4D flow MRI with voxel size 0.5 mm 3 and CS = 2.5 was compared to all other combinations of CS factors and voxel sizes. In this case, the difference between medians ranged from −20.82% (voxel size = 1.5 mm 3 , CS = 6.5) to −3.09 % (voxel size = 1 mm 3 , CS = 2.5). Overall, the median velocity obtained with voxel size 0.5 and CS = 2.5 was underestimated by average on −11.42 ± 6.1% in comparison to all other combinations of voxel size and CS factors used.

Discussion
In previous studies, 4D flow MRI has been compared with numerous invasive and noninvasive methods [22,[24][25][26][27]. However, the studies systematically addressing the impact of 4D flow MRI acquisition parameters on flow quantification are limited. In a recent review [41], the authors performed a systematic analysis on the comparability of the flow volume and peak velocity measured with 4D flow to 2D flow MRI in heart and great vessels. As a result, they found that 4D flow results are in the agreement with 2D flow in only 55% of the manuscript considered. The authors emphasize that 4D flow acquisition parameters have a great impact on flow accuracy, and each center needs to perform its quality assurance before using 4D flow MRI for clinical needs. In our work, we wanted to contribute to this awareness too. Thus, we aimed to assess the effect of spatial resolution and CS acceleration factors on flow quantification in small vessels. Firstly, a strong correlation (rho > 0.97) and linear relationship (R2 > 0.96) were observed between 2D and 4D flow MRI for all the various CS acceleration factors and spatial resolutions considered in the silicone tubes. This finding is in agreement with previously reported values. When comparing 4D flow with 2D flow MRI, transonic flow-sensor or Doppler US, a strong correlation was reported in assessing flow in a pulsatile phantom (rho = 0.96) [22], in healthy human volunteers (rho > 0.95) [42], in a swine aorta (rho = 0.86) [43], and in canine arteries (rho = 0.95 [28] and rho = 0.89 [24]).

Effect of Spatial Resolution and MR Acceleration on the Flow in Silicone Tubes
It is well known that the accuracy of phase-contrast imaging in small vessels is lower than in large ones because the relative number of voxels at the vessel's boundaries increases with decreasing vessel diameter [44]. Overall, flow values between 2D and 4D flow MRI were close, and a linear slope ranged from 0.93 to 1.09 for the smaller voxel sizes of 0.5 and 1 mm 3 . However, flow was overestimated in comparison to 2D flow MRI for a greater voxel size (1.5 mm 3 ), where the linear slope varied from 1.05 to 1.29. The flow was overestimated in all experiments by an average of 1.33 ± 8.31% and 0.24 ± 10.73% in comparison to 2D flow MRI and the US sensor, respectively, for all voxel sizes and CS acceleration factors used for 4D flow MRI. Flow was overestimated by 16%, as reported previously in a swine study (ID at different segments = 11.6, 16.55, and 18.1 mm, voxel size: 2.1 × 1.7 × 2.8 mm 3 , nROI = 5, 7, and 8) [43]. In another study, the flow was overestimated by 5% in a pulsatile phantom (ID = 20 mm) [22] compared to the transonic sensor. In some cases, the difference was eliminated when a high resolution was used (voxel size: 1.9 and 1.5 mm 3 instead of 2.8 × 2.8 × 2.2 mm 3 , nROI = 10 and 13 instead of 7, respectively) [22]. However, another study reports that no statistical difference was found for flow estimated with 2D flow MRI obtained with 1 × 1 mm 2 and 2 × 2 mm 2 (single-slice acquisition) compared to the flow sensor (ID = 3.5 ± 0.7 mm, nROI = 2-4 instead of 1-2) [28]. The spatial resolution itself without a context on vessel diameters does not allow to estimate the flow quantification error. It might be no effect of spatial resolution seen because it was already sufficient resolution and its increase did not change anything, or it might not be sufficient in both cases, and the flow is calculated with error in both cases. In the recent review [41], the authors showed that the spatial resolution did not have a statistical impact on the agreement between 2D flow and 4D flow MRI; likely some other parameters had a stronger impact, or the spatial resolution was sufficient to eliminate an effect on the flow quantification error.
The flow estimation error is the greatest on the boundary voxels and, therefore, it depends on the number of boundary voxels N boundary ∼ = π·D/voxel size [44], where D corresponds to the diameter of the ROI. In other words, the error depends linearly on the number of voxels per vessel diameter (nROI). If this value is small, then it will limit the accuracy of the flow determination. We observed that in areas with less than six voxels per diameter, the RMS value decreases as the number of voxels increases. However, if more than seven voxels are used, the RMS value tends to plateau. The difference between median flow values obtained with 2D and 4D flow MRI (CS =2.5) was up to 20% for nROI in the range 2-5 voxels per diameter and less than 5% for nROI in the range 7-10 voxels per diameter. Previous studies have shown that an insufficient number of voxels per vessel diameter might lead to blood flow overestimation. Aristova et al. [45] investigated the accuracy of flow measurement with dual-venc 4D flow MRI accelerated with PEAK-GRAPPA (R = 2, 3, and 5) in comparison to 4D flow accelerated with GRAPPA (R = 2) using in vitro channels (ID = 4, 6, and 8 mm). In that study, flow measurement accuracy was 10-15% when nROI = 5, and less when more than five voxels were used, which is in accordance with our results. Another study reported simulated flow error to be less than 4% (nROI = 10), 20% (nROI = 4), and up to 60% (nROI = 2) [44]. In the same study, experimental MRI flow error was less than 9% (nROI = 10) and 20% (nROI = 5). In another study performed by Bouillot et al. [46], flow was overestimated by up to 120% (nROI = 1.3) and 60% (nROI = 2). However, in a different study, a weak trend of decreasing the flow error was apparent when more than four voxels per diameter were used [28]. Fewer voxels per diameter render flow determination less accurate. This effect must be considered when estimating the flow in small vessels. In addition to the MRI spatial resolution, it is desirable to report the size of the vessel diameter, too, or the number of voxels per diameter used for calculating flow.
Increasing spatial resolution might offer one solution for providing more voxels and better flow estimation. However, that cannot always be achieved easily and is compromised by longer measurement times and reduced signal-to-noise ratio (SNR). Using post-processing pipelines has been proposed as an alternative solution to decrease the error associated with boundary voxels. Recently, Boulleot et al. [46] suggested combining velocity information measured by 4D flow MRI with vessel geometry data measured with TOF or 3D rotational angiography to correct for insufficient spatial resolution. Globally, such correction led to 10-15% lower flow rates, and low sensitivity to the spatial resolution was observed when >2 voxels per diameter were used [46]. Another post-processing algorithm is based on the assumption that low flow at the boundaries induced small-phase shifts [47]. This algorithm reduced error from 36% without correction to less than 10% with correction. One more algorithm based on the predefined model-based method assumes a parabolic flow profile and cylindrical vessel geometry. It results in the reduction of errors from more than 50% to less than 10% (nROI > 5) [48]. Still, it is worth noting that even small uncertainties in flow definition can lead to a large variation in the determination of WSS when, for example, PC MRI flow is used to estimate aneurysm rupture risk with patient-specific CFD simulations [49].
Additionally, advanced accelerating imaging techniques can improve spatial resolution by reducing the examination time. Therefore, we also assessed the influence of the acceleration factor on flow measurements. No clear effect of the acceleration factor was found on the correlation coefficient, the linear slope. However, a trend of decreasing flow difference error and increasing RMS with increased CS factors was observed, although a small acceleration factor peak flow is better assessed for larger voxel sizes (1.0 and 1.5 mm 3 ). Moreover, the smallest linear slope (1.05) was yielded when 1.5 mm 3 voxels were used. Additionally, a higher acceleration factor (CS = 6.5) resulted in statistically higher RMS values than the ones obtained with a smaller acceleration factor (CS = 2.5); however, no statistical difference was found between CS = 2.5 and CS = 4.5. Additionally, the dependence of flow error from nROI was observed only for CS = 2.5 and 4.5; for CS = 6.5 and 13, some other factors likely play the dominant role. Previously, the feasibility of using acceleration factors of 4 to 10 [50] and even 10 to 30 [51] was demonstrated for in vivo aortic applications. Therefore, the use of higher CS factors might be appropriate. Yet, in the present study, an acceleration factor of 13 resulted in severe image distortions. No significant impact of CS = 13 on flow quantification was seen here likely due to the consistent placement of ROIs used for flow quantification. Still, peak flow was previously underestimated in the order of 25% [51] for an acceleration factor of 20, and 4% and 5% for acceleration factors of 8 and 13, respectively [42]. However, previously, Aristova et al. [45] showed no significant relationship between acceleration factor and flow measurement error for dual-venc 4D flow MRI accelerated with PEAK-GRAPPA. Overall, the number of studies considering the effect of acceleration of flow quantification even in the big vessels is limited. Moreover, it is worth noting that according to the review on 4D flow MRI in the heart and great vessels [41], most of the clinical centers do not use any acceleration technique at all. Recently, Pathrose et al. [52] studied CS 4D flow acquisitions in patients with aortic disease in comparison to 4D flow accelerated with GRAPPA. On average, the net flow was underestimated in the aortic arch by −7.1 ± 10.5%, −6.8 ± 13.8%, and −5.0 ± 10.7% for CS = 5.7, 7.7, and 10.2, respectively. However, at descending aorta, it was underestimated by −5.3 ± 11.8% only when CS = 10.2 was used. Moreover, at the ascending aorta, no significant impact of CS factor was found. Moreover, for some patients, the flow values were in good agreement between all 4D flow acquisition protocols. That being said, likely some other factors lead to the flow quantification error in combination with acceleration. In our study, on the contrary, flow overestimation was observed with increased CS factor. We hypothesize that the accelerated 4D flow acquisition is more capable of measuring more accurately peak flow values, and thus, better agreement between flow measured with CS = 6.5 and flow sensor was observed in Figure 2.

Velocity in an Aneurysm Model
As a neurovascular application of 4D flow MRI, we assessed the velocity field in the intracranial aneurysm model. Aneurysms may rupture, and a potentially fatal subarachnoid hemorrhage can develop [53]. Thus, aneurysms at risk of rupture need to be identified and treated. Some current options for treating aneurysms are based on flow disruption. Four-dimensional flow MRI has already been applied after aneurysm treatment to assess the blood flow changes [12,13]. Therefore, accurate assessment of aneurysm hemodynamics in vivo is desirable. Several studies compared 4D flow MRI velocity maps in the aneurysm to CFD and PIV [14,27,54]. However, studies addressing the impact of 4D flow MRI parameters on aneurysm flow are limited. In the present study, a good qualitative reproduction of time-averaged velocity magnitude maps calculated by CFD was achieved by all values of voxel sizes and acceleration factors. The main features of flow patterns, such as stagnation and high-velocity zones, were revealed with all voxel sizes and CS factors. Moreover, the velocity datasets measured with different voxel sizes yielded very close values of median velocity over the aneurysm sac. However, the larger voxel sizes resulted in smaller velocity values. This agrees with previous data by Roloff et al. [27], who reported underestimation of averaged velocity magnitude on 7.5 cm/s by 4D flow MRI (isotropic voxel size 0.57 mm 3 ) compared to tomographic PIV (isotropic voxel size 0.23 mm 3 ). However, when the PIV was downsampled to 4D flow MRI voxel size, the difference was reduced to 0.02 cm/s. Furthermore, further downsampling of the PIV data (isotropic voxel size 0.83 mm 3 ) resulted in an average velocity magnitude reduction of 9.8 cm/s compared to the original data. The large voxel size likely represents a major factor in underestimating the velocity magnitude, which might be independent from the measurement technique. In another study, significant variation between velocity fields was observed in a comparison of 4D flow MRI (isotropic voxel size 0.5 mm 3 ) and PIV (isotropic voxel size 0.17 mm 3 ), which was eliminated after downsampling the spatial resolution of PIV [54].
The overall difference in velocity maps and median values acquired with various voxel sizes was relatively small in our study. Still, owing to the high-resolution secondary flow features, such as recirculation flow areas, can be captured, as discussed in Medero et al. [54]. Additionally, velocity-derived parameters such as WSS and OSI might also benefit from high resolution [55], as well as pressure [56]. The aneurysm considered in this study had a relatively simple shape, with no bulbs or aneurysm daughters. Therefore, only one vortex flow in the center of the aneurysm sac was present and well captured even by the coarse voxel size. The aneurysm investigated in our study was relatively large at 17.5 mm in the maximum dimension, resulting in 35, 18, and 12 voxels per maximum dimension using voxel sizes of 0.5, 1, and 1.5 mm 3 , respectively. Thus, spatial resolution was likely already sufficient for all voxel sizes. For smaller aneurysms with additional morphological features, a higher resolution is recommended.
The varying acceleration factor resulted in very close, but significantly different median velocity values for almost all cases. Additionally, the acceleration factor has a greater effect on velocity distribution than voxel size, likely owing to the relatively large aneurysm size. Moreover, the increased acceleration factor tended to decrease velocity. However, the changes did not exceed 11% for CS = 4.5 and 17% for CS = 6.5 in comparison to CS = 2.5, and no statistical difference was reported between the velocity field obtained with 4D flow MRI using another acceleration scheme [57]. The reduction in examination time gained by using a high acceleration factor can be used to improve spatial and temporal resolution. Here, we did not vary the temporal resolution; however, it may impact peak flow and velocity-derived parameters such as OSI [58,59].

Limitation
A limited number of models were considered in the present study: silicone tubes and one patient-specific aneurysm model. Other configurations of aneurysms or vessels might yield different results. Furthermore, this is an in vitro study only. It would be beneficial to vary spatial resolution and acceleration factor in vivo as well. In addition, some 4D flow MRI data were subjected to manual phase-wrapping correction, and some data could not be fully recovered, which may introduce errors in the flow estimation. The flow measured with 4D flow MRI data was compared to that measured with 2D flow MRI, representing another limitation of the study. Indeed, 2D flow MRI does not constitute a ground truth method, as it uses non-isotropic voxel sizes and non-real-time data collection (cine approach). Still, flow measured with 2D flow MRI was in good agreement with the flow values obtained with a flow sensor.

Conclusions
The spatial resolution affects the flow rate estimated by using 4D flow MRI. By coarse voxel size, flow is overestimated. The use of fewer than seven voxels per vessel diameter overestimated flow in more than 5% of flow measured with 2D flow MRI. Furthermore, we encourage reporting the inner diameter of the investigated vessel in addition to the spatial resolution of 4D flow MRI. This will make it possible to estimate the possible flow error. The velocity error in the aneurysm sac did not exceed 11% for CS = 4.5 in comparison to the velocity obtained using 4D flow MRI with CS = 2.5. Therefore, CS factors from 2.5-4.5 can be used to improve the spatial or temporal resolution for accurate quantification of flow rate and velocity.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/tomography8010038/s1, Figure S1: Time-resolved flow rate curves measured with US-sensor three times each 90 minutes. Slight variations are present in the curves indicating good temporal stability of the circulating system.; Figure S2: Scatter plots and Spearman-rank correlation coefficient between time-resolved net flow rate values obtained by 2D flow and US sensor. A correlation coefficient was 0.96.l; Figure S3: Scatter plots and Spearman-rank correlation coefficient between time-resolved net flow rate values obtained by 4D flow US sensor in silicone tubes. A strong correlation was observed for all voxel sizes and acceleration factors (rho > 0.94).; Figure S4: Box plots representing flow values obtained with 4D flow MRI and flow differences calculated between 4D flow MRI and US sensor.; Figure S5: Scatter and box plots for experiments performed on silicone tubes.; Figure S6: Time-resolved flow rate curves measured with 2D flow and 4D flow MRI in patient-specific aneurysm model.; Figure S7: Time-resolved flow rate curves measured with 2D flow MRI at ROI 1-3 (a top) of an aneurysm model.; Table S1: Results of Shapiro-Wilk test for normality of experimental data: flow rate values.; Table S2: Results of Shapiro-Wilk test for normality of experimental data: flow differences.