# PlotLine3D

VTKExamples/Cxx/Plotting/PlotLine3D

### Description

This example illustrates plotting 3D points using vtkPlotLine3D. It plots the solution to the Lorenz Attractor.

Warning

The current vtkPlotLine3D produces the warning, vtkOpenGLContextDevice3D (): a line width has been requested that is larger than your system supports" which appears to have no effect on the results.

### Code

PlotLine3D.cxx

    #include <vtkNew.h>

    #include <vtkChartXYZ.h>
    #include <vtkPen.h>
    #include <vtkPlotLine3D.h>

    #include <vtkNamedColors.h>
    #include <vtkRenderWindow.h>
    #include <vtkTable.h>
    #include <vtkFloatArray.h>
    #include <vtkContextView.h>
    #include <vtkContextScene.h>
    #include <vtkContext3D.h>
    #include <vtkRenderWindowInteractor.h>
    #include <vtkRenderer.h>
    #include <vtkCamera.h>

    // Plot the solution to the Lorenz attractor.
    // http://en.wikipedia.org/wiki/Lorenz_system
    namespace
    {
    void lorenz(const float * varX, float * varXDerivative)
    {
      const float sigma = 10.f;
      const float rho = 28.f;
      const float beta = 2.66666666666f;

      varXDerivative[0] = sigma * (varX[1] - varX[0]);
      varXDerivative[1] = varX[0] * (rho - varX[2]) - varX[1];
      varXDerivative[2] = varX[0] * varX[1] - beta * varX[2];
    }
    } // end anonymous namespace

    //----------------------------------------------------------------------------
    int main(int, char * [])
    {
      vtkNew<vtkNamedColors> colors;

      // Create the data.
      vtkNew<vtkTable> varXSolution;

      vtkNew<vtkFloatArray> arrX0;
      arrX0->SetName("X");
      varXSolution->AddColumn(arrX0);

      vtkNew<vtkFloatArray> arrX1;
      arrX1->SetName("Y");
      varXSolution->AddColumn(arrX1);

      vtkNew<vtkFloatArray> arrX2;
      arrX2->SetName("Z");
      varXSolution->AddColumn(arrX2);

      const unsigned int numberOfTimePoints = 1000;
      varXSolution->SetNumberOfRows(numberOfTimePoints);
      float varX[3];
      varX[0] = 0.0f;
      varX[1] = 1.0f;
      varX[2] = 1.05f;

      float varXDerivative[3];
      const float deltaT = 0.01f;
      for (unsigned int ii = 0; ii < numberOfTimePoints; ++ii)
      {
        varXSolution->SetValue(ii, 0, varX[0]);
        varXSolution->SetValue(ii, 1, varX[1]);
        varXSolution->SetValue(ii, 2, varX[2]);
        lorenz(varX, varXDerivative);
        varX[0] += varXDerivative[0] * deltaT;
        varX[1] += varXDerivative[1] * deltaT;
        varX[2] += varXDerivative[2] * deltaT;
      }

      // Set up a 3D scene and add an XYZ chart to it.
      vtkNew<vtkContextView> view;
      view->GetRenderWindow()->SetSize(640, 480);

      vtkNew<vtkChartXYZ> chart;
      chart->SetGeometry(vtkRectf(5.0, 5.0, 635.0, 475.0));
      view->GetScene()->AddItem(chart);

      // Add a line plot.
      vtkNew<vtkPlotLine3D> plot;
      plot->SetInputData(varXSolution);
      plot->GetPen()->SetColorF(colors->GetColor3d("Tomato").GetData());
      view->GetRenderWindow()->SetMultiSamples(0);
      plot->GetPen()->SetWidth(6.0);
      chart->AddPlot(plot);

      // Finally render the scene.
      view->GetRenderer()->SetBackground(colors->GetColor3d("Mint").GetData());
      view->GetRenderWindow()->Render();
      view->GetInteractor()->Initialize();
      view->GetInteractor()->Start();

      return EXIT_SUCCESS;
    }

### CMakeLists.txt

    cmake_minimum_required(VERSION 3.3 FATAL_ERROR)

    project(PlotLine3D)

    find_package(VTK COMPONENTS
      vtkChartsCore
      vtkCommonColor
      vtkCommonCore
      vtkCommonDataModel
      vtkInteractionStyle
      vtkRenderingContext2D
      vtkRenderingContextOpenGL2
      vtkRenderingCore
      vtkRenderingFreeType
      vtkRenderingGL2PSOpenGL2
      vtkRenderingOpenGL2
      vtkViewsContext2D QUIET)
    if (NOT VTK_FOUND)
      message("Skipping PlotLine3D: ${VTK_NOT_FOUND_MESSAGE}")
      return ()
    endif()
    message (STATUS "VTK_VERSION: ${VTK_VERSION}")
    if (VTK_VERSION VERSION_LESS "8.90.0")
      # old system
      include(${VTK_USE_FILE})
      add_executable(PlotLine3D MACOSX_BUNDLE PlotLine3D.cxx )
      target_link_libraries(PlotLine3D PRIVATE ${VTK_LIBRARIES})
    else ()
      # include all components
      add_executable(PlotLine3D MACOSX_BUNDLE PlotLine3D.cxx )
      target_link_libraries(PlotLine3D PRIVATE ${VTK_LIBRARIES})
      # vtk_module_autoinit is needed
      vtk_module_autoinit(
        TARGETS PlotLine3D
        MODULES ${VTK_LIBRARIES}
        )
    endif ()

### Download and Build PlotLine3D

Click here to download PlotLine3D and its CMakeLists.txt file. Once the tarball PlotLine3D.tar has been downloaded and extracted,

    cd PlotLine3D/build

If VTK is installed:

    cmake ..

If VTK is not installed but compiled on your system, you will need to specify the path to your VTK build:

    cmake -DVTK_DIR:PATH=/home/me/vtk_build ..

Build the project:

    make

and run it:

    ./PlotLine3D

WINDOWS USERS

Be sure to add the VTK bin directory to your path. This will resolve the VTK dll's at run time.