Get Voxel array

[rivuletstudio]

Get Voxel array

[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.