# PlotLine3D

vtk-examples/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.

Question

### Code¶

PlotLine3D.cxx

#include <vtkCamera.h>
#include <vtkChartXYZ.h>
#include <vtkContext3D.h>
#include <vtkContextScene.h>
#include <vtkContextView.h>
#include <vtkFloatArray.h>
#include <vtkNamedColors.h>
#include <vtkNew.h>
#include <vtkPen.h>
#include <vtkPlotLine3D.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkTable.h>
#include <vtkVersion.h>

#if VTK_VERSION_NUMBER >= 90220220630ULL
#define VTK_HAS_SETCOLORF 1
#endif

// 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];
}
} // namespace

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

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

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

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

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

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);
view->GetRenderWindow()->SetWindowName("PlotLine3D");

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

vtkNew<vtkPlotLine3D> plot;
plot->SetInputData(varXSolution);
#if VTK_HAS_SETCOLORF
plot->GetPen()->SetColorF(colors->GetColor3d("LightCoral").GetData());
#else
plot->GetPen()->SetColor(colors->GetColor3d("LightCoral").GetData());
#endif
view->GetRenderWindow()->SetMultiSamples(0);
plot->GetPen()->SetWidth(2.0);

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

return EXIT_SUCCESS;
}


### CMakeLists.txt¶

cmake_minimum_required(VERSION 3.12 FATAL_ERROR)

project(PlotLine3D)

find_package(VTK COMPONENTS
ChartsCore
CommonColor
CommonCore
CommonDataModel
InteractionStyle
RenderingContext2D
RenderingContextOpenGL2
RenderingCore
RenderingFreeType
RenderingGL2PSOpenGL2
RenderingOpenGL2
ViewsContext2D
)

if (NOT VTK_FOUND)
message(FATAL_ERROR "PlotLine3D: Unable to find the VTK build folder.")
endif()

# Prevent a "command line is too long" failure in Windows.
set(CMAKE_NINJA_FORCE_RESPONSE_FILE "ON" CACHE BOOL "Force Ninja to use response files.")
target_link_libraries(PlotLine3D PRIVATE ${VTK_LIBRARIES} ) # vtk_module_autoinit is needed vtk_module_autoinit( TARGETS PlotLine3D MODULES${VTK_LIBRARIES}
)


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.