Filter Gaussian Curvature get wrong results
To reproduce:
- Sources -> Box
- Filter -> Triangulate
- Filter -> Curvature
Gaussian Curvature of points in the corner of the cube should be PI / 2, it gives 28.274334
If we set properties(X/Y/Z Length) of the Box to 0.01, result becomes 282743.351463
I have checked the source in vtkCurvatures.cxx line 289, curvature assignment is gaussCurvatureData[v] = 3 * K[v] / dA[v] in which dA is total area of the neighboring triangles. This devision makes no sense in gaussian curvature computation (2 * PI - total angle).