[vtkusers] Problem with "Compute Gaussian, Mean, Minm and Max Curvatures" example

[Sky77]

It seems that my version of VTK is to old. The header file
vtkBooleanOperationPolyDataFilter.h can't be found. For this purpose, i have
try something else to get an object with negative Gaussian curvature. A
torus have positive Gaussian curvature on the outside and negative Gaussian
curvature on the inside. I have created one with blender and import it as
obj file. 

This ist the input file:  torus.obj
<http://vtk.1045678.n5.nabble.com/file/n5728401/torus.obj>   

Here <http://vtk.1045678.n5.nabble.com/file/n5728401/result.png>   you can
see my result. Also in this case no negative curvature is computed.

Has someone of you experience with the curvature computation with vtk? I
mean has someone used this functionality with success?

#include <vtkSmartPointer.h>
#include <vtkCurvatures.h>
#include <vtkXMLPolyDataReader.h>
#include <vtkLookupTable.h>
#include <vtkColorTransferFunction.h>
#include <vtkColorSeries.h>
#include <vtkPointData.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkScalarBarActor.h>
#include <vtkRenderWindow.h>
#include <vtkRenderer.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkOBJReader.h>
#include <vtkUnstructuredGrid.h>
#include <vtkCell.h>
#include <vtkCellArray.h>
#include <vtkIdList.h>
#include <vtkUnsignedCharArray.h>
#include <string>

int main(int argc, char *argv[])
{
    // Parse command line arguments
    if(argc != 2)
    {
        std::cout << "Usage: " << argv[0] << " Filename(.obj)" << std::endl;
        return EXIT_FAILURE;
    }

    std::string filename = argv[1];
    vtkSmartPointer<vtkOBJReader> reader =
vtkSmartPointer<vtkOBJReader>::New();
    reader->SetFileName(filename.c_str());
    reader->Update();

    vtkSmartPointer<vtkCurvatures> curvaturesFilter =
vtkSmartPointer<vtkCurvatures>::New();
    curvaturesFilter->SetInputConnection(reader->GetOutputPort());
    curvaturesFilter->SetCurvatureTypeToGaussian();
    curvaturesFilter->Update();

    // Get scalar range from command line if present, otherwise use
    // range of computed curvature
    double scalarRange[2];
    if (argc >= 4)
    {
      scalarRange[0] = atof(argv[2]);
      scalarRange[1] = atof(argv[3]);
    }
    else
    {
      curvaturesFilter->GetOutput()->GetScalarRange(scalarRange);
    }

    int scheme = 16;
    if (argc >= 5)
    {
      scheme = atoi(argv[4]);
    }

    // Build a lookup table
    vtkSmartPointer<vtkColorSeries> colorSeries =
vtkSmartPointer<vtkColorSeries>::New();
    colorSeries->SetColorScheme(scheme);

    vtkSmartPointer<vtkColorTransferFunction> lut =
vtkSmartPointer<vtkColorTransferFunction>::New();
    lut->SetColorSpaceToHSV();

    // Use a color series to create a transfer function
    int numColors = colorSeries->GetNumberOfColors();
    for (int i = 0; i < numColors; i++)
    {
      vtkColor3ub color = colorSeries->GetColor(i);
      double dColor[3];
      dColor[0] = static_cast<double> (color[0]) / 255.0;
      dColor[1] = static_cast<double> (color[1]) / 255.0;
      dColor[2] = static_cast<double> (color[2]) / 255.0;
      double t = scalarRange[0] + (scalarRange[1] - scalarRange[0]) /
(numColors - 1) * i;
      lut->AddRGBPoint(t, dColor[0], dColor[1], dColor[2]);
    }

    // Create a mapper and actor
    vtkSmartPointer<vtkPolyDataMapper> mapper =
vtkSmartPointer<vtkPolyDataMapper>::New();
    mapper->SetInputConnection(curvaturesFilter->GetOutputPort());
    //    mapper->SetLookupTable(lut);
    mapper->SetScalarRange(scalarRange);

    vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
    actor->SetMapper(mapper);

    // Create a scalar bar
    vtkSmartPointer<vtkScalarBarActor> scalarBar =
vtkSmartPointer<vtkScalarBarActor>::New();
    scalarBar->SetLookupTable(mapper->GetLookupTable());
   
scalarBar->SetTitle(curvaturesFilter->GetOutput()->GetPointData()->GetScalars()->GetName());
    scalarBar->SetNumberOfLabels(5);

    // Create a renderer, render window, and interactor
    vtkSmartPointer<vtkRenderer> renderer =
vtkSmartPointer<vtkRenderer>::New();
    vtkSmartPointer<vtkRenderWindow> renderWindow =
vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);
    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
vtkSmartPointer<vtkRenderWindowInteractor>::New();
    renderWindowInteractor->SetRenderWindow(renderWindow);

    // Add the actors to the scene
    renderer->AddActor(actor);
    renderer->AddActor2D(scalarBar);

    renderer->SetBackground(.1, .2, .3); // Background color blue

    // Render and interact
    renderWindow->Render();
    renderWindowInteractor->Start();

    return EXIT_SUCCESS;

}



--
View this message in context: http://vtk.1045678.n5.nabble.com/Problem-with-Compute-Gaussian-Mean-Minm-and-Max-Curvatures-example-tp5728377p5728401.html
Sent from the VTK - Users mailing list archive at Nabble.com.