Vessel reconstruction from centerlines/radius

Support
1

Hello everyone,

I am working on segmenting small cerebral arteries for future CFD studies.

After segmenting the small arteries that interest me, I am looking for a way to reconstruct an ideal (cylindrical) geometry from the average radius obtained for each segment/centerline with a smoothed intersection (blending type).

I have tried to create cylinders around each centerline with VTK librairy, but the intersections between cylinders are bad.

topic_1
topic_11140Ă—674 191 KB

Is there a method already implemented in 3D Slicer that could address this?

Thank you very much.

Gred

2

I’d suggest making a distance transform from the centerlines and then thresholding that to make the surface.

3

Hi, there is a vtk filter that should help a bit in this situation:

1 Like
4

Hi Pieper, by making a distance transform from the centerlines, do you mean somethink like polyball modelling along the centerlines and get the isosurface?

5

I’ve never used the polyball feature of vmtk, but maybe it would work. I was thinking you can just fill a volume with the distance to the closest vessel segment and threshold that. There are probably many ways of doing this depending on how efficient you need it to be.

6

Could you describe in more detail what you mean by “thresholding”? Perhaps a use case? I don’t quite understand what you mean, and I have a similar problem to the thread’s creator.

7

I mean using the threshold tool in the segment editor. Or for creating a mesh from a distance transform you could use the GrayscaleModelMaker directly. Just set a threshold for the isosurface distance you want.

8

I’m still stuck on 3D reconstruction with smoothing and physiological connections. I’ll take a look at greyscalemodel.

If you have any other ideas, please feel free to share them.

Thank you very much for your feedback.

9

Hello there,

I have finaly managed to obtain what I wanted by doing the following:

  1. creating cylinders around each centerline curve with a mean radius taken from the centerline quantification table
  2. converting to segment node and fill holes
  3. Finally converting it to model.

It is not the best and most obvious way to do it but the result is great.

2 Likes