Hierarchy Aware CT-Driven Constrained Constructive Optimization: A Framework for Coronary Network Generation

Subject-specific three-dimensional models of the coronary arterial tree are widely used for simulations of coronary flow and myocardial perfusion. Current image-based reconstructions, however, cannot capture the full microvascular hierarchy, whereas existing synthetic tree generation methods either operate in idealized left ventricular geometries or lack consistency with experimental vessel hierarchy when constrained by anatomy.

To address these issues, we present a hierarchy-aware CT-driven constrained constructive optimization (HCT-CCO) framework that generates coronary arterial trees as extensions from segmented subject-specific base epicardial coronary networks. The method combines a distance field that restricts candidate segments to the segmented myocardial wall, applies Kamiya-type optimization seeded with HK-type daughter radii and a flow-weighted bifurcation point, and a hierarchy-aware, stage-wise branching schedule that progressively shifts growth toward distal segments as the domain fills. The resulting networks are assigned diameter-defined Strahler orders and analyzed in terms of diameter-order scaling, order-order connectivity, and transmural segment distributions.

We show that the HCT-CCO algorithm is more stable across stochastic realizations and produces connectivity matrices and diameter-order curves that show improved agreement with experimental benchmarks compared to the conventional CCO algorithm. Application to one synthetic case and five subject-specific CT datasets yielded coronary trees spanning diameter-defined Strahler orders 1-11 with monotonic diameter-order trends and mean diameters consistent with measured porcine coronary morphometry and transmural segment distributions across cases.

The HCT-CCO algorithm thus provides a hierarchy-aware, image-constrained approach for generating subject-specific coronary arterial networks that jointly capture experimentally supported diameter-order relations, hierarchical connectivity, and transmural vessel distributions, and offers anatomically grounded substrates for coronary hemodynamic and perfusion modeling.