[vtk-developers] problem using vtkcurvatures

[akashlad]

Hi everybody! 
i'm new to VTK and i'm trying to segment the brain sulci and gyri. i am
using vtk on 3D slicer.
the code that i wrote so far is : 

node=slicer.util.getNode('grey')
vol=node.GetPolydata()
curve=vtk.vtkCurvatures()
curve.SetCurvatureTypeToMean()
curve.SetInputData(vol)
curve.Update()
curved_vol=curve.GetOutput()

th=vtk.vtkThreshold()
th.SetInputData(curved_vol)
th.thresholdByupper(0)
th.Update()

surf=vtk.vtkDataSetSurfaceFilter()
surf.SetInputData(th.GetOutput())
surf.Update()
output_volume=surf.GetOutput()

after applying the curvature my idea was to use a threshold filter to cut
out the sulci from the model and obtain only the gyri. 
the problem is that the curvatures filter does not work very well with my
model and i don't know why. the result that i get is this 
<http://vtk.1045678.n5.nabble.com/file/n5738159/13275218_10208217634112196_787516635_o.png> 

the problem is that the filter seems to not recognize well the gyri in
respect to the sulci
here you can find my model:
https://www.dropbox.com/s/z99x7s7tj5is01a/grey.ply?dl=0
So, basically, now the idea is to make the filter work on a major number of
triangles of the model mesh, so that it won't recognise the small dips on
the gyri. How could I do that?
Thanks for any help! 

best regards

Akash Lad



--
View this message in context: http://vtk.1045678.n5.nabble.com/problem-using-vtkcurvatures-tp5738159.html
Sent from the VTK - Dev mailing list archive at Nabble.com.