Mathematical Reconstruction of Patient-Specific Vascular Networks Based on Clinical Images and Global Optimization


We present a hybrid reconstruction framework that integrates medical image information with angiogenesis-based optimization to generate a complete, 3-D, patient-specific vascular network of the human brain. In particular, we use segmentation techniques to obtain a coarse structure of the brain vascular network and then search for a refined configuration with minimum network material cost and power cost using global constructive optimization.

Submitted to Computers in Biology and Medicine

The presentation slide can be found here.


Cancer is a major cause of death worldwide and it becomes particularly threatening once it begins to metastasize. During metastasis, the cancerous cells of the primary tumor start to spread in the body and form secondary tumors. The blood vessels serve as pathways for this transportation and hence are crucial for understanding and monitoring cancer growth. Existing medical imaging modalities, such as computed tomography (CT) and magnetic resonance angiography (MRA), are able to provide 3-D contrast images of the vascular tissues, but the data acquired are often incomplete and lack essential details. A much-needed tool for studying blood vessels is one that could reconstruct patient-specific vascular network models based on incomplete data obtained from clinical images.


To this end, we developed a computational framework that takes a medical image as input and reconstructs a complete, patient-specific vascular network model using a mathematical optimization procedure. Our framework extracts major vessel segments from the provided image and uses the organ geometry to select vessel termination points. Then, it generates the remainder vessels based on physiological optimality principles.


Using the framework, we generated a set of vascular network models with over 3000 terminal segments from a patient’s brain MRA scan. The resulted networks were validated in two ways: (1) We analyzed the distributions of the Strahler orders, vessel radii, and branch lengths of the vascular models. These morphometric properties match with actual human data. (2) We performed fluid dynamics simulation inside the reconstructed vessels and showed that the pressure and wall shear stress distributions agree with in vivo measurements in the existing literature.

Conclusion and significance

The qualitative and quantitative agreements in morphometric and hemodynamic properties of the generated vasculatures demonstrate that the proposed framework is effective in reconstructing brain vascular network structures. It bridges the gap between image-based vessel models, which are limited by the resolution of the clinical images, and hypothetical models, which cannot be used in simulating blood flow and planning treatment for individual patients. This is an important step towards developing novel diagnostic as well as prognostic tools in cancer research.

Junhong Shen
Junhong Shen
Undergraduate in Math. of Comp.

My research interests include theories and applications of reinforcement learning and machine learning.