Get Voxel array

[River-zy]

Hi ,I'm try Rivulet2 on other tree structures from 3D medical images.Finally I also got .swc file and .vtk file.I opened the .vtk file, and the 3D model is the same as Lung Airway Tracing example mentioned in the Readme, but I don’t know how to get the voxel array of the model.
Thanks for your suggestions

[lsqshr]

Thanks for your interests.

Not sure if I understood your question correctly. Do you mean to rasterize the vtk / swc result lines in a 3D numpy array? so that the centerline voxels are 1 and the rest are 0. This is not included in rivuletpy since the expected output here are lines defined in world coordinates (mm) in vtk and image coordinates in swc. However you could check out the make_vol function in https://github.com/RivuletStudio/Tuby/blob/master/tuby.py to get the idea of how to rasterize the 3D curves into a 3D voxel volume.

Pls let me know if I made sense.