/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVTestExtractVOI(String [] argv)
{
//Prefix Content is: ""
// to mark the origin[]
sphere = new vtkSphereSource();
sphere.SetRadius((double)2.0);
sphereMapper = vtkPolyDataMapper.New();
sphereMapper.SetInputConnection((vtkAlgorithmOutput)sphere.GetOutputPort());
sphereMapper.ImmediateModeRenderingOn();
sphereActor = new vtkActor();
sphereActor.SetMapper((vtkMapper)sphereMapper);
rt = new vtkRTAnalyticSource();
rt.SetWholeExtent((int)-50,(int)50,(int)-50,(int)50,(int)0,(int)0);
voi = new vtkExtractVOI();
voi.SetInputConnection((vtkAlgorithmOutput)rt.GetOutputPort());
voi.SetVOI((int)-11,(int)39,(int)5,(int)45,(int)0,(int)0);
voi.SetSampleRate((int)5,(int)5,(int)1);
// Get rid ambiguous triagulation issues.[]
surf = new vtkDataSetSurfaceFilter();
surf.SetInputConnection((vtkAlgorithmOutput)voi.GetOutputPort());
tris = new vtkTriangleFilter();
tris.SetInputConnection((vtkAlgorithmOutput)surf.GetOutputPort());
mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection((vtkAlgorithmOutput)tris.GetOutputPort());
mapper.ImmediateModeRenderingOn();
mapper.SetScalarRange((double)130,(double)280);
actor = new vtkActor();
actor.SetMapper((vtkMapper)mapper);
ren = vtkRenderer.New();
ren.AddActor((vtkProp)actor);
ren.AddActor((vtkProp)sphereActor);
ren.ResetCamera();
camera = ren.GetActiveCamera();
//$camera SetPosition 68.1939 -23.4323 12.6465[]
//$camera SetViewUp 0.46563 0.882375 0.0678508 []
//$camera SetFocalPoint 3.65707 11.4552 1.83509 []
//$camera SetClippingRange 59.2626 101.825 []
renWin = vtkRenderWindow.New();
renWin.AddRenderer((vtkRenderer)ren);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
iren.Initialize();
//deleteAllVTKObjects();
}
private void ReadPDB()
{
// 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\caffeine.pdb");
vtkPDBReader pdb = vtkPDBReader.New();
pdb.SetFileName(filePath);
pdb.SetHBScale(1.0);
pdb.SetBScale(1.0);
pdb.Update();
Debug.WriteLine("# of atoms is: " + pdb.GetNumberOfAtoms());
// if molecule contains a lot of atoms, reduce the resolution of the sphere (represents an atom) for faster rendering
int resolution = (int)Math.Floor(Math.Sqrt(300000.0 / pdb.GetNumberOfAtoms())); // 300000.0 is an empriric value
if (resolution > 20)
{
resolution = 20;
}
else if (resolution < 4)
{
resolution = 4;
}
Debug.WriteLine("Resolution is: " + resolution);
vtkSphereSource sphere = vtkSphereSource.New();
sphere.SetCenter(0, 0, 0);
sphere.SetRadius(1);
sphere.SetThetaResolution(resolution);
sphere.SetPhiResolution(resolution);
vtkGlyph3D glyph = vtkGlyph3D.New();
glyph.SetInputConnection(pdb.GetOutputPort());
glyph.SetOrient(1);
glyph.SetColorMode(1);
// glyph.ScalingOn();
glyph.SetScaleMode(2);
glyph.SetScaleFactor(.25);
glyph.SetSourceConnection(sphere.GetOutputPort());
vtkPolyDataMapper atomMapper = vtkPolyDataMapper.New();
atomMapper.SetInputConnection(glyph.GetOutputPort());
atomMapper.ImmediateModeRenderingOn();
atomMapper.UseLookupTableScalarRangeOff();
atomMapper.ScalarVisibilityOn();
atomMapper.SetScalarModeToDefault();
vtkLODActor atom = vtkLODActor.New();
atom.SetMapper(atomMapper);
atom.GetProperty().SetRepresentationToSurface();
atom.GetProperty().SetInterpolationToGouraud();
atom.GetProperty().SetAmbient(0.15);
atom.GetProperty().SetDiffuse(0.85);
atom.GetProperty().SetSpecular(0.1);
atom.GetProperty().SetSpecularPower(30);
atom.GetProperty().SetSpecularColor(1, 1, 1);
atom.SetNumberOfCloudPoints(30000);
vtkTubeFilter tube = vtkTubeFilter.New();
tube.SetInputConnection(pdb.GetOutputPort());
tube.SetNumberOfSides(resolution);
tube.CappingOff();
tube.SetRadius(0.2);
// turn off variation of tube radius with scalar values
tube.SetVaryRadius(0);
tube.SetRadiusFactor(10);
vtkPolyDataMapper bondMapper = vtkPolyDataMapper.New();
bondMapper.SetInputConnection(tube.GetOutputPort());
bondMapper.ImmediateModeRenderingOn();
bondMapper.UseLookupTableScalarRangeOff();
bondMapper.ScalarVisibilityOff();
bondMapper.SetScalarModeToDefault();
vtkLODActor bond = vtkLODActor.New();
bond.SetMapper(bondMapper);
bond.GetProperty().SetRepresentationToSurface();
bond.GetProperty().SetInterpolationToGouraud();
bond.GetProperty().SetAmbient(0.15);
bond.GetProperty().SetDiffuse(0.85);
bond.GetProperty().SetSpecular(0.1);
bond.GetProperty().SetSpecularPower(30);
bond.GetProperty().SetSpecularColor(1, 1, 1);
bond.GetProperty().SetDiffuseColor(1.0000, 0.8941, 0.70981);
// 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(atom);
renderer.AddActor(bond);
}
/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVStreamPolyData(String [] argv)
{
//Prefix Content is: ""
NUMBER_OF_PIECES = 5;
// Generate implicit model of a sphere[]
//[]
// Create renderer stuff[]
//[]
ren1 = vtkRenderer.New();
renWin = vtkRenderWindow.New();
renWin.AddRenderer((vtkRenderer)ren1);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
// create pipeline that handles ghost cells[]
sphere = new vtkSphereSource();
sphere.SetRadius((double)3);
sphere.SetPhiResolution((int)100);
sphere.SetThetaResolution((int)150);
// sphere AddObserver StartEvent {tk_messageBox -message "Executing with piece [[sphere GetOutput] GetUpdatePiece]"}[]
// Just playing with an alternative that is not currently used.[]
//method moved
// Just playing with an alternative that is not currently used.[]
deci = new vtkDecimatePro();
deci.SetInputConnection((vtkAlgorithmOutput)sphere.GetOutputPort());
// this did not remove seams as I thought it would[]
deci.BoundaryVertexDeletionOff();
//deci PreserveTopologyOn[]
// Since quadric Clustering does not handle borders properly yet,[]
// the pieces will have dramatic "eams"[]
q = new vtkQuadricClustering();
q.SetInputConnection((vtkAlgorithmOutput)sphere.GetOutputPort());
q.SetNumberOfXDivisions((int)5);
q.SetNumberOfYDivisions((int)5);
q.SetNumberOfZDivisions((int)10);
q.UseInputPointsOn();
streamer = new vtkPolyDataStreamer();
//streamer SetInputConnection [deci GetOutputPort][]
streamer.SetInputConnection((vtkAlgorithmOutput)q.GetOutputPort());
//streamer SetInputConnection [pdn GetOutputPort][]
streamer.SetNumberOfStreamDivisions((int)NUMBER_OF_PIECES);
mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection((vtkAlgorithmOutput)streamer.GetOutputPort());
mapper.ScalarVisibilityOff();
mapper.SetPiece((int)0);
mapper.SetNumberOfPieces((int)2);
mapper.ImmediateModeRenderingOn();
actor = new vtkActor();
actor.SetMapper((vtkMapper)mapper);
actor.GetProperty().SetColor((double) 0.8300, 0.2400, 0.1000 );
// Add the actors to the renderer, set the background and size[]
//[]
ren1.GetActiveCamera().SetPosition((double)5,(double)5,(double)10);
ren1.GetActiveCamera().SetFocalPoint((double)0,(double)0,(double)0);
ren1.AddActor((vtkProp)actor);
ren1.SetBackground((double)1,(double)1,(double)1);
renWin.SetSize((int)300,(int)300);
iren.Initialize();
// render the image[]
//[]
// 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 AVTestGridSynchronizedTemplates3D(String [] argv)
{
//Prefix Content is: ""
// cut data[]
pl3d = new vtkMultiBlockPLOT3DReader();
pl3d.SetXYZFileName((string)"" + (VTK_DATA_ROOT.ToString()) + "/Data/combxyz.bin");
pl3d.SetQFileName((string)"" + (VTK_DATA_ROOT.ToString()) + "/Data/combq.bin");
pl3d.SetScalarFunctionNumber((int)100);
pl3d.SetVectorFunctionNumber((int)202);
pl3d.Update();
range = ((vtkDataSet)pl3d.GetOutput().GetBlock(0)).GetPointData().GetScalars().GetRange();
min = (double)(lindex(range,0));
max = (double)(lindex(range,1));
value = (min+max)/2.0;
//vtkGridSynchronizedTemplates3D cf[]
cf = new vtkContourFilter();
cf.SetInputData((vtkDataSet)pl3d.GetOutput().GetBlock(0));
cf.SetValue((int)0,(double)value);
//cf ComputeNormalsOff[]
cfMapper = vtkPolyDataMapper.New();
cfMapper.ImmediateModeRenderingOn();
cfMapper.SetInputConnection((vtkAlgorithmOutput)cf.GetOutputPort());
cfMapper.SetScalarRange((double)((vtkDataSet)pl3d.GetOutput().GetBlock(0)).GetPointData().GetScalars().GetRange()[0], (double)((vtkDataSet)pl3d.GetOutput().GetBlock(0)).GetPointData().GetScalars().GetRange()[1]);
cfActor = new vtkActor();
cfActor.SetMapper((vtkMapper)cfMapper);
//outline[]
outline = new vtkStructuredGridOutlineFilter();
outline.SetInputData((vtkDataSet)pl3d.GetOutput().GetBlock(0));
outlineMapper = vtkPolyDataMapper.New();
outlineMapper.SetInputConnection((vtkAlgorithmOutput)outline.GetOutputPort());
outlineActor = new vtkActor();
outlineActor.SetMapper((vtkMapper)outlineMapper);
outlineActor.GetProperty().SetColor((double)0,(double)0,(double)0);
//# Graphics stuff[]
// Create the RenderWindow, Renderer and both Actors[]
//[]
ren1 = vtkRenderer.New();
renWin = vtkRenderWindow.New();
renWin.SetMultiSamples(0);
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)outlineActor);
ren1.AddActor((vtkProp)cfActor);
ren1.SetBackground((double)1,(double)1,(double)1);
renWin.SetSize((int)400,(int)400);
cam1 = ren1.GetActiveCamera();
cam1.SetClippingRange((double)3.95297,(double)50);
cam1.SetFocalPoint((double)9.71821,(double)0.458166,(double)29.3999);
cam1.SetPosition((double)2.7439,(double)-37.3196,(double)38.7167);
cam1.SetViewUp((double)-0.16123,(double)0.264271,(double)0.950876);
iren.Initialize();
// render the image[]
//[]
// loop over surfaces[]
i = 0;
while((i) < 17)
{
cf.SetValue((int)0,(double)min+(i/16.0)*(max-min));
renWin.Render();
i = i + 1;
}
cf.SetValue((int)0,(double)min+(0.2)*(max-min));
renWin.Render();
// prevent the tk window from showing up then start the event loop[]
//deleteAllVTKObjects();
}
