A framework for automated and objective modification of tubular structures: Application to the internal carotid artery

A parametric study of the effect of 3D plaque shape on local hemodynamics and implications for plaque instability

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier-Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

On the numerical treatment of viscous and convective effects in relative pressure reconstruction methods

The mechanism of many cardiovascular diseases can be understood by studying the pressure distribution in blood vessels. Direct pressure measurements, however, require invasive probing and provide only single-point data. Alternatively, relative pressure fields can be reconstructed from imaging-based velocity measurements by considering viscous and inertial forces. Both contributions can be potential troublemakers in pressure reconstruction: the former due to its higher-order derivatives, and the latter because of the quadratic nonlinearity in the convective acceleration. Viscous and convective terms can be treated in various forms, which, although equivalent for ideal measurements, can perform differently in practice. In fact, multiple versions are often used in literature, with no apparent consensus on the more suitable variants. In this context, the present work investigates the performance of different versions of relative pressure estimators. For viscous effects, in particular, two new modified estimators are presented to circumvent second-order differentiation without requiring surface integrals. In-silico and in-vitro data in the typical regime of cerebrovascular flows are considered, allowing a systematic noise sensitivity study. Convective terms are shown to be the main source of error, even for flows with pronounced viscous component. Moreover, the conservation (often integrated) form of convection exhibits higher noise sensitivity than the standard convective description, in all three families of estimators considered here. For the classical pressure Poisson estimator, the present modified version of the viscous term tends to yield better accuracy than the (recently introduced) integrated form, although this effect is in most cases negligible when compared to convection-related errors.

Automated landmarking of bends in vascular structures: a comparative study with application to the internal carotid artery

Automated tools for landmarking the internal carotid artery (ICA) bends have the potential for efficient and objective medical image-based morphometric analysis. The two existing algorithms rely on numerical approximations of curvature and torsion of the centerline. However, input parameters, original source code, comparability, and robustness of the algorithms remain unknown. To address the former two, we have re-implemented the algorithms, followed by sensitivity analyses. Of the input parameters, the centerline smoothing had the least impact resulting in 6-7 bends, which is anatomically realistic. In contrast, centerline resolution showed to completely over- and underestimated the number of bends varying from 3 to 33. Applying the algorithms to the same cohort revealed a variability that makes comparison of results between previous studies questionable. Assessment of robustness revealed how one algorithm is vulnerable to model smoothness and noise, but conceptually independent of application. In contrast, the other algorithm is robust and consistent, but with limited general applicability. In conclusion, both algorithms are equally valid albeit they produce vastly different results. We have provided a well-documented open-source implementation of the algorithms. Finally, we have successfully performed this study on the ICA, but application to other vascular regions should be performed with caution.

A framework for automated and objective modification of tubular structures: Application to the internal carotid artery

Patient-specific medical image-based computational fluid dynamics has been widely used to reveal fundamental insight into mechanisms of cardiovascular disease, for instance, correlating morphology to adverse vascular remodeling. However, segmentation of medical images is laborious, error-prone, and a bottleneck in the development of large databases that are needed to capture the natural variability in morphology. Instead, idealized models, where morphological features are parameterized, have been used to investigate the correlation with flow features, but at the cost of limited understanding of the complexity of cardiovascular flows. To combine the advantages of both approaches, we developed a tool that preserves the patient-specificness inherent in medical images while allowing for parametric alteration of the morphology. In our open-source framework morphMan we convert the segmented surface to a Voronoi diagram, modify the diagram to change the morphological features of interest, and then convert back to a new surface. In this paper, we present algorithms for modifying bifurcation angles, location of branches, cross-sectional area, vessel curvature, shape of bends, and surface roughness. We show qualitative and quantitative validation of the algorithms, performing with an accuracy exceeding 97% in general, and proof-of-concept on combining the tool with computational fluid dynamics. By combining morphMan with appropriate clinical measurements, one could explore the morphological parameter space and resulting hemodynamic response using only a handful of segmented surfaces, effectively minimizing the main bottleneck in image-based computational fluid dynamics.