public void WriteStructuredGrid(string filePath, List <Vector3d> point_data, Dictionary <String, vtkDataArray> scalar_dataArray, List <String> arrayNames, int numberOfPoints, int[] dimensions)
{
vtkPoints points = vtkPoints.New();
for (int i = 0; i < numberOfPoints; i++)
{
points.InsertNextPoint(point_data[i].X, point_data[i].Y, point_data[i].Z);
}
vtkStructuredGrid structuredGrid = vtkStructuredGrid.New();
structuredGrid.SetDimensions(dimensions[0], dimensions[1], dimensions[2]);
structuredGrid.SetPoints(points);
int numberOfScalarData = scalar_dataArray.Count();
for (int i = 0; i < numberOfScalarData; i++)
{
scalar_dataArray.TryGetValue(arrayNames[i], out vtkDataArray scalars);
structuredGrid.GetPointData().AddArray(scalars);
}
// add file ending if it is not existent
string suffix = (".vts");
if (!(filePath.EndsWith(suffix)))
{
filePath += suffix;
}
// Write file
vtkXMLStructuredGridWriter writer = vtkXMLStructuredGridWriter.New();
writer.SetFileName(filePath);
writer.SetInput(structuredGrid);
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("VTU file was writen and is saved at {0}", filePath);
}
public void WriteSimpleStructuredGrid(string filePath)
{
vtkPoints points = vtkPoints.New();
double x, y, z;
x = 0.0;
y = 0.0;
z = 0.0;
for (uint k = 0; k < 2; k++)
{
z += 2.0;
for (uint j = 0; j < 3; j++)
{
y += 1.0;
for (uint i = 0; i < 2; i++)
{
x += .5;
points.InsertNextPoint(x, y, z);
}
}
}
vtkDataArray dataArrayVx;
vtkDataArray dataArrayVy;
vtkDataArray dataArrayVz;
vtkDoubleArray Vx = new vtkDoubleArray();
vtkDoubleArray Vy = new vtkDoubleArray();
vtkDoubleArray Vz = new vtkDoubleArray();
for (int i = 0; i < 12; i++)
{
Vx.InsertNextValue(-1 + i);
Vz.InsertNextValue(2 - i);
Vy.InsertNextValue(-4 + i);
}
dataArrayVx = Vx;
dataArrayVx.SetName("Vx");
dataArrayVy = Vy;
dataArrayVy.SetName("Vy");
dataArrayVz = Vz;
dataArrayVz.SetName("Vz");
vtkStructuredGrid structuredGrid = vtkStructuredGrid.New();
structuredGrid.SetDimensions(2, 3, 2);
structuredGrid.SetPoints(points);
structuredGrid.GetPointData().AddArray(dataArrayVx);
structuredGrid.GetPointData().AddArray(dataArrayVy);
structuredGrid.GetPointData().AddArray(dataArrayVz);
// add file ending if it is not existent
string suffix = (".vts");
if (!(filePath.EndsWith(suffix)))
{
filePath += suffix;
}
// Write file
vtkXMLStructuredGridWriter writer = vtkXMLStructuredGridWriter.New();
writer.SetFileName(filePath);
writer.SetInput(structuredGrid);
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("VTU file was writen and is saved at {0}", filePath);
}
private void XMLStructuredGridWriter()
{
// 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\structuredgrid_test.vts");
// Create a grid
vtkStructuredGrid structuredGrid = vtkStructuredGrid.New();
vtkPoints points = vtkPoints.New();
points.InsertNextPoint(0, 0, 0);
points.InsertNextPoint(1, 0, 0);
points.InsertNextPoint(0, 1, 0);
points.InsertNextPoint(1, 1, 0);
points.InsertNextPoint(0, 2, 0);
points.InsertNextPoint(1, 2, 1);
// Specify the dimensions of the grid
structuredGrid.SetDimensions(2, 3, 1);
structuredGrid.SetPoints(points);
// Write file
vtkXMLStructuredGridWriter writer = vtkXMLStructuredGridWriter.New();
writer.SetFileName(filePath);
writer.SetInput(structuredGrid);
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;
}
reader.SetFileName(filePath);
reader.Update(); // here we read the file actually
vtkStructuredGridGeometryFilter geometryFilter = vtkStructuredGridGeometryFilter.New();
geometryFilter.SetInputConnection(reader.GetOutputPort());
geometryFilter.Update();
// Visualize
vtkPolyDataMapper mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection(geometryFilter.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);
}
