/// <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();
}
private void WarpVector()
{
vtkPoints points = vtkPoints.New();
points.InsertNextPoint(0.0, 0.0, 0.0);
points.InsertNextPoint(1.0, 0.0, 0.0);
points.InsertNextPoint(2.0, 0.0, 0.0);
points.InsertNextPoint(3.0, 0.0, 0.0);
points.InsertNextPoint(4.0, 0.0, 0.0);
vtkCellArray lines = vtkCellArray.New();
vtkLine line = vtkLine.New();
line.GetPointIds().SetId(0, 0);
line.GetPointIds().SetId(1, 1);
lines.InsertNextCell(line);
line.GetPointIds().SetId(0, 1);
line.GetPointIds().SetId(1, 2);
lines.InsertNextCell(line);
line.GetPointIds().SetId(0, 2);
line.GetPointIds().SetId(1, 3);
lines.InsertNextCell(line);
line.GetPointIds().SetId(0, 3);
line.GetPointIds().SetId(1, 4);
lines.InsertNextCell(line);
vtkDoubleArray warpData = vtkDoubleArray.New();
warpData.SetNumberOfComponents(3);
warpData.SetName("warpData");
double[] warp = new double[] { 0.0, 0.0, 0.0 };
warp[1] = 0.0;
warpData.InsertNextTuple3(warp[0], warp[1], warp[2]);
warp[1] = 0.1;
warpData.InsertNextTuple3(warp[0], warp[1], warp[2]);
warp[1] = 0.3;
warpData.InsertNextTuple3(warp[0], warp[1], warp[2]);
warp[1] = 0.0;
warpData.InsertNextTuple3(warp[0], warp[1], warp[2]);
warp[1] = 0.1;
warpData.InsertNextTuple3(warp[0], warp[1], warp[2]);
vtkPolyData polydata = vtkPolyData.New();
polydata.SetPoints(points);
polydata.SetLines(lines);
polydata.GetPointData().AddArray(warpData);
polydata.GetPointData().SetActiveVectors(warpData.GetName());
//WarpVector will use the array marked as active vector in polydata
//it has to be a 3 component array
//with the same number of tuples as points in polydata
vtkWarpVector warpVector = vtkWarpVector.New();
#if VTK_MAJOR_VERSION_5
warpVector.SetInput(polydata);
#else
warpVector.SetInputData(polydata);
#endif
warpVector.Update();
vtkPolyDataMapper mapper = vtkPolyDataMapper.New();
#if VTK_MAJOR_VERSION_5
mapper.SetInput(warpVector.GetPolyDataOutput());
#else
mapper.SetInputData(warpVector.GetPolyDataOutput());
#endif
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(.2, .6, .3);
renderer.AddActor(actor);
}
