vv/vvSlicerManagerCommand.cxx

   1 /*=========================================================================
   2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
   3
   4   Authors belong to:
   5   - University of LYON              http://www.universite-lyon.fr/
   6   - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
   7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
   8
   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.
  12
  13   It is distributed under dual licence
  14
  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"
  20
  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>
  29 #include <vtkMath.h>
  30 #include <vtkAbstractPropPicker.h>
  31 #include <vtkAssemblyPath.h>
  32 #include <vtkCornerAnnotation.h>
  33 #include <vtkRenderWindow.h>
  34 #include <vvImage.h>
  35 #include <vtkImageReslice.h>
  36 #include "vvSlicer.h"
  37 #include "vvInteractorStyleNavigator.h"
  38 #include "vtkTransform.h"
  39
  40 #include <cmath>
  41
  42 //------------------------------------------------------------------------------
  43 vvSlicerManagerCommand::vvSlicerManagerCommand()
  44 {
  45   mStartSlicer = -1;
  46   mSlicerNumber=-1;
  47 }
  48 //------------------------------------------------------------------------------
  49
  50
  51 //------------------------------------------------------------------------------
  52 //return the num of the current slicer if visible (-1 else)
  53 int vvSlicerManagerCommand::FindSlicerNumber(vtkRenderWindow* renwin)
  54 {
  55   int rvalue;
  56   if (renwin != SM->GetSlicer(mSlicerNumber)->GetRenderWindow() ||
  57       !SM->GetSlicer(mSlicerNumber)->GetRenderer()->GetDraw())
  58     rvalue = -1;
  59   else rvalue = mSlicerNumber;
  60   //std::cerr << this << ":" << mSlicerNumber << ": " << rvalue << endl;
  61   return rvalue;
  62 }
  63 //------------------------------------------------------------------------------
  64
  65
  66 //------------------------------------------------------------------------------
  67 void vvSlicerManagerCommand::Execute(vtkObject *caller,
  68                                      unsigned long event,
  69                                      void *vtkNotUsed(callData))
  70 {
  71   //KeyPress event
  72   vvInteractorStyleNavigator *isi =
  73     dynamic_cast<vvInteractorStyleNavigator *>(caller);
  74   if (isi) {
  75     double x = isi->GetInteractor()->GetEventPosition()[0];
  76     double y = isi->GetInteractor()->GetEventPosition()[1];
  77     double z;
  78
  79     int VisibleInWindow = this->FindSlicerNumber(isi->GetInteractor()->GetRenderWindow());
  80     vtkRenderer* renderer=NULL;
  81     if (VisibleInWindow>-1)
  82       renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
  83     newLandmark = false;
  84
  85
  86
  87     if (event == vtkCommand::StartPickEvent && VisibleInWindow == -1) {
  88       for (int i = 0; i < this->SM->NumberOfSlicers(); i++) {
  89         if (this->SM->GetSlicer(i)->GetCursorVisibility()) {
  90           this->SM->GetSlicer(i)->SetCursorVisibility(0);
  91           this->SM->GetSlicer(i)->Render();
  92         }
  93       }
  94     }
  95     if ( VisibleInWindow > -1 ) {
  96       if (event == vtkCommand::KeyPressEvent) {
  97         std::string KeyPress = isi->GetInteractor()->GetKeySym();
  98         if (KeyPress == "Tab") {
  99           if(isi->GetInteractor()->GetShiftKey())
 100             this->SM->PrevImage(VisibleInWindow);
 101           else
 102             this->SM->NextImage(VisibleInWindow);
 103           return;
 104         }
 105         if (KeyPress == "f" || KeyPress == "F") {
 106           FlyToPosition(isi->GetInteractor(),this->SM->GetSlicer(VisibleInWindow));
 107           return;
 108         }
 109         if (KeyPress == "w") {
 110           this->SM->SetLocalColorWindowing(VisibleInWindow);
 111           return;
 112         }
 113         if (KeyPress == "0") {
 114           this->SM->SetPreset(0);
 115           this->SM->UpdateWindowLevel();
 116           return;
 117         }
 118         if (KeyPress == "1") {
 119           this->SM->SetPreset(1);
 120           this->SM->UpdateWindowLevel();
 121           return;
 122         }
 123         if (KeyPress == "2") {
 124           this->SM->SetPreset(2);
 125           this->SM->UpdateWindowLevel();
 126
 127           return;
 128         }
 129         if (KeyPress == "3") {
 130           this->SM->SetPreset(3);
 131           this->SM->UpdateWindowLevel();
 132           return;
 133         }
 134         if (KeyPress == "4") {
 135           this->SM->SetPreset(4);
 136           this->SM->UpdateWindowLevel();
 137           return;
 138         }
 139         if (KeyPress == "5") {
 140           this->SM->SetPreset(5);
 141           this->SM->UpdateWindowLevel();
 142           return;
 143         }
 144         if (KeyPress == "6") {
 145           this->SM->SetColorMap(0);
 146           this->SM->UpdateWindowLevel();
 147           return;
 148         }
 149         if (KeyPress == "7") {
 150           this->SM->SetColorMap(1);
 151           this->SM->UpdateWindowLevel();
 152           return;
 153         }
 154         if (KeyPress == "8") {
 155           this->SM->SetColorMap(2);
 156           this->SM->UpdateWindowLevel();
 157           return;
 158         }
 159         if (KeyPress == "9") {
 160           this->SM->SetColorMap(3);
 161           this->SM->UpdateWindowLevel();
 162           return;
 163         }
 164         if (KeyPress == "equal") { //keycodes are in vtkWin32RenderWindowInteractor
 165           this->SM->SetPreset(7);
 166           //this->SM->SetColorMap(1);
 167           this->SM->UpdateWindowLevel();
 168           return;
 169         }
 170         if (KeyPress == "minus") {
 171           this->SM->SetColorWindow(-this->SM->GetColorWindow());
 172           this->SM->SetColorMap(-1);
 173           this->SM->UpdateWindowLevel();
 174           return;
 175         }
 176         if (KeyPress == "c") {
 177           this->SM->ToggleContourSuperposition();
 178           this->SM->Render();
 179           return;
 180         }
 181         if (KeyPress == "l") {
 182           this->SM->ToggleInterpolation();
 183           this->SM->Render();
 184           return;
 185         }
 186         if (KeyPress == "h") {
 187           this->SM->SetCursorVisibility(0);
 188           this->SM->Render();
 189           return;
 190         }
 191         if (KeyPress == "u") {
 192           this->SM->Reload();
 193           this->SM->Render();
 194           return;
 195         }
 196         if (KeyPress == "r" || KeyPress=="R") {
 197           this->SM->GetSlicer(VisibleInWindow)->ResetCamera();
 198           this->SM->GetSlicer(VisibleInWindow)->Render();
 199           return;
 200         }
 201         if (KeyPress == "g") {
 202           double* cursorPos = this->SM->GetSlicer(VisibleInWindow)->GetCursorPosition();
 203           this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(
 204             cursorPos[0],cursorPos[1],cursorPos[2],cursorPos[3]);
 205           this->SM->UpdateViews(1,VisibleInWindow);
 206           this->SM->UpdateLinked(VisibleInWindow);
 207           return;
 208         }
 209         if (KeyPress == "F5") {
 210           this->SM->GetSlicer(VisibleInWindow)->FlipHorizontalView();
 211           this->SM->GetSlicer(VisibleInWindow)->Render();
 212           return;
 213         }
 214         if (KeyPress == "F6") {
 215           this->SM->GetSlicer(VisibleInWindow)->FlipVerticalView();
 216           this->SM->GetSlicer(VisibleInWindow)->Render();
 217           return;
 218         }
 219         if (KeyPress == "Up") {
 220           DD("------------ up");
 221           this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()+1);
 222           DD("------------ after set slice");
 223           this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 224           DD("------------ after UpdateInfoOnCursorPosition");
 225           this->SM->UpdateSlice(VisibleInWindow);
 226           DD("------------ after updateslice");
 227         }
 228         if (KeyPress == "Down") {
 229           this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()-1);
 230           this->SM->UpdateSlice(VisibleInWindow);
 231           this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 232         }
 233         if (KeyPress == "space") {
 234           newLandmark = true;
 235         }
 236         if (KeyPress == "Left")
 237           this->SM->SetPreviousTSlice(VisibleInWindow);
 238         if (KeyPress == "Right")
 239           this->SM->SetNextTSlice(VisibleInWindow);
 240
 241         if (KeyPress == "F2") {
 242           this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Sagital\n<slice>");
 243           this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(0);
 244           this->SM->UpdateSliceRange(VisibleInWindow);
 245           this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 246         }
 247         if (KeyPress == "F3") {
 248           this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Coronal\n<slice>");
 249           this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(1);
 250           this->SM->UpdateSliceRange(VisibleInWindow);
 251           this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 252         }
 253         if (KeyPress == "F4") {
 254           this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetText(2,"Axial\n<slice>");
 255           this->SM->GetSlicer(VisibleInWindow)->SetSliceOrientation(2);
 256           this->SM->UpdateSliceRange(VisibleInWindow);
 257           this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 258         }
 259
 260       }
 261
 262       //All type of mouse events
 263       if (event == vtkCommand::LeaveEvent) {
 264         this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
 265             -VTK_DOUBLE_MAX,this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
 266         this->SM->GetSlicer(VisibleInWindow)->Render();
 267         return;
 268       }
 269
 270       //DD(event);
 271       // Mouse release HERE
 272       if (event == vtkCommand::EndPickEvent) {
 273         //            DD(VisibleInWindow);
 274         if (VisibleInWindow > -1)
 275           this->SM->LeftButtonReleaseEvent(VisibleInWindow);
 276         return; // no return !!!! ???
 277       }
 278
 279       if (event == vtkCommand::StartWindowLevelEvent) {
 280         mStartSlicer = -1;
 281         this->InitialWindow = this->SM->GetColorWindow();
 282         this->InitialLevel = this->SM->GetColorLevel();
 283
 284         if (VisibleInWindow > -1) {
 285           mStartSlicer = VisibleInWindow;
 286         }
 287         return;
 288       }
 289
 290       if (event == vtkCommand::EndWindowLevelEvent) {
 291         mStartSlicer = -1;
 292       }
 293
 294     }
 295     if (VisibleInWindow > -1) {
 296       this->SM->Activated();
 297       //if(!this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->GetVisibility())
 298       this->SM->GetSlicer(VisibleInWindow)->GetAnnotation()->SetVisibility(1);
 299
 300       if (event == vtkCommand::MouseWheelForwardEvent && !isi->GetInteractor()->GetControlKey()) {
 301         this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()+1);
 302         this->SM->UpdateSlice(VisibleInWindow);
 303         this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 304       } else if (event == vtkCommand::MouseWheelForwardEvent && isi->GetInteractor()->GetControlKey()) {
 305         double factor = 2;
 306         this->Dolly(pow((double)1.1, factor),isi->GetInteractor());
 307       } else if (event == vtkCommand::MouseWheelBackwardEvent && !isi->GetInteractor()->GetControlKey()) {
 308         this->SM->GetSlicer(VisibleInWindow)->SetSlice(this->SM->GetSlicer(VisibleInWindow)->GetSlice()-1);
 309         this->SM->UpdateSlice(VisibleInWindow);
 310         this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 311       } else if (event == vtkCommand::MouseWheelBackwardEvent && isi->GetInteractor()->GetControlKey()) {
 312         double factor = -2;
 313         this->Dolly(pow((double)1.1, factor),isi->GetInteractor());
 314       }
 315       double xWorld=0;
 316       double yWorld=0;
 317       double zWorld=0;
 318
 319       //Move into World Coordinate
 320       renderer->DisplayToNormalizedDisplay(x,y);
 321       renderer->NormalizedDisplayToViewport(x,y);
 322       renderer->ViewportToNormalizedViewport(x,y);
 323       renderer->NormalizedViewportToView(x,y,z);
 324       renderer->ViewToWorld(x,y,z);
 325       switch (this->SM->GetSlicer(VisibleInWindow)->GetSliceOrientation()) {
 326       case vtkImageViewer2::SLICE_ORIENTATION_XY:
 327         xWorld = x;
 328         yWorld = y;
 329         zWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
 330                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[2] +
 331                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[2];
 332         break;
 333
 334       case vtkImageViewer2::SLICE_ORIENTATION_XZ:
 335         xWorld = x;
 336         yWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
 337                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[1] +
 338                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[1];
 339         zWorld = z;
 340         break;
 341
 342       case vtkImageViewer2::SLICE_ORIENTATION_YZ:
 343         xWorld = this->SM->GetSlicer(VisibleInWindow)->GetSlice()*
 344                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetSpacing()[0] +
 345                  this->SM->GetSlicer(VisibleInWindow)->GetInput()->GetOrigin()[0];
 346         yWorld = y;
 347         zWorld = z;
 348         break;
 349       }
 350       this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(xWorld,yWorld,zWorld,
 351           this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
 352       if (newLandmark) {
 353         this->SM->AddLandmark(xWorld,yWorld,zWorld,
 354                               this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
 355         this->SM->GetSlicer(VisibleInWindow)->UpdateLandmarks();
 356         this->SM->Render();
 357       }
 358       if (event == vtkCommand::PickEvent || event == vtkCommand::StartPickEvent) {
 359         this->SM->UpdateViews(1,VisibleInWindow);
 360         this->SM->UpdateLinked(VisibleInWindow);
 361         this->SM->UpdateInfoOnCursorPosition(VisibleInWindow);
 362       } else {
 363         this->SM->GetSlicer(VisibleInWindow)->Render();
 364       }
 365       //this->SM->GetSlicer(VisibleInWindow)->SetCurrentPosition(-VTK_DOUBLE_MAX,-VTK_DOUBLE_MAX,
 366       //-VTK_DOUBLE_MAX,this->SM->GetSlicer(VisibleInWindow)->GetTSlice());
 367       //this->SM->GetSlicer(VisibleInWindow)->Render();
 368     }
 369
 370     if (event == vtkCommand::WindowLevelEvent && mStartSlicer > -1) {
 371       this->SM->GetSlicer(mStartSlicer)->GetAnnotation()->SetVisibility(1);
 372       // Adjust the window level here
 373       int *size = isi->GetInteractor()->GetRenderWindow()->GetSize();
 374       double window = this->InitialWindow;
 375       double level = this->InitialLevel;
 376       double range[2];
 377       this->SM->GetImage()->GetScalarRange(range);
 378
 379       // Compute normalized delta
 380       double dx = static_cast<double>(isi->GetWindowLevelCurrentPosition()[0] -
 381                                       isi->GetWindowLevelStartPosition()[0]) / size[0];
 382       double dy = static_cast<double>(isi->GetWindowLevelStartPosition()[1] -
 383                                       isi->GetWindowLevelCurrentPosition()[1]) / size[1];
 384       //Window is exponential in nature, use exponential to avoid falling into negative numbers
 385       dx = std::exp(1.0 * (dx*fabs(dx) + dx)) ; //Quadratic behavior for more reactive interface
 386       dy = 0.15 * (dy*fabs(dy) + dy) * (range[1]-range[0]);//Quadratic behavior for more reactive interface
 387
 388       this->SM->SetColorWindow(window*dx);
 389       this->SM->SetColorLevel(level-dy);
 390       this->SM->SetPreset(6);
 391       this->SM->UpdateWindowLevel();
 392       this->SM->Render();
 393       return;
 394     }
 395   }
 396 }
 397 //------------------------------------------------------------------------------
 398
 399
 400 //------------------------------------------------------------------------------
 401 void vvSlicerManagerCommand::Dolly(double factor, vtkRenderWindowInteractor *interactor)
 402 {
 403   int VisibleInWindow = this->FindSlicerNumber(interactor->GetRenderWindow());
 404   vtkRenderer* renderer;
 405   if (VisibleInWindow>-1)
 406     renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
 407   else {
 408     return;
 409   }
 410
 411   double viewFocus[4],viewPoint[4],motionVector[3], focalDepth;
 412   double oldPos[3], newPos[3], distance[2];
 413   vtkCamera *camera = renderer->GetActiveCamera();
 414   camera->GetFocalPoint(viewFocus);
 415
 416   renderer->SetWorldPoint(viewFocus[0], viewFocus[0], viewFocus[0], 1.0);
 417   renderer->WorldToDisplay();
 418   renderer->GetDisplayPoint(viewFocus);
 419
 420   focalDepth = viewFocus[2];
 421
 422   oldPos[0] = renderer->GetCenter()[0];
 423   oldPos[1] = renderer->GetCenter()[1];
 424   oldPos[2] = focalDepth;
 425
 426   distance[0] = 1/factor*
 427                 (interactor->GetEventPosition()[0]-renderer->GetCenter()[0]);
 428   distance[1] = 1/factor*
 429                 (interactor->GetEventPosition()[1]-renderer->GetCenter()[1]);
 430
 431   newPos[0] = interactor->GetEventPosition()[0] - distance[0];
 432   newPos[1] = interactor->GetEventPosition()[1] - distance[1];
 433   newPos[2] = focalDepth;
 434
 435   renderer->DisplayToNormalizedDisplay(oldPos[0],oldPos[1]);
 436   renderer->NormalizedDisplayToViewport(oldPos[0],oldPos[1]);
 437   renderer->ViewportToNormalizedViewport(oldPos[0],oldPos[1]);
 438   renderer->NormalizedViewportToView(oldPos[0],oldPos[1],oldPos[2]);
 439   renderer->ViewToWorld(oldPos[0],oldPos[1],oldPos[2]);
 440
 441   renderer->DisplayToNormalizedDisplay(newPos[0],newPos[1]);
 442   renderer->NormalizedDisplayToViewport(newPos[0],newPos[1]);
 443   renderer->ViewportToNormalizedViewport(newPos[0],newPos[1]);
 444   renderer->NormalizedViewportToView(newPos[0],newPos[1],newPos[2]);
 445   renderer->ViewToWorld(newPos[0],newPos[1],newPos[2]);
 446
 447   motionVector[0] = newPos[0] - oldPos[0];
 448   motionVector[1] = newPos[1] - oldPos[1];
 449   motionVector[2] = newPos[2] - oldPos[2];
 450
 451   camera->GetFocalPoint(viewFocus);
 452   camera->GetPosition(viewPoint);
 453   camera->SetFocalPoint(motionVector[0] + viewFocus[0],
 454                         motionVector[1] + viewFocus[1],
 455                         motionVector[2] + viewFocus[2]);
 456
 457   camera->SetPosition(motionVector[0] + viewPoint[0],
 458                       motionVector[1] + viewPoint[1],
 459                       motionVector[2] + viewPoint[2]);
 460
 461   if (camera->GetParallelProjection()) {
 462     camera->SetParallelScale(camera->GetParallelScale() / factor);
 463   } else {
 464     camera->Dolly(factor);
 465   }
 466
 467   if (interactor->GetLightFollowCamera()) {
 468     renderer->UpdateLightsGeometryToFollowCamera();
 469   }
 470   renderer->ResetCameraClippingRange();
 471   //interactor->Render();
 472 }
 473 //------------------------------------------------------------------------------
 474
 475
 476 //------------------------------------------------------------------------------
 477 void vvSlicerManagerCommand::FlyToPosition(vtkRenderWindowInteractor *interactor,vvSlicer* slicer)
 478 {
 479   double flyFrom[3], flyTo[3];
 480   double d[3], focalPt[3], position[3], positionFrom[3];
 481   int i, j;
 482   int VisibleInWindow = this->FindSlicerNumber(interactor->GetRenderWindow());
 483   vtkRenderer* renderer=NULL;
 484   if (VisibleInWindow>-1)
 485     renderer=this->SM->GetSlicer(VisibleInWindow)->GetRenderer();
 486   else
 487     return;
 488
 489   interactor->GetPicker()->Pick(interactor->GetEventPosition()[0],
 490                                 interactor->GetEventPosition()[1], 0.0,
 491                                 renderer);
 492
 493   vtkAssemblyPath *path=NULL;
 494   vtkAbstractPropPicker *picker;
 495   if ( (picker=vtkAbstractPropPicker::SafeDownCast(interactor->GetPicker()))) {
 496     path = picker->GetPath();
 497   }
 498   if ( path != NULL ) {
 499     flyTo[0] = picker->GetPickPosition()[0];
 500     flyTo[1] = picker->GetPickPosition()[1];
 501     flyTo[2] = picker->GetPickPosition()[2];
 502     renderer->GetActiveCamera()->GetFocalPoint(flyFrom);
 503     renderer->GetActiveCamera()->GetPosition(positionFrom);
 504
 505     switch (slicer->GetSliceOrientation()) {
 506     case vtkImageViewer2::SLICE_ORIENTATION_XY:
 507       flyTo[2] = flyFrom[2];
 508       break;
 509
 510     case vtkImageViewer2::SLICE_ORIENTATION_XZ:
 511       flyTo[1] = flyFrom[1];
 512       break;
 513
 514     case vtkImageViewer2::SLICE_ORIENTATION_YZ:
 515       flyTo[0] = flyFrom[0];
 516       break;
 517     }
 518
 519
 520     for (i=0; i<3; i++) {
 521       d[i] = flyTo[i] - flyFrom[i];
 522     }
 523     double distance = vtkMath::Normalize(d);
 524     double delta = distance/15;
 525
 526     for (i=1; i<=15; i++) {
 527       for (j=0; j<3; j++) {
 528         focalPt[j] = flyFrom[j] + d[j]*i*delta;
 529         position[j] = positionFrom[j] + d[j]*i*delta;
 530       }
 531       renderer->GetActiveCamera()->SetFocalPoint(focalPt);
 532       renderer->GetActiveCamera()->SetPosition(position);
 533       renderer->GetActiveCamera()->Dolly(0.3/15 + 1.0);
 534       renderer->ResetCameraClippingRange();
 535       interactor->Render();
 536     }
 537   }
 538 }
 539 //------------------------------------------------------------------------------