1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
5 - University of LYON http://www.universite-lyon.fr/
6 - Léon Bérard cancer center http://www.centreleonberard.fr
7 - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
9 This software is distributed WITHOUT ANY WARRANTY; without even
10 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
11 PURPOSE. See the copyright notices for more information.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ===========================================================================**/
18 #include "vvSlicerManagerCommand.h"
19 #include "vvSlicerManager.h"
21 #include <vtkTextProperty.h>
22 #include <vtkRenderer.h>
23 #include <vtkImageActor.h>
24 #include <vtkRenderWindowInteractor.h>
25 #include <vtkPropPicker.h>
26 #include <vtkCamera.h>
27 #include <vtkImageMapToWindowLevelColors.h>
28 #include <vtkLookupTable.h>
30 #include <vtkAbstractPropPicker.h>
31 #include <vtkAssemblyPath.h>
32 #include <vtkCornerAnnotation.h>
33 #include <vtkRenderWindow.h>
35 #include <vtkImageReslice.h>
37 #include "vvInteractorStyleNavigator.h"
38 #include "vtkTransform.h"
42 //------------------------------------------------------------------------------
43 vvSlicerManagerCommand::vvSlicerManagerCommand()
48 //------------------------------------------------------------------------------
51 //------------------------------------------------------------------------------
52 //return the num of the current slicer if visible (-1 else)
53 int vvSlicerManagerCommand::FindSlicerNumber(vtkRenderWindow* renwin)
56 if (renwin != SM->GetSlicer(mSlicerNumber)->GetRenderWindow() ||
57 !SM->GetSlicer(mSlicerNumber)->GetRenderer()->GetDraw())
59 else rvalue = mSlicerNumber;
60 //std::cerr << this << ":" << mSlicerNumber << ": " << rvalue << endl;
63 //------------------------------------------------------------------------------
66 //------------------------------------------------------------------------------
67 void vvSlicerManagerCommand::Execute(vtkObject *caller,
69 void *vtkNotUsed(callData))
72 vvInteractorStyleNavigator *isi =
73 dynamic_cast<vvInteractorStyleNavigator *>(caller);
75 double x = isi->GetInteractor()->GetEventPosition()[0];
76 double y = isi->GetInteractor()->GetEventPosition()[1];
79 int VisibleInWindow = this->FindSlicerNumber(isi->GetInteractor()->GetRenderWindow());
80 vtkRenderer* renderer=NULL;
81 if (VisibleInWindow>-1)
82 renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
87 if (event == vtkCommand::StartPickEvent && VisibleInWindow == -1) {
88 for (int i = 0; i < this->SM->GetNumberOfSlicers(); i++) {
89 if (this->SM->GetSlicer(i)->GetCursorVisibility() && !this->SM->IsLinked()) {
90 this->SM->GetSlicer(i)->SetCursorVisibility(0);
91 this->SM->GetSlicer(i)->SetCornerAnnotationVisibility(0);
92 this->SM->GetSlicer(i)->Render();
96 if ( VisibleInWindow > -1 ) {
97 if (event == vtkCommand::KeyPressEvent) {
98 std::string KeyPress = isi->GetInteractor()->GetKeySym();
99 bool bCtrlKey = isi->GetInteractor()->GetControlKey();
100 if (KeyPress == "Tab") {
101 if(isi->GetInteractor()->GetShiftKey())
102 this->SM->PrevImage(VisibleInWindow);
104 this->SM->NextImage(VisibleInWindow);
107 if (KeyPress == "f" || KeyPress == "F") {
108 FlyToPosition(isi->GetInteractor(),this->SM->GetSlicer(VisibleInWindow));
111 if (KeyPress == "w") {
112 this->SM->SetLocalColorWindowing(VisibleInWindow);
115 if (KeyPress == "0") {
116 this->SM->SetPreset(0);
117 this->SM->UpdateWindowLevel();
120 if (KeyPress == "1") {
121 this->SM->SetPreset(1);
122 this->SM->UpdateWindowLevel();
125 if (KeyPress == "2") {
126 this->SM->SetPreset(2);
127 this->SM->UpdateWindowLevel();
131 if (KeyPress == "3") {
132 this->SM->SetPreset(3);
133 this->SM->UpdateWindowLevel();
136 if (KeyPress == "4") {
137 this->SM->SetPreset(4);
138 this->SM->UpdateWindowLevel();
141 if (KeyPress == "5") {
142 this->SM->SetPreset(5);
143 this->SM->UpdateWindowLevel();
146 if (KeyPress == "6") {
147 this->SM->SetColorMap(0);
148 this->SM->UpdateWindowLevel();
151 if (KeyPress == "7") {
152 this->SM->SetColorMap(1);
153 this->SM->UpdateWindowLevel();
156 if (KeyPress == "8") {
157 this->SM->SetColorMap(2);
158 this->SM->UpdateWindowLevel();
161 if (KeyPress == "9") {
162 this->SM->SetColorMap(3);
163 this->SM->UpdateWindowLevel();
166 if (KeyPress == "equal") { //keycodes are in vtkWin32RenderWindowInteractor
167 this->SM->SetPreset(7);
168 //this->SM->SetColorMap(1);
169 this->SM->UpdateWindowLevel();
172 if (KeyPress == "minus") {
173 this->SM->SetColorWindow(-this->SM->GetColorWindow());
174 this->SM->SetColorMap(-1);
175 this->SM->UpdateWindowLevel();
178 if (KeyPress == "c") {
179 this->SM->ToggleContourSuperposition();
183 if (KeyPress == "l") {
184 this->SM->ToggleInterpolation();
188 if (KeyPress == "h") {
189 std::cout << "KeyPress == \"h\"\n";
190 this->SM->SetCursorAndCornerAnnotationVisibility(0);
194 if (KeyPress == "x") {
195 for(int i=0; i<SM->GetNumberOfSlicers(); i++) {
196 SM->RemoveActor("overlay",0);
202 if (KeyPress == "u") {
207 if (KeyPress == "r" || KeyPress=="R") {
208 this->SM->GetSlicer(VisibleInWindow)->ResetCamera();
209 this->SM->GetSlicer(VisibleInWindow)->Render();
210 this->SM->UpdateLinkedNavigation( this->SM->GetSlicer(VisibleInWindow) );
213 if (KeyPress == "g") {
215 this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(0,0,0,0);
217 double* cursorPos = this->SM->GetSlicer(VisibleInWindow)->GetCursorPosition();
218 this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(
219 cursorPos[0],cursorPos[1],cursorPos[2],cursorPos[3]);
221 this->SM->UpdateViews(1,VisibleInWindow);
222 this->SM->UpdateLinked(VisibleInWindow);
225 if (KeyPress == "F5") {
226 this->SM->GetSlicer(VisibleInWindow)->FlipHorizontalView();
227 this->SM->GetSlicer(VisibleInWindow)->Render();
230 if (KeyPress == "F6") {
231 this->SM->GetSlicer(VisibleInWindow)->FlipVerticalView();
232 this->SM->GetSlicer(VisibleInWindow)->Render();
235 if (KeyPress == "Up") {
236 this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()+1);
237 this->SM->UpdateSlice(VisibleInWindow);
239 if (KeyPress == "Down") {
240 this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()-1);
241 this->SM->UpdateSlice(VisibleInWindow);
243 if (KeyPress == "space") {
246 if (KeyPress == "Left")
247 this->SM->SetPreviousTSlice(VisibleInWindow);
248 if (KeyPress == "Right")
249 this->SM->SetNextTSlice(VisibleInWindow);
251 if (KeyPress == "F2") {
252 this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Sagital\n<slice>");
253 //this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(0);
254 this->SM->SetSliceOrientation(VisibleInWindow, 0);
255 this->SM->UpdateSliceRange(VisibleInWindow);
257 if (KeyPress == "F3") {
258 this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Coronal\n<slice>");
259 //this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(1);
260 this->SM->SetSliceOrientation(VisibleInWindow, 1);
261 this->SM->UpdateSliceRange(VisibleInWindow);
263 if (KeyPress == "F4") {
264 this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Axial\n<slice>");
265 //this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(2);
266 this->SM->SetSliceOrientation(VisibleInWindow, 2);
267 this->SM->UpdateSliceRange(VisibleInWindow);
272 //All type of mouse events
273 if (event == vtkCommand::LeaveEvent) {
274 this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
275 -VTK_DOUBLE_MAX,this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
276 this->SM->GetSlicer(VisibleInWindow)->Render();
281 // Mouse release HERE
282 if (event == vtkCommand::EndPickEvent) {
283 // DD(VisibleInWindow);
284 if (VisibleInWindow > -1)
285 this->SM->LeftButtonReleaseEvent(VisibleInWindow);
286 return; // no return !!!! ???
289 if (event == vtkCommand::StartWindowLevelEvent) {
291 this->InitialWindow = this->SM->GetColorWindow();
292 this->InitialLevel = this->SM->GetColorLevel();
294 if (VisibleInWindow > -1) {
295 mStartSlicer = VisibleInWindow;
300 if (event == vtkCommand::EndWindowLevelEvent) {
304 if (event == vtkCommand::EndInteractionEvent) {
306 this->SM->UpdateLinkedNavigation(this->SM->GetSlicer(VisibleInWindow),true);
310 if (VisibleInWindow > -1) {
311 this->SM->Activated();
312 this->SM->GetSlicer(VisibleInWindow)->SetCornerAnnotationVisibility(1);
314 if (event == vtkCommand::MouseWheelForwardEvent && !isi->GetInteractor()->GetControlKey()) {
316 this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()+1);
317 this->SM->UpdateSlice(VisibleInWindow);
318 } else if (event == vtkCommand::MouseWheelForwardEvent && isi->GetInteractor()->GetControlKey()) {
320 this->Dolly(pow((double)1.1, factor),isi->GetInteractor());
321 Execute(caller, vtkCommand::EndInteractionEvent, NULL);
322 } else if (event == vtkCommand::MouseWheelBackwardEvent && !isi->GetInteractor()->GetControlKey()) {
324 this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()-1);
325 this->SM->UpdateSlice(VisibleInWindow);
326 } else if (event == vtkCommand::MouseWheelBackwardEvent && isi->GetInteractor()->GetControlKey()) {
328 this->Dolly(pow((double)1.1, factor),isi->GetInteractor());
329 Execute(caller, vtkCommand::EndInteractionEvent, NULL);
335 //Move into World Coordinate
336 renderer->DisplayToNormalizedDisplay(x,y);
337 renderer->NormalizedDisplayToViewport(x,y);
338 renderer->ViewportToNormalizedViewport(x,y);
339 renderer->NormalizedViewportToView(x,y,z);
340 renderer->ViewToWorld(x,y,z);
341 switch (this->SM->GetSlicer(VisibleInWindow)->GetSliceOrientation()) {
342 case vtkImageViewer2::SLICE_ORIENTATION_XY:
345 zWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
346 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[2] +
347 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[2];
350 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
352 yWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
353 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[1] +
354 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[1];
358 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
359 xWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
360 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[0] +
361 this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[0];
366 this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(xWorld,yWorld,zWorld,
367 this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
369 this->SM->AddLandmark(xWorld,yWorld,zWorld,
370 this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
371 this->SM->GetSlicer(VisibleInWindow)->UpdateLandmarks();
374 if (event == vtkCommand::PickEvent || event == vtkCommand::StartPickEvent) {
376 this->SM->UpdateViews(1,VisibleInWindow);
377 this->SM->UpdateLinked(VisibleInWindow);
378 this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
380 this->SM->GetSlicer(VisibleInWindow)->Render();
382 //this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
383 //-VTK_DOUBLE_MAX,this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
384 //this->SM->GetSlicer(VisibleInWindow)->Render();
387 if (event == vtkCommand::WindowLevelEvent && mStartSlicer > -1) {
388 this->SM->GetSlicer(mStartSlicer)->SetCornerAnnotationVisibility(1);
389 // Adjust the window level here
390 int *size = isi->GetInteractor()->GetRenderWindow()->GetSize();
391 double window = this->InitialWindow;
392 double level = this->InitialLevel;
394 this->SM->GetImage()->GetScalarRange(range);
396 // Compute normalized delta
397 double dx = static_cast<double>(isi->GetWindowLevelCurrentPosition()[0] -
398 isi->GetWindowLevelStartPosition()[0]) / size[0];
399 double dy = static_cast<double>(isi->GetWindowLevelStartPosition()[1] -
400 isi->GetWindowLevelCurrentPosition()[1]) / size[1];
401 //Window is exponential in nature, use exponential to avoid falling into negative numbers
402 dx = std::exp(1.0 * (dx*fabs(dx) + dx)) ; //Quadratic behavior for more reactive interface
403 dy = 0.15 * (dy*fabs(dy) + dy) * (range[1]-range[0]);//Quadratic behavior for more reactive interface
405 this->SM->SetColorWindow(window*dx);
406 this->SM->SetColorLevel(level-dy);
407 this->SM->SetPreset(6);
408 this->SM->UpdateWindowLevel();
414 //------------------------------------------------------------------------------
417 //------------------------------------------------------------------------------
418 void vvSlicerManagerCommand::Dolly(double factor, vtkRenderWindowInteractor *interactor)
420 int VisibleInWindow = this->FindSlicerNumber(interactor->GetRenderWindow());
421 vtkRenderer* renderer;
422 if (VisibleInWindow>-1)
423 renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
428 double viewFocus[4],viewPoint[4],motionVector[3], focalDepth;
429 double oldPos[3], newPos[3], distance[2];
430 vtkCamera *camera = renderer->GetActiveCamera();
431 camera->GetFocalPoint(viewFocus);
433 renderer->SetWorldPoint(viewFocus[0], viewFocus[0], viewFocus[0], 1.0);
434 renderer->WorldToDisplay();
435 renderer->GetDisplayPoint(viewFocus);
437 focalDepth = viewFocus[2];
439 oldPos[0] = renderer->GetCenter()[0];
440 oldPos[1] = renderer->GetCenter()[1];
441 oldPos[2] = focalDepth;
443 distance[0] = 1/factor*
444 (interactor->GetEventPosition()[0]-renderer->GetCenter()[0]);
445 distance[1] = 1/factor*
446 (interactor->GetEventPosition()[1]-renderer->GetCenter()[1]);
448 newPos[0] = interactor->GetEventPosition()[0] - distance[0];
449 newPos[1] = interactor->GetEventPosition()[1] - distance[1];
450 newPos[2] = focalDepth;
452 renderer->DisplayToNormalizedDisplay(oldPos[0],oldPos[1]);
453 renderer->NormalizedDisplayToViewport(oldPos[0],oldPos[1]);
454 renderer->ViewportToNormalizedViewport(oldPos[0],oldPos[1]);
455 renderer->NormalizedViewportToView(oldPos[0],oldPos[1],oldPos[2]);
456 renderer->ViewToWorld(oldPos[0],oldPos[1],oldPos[2]);
458 renderer->DisplayToNormalizedDisplay(newPos[0],newPos[1]);
459 renderer->NormalizedDisplayToViewport(newPos[0],newPos[1]);
460 renderer->ViewportToNormalizedViewport(newPos[0],newPos[1]);
461 renderer->NormalizedViewportToView(newPos[0],newPos[1],newPos[2]);
462 renderer->ViewToWorld(newPos[0],newPos[1],newPos[2]);
464 motionVector[0] = newPos[0] - oldPos[0];
465 motionVector[1] = newPos[1] - oldPos[1];
466 motionVector[2] = newPos[2] - oldPos[2];
468 camera->GetFocalPoint(viewFocus);
469 camera->GetPosition(viewPoint);
470 camera->SetFocalPoint(motionVector[0] + viewFocus[0],
471 motionVector[1] + viewFocus[1],
472 motionVector[2] + viewFocus[2]);
474 camera->SetPosition(motionVector[0] + viewPoint[0],
475 motionVector[1] + viewPoint[1],
476 motionVector[2] + viewPoint[2]);
478 if (camera->GetParallelProjection()) {
479 camera->SetParallelScale(camera->GetParallelScale() / factor);
481 camera->Dolly(factor);
484 if (interactor->GetLightFollowCamera()) {
485 renderer->UpdateLightsGeometryToFollowCamera();
487 renderer->ResetCameraClippingRange();
488 //interactor->Render();
490 //------------------------------------------------------------------------------
493 //------------------------------------------------------------------------------
494 void vvSlicerManagerCommand::FlyToPosition(vtkRenderWindowInteractor *interactor,vvSlicer* slicer)
496 double flyFrom[3], flyTo[3];
497 double d[3], focalPt[3], position[3], positionFrom[3];
499 int VisibleInWindow = this->FindSlicerNumber(interactor->GetRenderWindow());
500 vtkRenderer* renderer=NULL;
501 if (VisibleInWindow>-1)
502 renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
506 interactor->GetPicker()->Pick(interactor->GetEventPosition()[0],
507 interactor->GetEventPosition()[1], 0.0,
510 vtkAssemblyPath *path=NULL;
511 vtkAbstractPropPicker *picker;
512 if ( (picker=vtkAbstractPropPicker::SafeDownCast(interactor->GetPicker()))) {
513 path = picker->GetPath();
515 if ( path != NULL ) {
516 flyTo[0] = picker->GetPickPosition()[0];
517 flyTo[1] = picker->GetPickPosition()[1];
518 flyTo[2] = picker->GetPickPosition()[2];
519 renderer->GetActiveCamera()->GetFocalPoint(flyFrom);
520 renderer->GetActiveCamera()->GetPosition(positionFrom);
522 switch (slicer->GetSliceOrientation()) {
523 case vtkImageViewer2::SLICE_ORIENTATION_XY:
524 flyTo[2] = flyFrom[2];
527 case vtkImageViewer2::SLICE_ORIENTATION_XZ:
528 flyTo[1] = flyFrom[1];
531 case vtkImageViewer2::SLICE_ORIENTATION_YZ:
532 flyTo[0] = flyFrom[0];
537 for (i=0; i<3; i++) {
538 d[i] = flyTo[i] - flyFrom[i];
540 double distance = vtkMath::Normalize(d);
541 double delta = distance/15;
543 for (i=1; i<=15; i++) {
544 for (j=0; j<3; j++) {
545 focalPt[j] = flyFrom[j] + d[j]*i*delta;
546 position[j] = positionFrom[j] + d[j]*i*delta;
548 renderer->GetActiveCamera()->SetFocalPoint(focalPt);
549 renderer->GetActiveCamera()->SetPosition(position);
550 renderer->GetActiveCamera()->Dolly(0.3/15 + 1.0);
551 renderer->ResetCameraClippingRange();
552 interactor->Render();
556 //------------------------------------------------------------------------------