This doesn't appear to be working for me. I have vtkPolyData that is the output of a vtkContourWidget so it only exists in 2 planes. Here's the code that I'm currently running<div><br></div><div><div> vtkCurvatures curvatures = new vtk.vtkCurvatures();</div>
<div> curvatures.SetInput(currentPolyData);</div><div> curvatures.SetCurvatureTypeToGaussian();</div><div> curvatures.Update();</div><div> vtkPolyData curvatureData = curvatures.GetOutput();</div>
</div><div> for (int j=0; j<curvatureData.GetNumberOfPoints(); j++) {</div><div> double curv = curvatureData.GetPointData().GetArray(0).GetTuple1(j);</div><div> ...</div>
<div>However, an exception is being thrown because the curvatureData.GetPointData().GetNumberOfArrays() method is returning 0. In your example I'm not sure what the meanCurvature array is. I don't want to visualize the curvature data, just get the curvature at each point and iterate over the values and compare them. The full goal is to get the point index for the n points in currentPolyData with the highest curvature.</div>
<div><br><br><div class="gmail_quote">On Tue, Oct 19, 2010 at 7:25 PM, David Doria <span dir="ltr"><<a href="mailto:daviddoria@gmail.com">daviddoria@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div></div><div class="h5">On Tue, Oct 19, 2010 at 10:04 PM, Jonathan Morra <<a href="mailto:jonmorra@gmail.com">jonmorra@gmail.com</a>> wrote:<br>
> I have a vtkPolyData and I would like to get the curvature at every point.<br>
> I noticed that there is a vtkCurvatures class. Here's the code that I have<br>
> written so far<br>
> vtkCurvatures curvatures = new vtk.vtkCurvatures();<br>
> curvatures.SetInput(currentPolyData);<br>
> curvatures.SetCurvatureTypeToGaussian();<br>
> curvatures.Update();<br>
> vtkPolyData curvatureData = curvatures.GetOutput();<br>
> However, I don't know what to do with the curvatureData vtkPolyData, and how<br>
> to extract the curvature at every point.<br>
> Any help would be greatly appreciated.<br>
> Thanks<br>
<br>
</div></div>Does this example help?<br>
<br>
<a href="http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/Curvatures" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/Curvatures</a><br>
<br>
If you look in Graphics/vtkCurvatures.cxx, you will see the names of<br>
the arrays that you would want to get: "Minimum_Curvature",<br>
"Maximum_Curvature", etc. And you will see that they are stored in the<br>
PointData of the PolyData object:<br>
<br>
mesh->GetPointData()->AddArray(meanCurvature);<br>
mesh->GetPointData()->SetActiveScalars("Mean_Curvature");<br>
<font color="#888888"><br>
David<br>
</font></blockquote></div><br></div>