Визуализация нормальных векторов - PullRequest
Новая
       39

Визуализация нормальных векторов

0 голосов
/ 15 декабря 2018

У меня есть несколько файлов облаков точек в формате .xyz.Каждый из этих файлов содержит приблизительно 1 миллион точек и соответствующие нормальные векторы для каждой точки.Каждая строка файлов .xyz принадлежит одной точке, и каждая строка содержит декартовы координаты (x, y и z) и нормальный вектор (nx, ny и nz).Моя цель - показать точки и их нормальные векторы.Я могу успешно читать файлы .xyz и успешно показывать очки.Тем не менее, я не могу показать нормальные векторы успешно.Я пытался решить проблему с помощью класса vtkGlyph3D, как показано в некоторых примерах Kitware, но каждая попытка была неудачной для меня.Мой код показан ниже.Не могли бы вы помочь мне решить эту проблему?Я с нетерпением жду ответа от вас.Заранее спасибо 

#include <vtkVersion.h>
#include <vtkSmartPointer.h>
#include <vtkProperty.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkDelimitedTextReader.h>
#include <vtkDoubleArray.h>
#include <vtkTable.h>
#include <vtkPointData.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkVertexGlyphFilter.h>

int main(int argc, char* argv[])
{
    vtkSmartPointer<vtkDelimitedTextReader> reader = vtkSmartPointer<vtkDelimitedTextReader>::New();
    reader->SetFileName("points_and_normals.xyz");
    reader->DetectNumericColumnsOn();
    reader->SetFieldDelimiterCharacters(" ");
    reader->Update();

    vtkTable* table = reader->GetOutput();
    vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
    vtkSmartPointer<vtkDoubleArray> normals = vtkSmartPointer<vtkDoubleArray>::New();
    normals->SetNumberOfComponents(3); //3d normals (ie x,y,z)

    std::cout << "Table has " << table->GetNumberOfRows()
            << " rows." << std::endl;
    std::cout << "Table has " << table->GetNumberOfColumns()
            << " columns." << std::endl;

    for (vtkIdType i = 0; i < table->GetNumberOfRows(); i++)
    {

            points->InsertNextPoint((table->GetValue(i, 0)).ToDouble(),
                    (table->GetValue(i, 1)).ToDouble(),
                    (table->GetValue(i, 2)).ToDouble());

            double n[3];
            n[0] = (table->GetValue(i, 3)).ToDouble();
            n[1] = (table->GetValue(i, 4)).ToDouble();
            n[2] = (table->GetValue(i, 5)).ToDouble();
            normals->InsertNextTuple(n);
    }

    std::cout << "There are " << points->GetNumberOfPoints()
            << " points." << std::endl;

    vtkSmartPointer<vtkPolyData> polydata = vtkSmartPointer<vtkPolyData>::New();
    polydata->SetPoints(points);
    polydata->GetPointData()->SetNormals(normals);
    vtkSmartPointer<vtkVertexGlyphFilter> glyphFilter = vtkSmartPointer<vtkVertexGlyphFilter>::New();
#if VTK_MAJOR_VERSION <= 5
        glyphFilter->SetInputConnection(polydata->GetProducerPort());
#else
        glyphFilter->SetInputData(polydata);
#endif
    glyphFilter->Update();
    // Visualize
    vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
    mapper->SetInputConnection(glyphFilter->GetOutputPort());
    vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
    actor->SetMapper(mapper);
    actor->GetProperty()->SetPointSize(3);
    actor->GetProperty()->SetColor(0, 0, 1);
    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);
    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    renderWindowInteractor->SetRenderWindow(renderWindow);
    renderer->AddActor(actor);
    renderer->SetBackground(.5, .5, .5);
    renderWindow->Render();
    renderWindowInteractor->Start();
    return EXIT_SUCCESS;
}
...