#######################################################################
### GET THE CENTERLINE

### This first part, we use only if we can not compute the centerline fully but only a part of it
# vmtksurfacesubdivision -method loop -ifile \
# "...\P001Bs\CL\SVR.vtk"\
# --pipe vmtksurfacedecimation -reduction 0.5


vmtksurfacereader -ifile \
"...\P001Bs\CL\SVR.vtk"\
--pipe vmtkcenterlines --pipe vmtkrenderer --pipe vmtksurfaceviewer -opacity 0.25 --pipe vmtksurfaceviewer -i @vmtkcenterlines.o -array MaximumInscribedSphereRadius -ofile \
"...\P001Bs\CL\CVR.vtp"


### GET THE ATTRIBUTES

vmtkcenterlineattributes -ifile \
"...\P001Bs\CL\CVR.vtp" \
--pipe vmtkcenterlinegeometry -smoothing 1 -ofile \
"...\P001Bs\CL\CVR.dat"