private void Create3DDRC(cDRC DRCToDraw, cDRC_Region AssociatedRegion, Color Color, double Min, double Max)
{
if (DRCToDraw.ResultFit == null) return;
Position = new cPoint3D(AssociatedRegion.PosXMin + 0.5, AssociatedRegion.PosYMin + 0.2, 0);
this.Colour = Color;
vtkPoints points = vtkPoints.New();
vtkUnsignedCharArray colors = vtkUnsignedCharArray.New();
colors.SetName("Colors");
colors.SetNumberOfComponents(3);
colors.SetNumberOfTuples(AssociatedRegion.NumConcentrations);
for (int i = 0; i < AssociatedRegion.NumConcentrations; i++)
{
if (i >= DRCToDraw.ResultFit.Y_Estimated.Count) continue;
double PosZ = 8 - ((DRCToDraw.ResultFit.GetNormalizedY_Estimated()[i]) * 8);
points.InsertPoint(i, i, 0, PosZ);
colors.InsertTuple3(i / AssociatedRegion.NumConcentrations, i / AssociatedRegion.NumConcentrations, 255, i / AssociatedRegion.NumConcentrations);
}
Spline = vtkParametricSpline.New();
Spline.SetPoints(points);
Spline.ClosedOff();
vtkParametricFunctionSource SplineSource = vtkParametricFunctionSource.New();
SplineSource.SetParametricFunction(Spline);
// SplineSource.GetPolyDataInput(0).GetPointData().AddArray(colors);
vtkLinearExtrusionFilter extrude = vtkLinearExtrusionFilter.New();
extrude.SetInputConnection(SplineSource.GetOutputPort());
//extrude.GetPolyDataInput(0).GetPointData().AddArray(colors);
extrude.SetScaleFactor(AssociatedRegion.NumReplicate - 0.2);
//extrude.SetExtrusionTypeToNormalExtrusion();
extrude.SetExtrusionTypeToVectorExtrusion();
extrude.SetVector(0, 1, 0);
vtk_PolyDataMapper = vtkPolyDataMapper.New();
vtk_PolyDataMapper.SetInputConnection(extrude.GetOutputPort()/*SplineSource.GetOutputPort()*/);
vtk_PolyDataMapper.GetInput().GetPointData().AddArray(colors);
vtk_PolyDataMapper.ScalarVisibilityOn();
vtk_PolyDataMapper.SetScalarModeToUsePointFieldData();
vtk_PolyDataMapper.SelectColorArray("Colors");
CreateVTK3DObject(3);
}
private void ParametricObjectsDemo()
{
// we create a matrix of 4x4 renderer in our renderwindow
// each renderer can be interacted with independently from one another
int rendererSize = 189; // width per renderer
int gridDimensions = 4;
this.Size = new System.Drawing.Size(756, 756);
Random rnd = new Random(2);
List <vtkParametricFunction> parametricObjects = new List <vtkParametricFunction>();
parametricObjects.Add(vtkParametricBoy.New());
parametricObjects.Add(vtkParametricConicSpiral.New());
parametricObjects.Add(vtkParametricCrossCap.New());
parametricObjects.Add(vtkParametricDini.New());
vtkParametricEllipsoid ellipsoid = vtkParametricEllipsoid.New();
ellipsoid.SetXRadius(0.5);
ellipsoid.SetYRadius(2.0);
parametricObjects.Add(ellipsoid);
parametricObjects.Add(vtkParametricEnneper.New());
parametricObjects.Add(vtkParametricFigure8Klein.New());
parametricObjects.Add(vtkParametricKlein.New());
parametricObjects.Add(vtkParametricMobius.New());
vtkParametricRandomHills randomHills = vtkParametricRandomHills.New();
randomHills.AllowRandomGenerationOff();
parametricObjects.Add(randomHills);
parametricObjects.Add(vtkParametricRoman.New());
vtkParametricSuperEllipsoid superEllipsoid = vtkParametricSuperEllipsoid.New();
superEllipsoid.SetN1(.50);
superEllipsoid.SetN2(.1);
parametricObjects.Add(superEllipsoid);
vtkParametricSuperToroid superToroid = vtkParametricSuperToroid.New();
superToroid.SetN1(0.2);
superToroid.SetN2(3.0);
parametricObjects.Add(superToroid);
parametricObjects.Add(vtkParametricTorus.New());
vtkParametricSpline spline = vtkParametricSpline.New();
vtkPoints inputPoints = vtkPoints.New();
vtkMath.RandomSeed(8775070);
for (int p = 0; p < 10; p++)
{
double x = vtkMath.Random(0.0, 1.0);
double y = vtkMath.Random(0.0, 1.0);
double z = vtkMath.Random(0.0, 1.0);
inputPoints.InsertNextPoint(x, y, z);
}
spline.SetPoints(inputPoints);
parametricObjects.Add(spline);
List <vtkParametricFunctionSource> parametricFunctionSources = new List <vtkParametricFunctionSource>();
List <vtkRenderer> renderers = new List <vtkRenderer>();
List <vtkPolyDataMapper> mappers = new List <vtkPolyDataMapper>();
List <vtkActor> actors = new List <vtkActor>();
List <vtkTextMapper> textMappers = new List <vtkTextMapper>();
List <vtkActor2D> textActors = new List <vtkActor2D>();
// Create one text property for all
vtkTextProperty textProperty = vtkTextProperty.New();
textProperty.SetFontSize(12);
textProperty.SetJustificationToCentered();
// Create a source, renderer, mapper, and actor
// for each object
for (int i = 0; i < parametricObjects.Count; i++)
{
parametricFunctionSources.Add(vtkParametricFunctionSource.New());
parametricFunctionSources[i].SetParametricFunction(parametricObjects[i]);
parametricFunctionSources[i].Update();
mappers.Add(vtkPolyDataMapper.New());
mappers[i].SetInputConnection(parametricFunctionSources[i].GetOutputPort());
actors.Add(vtkActor.New());
actors[i].SetMapper(mappers[i]);
textMappers.Add(vtkTextMapper.New());
textMappers[i].SetInput(parametricObjects[i].GetClassName());
textMappers[i].SetTextProperty(textProperty);
textActors.Add(vtkActor2D.New());
textActors[i].SetMapper(textMappers[i]);
textActors[i].SetPosition(rendererSize / 2, 16);
renderers.Add(vtkRenderer.New());
}
// Need a renderer even if there is no actor
for (int i = parametricObjects.Count; i < gridDimensions * gridDimensions; i++)
{
renderers.Add(vtkRenderer.New());
}
vtkRenderWindow renderWindow = renderWindowControl1.RenderWindow;
renderWindow.SetSize(rendererSize * gridDimensions, rendererSize * gridDimensions);
for (int row = 0; row < gridDimensions; row++)
{
for (int col = 0; col < gridDimensions; col++)
{
int index = row * gridDimensions + col;
// (xmin, ymin, xmax, ymax)
double[] viewport = new double[] {
(col) * rendererSize / (double)(gridDimensions * rendererSize),
(gridDimensions - (row + 1)) * rendererSize / (double)(gridDimensions * rendererSize),
(col + 1) * rendererSize / (double)(gridDimensions * rendererSize),
(gridDimensions - row) * rendererSize / (double)(gridDimensions * rendererSize)
};
//Debug.WriteLine(viewport[0] + " " + viewport[1] + " " + viewport[2] + " " + viewport[3]);
renderWindow.AddRenderer(renderers[index]);
IntPtr pViewport = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(double)) * 4);
Marshal.Copy(viewport, 0, pViewport, 4);
renderers[index].SetViewport(pViewport);
Marshal.FreeHGlobal(pViewport);
if (index > parametricObjects.Count - 1)
{
continue;
}
renderers[index].AddActor(actors[index]);
renderers[index].AddActor(textActors[index]);
renderers[index].SetBackground(.2 + rnd.NextDouble() / 8, .3 + rnd.NextDouble() / 8, .4 + rnd.NextDouble() / 8);
renderers[index].ResetCamera();
renderers[index].GetActiveCamera().Azimuth(30);
renderers[index].GetActiveCamera().Elevation(-50);
renderers[index].GetActiveCamera().Pitch(-2);
renderers[index].ResetCameraClippingRange();
}
}
}
private void DrawParametricObjects(VTKParametric type)
{
// Select one of the following (matching the selection above)
vtkParametricFunctionSource parametricFunctionSource = vtkParametricFunctionSource.New();
switch (type)
{
case VTKParametric.Torus:
vtkParametricTorus Torus = vtkParametricTorus.New();
parametricFunctionSource.SetParametricFunction(Torus);
break;
case VTKParametric.Boy:
vtkParametricBoy Boy = vtkParametricBoy.New();
parametricFunctionSource.SetParametricFunction(Boy);
break;
case VTKParametric.ConicSpiral:
vtkParametricConicSpiral ConicSpiral = vtkParametricConicSpiral.New();
parametricFunctionSource.SetParametricFunction(ConicSpiral);
break;
case VTKParametric.CrossCap:
vtkParametricCrossCap CrossCap = vtkParametricCrossCap.New();
parametricFunctionSource.SetParametricFunction(CrossCap);
break;
case VTKParametric.Dini:
vtkParametricDini Dini = vtkParametricDini.New();
parametricFunctionSource.SetParametricFunction(Dini);
break;
case VTKParametric.Ellipsoid:
vtkParametricEllipsoid Ellipsoid = vtkParametricEllipsoid.New();
Ellipsoid.SetXRadius(0.5);
Ellipsoid.SetYRadius(2.0);
parametricFunctionSource.SetParametricFunction(Ellipsoid);
break;
case VTKParametric.Enneper:
vtkParametricEnneper Enneper = vtkParametricEnneper.New();
parametricFunctionSource.SetParametricFunction(Enneper);
break;
case VTKParametric.Figure8Klein:
vtkParametricFigure8Klein Figure8Klein = vtkParametricFigure8Klein.New();
parametricFunctionSource.SetParametricFunction(Figure8Klein);
break;
case VTKParametric.Klein:
vtkParametricKlein Klein = vtkParametricKlein.New();
parametricFunctionSource.SetParametricFunction(Klein);
break;
case VTKParametric.Mobius:
vtkParametricMobius Mobius = vtkParametricMobius.New();
parametricFunctionSource.SetParametricFunction(Mobius);
break;
case VTKParametric.RandomHills:
vtkParametricRandomHills RandomHills = vtkParametricRandomHills.New();
parametricFunctionSource.SetParametricFunction(RandomHills);
break;
case VTKParametric.Roman:
vtkParametricRoman Roman = vtkParametricRoman.New();
parametricFunctionSource.SetParametricFunction(Roman);
break;
case VTKParametric.Spline:
vtkParametricSpline spline = vtkParametricSpline.New();
vtkPoints inputPoints = vtkPoints.New();
vtkMath.RandomSeed(8775070);
for (int p = 0; p < 10; p++)
{
double x = vtkMath.Random(0.0, 1.0);
double y = vtkMath.Random(0.0, 1.0);
double z = vtkMath.Random(0.0, 1.0);
inputPoints.InsertNextPoint(x, y, z);
}
spline.SetPoints(inputPoints);
parametricFunctionSource.SetParametricFunction(spline);
break;
case VTKParametric.SuperEllipsoid:
vtkParametricSuperEllipsoid superEllipsoid = vtkParametricSuperEllipsoid.New();
superEllipsoid.SetN1(.50);
superEllipsoid.SetN2(.1);
parametricFunctionSource.SetParametricFunction(superEllipsoid);
break;
case VTKParametric.SuperToroid:
vtkParametricSuperToroid SuperToroid = vtkParametricSuperToroid.New();
SuperToroid.SetN1(0.5);
SuperToroid.SetN2(0.1);
parametricFunctionSource.SetParametricFunction(SuperToroid);
break;
}
parametricFunctionSource.Update();
// Setup mapper and actor
vtkPolyDataMapper mapper = vtkPolyDataMapper.New();
mapper.SetInputConnection(parametricFunctionSource.GetOutputPort());
vtkActor actor = vtkActor.New();
actor.SetMapper(mapper);
// Visualize
RenderAddActor(actor);
}