vmtkcenterlinesections doesn't work sometime


I am using vmtkcenterlinesections to calculate the cross-sectional area of my vessel. But sometimes the function brokes for certain vessel geometry… It says ‘executing vmtkcenterlinesections’, then directly exit without returning anything. I increase vmtkcenterlineresampling parameter from 1 to, for example, to 3. Then it works. However, this may cause to miss the minimum cross-sectional area due to the too sparse sampling on the centerline. Does anyone has any idea to solve the situation? Like, could it be the roughness of the surface that broke the function? (Cause mine vessel wall is pretty bumpy)
Thank you!

1 Like

Absolutely. If the vessel polydata quality is bad then anything can happen. All computational geometry algorithms are inherently unstable, so if you provide bad input then anything can happen.

Could you share (upload somewhere and post the link) your problematic mesh files? Or at least post a few screenshots?