Пример трехмерного контурного графика С++ с использованием VTK?

Некоторое время назад я хотел сделать то же самое и провел обширный поиск. Это не так уж много, поэтому я хочу поделиться с вами своим решением. Он содержит лишний код, но работает. http://pastie.org/4721543

Я только что заново открыл этот вопрос, потому что в VTK 6.0 будет функция plotsurface. Что, скорее всего, будет гораздо более производительным, чем мой код. Я попробую позже и обновлю этот пост.

Небольшое объяснение моего кода может быть полезным. Он делает следующее:

1. создайте плоскость размером с вашу 2D-матрицу, которая содержит ваши значения z (столбцы и строки представляют значения x и y).
2. создать xyz-точки из входной матрицы
3. сохранить значения z как скаляры для каждой точки
4. деформируйте эту плоскость в соответствии со значениями z матрицы и примените цвета (mapper->SetScalarRange(..)
5. нарисуй свой сюжет

Он также добавляет таймер и команду обновления, выводит статистику рендеринга и т. д., как я уже сказал, многое можно опустить. В моем сообщении в блоге рассказывается немного больше на этом, но не так много.

Я считаю, что грубый набросок конвейера VTK, который может вам подойти, выглядит следующим образом. Если вы создадите несколько точек x и y как vtkPoints, а затем переместите их в тип данных StructuredGrid, а затем добавите ось Z с помощью GeometryFilter, вы должны создать поверхность, пройдя через vtkWarpScalar и обычный конвейер Mapper и Actor.

Следующий пример кода может помочь /Examples/DataManipulation/Cxx/SGrid.cxx в вашей сборке vtk. Поверхностное построение также может быть выполнено, если у вас есть полный объем и извлеченная изоповерхность. Если это не так, и вы строите 3D-функцию, вероятно, лучше всего использовать структурированную сетку.

С помощью кода Кайкухна я смог создать простой поверхностный график. Надеюсь, это поможет кому-то, кто ищет подобное решение.

#include "stdafx.h"


#include <vtkRenderer.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderWindow.h>
#include <vtkSmartPointer.h>
#include <vtkPoints.h>
#include <vtkPointData.h>
#include <vtkPolyData.h>
#include <vtkPolyDataMapper.h>
#include <vtkDataSetMapper.h>
#include <vtkProperty.h>
#include <vtkCubeAxesActor2D.h>
#include <vtkInteractorStyleTrackballCamera.h>
#include <vtkInteractorStyleTrackball.h>
#include <vtkSimplePointsReader.h>
#include <vtkWarpScalar.h>
#include <vtkAxisActor2D.h>


int mains(int, char *[])
{


// Read the file
  vtkSmartPointer<vtkSimplePointsReader> reader =vtkSmartPointer<vtkSimplePointsReader>::New();
  reader->SetFileName ( "simple.xyz" );
  reader->Update();

  vtkSmartPointer<vtkPolyData> inputPolyData = vtkSmartPointer<vtkPolyData>::New();
  inputPolyData ->CopyStructure(reader->GetOutput());


  // warp plane
  vtkSmartPointer<vtkWarpScalar> warp = vtkSmartPointer<vtkWarpScalar>::New();
  warp->SetInput(inputPolyData);
  warp->SetScaleFactor(0.0);

  // Visualize
  vtkSmartPointer<vtkDataSetMapper> mapper = vtkSmartPointer<vtkDataSetMapper>::New();
  mapper->SetInputConnection(warp->GetOutputPort());



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

  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(.3, .6, .3);
  renderWindow->Render();

  vtkSmartPointer<vtkInteractorStyleTrackballCamera> style = vtkSmartPointer<vtkInteractorStyleTrackballCamera>::New();
  renderWindowInteractor->SetInteractorStyle(style);

  // add & render CubeAxes
  vtkSmartPointer<vtkCubeAxesActor2D> axes = vtkSmartPointer<vtkCubeAxesActor2D>::New();
  axes->SetInput(warp->GetOutput());
  axes->SetFontFactor(3.0);
  axes->SetFlyModeToNone();
  axes->SetCamera(renderer->GetActiveCamera());

  vtkSmartPointer<vtkAxisActor2D> xAxis = axes->GetXAxisActor2D();
  xAxis->SetAdjustLabels(1);

  renderer->AddViewProp(axes);
  renderWindowInteractor->Start();

  return EXIT_SUCCESS;
}

простой .xyz;

0.0 0.0 3.0
1.0 4.0 0.0
0.0 1.0 0.0
4.0 0.0 3.0
1.0 4.0 7.0
0.0 6.0 0.0

В моей системе код работал, когда SetInput() заменялся на SetInputData() для «осей» и «деформации».