A porosity model for medical image segmentation of vessels

A physics-based medical image segmentation method is developed. Specifically, the image greyscale intensity is used to infer the voxel partial volumes and subsequently formulate a porous medium analogy. The method involves first translating the medical image volumetric data into a three-dimensional computational domain of a porous material. A velocity field is then obtained from numerical simulations of incompressible fluid flow in the porous material, and finally a velocity iso-surface provides the surface description of the target object. The approach is tested on CT images of eight patient-specific cases, where cerebral aneurysms, nasal cavities (NC), and an aortic arch (AA) are the objects of interest. In the aneurysm cases, the results are compared against constant greyscale thresholding and manual segmentation. The manual segmentations of the aneurysms are validated by a clinical practitioner. Only a qualitative comparison is available for the NC, and the AA geometries. The results show that the proposed method is effective and capable of extracting the target object in a noisy domain. A sensitivity study is carried out to verify the method's performance with respect to modelling or user choices. The segmentation by the proposed method is also evaluated by performing computational fluid dynamics simulation, including a near-wall flow analysis, to ensure that the segmented geometry and the resulting computed solution are representative and meaningful.