public void WriteVTP(string filePath, List <Vector3d> point_data, List <Vector3> cellPointsList, Dictionary <String, vtkDataArray> scalar_dataArray, List <String> arrayNames)
{
vtkPoints points = vtkPoints.New();
for (int i = 0; i < point_data.Count(); i++)
{
points.InsertNextPoint(point_data[i].X, point_data[i].Y, point_data[i].Z);
}
vtkPolyData polyData = vtkPolyData.New();
polyData.SetPoints(points);
vtkCellArray cellArrayOne = vtkCellArray.New();
for (int i = 0; i < cellPointsList.Count(); i++)
{
vtkTetra tetra = vtkTetra.New();
tetra.GetPointIds().SetId(0, (long)cellPointsList[i][0]);
tetra.GetPointIds().SetId(1, (long)cellPointsList[i][1]);
tetra.GetPointIds().SetId(2, (long)cellPointsList[i][2]);
tetra.GetPointIds().SetId(3, (long)cellPointsList[i][2]);
cellArrayOne.InsertNextCell(tetra);
}
polyData.SetPolys(cellArrayOne);
int numberOfScalarData = scalar_dataArray.Count();
for (int i = 0; i < numberOfScalarData; i++)
{
scalar_dataArray.TryGetValue(arrayNames[i], out vtkDataArray scalars);
polyData.GetPointData().AddArray(scalars);
}
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
// add file ending if it is not existent
string suffix = (".vtp");
if (!(filePath.EndsWith(suffix)))
{
filePath += suffix;
}
writer.SetFileName(filePath);
writer.SetInput(polyData);
writer.Write();
// Read and display file for verification that it was written correctly
vtkXMLStructuredGridReader reader = vtkXMLStructuredGridReader.New();
if (reader.CanReadFile(filePath) == 0)
{
//MessageBox.Show("Cannot read file \"" + filePath + "\"", "Error", MessageBoxButtons.OK);
return;
}
Console.WriteLine("VTP file was writen and is saved at {0}", filePath);
}
internal void ApplyTransform(ref vtkTransform transform, string filename)
{
vtkPolyData polydata = vtkPolyData.New();
CreatePolydata(ref polydata);
vtkTransformFilter transformFilter = vtkTransformFilter.New();
#if VTK_MAJOR_VERSION_5
transformFilter.SetInputConnection(polydata.GetProducerPort());
#else
transformFilter.SetInputData(polydata);
#endif
transformFilter.SetTransform(transform);
transformFilter.Update();
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetFileName(filename);
#if VTK_MAJOR_VERSION_5
writer.SetInputConnection(transformFilter.GetOutputPort());
#else
writer.SetInputData(transformFilter);
#endif
writer.Write();
}
private void MatrixMathFilter()
{
// Path to vtk data must be set as an environment variable
// VTK_DATA_ROOT = "C:\VTK\vtkdata-5.8.0"
vtkTesting test = vtkTesting.New();
string root = test.GetDataRoot();
string filePath = System.IO.Path.Combine(root, @"Data\tensors.vtk");
vtkUnstructuredGridReader reader = vtkUnstructuredGridReader.New();
reader.SetFileName(filePath);
reader.Update();
vtkDataSetSurfaceFilter surfaceFilter = vtkDataSetSurfaceFilter.New();
surfaceFilter.SetInputConnection(reader.GetOutputPort());
surfaceFilter.Update();
vtkMatrixMathFilter matrixMathFilter = vtkMatrixMathFilter.New();
//matrixMathFilter.SetOperationToDeterminant();
matrixMathFilter.SetOperationToEigenvalue();
matrixMathFilter.SetInputConnection(surfaceFilter.GetOutputPort());
matrixMathFilter.Update();
matrixMathFilter.GetOutput().GetPointData().SetActiveScalars("Eigenvalue");
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetInputConnection(matrixMathFilter.GetOutputPort());
writer.SetFileName(System.IO.Path.Combine(root, @"Data\output.vtp"));
writer.Write();
vtkPolyDataMapper mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection(matrixMathFilter.GetOutputPort());
// actor
vtkActor actor = vtkActor.New();
actor.SetMapper(mapper);
// get a reference to the renderwindow of our renderWindowControl1
vtkRenderWindow renderWindow = renderWindowControl1.RenderWindow;
// renderer
vtkRenderer renderer = renderWindow.GetRenderers().GetFirstRenderer();
// set background color
renderer.SetBackground(0.2, 0.3, 0.4);
// add our actor to the renderer
renderer.AddActor(actor);
}
internal void Write(string filename)
{
vtkPolyData polydata = vtkPolyData.New();
CreatePolydata(ref polydata);
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetFileName(filename);
#if VTK_MAJOR_VERSION_5
writer.SetInputConnection(polydata.GetProducerPort());
#else
writer.SetInputData(polydata);
#endif
writer.Write();
}
private void Curvature()
{
// Path to vtk data must be set as an environment variable
// VTK_DATA_ROOT = "C:\VTK\vtkdata-5.8.0"
vtkTesting test = vtkTesting.New();
string root = test.GetDataRoot();
// note: one should use a coarse surface mesh to achieve a good visual effect
string filePath = System.IO.Path.Combine(root, @"Data\bunny.vtp");
// Create a polydata
vtkXMLPolyDataReader reader = vtkXMLPolyDataReader.New();
reader.SetFileName(filePath);
vtkCurvatures curvaturesFilter = vtkCurvatures.New();
curvaturesFilter.SetInputConnection(reader.GetOutputPort());
//curvaturesFilter.SetCurvatureTypeToGaussian();
//curvaturesFilter.SetCurvatureTypeToMean();
//curvaturesFilter.SetCurvatureTypeToMaximum();
curvaturesFilter.SetCurvatureTypeToMinimum();
// To inspect more closely, if required
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetInputConnection(curvaturesFilter.GetOutputPort());
writer.SetFileName(System.IO.Path.Combine(root, @"Data\gauss.vtp"));
writer.Write();
// Create a mapper and actor
vtkPolyDataMapper mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection(curvaturesFilter.GetOutputPort());
double[] range = curvaturesFilter.GetOutput().GetScalarRange();
mapper.SetScalarRange(range[0], range[1]);
vtkActor actor = vtkActor.New();
actor.SetMapper(mapper);
// get a reference to the renderwindow of our renderWindowControl1
vtkRenderWindow renderWindow = renderWindowControl1.RenderWindow;
// renderer
vtkRenderer renderer = renderWindow.GetRenderers().GetFirstRenderer();
// set background color
renderer.SetBackground(0.0, 0.0, 0.0);
// add our actor to the renderer
renderer.AddActor(actor);
}
private static void FindAllArrayNames(string filePath)
{
vtkPolyData polydata = vtkPolyData.New();
if (filePath == null)
{
vtkSphereSource sphereSource = vtkSphereSource.New();
sphereSource.Update();
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetFileName(@"c:\vtk\vtkdata-5.8.0\Data\testFindAllArrayNames.vtp");
writer.SetInputConnection(sphereSource.GetOutputPort());
writer.Write();
polydata = sphereSource.GetOutput();
}
else
{
vtkXMLPolyDataReader reader = vtkXMLPolyDataReader.New();
reader.SetFileName(filePath);
reader.Update();
polydata = reader.GetOutput();
}
FindAllData(ref polydata);
}
private void WritePolyData()
{
// Path to vtk data must be set as an environment variable
// VTK_DATA_ROOT = "C:\VTK\vtkdata-5.8.0"
vtkTesting test = vtkTesting.New();
string root = test.GetDataRoot();
string filePath = System.IO.Path.Combine(root, @"Data\poly_test.vtp");
// Create 4 points for a tetrahedron
vtkPoints points = vtkPoints.New();
points.InsertNextPoint(0, 0, 0);
points.InsertNextPoint(0, 0, 1);
points.InsertNextPoint(1, 0, 0);
points.InsertNextPoint(0, 1, 0);
// Create a polydata object and add the points to it.
vtkPolyData polydata = vtkPolyData.New();
polydata.SetPoints(points);
// it's not enough only to define points
// we need to define faces too
// (must be defined in counter clockwise order as viewed from the outside)
vtkTriangle face0 = vtkTriangle.New();
face0.GetPointIds().SetId(0, 0);
face0.GetPointIds().SetId(1, 2);
face0.GetPointIds().SetId(2, 1);
vtkTriangle face1 = vtkTriangle.New();
face1.GetPointIds().SetId(0, 0);
face1.GetPointIds().SetId(1, 3);
face1.GetPointIds().SetId(2, 2);
vtkTriangle face2 = vtkTriangle.New();
face2.GetPointIds().SetId(0, 0);
face2.GetPointIds().SetId(1, 1);
face2.GetPointIds().SetId(2, 3);
vtkTriangle face3 = vtkTriangle.New();
face3.GetPointIds().SetId(0, 1);
face3.GetPointIds().SetId(1, 2);
face3.GetPointIds().SetId(2, 3);
vtkCellArray faces = vtkCellArray.New();
faces.InsertNextCell(face0);
faces.InsertNextCell(face1);
faces.InsertNextCell(face2);
faces.InsertNextCell(face3);
polydata.SetPolys(faces);
// Write the file
vtkXMLPolyDataWriter writer = vtkXMLPolyDataWriter.New();
writer.SetFileName(filePath);
writer.SetInput(polydata);
// Optional - set the mode. The default is binary.
//writer.SetDataModeToBinary();
writer.SetDataModeToAscii();
writer.Write();
// Read and display file for verification that it was written correctly
vtkXMLPolyDataReader reader = vtkXMLPolyDataReader.New();
if (reader.CanReadFile(filePath) == 0)
{
MessageBox.Show("Cannot read file \"" + filePath + "\"", "Error", MessageBoxButtons.OK);
return;
}
reader.SetFileName(filePath);
reader.Update();
vtkDataSetMapper mapper = vtkDataSetMapper.New();
mapper.SetInputConnection(reader.GetOutputPort());
// actor
vtkActor actor = vtkActor.New();
actor.SetMapper(mapper);
// get a reference to the renderwindow of our renderWindowControl1
vtkRenderWindow renderWindow = renderWindowControl1.RenderWindow;
// renderer
vtkRenderer renderer = renderWindow.GetRenderers().GetFirstRenderer();
// set background color
renderer.SetBackground(0.2, 0.3, 0.4);
// add our actor to the renderer
renderer.AddActor(actor);
}
