3 #include <creaVtkBasicSlicer.h>
5 #include "vtkRenderWindow.h"
6 #include "vtkRenderWindowInteractor.h"
8 #include "vtkOutlineFilter.h"
9 #include "vtkPolyDataMapper.h"
11 #include "vtkImagePlaneWidget.h"
12 #include "vtkCellPicker.h"
13 #include "vtkProperty.h"
14 #include "vtkRenderer.h"
15 #include "vtkCamera.h"
17 #include "vtkInteractorStyleTrackballCamera.h"
22 void VtkBasicSlicer( vtkImageData* I )
25 int xMin, xMax, yMin, yMax, zMin, zMax;
26 I->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
29 I->GetSpacing(sx, sy, sz);
32 I->GetOrigin(ox,oy,oz);
35 // An outline is shown for context.
37 vtkOutlineFilter* outline = vtkOutlineFilter::New();
40 vtkPolyDataMapper* outlineMapper = vtkPolyDataMapper::New();
41 outlineMapper->SetInputConnection(outline->GetOutputPort());
43 vtkActor* outlineActor = vtkActor::New();
44 outlineActor->SetMapper(outlineMapper);
47 // The shared picker enables us to use 3 planes at one time
48 // and gets the picking order right
49 vtkCellPicker* picker = vtkCellPicker::New();
50 picker->SetTolerance(0.005);
52 // The 3 image plane widgets are used to probe the dataset.
53 vtkImagePlaneWidget* planeWidgetX = vtkImagePlaneWidget::New();
54 planeWidgetX->DisplayTextOn();
55 planeWidgetX->SetInput(I);
56 planeWidgetX->SetPlaneOrientationToXAxes();
57 planeWidgetX->SetSlicePosition((xMax+xMin)/2.);
58 planeWidgetX->SetPicker(picker);
59 planeWidgetX->SetKeyPressActivationValue('x');
60 vtkProperty* prop1 = planeWidgetX->GetPlaneProperty();
61 prop1->SetColor(1, 0, 0);
63 vtkImagePlaneWidget* planeWidgetY = vtkImagePlaneWidget::New();
64 planeWidgetY->DisplayTextOn();
65 planeWidgetY->SetInput(I);
66 planeWidgetY->SetPlaneOrientationToYAxes();
67 planeWidgetY->SetSlicePosition((yMax+yMin)/2.);
68 planeWidgetY->SetPicker(picker);
69 planeWidgetY->SetKeyPressActivationValue('y');
70 vtkProperty* prop2 = planeWidgetY->GetPlaneProperty();
71 prop2->SetColor(1, 1, 0);
72 planeWidgetY->SetLookupTable(planeWidgetX->GetLookupTable());
74 // for the z-slice, turn off texture interpolation:
75 // interpolation is now nearest neighbour, to demonstrate
76 // cross-hair cursor snapping to pixel centers
77 vtkImagePlaneWidget* planeWidgetZ = vtkImagePlaneWidget::New();
78 planeWidgetZ->DisplayTextOn();
79 planeWidgetZ->SetInput(I);
80 planeWidgetZ->SetPlaneOrientationToZAxes();
81 planeWidgetZ->SetSlicePosition((zMax+zMin)/2.);
82 planeWidgetZ->SetPicker(picker);
83 planeWidgetZ->SetKeyPressActivationValue('z');
84 vtkProperty* prop3 = planeWidgetZ->GetPlaneProperty();
85 prop3->SetColor(0, 0, 1);
86 planeWidgetZ->SetLookupTable(planeWidgetX->GetLookupTable());
88 planeWidgetZ->SetWindowLevel(512,256);
90 // Create the RenderWindow and Renderer
91 vtkRenderer* ren = vtkRenderer::New();
92 vtkRenderWindow* renWin = vtkRenderWindow::New();
93 renWin->AddRenderer(ren);
95 // Add the outline actor to the renderer, set the background color and size
96 ren->AddActor(outlineActor);
97 renWin->SetSize(600, 600);
98 ren->SetBackground(0.1, 0.1, 0.2);
101 vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
102 iren->SetInteractorStyle( vtkInteractorStyleTrackballCamera::New() );
103 iren->SetRenderWindow(renWin);
105 // Set the interactor for the widgets
106 planeWidgetX->SetInteractor(iren);
108 planeWidgetY->SetInteractor(iren);
110 planeWidgetZ->SetInteractor(iren);
113 // Create an initial interesting view
114 vtkCamera* cam1 = ren->GetActiveCamera();
115 cam1->SetFocalPoint(0, 0, 0);
116 cam1->SetPosition(0, 0, -1500);
117 //cam1->OrthogonalizeViewUp();
118 cam1->Elevation(110);
119 cam1->SetViewUp(0, 0, -1);
121 ren->ResetCameraClippingRange();
124 // render_widget->Render();
133 planeWidgetX->Delete();
134 planeWidgetY->Delete();
135 planeWidgetZ->Delete();
137 outlineMapper->Delete();
138 outlineActor->Delete();