- 1Department of Medical Physics, Faculty of Mathematics and Natural Sciences, University of Greifswald, Greifswald, Germany
- 2Department of Radiology, Northwestern University, Feinberg School of Medicine, Chicago, IL, United States
- 3Department of Neurology, Northwestern University, Feinberg School of Medicine, Chicacgo, IL, United States
- 4Department of Experimental Physics V, University of Wuerzburg, Wuerzburg, Germany
- 5Department of Diagnostic Radiology, University Hospital of Greifswald, Greifswald, Germany
Introduction: Intracranial 4D flow MRI enables quantitative assessment of hemodynamics in patients with intracranial atherosclerotic disease (ICAD). However, quantitative assessments are still challenging due to the time-consuming vessel segmentation, especially in the presence of stenoses, which can often result in user variability. To improve the reproducibility and robustness as well as to accelerate data analysis, we developed an accurate, fully automated segmentation for stenosed intracranial vessels using deep learning.
Methods: 154 dual-VENC 4D flow MRI scans (68 ICAD patients with stenosis, 86 healthy controls) were retrospectively selected. Manual segmentations were used as ground truth for training. For automated segmentation, deep learning was performed using a 3D U-Net. 20 randomly selected cases (10 controls, 10 patients) were separated and solely used for testing. Cross-sectional areas and flow parameters were determined in the Circle of Willis (CoW) and the sinuses. Furthermore, the flow conservation error was calculated. For statistical comparisons, Dice scores (DS), Hausdorff distance (HD), average symmetrical surface distance (ASSD), Bland-Altman analyses, and interclass correlations were computed using the manual segmentations from two independent observers as reference. Finally, three stenosis cases were analyzed in more detail by comparing the 4D flow-based segmentations with segmentations from black blood vessel wall imaging (VWI).
Results: Training of the network took approximately 10 h and the average automated segmentation time was 2.2 ± 1.0 s. No significant differences in segmentation performance relative to two independent observers were observed. For the controls, mean DS was 0.85 ± 0.03 for the CoW and 0.86 ± 0.06 for the sinuses. Mean HD was 7.2 ± 1.5 mm (CoW) and 6.6 ± 3.7 mm (sinuses). Mean ASSD was 0.15 ± 0.04 mm (CoW) and 0.22 ± 0.17 mm (sinuses). For the patients, the mean DS was 0.85 ± 0.04 (CoW) and 0.82 ± 0.07 (sinuses), the HD was 8.4 ± 3.1 mm (CoW) and 5.7 ± 1.9 mm (sinuses) and the mean ASSD was 0.22 ± 0.10 mm (CoW) and 0.22 ± 0.11 mm (sinuses). Small bias and limits of agreement were observed in both cohorts for the flow parameters. The assessment of the cross-sectional lumen areas in stenosed vessels revealed very good agreement (ICC: 0.93) with the VWI segmentation but a consistent overestimation (bias ± LOA: 28.1 ± 13.9%).
Discussion: Deep learning was successfully applied for fully automated segmentation of stenosed intracranial vasculatures using 4D flow MRI data. The statistical analysis of segmentation and flow metrics demonstrated very good agreement between the CNN and manual segmentation and good performance in stenosed vessels. To further improve the performance and generalization, more ICAD segmentations as well as other intracranial vascular pathologies will be considered in the future.
1 Introduction
Intracranial 4D flow magnetic resonance imaging (MRI) is a promising imaging modality enabling 3D visualization and quantification of blood flow values (1), and flow-related parameters (2). Previous studies already demonstrated that this phase-contrast (PC) technique can be successfully applied to a variety of pathologies, for example, to explore the hemodynamic alterations due to aneurysms (3), cerebral arteriovenous malformations (4) and intracranial atherosclerotic disease (ICAD) (5). The quantitative analysis of intracranial 4D flow MRI, however, still poses several practical challenges due to its complexity and time-consuming manual 3D segmentation required for quantification, especially in the presence of pathologies. Besides the complicated intracranial vessel geometry, structural and morphological changes due to atherosclerotic plaque formation can aggravate manual segmentation, thus leading to low reproducibility. For example, vascular stenoses can lead to flow artifacts and signal loss, hampering an accurate segmentation of the vessel. To improve the accuracy, reproducibility, and robustness of the analysis of hemodynamic parameters and to accelerate data analysis, an accurate, automated segmentation algorithm for stenosed intracranial vessels is required. A large variety of techniques have been developed to address the problem of semi or fully-automatic vessel segmentation (6). While previous segmentation approaches already drastically improve temporal efficiency in comparison to manual segmentations, they still often require manual labor and lack robustness and consistency, therefore often requiring user interactions from the technologist (7).
With the recent rise of deep learning and convolutional neural networks (CNN), new algorithms have been proposed, promising more reliable and less user-dependent vascular segmentation (8). In particular, the introduction of the U-NET and its variants (9, 10) led to a broad range of new techniques for vessel segmentation using clinical vessel imaging techniques, already achieving very good agreement in comparison to manual segmentations performed by radiologists (7, 11). For example, U-NET was successfully applied for segmentations of cerebral vessels in time-of-flight (TOF) magnetic resonance angiography (MRA) images in patients with cerebrovascular disease (12) and for digital subtraction angiography (DSA) images in patients with intracranial aneurysms (13).
However, up to now, most deep learning-based segmentation approaches for intracranial vessels use TOF or contrast-enhanced (CE) MRA images, as well as DSA, or computed tomography (CT) angiography images while there are no approaches based on 4D flow MRI to the knowledge of the authors at the time of writing this manuscript. For intracranial 4D flow MRI applications, non-deep learning algorithms such as centerline processing schemes (5, 14) or using the standard difference of mean velocity (15) were proposed. The use of 4D flow MRI for deep-learning-based segmentations, however, would have several advantages in comparison to other imaging and segmentation modalities: First, no registration is required to spatially match the segmentation of a different imaging modality with the 4D flow MRI measurement, which can be computationally expensive and prone to errors. Secondly, techniques such as dual-VENC-4D flow MRI enable the assessment of morphological and functional information of the complete vascular tree including both the arteries and the veins in a single measurement (1) while with TOF usually, two separate measurements (angio- and venograms) are necessary. Finally, the use of deep learning has the potential to reduce the aforementioned dependency on user interactions from the technologist.
Recently, Berhane et al. proposed a U-NET-based convolutional neural network technique for the automated segmentation of the aorta using phase-contrast angiography (PCMRA) images derived from aortic 4D flow MRI (16). Here, the automated segmentation achieved an excellent agreement with manual segmentations, however, its use for intracranial 4D flow MRI and its performance in pathology such as intracranial atherosclerotic stenosis still needs to be investigated.
Therefore, in this study, the neural network developed for aortic 4D flow measurements (16) was re-trained for the automated vessel segmentation of intracranial 4D flow MRI in healthy controls and intracranial atherosclerotic stenosis patients. To assess possible differences in segmentation performance between cases with and without disease, the results were compared with segmentations of healthy controls. For performance assessment, Dice score, Hausdorff distance, and average symmetrical surface distance were computed using two independent manual observers as reference. In addition, parameters such as peak velocity, flow rate, and flow conservation error were computed and compared to the manual analyses. Furthermore, the segmentation performance in stenosed vessels was analyzed and compared with segmentation results based on black blood vessel wall imaging (VWI).
2 Methods
2.1 Study cohort
As part of a clinical ICAD protocol at Northwestern Memorial Hospital, 4D flow MRI scans were acquired in ICAD patients. 35 cases expressed severe stenosis with >70% constriction, 25 cases expressed moderate stenosis with >50% and <70% constriction, and 5 cases had mild stenosis with <50% constriction. Additional 3 cases didn't have a significant stenosis. The data acquired between 2014 and 2022 were retrospectively selected (n = 68, n = 30 women) for this institutional review board (IRB) approved study. All ICAD-related stenoses were confirmed using the clinical electronic medical record, MRI/MRA, and MR vessel wall imaging review by two interventional neuroradiologists (RA, SAA).
In addition, 4D flow MRI data of healthy volunteers (n = 86, n = 43 women) was included in this study. Informed consent was obtained from all volunteers. An overview of all patients and volunteers can be found in Table 1 (White background: all cases. Shaded background: Testing cohorts only).
Table 1. Sex, median age, BMI, and heart rate (max and min values in brackets) for the control and ICAD group, respectively (top: all cases. Bottom: testing cases).
2.2 MRI
Patients were scanned using a clinical MRI protocol including 4D flow MRI and VWI with 3D-T1-SPACE (17) (sampling perfection with application-optimized contrasts using different flip angle evolutions). In addition, 4D flow MRI scans were acquired in healthy volunteers. The VWI parameters were: Spatial resolution 0.52 mm × 0.52 mm × 0.70 mm, TR = 800 ms, TE = 23 ms. The relevant scan parameters for the 4D flow scans can be found in Table 2. For both cohorts kt-GRAPPA accelerated (R = 5) dual-VENC 4D Flow MRI was utilized (1). The field of view (FOV) was positioned to cover the circle of Willis (CoW), including the basilar artery (BA), left and right internal carotid arteries (ICA), middle cerebral arteries (MCA), anterior cerebral arteries (ACA), the posterior cerebral arteries (PCA), the posterior communicating arteries (PCOM), the superior cerebral arteries (SCA), the vertebral arteries (VA) as well as the superior sagittal sinus (SSS), straight sinus (STR) and left and right transverse sinus (TS). All measurements were performed at 3 T using a Prisma Fit or Skyra (both Siemens Healthineers Inc., Erlangen, Germany).
2.3 Post-processing
A custom-built MATLAB (The MathWorks, Natick, USA) tool was used for Eddy current correction, noise masking, and anti-aliasing of the phase difference images (1, 18). Phase-contrast MRA images (PCMRA) were calculated using the pseudo-complex difference method (Equations 1, 2) (1, 3):
with
Here, denotes the magnitude images derived from the 4D flow measurement, i the index and N the total number of cardiac phases (1, 3).
For the creation of training and validation data, the PCMRA images were manually segmented in MIMICS (Mimics, Materialise, Belgium). This was achieved by applying a threshold to remove noisy voxels. Subsequently, the neurovascular architecture was identified using the region-growing tool in MIMICS to select areas of intracranial vessel voxels. Noisy voxels captured with the region-growing processes were manually removed from the segmentation. All cases were subsequently edited by a second investigator of more than 10 years of experience (PW). This second step was to achieve consensus segmentations so that these Observer 1 segmentations can serve as “ground truth”. Due to the lack of availability of the commercial software MIMICS, the second user changed to the open-source software 3DSlicer (Slicer 5.2.2., SlicerCommunity) ensuring repeatability. All initial segmentations were performed by operators with at least 2 years of experience (MA, JM, AR).
2.4 Automated segmentation using convolutional neural networks
For the automated segmentation, the CNN developed by Berhane et al. was used (16). The network consists of a 3D U-Net [(10), see Figure 1A]. The original convolution layers were replaced by dense blocks (19), as described previously (16). For each dense block, batch normalization, a linear rectifier unit (ReLu), a 3D convolution (3 × 3 × 3), and a dropout layer (dropout rate 0.1) were computed. For the training, the calculated PCMRA images [see Equation (1)] were center-cropped or padded to obtain a fixed dimension of 224 × 192 × 64 and used as input for the CNN. No patching of the data was used. Instead, center-cropping was applied to reduce dimensions, since the intracranial vessels are always located around the center. Furthermore, no data augmentation was performed since the intracranial 4D flow scans are always acquired in the same orientation. To increase efficiency, all prior feature maps were concatenated and used as inputs for the subsequent layers (16). In the encoding part of the U-Net (see left-hand side of Figure 1A), a max-pooling layer is applied for downsampling while transposed convolution is used for upsampling in the decoding part (see right-hand side of Figure 1A). In the final layer, a 1 × 1 × 1 convolution and a softmax function are applied. The last step generates a binary value for each voxel of the input image (0: background, 1: foreground). Segmentation masks were created by selecting the class with the highest probability per voxel. For the training, a composite loss function (softmax-cross entropy and Dice loss), a batch size of 1, a learning rate of 0.0001, and 300 epochs were used. All computations were performed in Python 3.6.13 (Python Software Foundation, Beaverton, OR) with Tensorflow 2.4.0 on a 13th Gen Intel Core i7-13700 (2,100 MHz, 16 Cores) CPU with an NVIDIA GeForce RTX 4070 Ti GPU with 16 GB VRAM. A cohort of 134 randomly selected cases (76 controls, 58 ICAD patients) were used for the training while the remaining 20 cases (10 controls, 10 ICAD patients: four with severe, four with moderate, one with mild, one without significant stenosis. See Table 1 for group statistics) were used for testing (see Figure 1B).
Figure 1. (A) Layer structures of the CNN. A symmetrical design is used based on the 3D U-Net architecture. Different from the original approach, dense blocks are implemented into each layer. Dense blocks enable the regulation of the growth of the CNN while efficiently applying feature maps extracted through the CNN by using concatenation after each convolution layer. Dense blocks consist of the serial application of batch normalization, activation with a linear rectifier unit (ReLu) a 3 × 3 × 3 convolution, and a dropout layer (dropout rate 0.1). (B) Chart illustrating the splitting of the patient and control data into the training and validation datasets.
2.5 CNN performance analysis
To compare the performance between the manual and automated segmentation, the Dice score (DS), Hausdorff distance (HD), and average symmetrical surface distance (ASSD) were calculated using (Equations 3–5) (20):
Here, X and Y are binary segmentation masks and d is the Euclidian distance between both segmentations. For the segmentation analysis, DS, HD, and ASSD were computed for:
(a) CoW and sinuses.
(b) Only the CoW.
(c) Only the sinuses.
The calculation of segmentation metrics was performed in MATLAB (DS) and Python (HD and ASSD). For visual presentations of the segmentation masks, Ensight 10.02 (CEI, inc., USA) was used.
2.6 Flow analysis
Magnitude, velocity, and segmentation data were imported into a semi-automatic MATLAB analysis tool (5). First, centerlines were created automatically and perpendicular analysis planes were placed equidistantly along the vessels at a 0.25 mm distance. Subsequently, lumen cross-sectional areas, peak velocity, and flow rates were extracted for all analysis planes. Planes close to branches and bifurcations were excluded to avoid systematic errors in the flow estimation. The same analysis planes were used for both the manual and the automated segmentations. For statistical comparisons between manual and automated segmentations, the median cross-sectional area, peak velocity, and temporally averaged flow rate values were calculated over all analysis planes for each vessel of interest. The vascular analysis was subdivided into:
(a) large arteries (BA, ICA, MCA),
(b) small arteries (ACA, PCA, PCOM, SCA, VA),
(c) sinuses (TS, STR, SSS).
In addition, the internal consistency of the flow rate was assessed by determining the flow conservation error (fce) for the arteries (Equation 6) (4):
2.7 Performance analysis in stenosed vessels
For an analysis of the segmentation performance in stenosed vessels, one ICAD patient with moderate stenosis (>50% constriction in the right MCA, male, 80 years old) and two patients with severe stenosis (>70% constriction in the right MCA, female, 68 years old, >70% constriction in the right ICA, female, 61 years old) were additionally validated with segmentations obtained from black blood VWI.
For comparison with VWI, rigid registration was applied using the SPM12 MATLAB tool box (21) to align the 3D-T1-SPACE images with the 4D flow images. Subsequently, 3D volume analysis of the black blood images was performed using a home-built 3D framework (22, 23). Based on an interactive specialized Dijkstra algorithm (24), the centerline, vessel volume and vessel wall were extracted and visualized. Beforehand, the black blood images were processed in multiple steps: First, a prior median filter was applied to reduce noise and to enhance the contrast for the Dijkstra-searching algorithm, which was introduced by manually set seed points. In the next step, the volume of the stenosed vessel was extracted along the centerline. In addition, a vertex model of the desired vessel structure was generated using a Marching Cube algorithm (25). Subsequently, the volume was imported to MATLAB. Since the SPM12 co-registration was not perfect and a few voxels off, a second rigid registration using the FLIRT (flexible image reconstruction toolbox) MATLAB tool box (26) was applied to co-register the VWI segmentation with the 4D flow segmentations.
Using the VWI segmentation as reference, DS, HD and ASSD were calculated for both manual segmentations and the automated segmentation. In addition, the lumen cross-sectional area profiles determined with the 4D flow segmentations were compared with segmentations obtained from the T1-SPACE images. Within a region of interest around the stenosis, the time-resolved median flow rates were determined for the automated and the Observer 1 and Observer 2 segmentations, respectively.
2.8 Interobserver study
For interobserver comparisons of standard manual segmentations, all datasets from the testing cohort (10 controls, 10 ICAD) were segmented by an additional observer (Observer 2) with medical background and one year of experience with 4D flow MRI, who was not part of the original segmentation instance (JW). The second observer was blinded to the original segmentations of the testing cases. All segmentations were performed in 3DSlicer and the segmentations were verified by a second investigator (PW) to avoid systematic errors with the segmentation process. Using the original manual (Observer 1) and the automated segmentation as a reference, DS, HD, and ASSD were computed and compared with the results from the automated segmentation vs. Observer 1 analysis. In addition, flow parameters and cross-sectional areas were compared using the same analysis planes as for the CNN segmentation and the segmentation performance was analyzed in the three stenosis cases described in section 2.7.
2.9 Statistical analysis
Cross-sectional areas, as well as flow metrics (peak velocity, flow rate, fce) derived from the automated segmentation and the second observer, were compared with results from the Observer 1 segmentation using correlation and Bland-Altman analysis. Using the Observer 1 segmentation as a reference, interclass correlation coefficients (ICC), relative bias, and limit of agreement (LOA) were assessed for all parameters. Normality was tested using a Shapiro–Wilk test. Depending on normality, a Mann–Whitney U-test or an unpaired t-test was utilized for statistical evaluations. A p-value <0.05 was considered statistically significant. All statistical analyses were performed in MATLAB.
3 Results
3.1 Performance of the automated segmentation
The training took approximately 10 h and the time required for a single CNN segmentation was 2.2 ± 1.0 s. Figure 2 displays exemplary results for Observer 1 (red) and CNN (blue) segmentations of a control case (Figure 2A) and an ICAD case (Figure 2B). Difference maps indicate regions of over- (blue) and underestimation (red) of the automated segmentation.
Figure 2. Exemplary results for manual (red) and CNN (blue) segmentations of a control case (A) DS = 0.89, HD = 5.09 mm, ASSD = 0.11 mm) and an ICAD case (B) DS = 0.84, HD = 3.77 mm, ASSD = 0.15 mm). Difference maps on the right indicate regions of over- (blue) and underestimation (red) of the automated segmentation. Orange arrows mark the location of a severe stenosis in the left MCA.
Table 3 displays the median values and the range of the DS, HD and ASSD for a comparison between the CNN segmentation and the original manual segmentation. No significant differences were observed when comparing the DS and HD values of the control group with the ICAD group (p ≥ 0.19). The average DS values were (mean ± STD) for CoW + sinuses: 0.86 ± 0.04 (controls) and 0.83 ± 0.04 (ICAD, p = 0.31). For CoW only: 0.85 ± 0.03 (controls) and 0.85 ± 0.04 (ICAD, p = 0.81). For sinuses only: 0.86 ± 0.06 (controls) and 0.82 ± 0.07 (ICAD, p = 0.19). Average HD values were for CoW + sinuses: 7.8 ± 1.6 mm (controls) and 7.9 ± 1.9 mm (ICAD, p = 0.89). For CoW only: 7.2 ± 1.5 mm (controls) and 8.4 ± 3.1 mm (ICAD, p = 0.28). For sinuses only: 6.6 ± 3.7 mm (controls) and 5.7 ± 1.9 mm (ICAD, p = 0.65).
For the ASSD values, significant differences between the control and ICAD group were observed for the CoW + sinus segmentation (0.14 ± 0.05 mm vs. 0.21 ± 0.08 mm, p = 0.02) and the CoW-only segmentation (0.14 ± 0.04 mm vs. 0.22 ± 0.10 mm, p < 0.05). However, for the sinuses-only segmentation, average ASSD values were 0.17 ± 0.08 mm (controls) and 0.22 ± 0.11 mm (ICAD, p = 0.23).
3.2 Flow metrics
Cross-sectional areas, peak velocity values, and flow rates were determined in the small and large arteries and the sinuses for both cohorts, respectively. Tables 4–6 display the median, range, relative bias, and limits of agreement of all of the above parameters for the automated segmentation vs. Observer 1 and the Observer 1 vs. Observer 2 comparisons. Furthermore, the p-values are shown for the Observer 1 vs. automated segmentations, the Observer 1 vs. Observer 2, and the automated segmentation vs. Observer 2 comparison. The top half of each table shows the results for the control group while the bottom half displays the results of the ICAD group. The correlation and Bland-Altman plots can be found in Supplementary Figures S1–S3 in the supplement. No significant differences for all parameters were observed when comparing the Observer 1 with the automated segmentation (p ≥ 0.40 for all vessels). The average cross-sectional areas determined with the automated segmentation were (mean ± STD) for the large arteries: 13.3 ± 4.4 mm2 (controls) vs. 12.1 ± 4.8 mm2 (ICAD, p = 0.20). For the small arteries: 5.7 ± 2.1 mm2 (controls) vs. 5.3 ± 1.1 mm2 (ICAD, p = 0.32). For the sinuses: 19.8 ± 11.6 mm2 (controls) vs. 22.1 ± 10.0 mm2 (ICAD, p = 0.48). The average peak velocity values were for the large arteries: 0.61 ± 0.15 m/s (controls) vs. 0.53 ± 0.16 m/s (ICAD, p = 0.06). For the small arteries: 0.45 ± 0.14 m/s (controls) vs. 0.41 ± 0.13 m/s (ICAD, p = 0.16). For the sinuses: 0.33 ± 0.14 m/s (controls) vs. 0.27 ± 0.12 m/s (ICAD, p = 0.15).
Table 4. Cross-sectional area values: comparison of the automated and the Observer 2 segmentation with the Observer 1 segmentation.
Table 5. Peak velocity values: comparison of the automated and the Observer 2 segmentation with the Observer 1 segmentation.
For the flow rates, a significant difference was observed between both cohorts in the large arteries: 3.9 ± 1.7 ml/s (controls) vs. 2.7 ± 1.2 ml/s (ICAD, p < 0.01). For the small arteries, the mean values were 1.4 ± 0.8 ml/s (controls) vs. 1.0 ± 0.5 ml/s (ICAD, p = 0.064). For the sinuses, the flow values were: 3.6 ± 2.6 ml/s (controls) and 3.1 ± 1.9 ml/s (ICAD, p = 0.58). When analyzing the Observer 1 segmentation, significant intergroup differences were observed for the flow rates in the large arteries (p < 0.01) and small arteries (p = 0.044) but not for the sinuses (p = 0.67). Using the flow rate values, the flow conservation error was determined for the manual and automated segmentation, respectively. For the control group, the fce was 0.16 ± 0.09 (manual) and 0.15 ± 0.10 (automated). For the ICAD group, the fce was 0.18 ± 0.12 (manual) and 0.20 ± 0.13 (automated). No significant differences were observed between the manual and CNN segmentation (p ≥ 0.80 for both groups).
3.3 Interobserver comparison
The segmentation time of an individual manual segmentation performed by Observer 2 was 1,103 ± 347 s. On the left of Table 7 the results for the DS, HD, and ASSD for the Observer 2 segmentation values are displayed using the Observer 1 segmentation as a reference. The right side of the table shows the performance analysis of the automated segmentation vs. Observer 2 comparison. No differences were observed when comparing the DS values from the automated segmentation vs. Observer 1 comparison in Table 3 with the Observer 1 vs. Observer 2 comparison (Controls: CoW + sinuses: 0.86 ± 0.03, CoW: 0.86 ± 0.04, sinuses: 0.85 ± 0.03. p-value for all vessels: ≥0.63. ICAD: CoW + sinuses: 0.81 ± 0.04. CoW: 0.84 ± 0.04. Sinuses: 0.79 ± 0.05. p-value for all vessels: ≥0.10).
The same applies for the automated segmentation vs. Observer 2 DS values (Controls: CoW + sinuses: 0.86 ± 0.02, CoW: 0.85 ± 0.03, sinuses: 0.86 ± 0.02. ICAD: CoW + sinuses: 0.82 ± 0.03, CoW: 0.83 ± 0.02, sinuses: 0.81 ± 0.04. p-value for all vessels: >0.05). When comparing the Dice scores of both cohorts, however, significantly lower DS values were observed in the ICAD group in the Observer 1 vs. Observer 2 comparison (CoW + sinuses: p < 0.01. Sinuses: p < 0.01) and in the automated segmentation vs. Observer 2 comparison (CoW + Sinuses: p < 0.01. CoW: p = 0.022. Sinuses: p = 0.045).
For the Observer 1 vs. Observer 2 HD values, no differences relative to the automated segmentation vs. Observer 1 comparison were noticeable in the control group (CoW + sinuses: 8.4 ± 2.0 mm. CoW: 7.0 ± 1.5 mm. Sinuses: 6.6 ± 3.2 mm. p-value for all vessels: ≥0.51) and the ICAD group (CoW + Sinuses: 9.8 ± 3.3 mm. CoW: 8.6 ± 3.0 mm. Sinuses: 7.7 ± 3.8 mm. p-value for all vessels: ≥0.14).
Similar HD values are also noticeable for the automated segmentation vs. Observer 2 comparison (Controls: CoW+ sinuses: 7.5 ± 1.3 mm, CoW: 7.2 ± 1.5 mm, sinuses: 5.0 ± 2.0 mm. ICAD: CoW + sinuses: 10.3 ± 3.9 mm, CoW: 8.0 ± 4.4 mm, sinuses: 8.6 ± 3.9 mm), with no differences compared to the Observer 1 vs. Observer 2 HD values (p ≥ 0.34). However, the Observer 2 vs. automated segmentation comparison of the ICAD group featured larger HD values in the sinuses compared to the automated segmentation vs. Observer 1 HD values (p = 0.047). Furthermore, the intergroup comparison revealed no differences for the Observer 1 vs. Observer 2 HD values (p-value for all vessels ≥0.14) but for the automated segmentation vs. Observer 2 HD values (CoW + sinuses: p < 0.01, sinuses: p < 0.01).
For the ASSD values, no differences between the Observer 1 vs. Observer 2 and the automated segmentation vs. Observer 1 comparison were detected (Controls: CoW + sinuses: 0.18 ± 0.05 mm. CoW: 0.18 ± 0.05 mm. Sinuses: 0.18 ± 0.05 mm. p-value for all vessels: ≥0.79. ICAD: CoW + Sinuses: 0.24 ± 0.08 mm. CoW: 0.22 ± 0.09 mm. Sinuses: 0.25 ± 0.08 mm. p-value for all vessels: ≥0.16).
Similar results were also observed for the automated segmentation vs. Observer 2 comparison (Controls: 0.15 ± 0.03 mm, CoW: 0.15 ± 0.05 mm, sinuses: 0.14 ± 0.02 mm. ICAD: CoW + sinuses: 0.21 ± 0.04 mm, CoW: 0.18 ± 0.06 mm, sinuses: 0.22 ± 0.06 mm. p-value for all vessels: p ≥ 0.10). However, in the Observer 1 vs. Observer 2 comparisons of the CoW, larger ASSD values are noticeable in the ICAD group relative to the control group (p = 0.022). Furthermore, in both comparisons presented in Table 7, larger sinusoidal ASSD values were observed (CoW + sinuses: p < 0.01, sinuses: p ≤ 0.045).
In the following, the flow metrics and cross-sectional area values determined with the Observer 2 segmentation were compared with the results from the original manual segmentation. Tables 4–6 display the bias, limits of agreement, the ICC values, and the p-values for the comparisons between the two human observers. The corresponding correlation and Bland-Altman plots can be found in the Supplementary Figures S4–S6 in the supplement. No significant differences were observed for the cross-sectional area values (p-value for all vessels: ≥0.19), the peak velocity values (p-value for all vessels: ≥0.64), and the flow rates (p-value for all vessels: ≥0.39). Furthermore, the fce analysis yielded no significant differences (p-value for both groups: ≥0.94). However, the comparison between the control and the ICAD groups revealed a significant difference for the cross-sectional areas in the large arteries (p = 0.016) but not for the small arteries (p = 0.068) and the sinuses (p = 0.75). In addition, significant differences were observed for the flow rates in the small (p = 0.011) and large arteries (p < 0.01). No significant differences were observed for the flow rates in the sinuses (p = 0.36) and for all peak velocity values (p ≥ 0.18).
The Observer 2 results were furthermore compared with the values obtained with the CNN segmentation. Similar to the comparison with the original human observer, no significant differences were observed for the cross-sectional area values (p-value for all vessels: ≥0.16), the peak velocity values (p-value for all vessels: ≥0.84), the flow rates (p-value for all vessels: ≥0.39), and the fce (p-value for both groups: ≥0.82).
3.4 Segmentation performance in stenosed vessels
For an analysis of the segmentation performance of the CNN in stenosed vessels, three intracranial stenosis cases were examined in further detail. Figures 3–5 display the segmentation results for a moderate (Figure 3) and a severe stenosis (Figure 4) in the right MCA and a severe stenosis in the right ICA (Figure 5). In Figures 3A, 4A, 5A the segmentation results of the stenosed vessel only are shown from Observer 1 (red), the automated segmentation (blue), Observer 2 (orange) and VWI (violet). Part B of the figures displays difference maps to illustrate the segmentation error between the automated segmentation and the Observer 1 (left) and between the automated segmentation and the VWI (right). The left side of part C of the figures shows profile plots of the lumen cross-sectional areas (for the Observer 1, automated, Observer 2 and VWI segmentations) as well as the peak velocity profiles. The gray shaded areas mark the location of the stenosis (see also the analysis planes in the segmentation plots in A and the orange arrow in the error maps in B). The plots on the right in part C of the figures show the time-resolved median flow rates around the stenosis (determined in the gray-shaded region above), assessed with the Observer 1 (red), the automated (blue), and the Observer 2 (orange) manual segmentation.
Figure 3. Analysis of segmentation performance in a moderate stenosis in the right MCA. (A) Isosurface renderings of Observer 1, Observer 2, automated and VWI segmentation. (B) Comparison between the automated and Observer 1 segmentation (left) and comparison between the automated and VWI segmentation (right). (C) Left: cross-sectional area and peak velocity profiles. Right: flow rates estimated in a region of interest around the stenosis (see grey shaded area in the profile plot) with the automated and Observer 1 and 2 segmentations. The stenosis is marked by analysis planes and orange arrows.
Figure 4. Analysis of segmentation performance in a severe stenosis in the right MCA. (A) Isosurface renderings of Observer 1, Observer 2, automated and VWI segmentation. (B) Comparison between the automated and Observer 1 segmentation (left) and comparison between the automated and VWI segmentation (right). (C) Left: cross-sectional area and peak velocity profiles. Right: flow rates estimated in a region of interest around the stenosis (see grey shaded area in the profile plot) with the automated and Observer 1 and 2 segmentations. The stenosis is marked by analysis planes and orange arrows.
Figure 5. Analysis of segmentation performance in a severe stenosis in the right ICA. (A) Isosurface renderings of Observer 1, Observer 2, automated and VWI segmentation. (B) Comparison between the automated and Observer 1 segmentation (left) and comparison between the automated and VWI segmentation (right). (C) Left: cross-sectional area and peak velocity profiles. Right: flow rates estimated in a region of interest around the stenosis (see grey shaded area in the profile plot) with the automated and Observer 1 and 2 segmentations. The stenosis is marked by analysis planes and orange arrows.
Table 8 illustrates the segmentation metrics of the stenosed vessels as well as statistical comparisons of the lumen area and flow rate values, determined around the stenosis (see gray shaded areas in the profile plots in Figures 3C–5C). Table 8A displays the DS, HD and ASSD values for the Observer 1 (O1) segmentation, the Observer 2 (O2) segmentation and the automated (CNN) segmentation using the VWI segmentation as reference. For calculation of segmentation metrics for this table, only the part of the vessel that was stenosed was considered. All three 4D flow MRI segmentations yield similar DS values (Observer 1: range 0.61–0.77. Observer 2: range 0.52–0.73. Automated segmentation: 0.59–0.75), similar HD values (Observer 1: range 2.1 mm–6.7 mm. Observer 2: range 2.3 mm–7.1 mm. Automated segmentation: range 2.6 mm–7.1 mm) and similar ASSD values (Observer 1: 0.20 mm–0.30 mm. Observer 2: 0.22 mm–0.45 mm. Automated segmentation: 0.23 mm–0.32 mm). Table 8B shows the lumen cross-sectional area values (median values and interquartile ranges) of the stenosed part of the vessels, estimated with the VWI segmentation, the Observer 1 (O1) segmentation, the Observer 2 (O2) segmentation and the automated (CNN) segmentation. In addition, bias and LOA values relative to the VWI segmentation are displayed. All three 4D flow segmentations feature a consistent overestimation relative to the segmentation obtained from the black blood images but good agreement (Observer 1: Overall bias ± LOA: 36.6 ± 16.1%, ICC: 0.93, p < 0.01. Observer 2: Overall bias ± LOA: 31.2 ± 42.4%, ICC: 0.84, p = 0.025. Automated segmentation: Overall bias ± LOA: 28.1 ± 13.9%, ICC: 0.93, p < 0.01). However, no differences were observed when comparing the automated segmentation with the manual segmentations (automated segmentation vs. Observer 1: Overall bias ± LOA: −8.8 ± 10.8%, ICC: 0.95, p = 0.05. Automated segmentation vs. Observer 2: Overall bias ± LOA: −3.2 ± 35.9%, ICC: 0.89, p = 0.15). Table 8C displays the flow rate values in the stenosed regions (median and interquartile ranges), estimated with both manual 4D flow segmentations and the automated segmentation. In addition, the table displays the bias and LOA values relative to the Observer 1 segmentation. No differences were observed when comparing the results from the automated segmentation with the manual segmentations of Observer 1 and 2 (automated segmentation vs. Observer 1: Overall bias ± LOA: −3.9 ± 4.3%, ICC > 0.99, p = 0.34. Automated segmentation vs. Observer 2: Overall bias ± LOA: −3.1 ± 13.9%, ICC: 0.98, p = 0.48).
4 Discussion
In this work, a CNN was trained using manually segmented intracranial 4D flow data and was successfully applied for the fully automated segmentation of stenosed intracranial vasculatures. In both a healthy control group and an ICAD patient cohort, similar segmentation performance could be achieved (comparison with Observer 1: median DS: ≥0.85 for controls, ≥0.84 for patients. Median HD: ≥6.3 mm for controls, ≥5.9 mm for patients. Median ASSD: ≥0.12 mm for controls, 0.21 mm for patients). In addition, no significant differences were observed when comparing the flow parameters and cross-sectional area values determined with the CNN segmentation with the original manual analysis and with the analysis performed by Observer 2. Interestingly, however, lower flow rates were observed in the CoW arteries in patients compared to controls regardless of the segmentation used for the analysis. These flow rate differences are likely due to the age difference of the two testing cohorts (see Table 1), as age related flow rate differences have already been reported by Wu et al. (27).
The automated segmentation required substantially less time than the manual segmentation performed by Observer 2 and was also significantly faster than reported for other automated segmentation techniques for intracranial 4D flow data (14). Furthermore, in contrast to manual segmentations, the CNN segmentation was not susceptible to often observed inter-observer variabilities (28).
Until recently, most deep learning-based segmentation networks of intracranial vasculature were based on TOF-MRA (29), CTA (30) or CTA in combination with DSA (31). Furthermore, the recent TopCoW challenge yielded very impressive results for CoW segmentations using CTA and MRA images (32). The large dataset available for this challenge may be used for transfer learning to further improve the segmentation performance of 4D flow images. A possible challenge, however, may be the large difference in voxel size between MRA [0.30 mm × 0.30 mm × 0.71 mm according to (32)] and 4D flow MRI (1 mm isotropic) and differences in image contrast, signal to noise ratio and artifacts due to the different sequence design. Furthermore, in contrast to the work presented in this paper, the TopCoW challenge only addressed the segmentation of the CoW arteries but not of veins such as the sinuses.
For 4D flow MRI, until recently, non-deep learning techniques such as threshold-based segmentation (33) or segmentation based on a centerline processing scheme (14) have been more common. Rothenberger et al. recently presented a post-processing technique using a standard difference of means, yielding average Dice score values of 0.76 for the CoW (15). However, in our study, we showed that using a deep learning approach, a larger mean Dice score of 0.85 could be achieved.
When comparing the Observer 2 segmentation with Observer 1 or the automated segmentation, no significant differences in segmentation performance were noticeable. Notable distinctions, however, were observed when comparing the control and stenosis cases. In both the Observer 1 vs. Observer 2 and the automated vs. Observer 2 comparisons, the ICAD group featured significantly lower Dice scores in the sinuses and overall larger ASSD values. Significantly larger ASSD values are also noticeable in the automated segmentations of the ICAD cases. One reason for the lower segmentation performance in the sinuses of the ICAD group may be the varying field of view size. In the control group, the number of slices varied between 40 and 44 while in the ICAD group, the number of slices was between 26 and 60. A too small FOV, however, may lead to incomplete coverage of the sinuses, which may exacerbate an accurate vessel segmentation. Variations in the number of slices may also partially explain the significantly larger ASSD values in the arteries, since a too small FOV may lead to insufficient coverage of the basilar artery and other vessels typically at the edges of the FOV such as the vertebral arteries. Another reason for the significantly larger ASSD values in the arteries may be the larger variability in vascular geometry noticeable in patients with ICAD due to the pathological changes caused by atherosclerosis and due to the significant differences in age between the two testing cohorts (see Table 1). Furthermore, especially in severely stenosed arteries, noticeable signal dropouts are often observed which aggravate accurate vessel segmentation.
The analysis of the cross-sectional area and flow metrics yielded no significant differences between the automated, the original manual, and the Observer 2 segmentation. In the large arteries, small bias and small limits of agreement are noticeable for both the control and the ICAD group. Slightly larger limits of agreement were observed in the small arteries, which may be attributed to partial volume effects due to the notably smaller vessel size. The reason for the larger variations observed in the sinuses may be explained again by the sometimes-incomplete coverage of the veins due to a too small FOV size. Also, exact segmentation of the sinuses is more challenging due to the much lower signal intensities as well as lower velocities in these vessels.
More detailed analyses of the segmentation performance in stenosed vessels revealed good performance of segmentation in stenosed areas that resulted in similar flow rates and cross-sectional area values between the automated segmentation and the two manual observers. Co-registration of 4D flow MRI with black blood vessel wall imaging confirmed correct segmentations of stenosed regions within the limits of spatial resolution of the 4D Flow MRI acquisition. As expected, an overestimation of the lumen areas relative to results obtained from black blood vessel wall imaging was noticeable in all 4D flow-based segmentations. Furthermore, due to the different sequence designs, the switching of the imaging gradients is different between the two imaging modalities. Thus, differences in image artifacts such as distortion, blurring and motion corruption are to be expected, exacerbating accurate co-registration. In this work, rigid co-registration was used to match the black blood segmentation with the 4D flow-based segmentations. For more accurate co-registration results, non-rigid co-registrations can be considered, however, this would increase the time investment of the post-processing and would be out of the scope of this study. In this work, we aimed to develop a CNN to automate intracranial vessel segmentation from 4D flow MRI data to ease the analysis of volumetric hemodynamic parameters. Our aim was not to use 4D flow MRI for the diagnosis of stenosis grade using luminal narrowing.
One limitation of this study is the small number of cases used for training and testing of the CNN architecture. In this work, 134 cases were used for training and 20 cases for testing. In contrast, 499 cases were used for training the CNN for automated aorta segmentation (16). The small number of training cases may be problematic since large variations in the CoW geometry have been reported (34). The limited number of training cases may also be a further explanation for the slightly worse segmentation performance of the sinuses in ICAD patients. Furthermore, the distribution of healthy and diseased training cases was skewed with 76 control cases but only 58 ICAD cases. However, we think that using as many training cases as possible was more important than an even distribution of healthy vs. ICAD training cases. In addition, due to the limited number of ICAD cases, we only had very few test cases. The random selection of the test cases caused an age-difference between the two testing cohorts, leading to significantly younger control cases compared to the ICAD patients (see Table 1). However, the focus of our study was to create a CNN for the segmentation of intracranial vessels in healthy controls as well as ICAD patients. This means that we selected as many as possible intracranial 4D flow MRI datasets while neglecting age matching. In a future more clinically focused study, quantitative results of the ICAD patients compared to an age- and gender matched healthy control cohort will be assessed. In addition, to further improve performance of the CNN, more stenosis cases as well as of other intracranial vascular diseases will be incorporated to further improve the segmentation performance and for a generalization of the automated segmentation.
5 Conclusion
In this work, a deep learning-based approach was presented for the fully automated vessel segmentation of intracranial 4D flow MRI data of healthy subjects and stenosis patients. The introduced CNN segmentation took only 2.2 s on average to complete. The automated segmentations of the intracranial arteries and veins are in very good agreement with the manual segmentations of two independent observers and the analysis of lumen cross-sectional areas and flow metrics yielded no significant differences between manual and automated segmentations. Furthermore, the accuracy of the automated segmentation of stenosed intracranial arteries could be verified by co-registered vessel wall imaging. The automation of intracranial vessel segmentation significantly reduces the analysis time and may improve the robustness of determining hemodynamic parameters with intracranial 4D flow MRI. This work could therefore be an integral factor in increasing its clinical application.
Data availability statement
The datasets presented in this article are not readily available because datasets from an IRB approved study involving patients were used. The source code for the convolutional neural network can be downloaded from https://github.com/drsol1986/AI-Segmentation-.git.
Ethics statement
The studies involving humans were approved by Northwestern University, Institutional Review Board Office. The studies were conducted in accordance with the local legislation and institutional requirements. The healthy controls gave their written informed consent. The stenosis cases were retrospectively selected.
Author contributions
PW: Conceptualization, Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing, Formal Analysis, Methodology. HB: Formal Analysis, Methodology, Software, Supervision, Writing – review & editing. JM: Data curation, Writing – review & editing. MA: Data curation, Writing – review & editing, Investigation, Methodology, Supervision. TR: Writing – review & editing, Methodology, Software. JW: Writing – review & editing, Data curation. AR: Writing – review & editing, Data curation. KJ: Formal Analysis, Writing – review & editing, Conceptualization, Methodology, Supervision. AP: Writing – review & editing, Data curation, Formal Analysis, Investigation, Validation. FC: Writing – review & editing, Project administration, Resources. RA: Formal Analysis, Investigation, Writing – review & editing, Data curation. SA: Project administration, Resources, Supervision, Writing – review & editing, Formal Analysis, Investigation, Validation. MM: Project administration, Resources, Supervision, Writing – review & editing, Funding acquisition. SS: Conceptualization, Formal Analysis, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Writing – review & editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article.
Research reported in this publication was supported by the National Institutes of Health (National Institute of Biomedical Imaging and Bioengineering: T32EB025766, National Heart, Lung and Blood Institute: 1R01HL149787, National Institute of Neurological Disorders and Stroke: 5R21NS122511; National Institute on Aging: P30AG059988 and K01AG080070). JW received funding from the Gerhard Domagk Nachwuchsförderprogramm.
Acknowledgments
The authors thank Dr. Patrick Vogel for his technical support with the black blood segmentations.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fradi.2024.1385424/full#supplementary-material
Supplementary Figure S1
Correlation and Bland-Altman analysis of the cross-sectional areas in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the automated segmentation.
Supplementary Figure S2
Correlation and Bland-Altman analysis of the peak velocity values in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the automated segmentation.
Supplementary Figure S3
Correlation and Bland-Altman analysis of the flow rates in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the automated segmentation.
Supplementary Figure S4
Correlation and Bland-Altman analysis of the cross-sectional areas in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the Observer 2 segmentation.
Supplementary Figure S5
Correlation and Bland-Altman analysis of the peak velocity values in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the Observer 2 segmentation.
Supplementary Figure S6
Correlation and Bland-Altman analysis of the flow rates in the large arteries (A), small arteries (B) and sinuses (C): Comparison between the manual ground truth and the Observer 2 segmentation.
References
1. Schnell S, Ansari SA, Wu C, Garcia J, Murphy IG, Rahman OA, et al. Accelerated dual-venc 4D flow MRI for neurovascular applications. J Magn Reson Imaging. (2017) 46(1):102–14. doi: 10.1002/jmri.25595
2. Marlevi D, Schollenberger J, Aristova M, Ferdian E, Ma Y, Young AA, et al. Noninvasive quantification of cerebrovascular pressure changes using 4D Flow MRI. Magn Reson Med. (2021) 86(6):3096–110. doi: 10.1002/mrm.28928
3. Schnell S, Ansari SA, Vakil P, Wasielewski M, Carr ML, Hurley MC, et al. Three-dimensional hemodynamics in intracranial aneurysms: Influence of size and morphology. J Magn Reson Imaging. (2014) 39(1):120–31. doi: 10.1002/jmri.24110
4. Aristova M, Vali A, Ansari SA, Shaibani A, Alden TD, Hurley MC, et al. Standardized evaluation of cerebral arteriovenous malformations using flow distribution network graphs and dual-venc 4D flow MRI. J Magn Reson Imaging. (2019) 50(6):1718–30. doi: 10.1002/jmri.26784
5. Vali A, Aristova M, Vakil P, Abdalla R, Prabhakaran S, Markl M, et al. Semi-automated analysis of 4D flow MRI to assess the hemodynamic impact of intracranial atherosclerotic disease. Magn Reson Med. (2019) 82(2):749–62. doi: 10.1002/mrm.27747
6. Moccia S, De Momi E, El Hadji S, Mattos LS. Blood vessel segmentation algorithms—review of methods, datasets and evaluation metrics. In: Computer Methods and Programs in Biomedicine, vol. 158. Elsevier Ireland Ltd. (2018). p. 71–91. doi: 10.1016/j.cmpb.2018.02.001
7. Goni MR, Ruhaiyem NIR, Mustapha M, Achuthan A, Che Mohd Nassir CMN. Brain vessel segmentation using deep learning—a review. IEEE Access. (2022) 10:111322–36. doi: 10.1109/ACCESS.2022.3214987
8. Phellan R, Peixinho A, Falçao A, Forkert ND. Vascular segmentation in TOF MRA images of the brain using a deep convolutional neural network. In: Cardoso MJ, Arbel T, Lee S-L, Cheplygina V, Balocco S, Mateus D, et al., editors. Intracranial Imaging and Computer Assisted Stenting, and Large-Scale Annotation of Biomedical Data and Expert Label Synthesis, Lecture Notes in Computer Science, vol. 10552. Cham: Springer International Publishing (2017). p. 39–46. doi: 10.1007/978-3-319-67534-3
9. Ronneberger O, Fischer P, Brox T. U-Net: convolutional networks for biomedical image segmentation. In: Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer Verlag (2015). p. 234–41. doi: 10.1007/978-3-319-24574-4_28
10. Çiçek Ö, Abdulkadir A, Lienkamp SS, Brox T, Ronneberger O. 3D U-Net: learning dense volumetric segmentation from sparse annotation. In: MICCAI 2016: Medical Image Computing and Computer-Assisted Intervention—MICCAI 2016, Lecture Notes in Computer Science, vol. 9901. (2016). p. 424–32. doi: 10.1007/978-3-319-46723-8
11. Siddique N, Paheding S, Elkin CP, Devabhaktuni V. U-Net and its variants for medical image segmentation: a review of theory and applications. IEEE Access. (2021) 9:82031–57. doi: 10.1109/ACCESS.2021.3086020
12. Livne M, Rieger J, Aydin OU, Taha AA, Akay EM, Kossen T, et al. A U-Net deep learning framework for high performance vessel segmentation in patients with cerebrovascular disease. Front Neurosci. (2019) 13(FEB):1–13. doi: 10.3389/fnins.2019.00097
13. Patel TR, Paliwal N, Jaiswal P, Waqas M, Mokin M, Siddiqui AH, et al. Multi-resolution CNN for brain vessel segmentation from cerebrovascular images of intracranial aneurysm: a comparison of U-Net and DeepMedic. Proceedings of SPIE Volume 11314, SPIE-Intl Soc Optical Eng; (2020). p. 101. doi: 10.1117/12.2549761
14. Schrauben E, Ambarki K, Spaak E, Malm J, Wieben O, Eklund A. Fast 4D flow MRI intracranial segmentation and quantification in tortuous arteries. J Magn Reson Imaging. (2015) 42(5):1458–64. doi: 10.1002/jmri.24900
15. Rothenberger SM, Patel NM, Zhang J, Schnell S, Craig BA, Ansari SA, et al. Automatic 4D flow MRI segmentation using the standardized difference of means velocity. IEEE Trans Med Imaging. (2023) 42(8):2360–73. doi: 10.1109/TMI.2023.3251734
16. Berhane H, Scott M, Elbaz M, Jarvis K, McCarthy P, Carr J, et al. Fully automated 3D aortic segmentation of 4D flow MRI for hemodynamic analysis using deep learning. Magn Reson Med. (2020) 84(4):2204–18. doi: 10.1002/mrm.28257
17. Komada T, Naganawa S, Ogawa H, Matsushima M, Kubota S, Kawai H, et al. Contrast-Enhanced MR Imaging of Metastatic Brain Tumor at 3 Tesla: Utility of T 1-weighted SPACE Compared with 2D Spin Echo and 3D Gradient Echo Sequence (2008).
18. Bock J, Kreher BW, Hennig J, Markl M. Optimized pre-processing of time-resolved 2D and 3D phase contrast MRI data. Proceedings of the 15th ISMRM Annual Meeting; May 2007; Berlin. p. 3138.
19. Huang G, Liu Z, Van Der Maaten L, Weinberger KQ. Densely connected convolutional networks. IEEE Conference on Computer Vision and Pattern Recognition; Las Vegas, Nevada, USA, (2016). p. 4700–08. Available online at: https://github.com/liuzhuang13/DenseNet26865468
20. Taha AA, Hanbury A. Metrics for evaluating 3D medical image segmentation: analysis, selection, and tool. BMC Med Imaging. (2015) 15(1):1–28. doi: 10.1186/s12880-015-0068-x
21. Ashburner J, Friston KJ. Unified segmentation. Neuroimage. (2005) 26(3):839–51. doi: 10.1016/j.neuroimage.2005.02.018
22. Guggenberger K, Krafft AJ, Ludwig U, Raithel E, Forman C, Meckel S, et al. Intracranial vessel wall imaging framework—data acquisition, processing, and visualization. Magn Reson Imaging. (2021) 83:114–24. doi: 10.1016/j.mri.2021.08.004
23. Vogel P, Kampf T, Guggenberger K, Raithel E, Forman C, Meckel S, et al. Robust centerline prediction for accurate vessel wall visualization of intracranial vessels in multi-contrast 3D MRI data. arXiv:2209.01422.
25. Lorensen WE, Cline HE. Marching cubes: a high resolution 3D surface construction algorithm (1987).
26. Fischer B, Modersitzki J. FLIRT: a flexible image registration toolbox. International Workshop on Biomedical Image Registration; (2003). p. 261–70.
27. Wu C, Honarmand AR, Schnell S, Kuhn R, Schoeneman SE, Ansari SA, et al. Age-related changes of normal cerebral and cardiac blood flow in children and adults aged 7 months to 61 years. J Am Heart Assoc. (2016) 5(1):1–13. doi: 10.1161/JAHA.115.002657
28. Juffermans JF, Westenberg JJM, van den Boogaard PJ, Roest AAW, van Assen HC, van der Palen RLF, et al. Reproducibility of aorta segmentation on 4D flow MRI in healthy volunteers. J Magn Reson Imaging. (2021) 53(4):1268–79. doi: 10.1002/jmri.27431
29. Hilbert A, Madai VI, Akay EM, Aydin OU, Behland J, Sobesky J, et al. BRAVE-NET: fully automated arterial brain vessel segmentation in patients with cerebrovascular disease. Front Artif Intell. (2020) 3:1–14. doi: 10.3389/frai.2020.552258
30. Chen Y, Fan S, Chen Y, Che C, Cao X, He X, et al. Vessel segmentation from volumetric images: a multi-scale double-pathway network with class-balanced loss at the voxel level. Med Phys. (2021) 48(7):3804–14. doi: 10.1002/mp.14934
31. Patel TR, Patel A, Veeturi SS, Shah M, Waqas M, Monteiro A, et al. Evaluating a 3D deep learning pipeline for cerebral vessel and intracranial aneurysm segmentation from computed tomography angiography–digital subtraction angiography image pairs. Neurosurg Focus. (2023) 54(6):1–11. doi: 10.3171/2023.3.FOCUS2374
32. Yang K, Musio F, Ma Y, Juchler N, Paetzold JC, Al-Maskari R, et al. TopCoW: benchmarking topology-aware anatomical segmentation of the circle of willis (CoW) for CTA and MRA. arXiv:2312:17670 (2023). Available online at: http://arxiv.org/abs/2312.17670 (Accessed May 17, 2024).
33. Roberts GS, Hoffman CA, Rivera-Rivera LA, Berman SE, Eisenmenger LB, Wieben O. Automated hemodynamic assessment for cranial 4D flow MRI. Magn Reson Imaging. (2023) 97:46–55. doi: 10.1016/j.mri.2022.12.016
34. Hakim A, Gralla J, Rozeik C, Mordasini P, Leidolt L, Piechowiak E, et al. Anomalies and normal variants of the cerebral arterial supply: a comprehensive pictorial review with a proposed workflow for classification and significance. J Neuroimaging. (2018) 28(1):14–35. Blackwell Publishing Inc. doi: 10.1111/jon.12475
Keywords: atherosclerosis, intracranial, segmentation, deep learning, stenoses, 4D flow, stenosis, convolutional neural network
Citation: Winter P, Berhane H, Moore JE, Aristova M, Reichl T, Wollenberg J, Richter A, Jarvis KB, Patel A, Caprio FZ, Abdalla RN, Ansari SA, Markl M and Schnell S (2024) Automated intracranial vessel segmentation of 4D flow MRI data in patients with atherosclerotic stenosis using a convolutional neural network. Front. Radiol. 4:1385424. doi: 10.3389/fradi.2024.1385424
Received: 12 February 2024; Accepted: 13 May 2024;
Published: 4 June 2024.
Edited by:
Alfonso Mastropietro, Istituto di Tecnologie Biomediche—Consiglio Nazionale delle Ricerche, ItalyReviewed by:
Jing Yuan, Zhejiang Normal University, ChinaChella Santhosh, K L University, India
Tatsat Rajendra Patel, University at Buffalo, United States
© 2024 Winter, Berhane, Moore, Aristova, Reichl, Wollenberg, Richter, Jarvis, Patel, Caprio, Abdalla, Ansari, Markl and Schnell. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Patrick Winter, cGF0cmljay53aW50ZXJAdW5pLWdyZWlmc3dhbGQuZGU=
‡ORCID:
Patrick Winter
orcid.org/0000-0003-1507-8524
Maria Aristova
orcid.org/0000-0002-9330-9415
Adam Richter
orcid.org/0000-0003-1904-4161
Kelly B. Jarvis
orcid.org/0000-0002-0419-5541
Ramez N. Abdalla
orcid.org/0000-0003-2611-2097
Sameer A. Ansari
orcid.org/0000-0001-8176-4616
Susanne Schnell
orcid.org/0000-0001-6009-4756