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);
}
/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVfieldToUGrid(String [] argv)
{
//Prefix Content is: ""
// Read a field representing unstructured grid and display it (similar to blow.tcl)[]
// create a reader and write out field daya[]
reader = new vtkUnstructuredGridReader();
reader.SetFileName((string)"" + (VTK_DATA_ROOT.ToString()) + "/Data/blow.vtk");
reader.SetScalarsName((string)"thickness9");
reader.SetVectorsName((string)"displacement9");
ds2do = new vtkDataSetToDataObjectFilter();
ds2do.SetInputConnection((vtkAlgorithmOutput)reader.GetOutputPort());
// we must be able to write here[]
try
{
channel = new StreamWriter("UGridField.vtk");
tryCatchError = "NOERROR";
}
catch(Exception)
{tryCatchError = "ERROR";}
if (tryCatchError.Equals("NOERROR"))
{
channel.Close();
write = new vtkDataObjectWriter();
write.SetInputConnection((vtkAlgorithmOutput)ds2do.GetOutputPort());
write.SetFileName((string)"UGridField.vtk");
write.Write();
// Read the field and convert to unstructured grid.[]
dor = new vtkDataObjectReader();
dor.SetFileName((string)"UGridField.vtk");
do2ds = new vtkDataObjectToDataSetFilter();
do2ds.SetInputConnection((vtkAlgorithmOutput)dor.GetOutputPort());
do2ds.SetDataSetTypeToUnstructuredGrid();
do2ds.SetPointComponent((int)0,(string)"Points",(int)0);
do2ds.SetPointComponent((int)1,(string)"Points",(int)1);
do2ds.SetPointComponent((int)2,(string)"Points",(int)2);
do2ds.SetCellTypeComponent((string)"CellTypes",(int)0);
do2ds.SetCellConnectivityComponent((string)"Cells",(int)0);
fd2ad = new vtkFieldDataToAttributeDataFilter();
fd2ad.SetInput((vtkDataObject)do2ds.GetUnstructuredGridOutput());
fd2ad.SetInputFieldToDataObjectField();
fd2ad.SetOutputAttributeDataToPointData();
fd2ad.SetVectorComponent((int)0,(string)"displacement9",(int)0);
fd2ad.SetVectorComponent((int)1,(string)"displacement9",(int)1);
fd2ad.SetVectorComponent((int)2,(string)"displacement9",(int)2);
fd2ad.SetScalarComponent((int)0,(string)"thickness9",(int)0);
// Now start visualizing[]
warp = new vtkWarpVector();
warp.SetInput((vtkDataObject)fd2ad.GetUnstructuredGridOutput());
// extract mold from mesh using connectivity[]
connect = new vtkConnectivityFilter();
connect.SetInputConnection((vtkAlgorithmOutput)warp.GetOutputPort());
connect.SetExtractionModeToSpecifiedRegions();
connect.AddSpecifiedRegion((int)0);
connect.AddSpecifiedRegion((int)1);
moldMapper = new vtkDataSetMapper();
moldMapper.SetInputConnection((vtkAlgorithmOutput)connect.GetOutputPort());
moldMapper.ScalarVisibilityOff();
moldActor = new vtkActor();
moldActor.SetMapper((vtkMapper)moldMapper);
moldActor.GetProperty().SetColor((double).2,(double).2,(double).2);
moldActor.GetProperty().SetRepresentationToWireframe();
// extract parison from mesh using connectivity[]
connect2 = new vtkConnectivityFilter();
connect2.SetInputConnection((vtkAlgorithmOutput)warp.GetOutputPort());
connect2.SetExtractionModeToSpecifiedRegions();
connect2.AddSpecifiedRegion((int)2);
parison = new vtkGeometryFilter();
parison.SetInputConnection((vtkAlgorithmOutput)connect2.GetOutputPort());
normals2 = new vtkPolyDataNormals();
normals2.SetInputConnection((vtkAlgorithmOutput)parison.GetOutputPort());
normals2.SetFeatureAngle((double)60);
lut = new vtkLookupTable();
lut.SetHueRange((double)0.0,(double)0.66667);
parisonMapper = vtkPolyDataMapper.New();
parisonMapper.SetInputConnection((vtkAlgorithmOutput)normals2.GetOutputPort());
parisonMapper.SetLookupTable((vtkScalarsToColors)lut);
parisonMapper.SetScalarRange((double)0.12,(double)1.0);
parisonActor = new vtkActor();
parisonActor.SetMapper((vtkMapper)parisonMapper);
cf = new vtkContourFilter();
cf.SetInputConnection((vtkAlgorithmOutput)connect2.GetOutputPort());
cf.SetValue((int)0,(double).5);
contourMapper = vtkPolyDataMapper.New();
contourMapper.SetInputConnection((vtkAlgorithmOutput)cf.GetOutputPort());
contours = new vtkActor();
contours.SetMapper((vtkMapper)contourMapper);
// Create graphics stuff[]
ren1 = vtkRenderer.New();
renWin = vtkRenderWindow.New();
renWin.AddRenderer((vtkRenderer)ren1);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
// Add the actors to the renderer, set the background and size[]
ren1.AddActor((vtkProp)moldActor);
ren1.AddActor((vtkProp)parisonActor);
ren1.AddActor((vtkProp)contours);
ren1.ResetCamera();
ren1.GetActiveCamera().Azimuth((double)60);
ren1.GetActiveCamera().Roll((double)-90);
ren1.GetActiveCamera().Dolly((double)2);
ren1.ResetCameraClippingRange();
ren1.SetBackground((double)1,(double)1,(double)1);
renWin.SetSize((int)375,(int)200);
iren.Initialize();
}
// prevent the tk window from showing up then start the event loop[]
//deleteAllVTKObjects();
}
/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVfieldToUGrid(String [] argv)
{
//Prefix Content is: ""
// Read a field representing unstructured grid and display it (similar to blow.tcl)[]
// create a reader and write out field daya[]
reader = new vtkUnstructuredGridReader();
reader.SetFileName((string)"" + (VTK_DATA_ROOT.ToString()) + "/Data/blow.vtk");
reader.SetScalarsName((string)"thickness9");
reader.SetVectorsName((string)"displacement9");
ds2do = new vtkDataSetToDataObjectFilter();
ds2do.SetInputConnection((vtkAlgorithmOutput)reader.GetOutputPort());
// we must be able to write here[]
try
{
channel = new StreamWriter("UGridField.vtk");
tryCatchError = "NOERROR";
}
catch (Exception)
{ tryCatchError = "ERROR"; }
if (tryCatchError.Equals("NOERROR"))
{
channel.Close();
write = new vtkDataObjectWriter();
write.SetInputConnection((vtkAlgorithmOutput)ds2do.GetOutputPort());
write.SetFileName((string)"UGridField.vtk");
write.Write();
// Read the field and convert to unstructured grid.[]
dor = new vtkDataObjectReader();
dor.SetFileName((string)"UGridField.vtk");
do2ds = new vtkDataObjectToDataSetFilter();
do2ds.SetInputConnection((vtkAlgorithmOutput)dor.GetOutputPort());
do2ds.SetDataSetTypeToUnstructuredGrid();
do2ds.SetPointComponent((int)0, (string)"Points", (int)0);
do2ds.SetPointComponent((int)1, (string)"Points", (int)1);
do2ds.SetPointComponent((int)2, (string)"Points", (int)2);
do2ds.SetCellTypeComponent((string)"CellTypes", (int)0);
do2ds.SetCellConnectivityComponent((string)"Cells", (int)0);
do2ds.Update();
fd2ad = new vtkFieldDataToAttributeDataFilter();
fd2ad.SetInputData((vtkDataObject)do2ds.GetUnstructuredGridOutput());
fd2ad.SetInputFieldToDataObjectField();
fd2ad.SetOutputAttributeDataToPointData();
fd2ad.SetVectorComponent((int)0, (string)"displacement9", (int)0);
fd2ad.SetVectorComponent((int)1, (string)"displacement9", (int)1);
fd2ad.SetVectorComponent((int)2, (string)"displacement9", (int)2);
fd2ad.SetScalarComponent((int)0, (string)"thickness9", (int)0);
fd2ad.Update();
// Now start visualizing[]
warp = new vtkWarpVector();
warp.SetInputData((vtkDataObject)fd2ad.GetUnstructuredGridOutput());
// extract mold from mesh using connectivity[]
connect = new vtkConnectivityFilter();
connect.SetInputConnection((vtkAlgorithmOutput)warp.GetOutputPort());
connect.SetExtractionModeToSpecifiedRegions();
connect.AddSpecifiedRegion((int)0);
connect.AddSpecifiedRegion((int)1);
moldMapper = new vtkDataSetMapper();
moldMapper.SetInputConnection((vtkAlgorithmOutput)connect.GetOutputPort());
moldMapper.ScalarVisibilityOff();
moldActor = new vtkActor();
moldActor.SetMapper((vtkMapper)moldMapper);
moldActor.GetProperty().SetColor((double).2, (double).2, (double).2);
moldActor.GetProperty().SetRepresentationToWireframe();
// extract parison from mesh using connectivity[]
connect2 = new vtkConnectivityFilter();
connect2.SetInputConnection((vtkAlgorithmOutput)warp.GetOutputPort());
connect2.SetExtractionModeToSpecifiedRegions();
connect2.AddSpecifiedRegion((int)2);
parison = new vtkGeometryFilter();
parison.SetInputConnection((vtkAlgorithmOutput)connect2.GetOutputPort());
normals2 = new vtkPolyDataNormals();
normals2.SetInputConnection((vtkAlgorithmOutput)parison.GetOutputPort());
normals2.SetFeatureAngle((double)60);
lut = new vtkLookupTable();
lut.SetHueRange((double)0.0, (double)0.66667);
parisonMapper = vtkPolyDataMapper.New();
parisonMapper.SetInputConnection((vtkAlgorithmOutput)normals2.GetOutputPort());
parisonMapper.SetLookupTable((vtkScalarsToColors)lut);
parisonMapper.SetScalarRange((double)0.12, (double)1.0);
parisonActor = new vtkActor();
parisonActor.SetMapper((vtkMapper)parisonMapper);
cf = new vtkContourFilter();
cf.SetInputConnection((vtkAlgorithmOutput)connect2.GetOutputPort());
cf.SetValue((int)0, (double).5);
contourMapper = vtkPolyDataMapper.New();
contourMapper.SetInputConnection((vtkAlgorithmOutput)cf.GetOutputPort());
contours = new vtkActor();
contours.SetMapper((vtkMapper)contourMapper);
// Create graphics stuff[]
ren1 = vtkRenderer.New();
renWin = vtkRenderWindow.New();
renWin.AddRenderer((vtkRenderer)ren1);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
// Add the actors to the renderer, set the background and size[]
ren1.AddActor((vtkProp)moldActor);
ren1.AddActor((vtkProp)parisonActor);
ren1.AddActor((vtkProp)contours);
ren1.ResetCamera();
ren1.GetActiveCamera().Azimuth((double)60);
ren1.GetActiveCamera().Roll((double)-90);
ren1.GetActiveCamera().Dolly((double)2);
ren1.ResetCameraClippingRange();
ren1.SetBackground((double)1, (double)1, (double)1);
renWin.SetSize((int)375, (int)200);
iren.Initialize();
}
// prevent the tk window from showing up then start the event loop[]
//deleteAllVTKObjects();
}