Marco Boegel1, Sonja Gehrisch2, Thomas Redel2, Christopher Rohkohl2, Philip Hoelter3, Arnd Doerfler3, Andreas Maier4, Markus Kowarschik2. 1. Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Martensstr. 3, 91058, Erlangen, Germany. marco.boegel@fau.de. 2. Siemens Healthcare GmbH, Siemensstr. 1, 91301, Forchheim, Germany. 3. Department of Neuroradiology, Universitätsklinikum Erlangen, 91058, Erlangen, Germany. 4. Pattern Recognition Lab, Department of Computer Science, Friedrich-Alexander Universität Erlangen-Nürnberg, Martensstr. 3, 91058, Erlangen, Germany.
Abstract
PURPOSE: Hemodynamic simulations are of increasing interest for the assessment of aneurysmal rupture risk and treatment planning. Achievement of accurate simulation results requires the usage of several patient-individual boundary conditions, such as a geometric model of the vasculature but also individualized inflow conditions. METHODS: We propose the automatic estimation of various parameters for boundary conditions for computational fluid dynamics (CFD) based on a single 3D rotational angiography scan, also showing contrast agent inflow. First the data are reconstructed, and a patient-specific vessel model can be generated in the usual way. For this work, we optimize the inflow waveform based on two parameters, the mean velocity and pulsatility. We use statistical analysis of the measurable velocity distribution in the vessel segment to estimate the mean velocity. An iterative optimization scheme based on CFD and virtual angiography is utilized to estimate the inflow pulsatility. Furthermore, we present methods to automatically determine the heart rate and synchronize the inflow waveform to the patient's heart beat, based on time-intensity curves extracted from the rotational angiogram. This will result in a patient-individualized inflow velocity curve. RESULTS: The proposed methods were evaluated on two clinical datasets. Based on the vascular geometries, synthetic rotational angiography data was generated to allow a quantitative validation of our approach against ground truth data. We observed an average error of approximately [Formula: see text] for the mean velocity, [Formula: see text] for the pulsatility. The heart rate was estimated very precisely with an average error of about [Formula: see text], which corresponds to about 6 ms error for the duration of one cardiac cycle. Furthermore, a qualitative comparison of measured time-intensity curves from the real data and patient-specific simulated ones shows an excellent match. CONCLUSION: The presented methods have the potential to accurately estimate patient-specific boundary conditions from a single dedicated rotational scan.
PURPOSE: Hemodynamic simulations are of increasing interest for the assessment of aneurysmal rupture risk and treatment planning. Achievement of accurate simulation results requires the usage of several patient-individual boundary conditions, such as a geometric model of the vasculature but also individualized inflow conditions. METHODS: We propose the automatic estimation of various parameters for boundary conditions for computational fluid dynamics (CFD) based on a single 3D rotational angiography scan, also showing contrast agent inflow. First the data are reconstructed, and a patient-specific vessel model can be generated in the usual way. For this work, we optimize the inflow waveform based on two parameters, the mean velocity and pulsatility. We use statistical analysis of the measurable velocity distribution in the vessel segment to estimate the mean velocity. An iterative optimization scheme based on CFD and virtual angiography is utilized to estimate the inflow pulsatility. Furthermore, we present methods to automatically determine the heart rate and synchronize the inflow waveform to the patient's heart beat, based on time-intensity curves extracted from the rotational angiogram. This will result in a patient-individualized inflow velocity curve. RESULTS: The proposed methods were evaluated on two clinical datasets. Based on the vascular geometries, synthetic rotational angiography data was generated to allow a quantitative validation of our approach against ground truth data. We observed an average error of approximately [Formula: see text] for the mean velocity, [Formula: see text] for the pulsatility. The heart rate was estimated very precisely with an average error of about [Formula: see text], which corresponds to about 6 ms error for the duration of one cardiac cycle. Furthermore, a qualitative comparison of measured time-intensity curves from the real data and patient-specific simulated ones shows an excellent match. CONCLUSION: The presented methods have the potential to accurately estimate patient-specific boundary conditions from a single dedicated rotational scan.
Authors: Marco Boegel; Philip Hoelter; Thomas Redel; Andreas Maier; Joachim Hornegger; Arnd Doerfler Journal: Conf Proc IEEE Eng Med Biol Soc Date: 2015
Authors: Odile Bonnefous; Vitor Mendes Pereira; Rafik Ouared; Olivier Brina; Hans Aerts; Roel Hermans; Fred van Nijnatten; Jean Stawiaski; Daniel Ruijters Journal: Med Phys Date: 2012-10 Impact factor: 4.071
Authors: B Davis; K Royalty; M Kowarschik; C Rohkohl; E Oberstar; B Aagaard-Kienitz; D Niemann; O Ozkan; C Strother; C Mistretta Journal: AJNR Am J Neuroradiol Date: 2013-04-25 Impact factor: 3.825
Authors: Matthew D Ford; Gordan R Stuhne; Hristo N Nikolov; Damiaan F Habets; Stephen P Lownie; David W Holdsworth; David A Steinman Journal: IEEE Trans Med Imaging Date: 2005-12 Impact factor: 10.048
Authors: Jürgen Endres; Markus Kowarschik; Thomas Redel; Puneet Sharma; Viorel Mihalef; Joachim Hornegger; Arnd Dörfler Journal: Comput Math Methods Med Date: 2012-11-04 Impact factor: 2.238
Authors: S Lang; P Hoelter; A I Birkhold; M Schmidt; J Endres; C Strother; A Doerfler; H Luecking Journal: AJNR Am J Neuroradiol Date: 2019-09 Impact factor: 3.825