lib/maracasVisuLib/src/interface/wxWindows/widgets/wxVtkMPR3DView.cxx

   1 /*# ---------------------------------------------------------------------
   2 #
   3 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
   4 #                        pour la Sant�)
   5 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
   6 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
   7 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
   8 #
   9 #  This software is governed by the CeCILL-B license under French law and
  10 #  abiding by the rules of distribution of free software. You can  use,
  11 #  modify and/ or redistribute the software under the terms of the CeCILL-B
  12 #  license as circulated by CEA, CNRS and INRIA at the following URL
  13 #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  14 #  or in the file LICENSE.txt.
  15 #
  16 #  As a counterpart to the access to the source code and  rights to copy,
  17 #  modify and redistribute granted by the license, users are provided only
  18 #  with a limited warranty  and the software's author,  the holder of the
  19 #  economic rights,  and the successive licensors  have only  limited
  20 #  liability.
  21 #
  22 #  The fact that you are presently reading this means that you have had
  23 #  knowledge of the CeCILL-B license and that you accept its terms.
  24 # ------------------------------------------------------------------------ */
  25
  26 #include "wxVtkMPR3DView.h"
  27
  28 #include "vtkDataSetMapper.h"
  29 #include "vtkPlaneSource.h"
  30 #include "vtkProbeFilter.h"
  31 #include "vtkStructuredPoints.h"
  32 #include "vtkMetaImageWriter.h"
  33 #include "vtkImageChangeInformation.h"
  34 #include "vtkPointData.h"
  35 #include "vtkWindowLevelLookupTable.h"
  36 #include "vtkPlane.h"
  37 #include "vtkCutter.h"
  38 #include "vtkTransform.h"
  39 #include "vtkmyPWCallback_3DPointWidget.h"
  40 #include "vtkProperty.h"
  41 #include "vtkCellPicker.h"
  42 #include "vtkInteractorStyleBaseView.h"
  43 #include <vtkImageReslice.h>
  44
  45
  46
  47 //----------------------------------------------------------------------------
  48 //----------------------------------------------------------------------------
  49 //----------------------------------------------------------------------------
  50 class vtkWindowLevelImagePlaneWidgetCallback : public vtkCommand
  51 {
  52          public:
  53            static vtkWindowLevelImagePlaneWidgetCallback *New()
  54            {
  55                    return new vtkWindowLevelImagePlaneWidgetCallback;
  56            }
  57
  58            void Execute( vtkObject *caller, unsigned long ev, void *callData )
  59                 {
  60                         bool needRefresh=false;
  61                         double slicePosition;
  62                         vtkImagePlaneWidget* ipw =dynamic_cast< vtkImagePlaneWidget* >( caller );
  63                     if ( ev == vtkCommand::WindowLevelEvent)
  64                         {
  65                                 double info[2];
  66                                 ipw->GetWindowLevel(info);
  67                                 _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorWindow(info[0]);
  68                                 _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorLevel(info[1]);
  69                                 needRefresh=true;
  70                         } // if ev
  71
  72                    if ( ev == vtkCommand::InteractionEvent)
  73                    {
  74                            needRefresh=true;
  75                            if (ipw->GetCursorDataStatus())
  76                            {
  77                                    double xyzv[4];
  78                                    ipw->GetCursorData(xyzv);
  79                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX(xyzv[0]);
  80                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY(xyzv[1]);
  81                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ(xyzv[2]);
  82                            } else {
  83                                    vtkImageData *image = _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
  84                                    double spc[3];
  85                                    if(image)
  86                                    {
  87                                            image->GetSpacing(spc);
  88                                    }
  89                                    double normal[3];
  90                                    slicePosition = ipw->GetSlicePosition();
  91                                    ipw->GetNormal(normal);
  92                                    // FreePlaneX
  93                                    if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
  94                                    {
  95                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX( slicePosition/spc[0] );
  96                                    }
  97                                    // FreePlaneY
  98                                    if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
  99                                    {
 100                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY( slicePosition/spc[1] );
 101                                    }
 102                                    // FreePlaneZ
 103                                    if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
 104                                    {
 105                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ( slicePosition/spc[2] );
 106                                    }
 107                            }// ipw->GetCursorDataStatus
 108                    }// ev
 109                    if (needRefresh)
 110                    {
 111                            vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)_wxvtkmpr3Dview->GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
 112                            isbv->SetParent_refresh_waiting();
 113                            isbv->EvaluateToRefresh();
 114                    }
 115            } // Execute
 116                 vtkWindowLevelImagePlaneWidgetCallback(){}
 117                 wxVtkMPR3DView *_wxvtkmpr3Dview;
 118 };
 119
 120 //-------------------------------------------------------------------
 121 //-------------------------------------------------------------------
 122 //-------------------------------------------------------------------
 123 wxVtkMPR3DView::wxVtkMPR3DView( wxVtk3DBaseView *wxvtk3Dbaseview )
 124 {
 125         _wxvtk3Dbaseview                        =       wxvtk3Dbaseview;
 126         _vtkmpr3Ddataviewer                     =       NULL;
 127         _wxvtkmpr3DviewCntrlPanel       =       NULL;
 128         //_myCallback                           =       NULL;
 129         _pointWidget                            =       NULL;
 130         _planeWidget                            =       NULL;
 131         _vtkplane                                       =       NULL;
 132         _probe                                          =       NULL;
 133         _contourMapper                          =       NULL;
 134         _planeWidgetX                           =       NULL;
 135         _planeWidgetY                           =       NULL;
 136         _planeWidgetZ                           =       NULL;
 137         //EED Fev 3 2010
 138         //Free Planes
 139 }
 140
 141 //-------------------------------------------------------------------
 142 wxVtkMPR3DView::~wxVtkMPR3DView()
 143 {
 144         VisiblePointWidget( false );
 145         VisiblePlaneWidget( false );
 146         if(_pointWidget!=NULL){
 147                 _pointWidget->Delete();
 148         }
 149         if(_planeWidget!=NULL){
 150                 _planeWidget->Delete();
 151         }
 152         if(_vtkplane!=NULL){
 153                 _vtkplane->Delete();
 154                 _probe->Delete();
 155                 _contourMapper->Delete();
 156                 _contourPlaneActor->Delete();
 157         }
 158 }
 159
 160 //-------------------------------------------------------------------
 161 void wxVtkMPR3DView::RemoveActor(vtkActor* actor)
 162 {
 163         _wxvtk3Dbaseview->GetRenderer()->RemoveActor(actor);
 164 }
 165 //-------------------------------------------------------------------
 166 void wxVtkMPR3DView::Configure()
 167 {
 168         vtkImageData *imageData = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 169         _wxvtk3Dbaseview->Configure();
 170         _wxvtkmpr3DviewCntrlPanel->UpdateControlPanel();
 171         // Actors are added to the renderer.
 172         vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
 173         _wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
 174
 175         vtkmyPWCallback_3DPointWidget *_myCallback = vtkmyPWCallback_3DPointWidget::New();
 176         _myCallback->SetWxVtkMPR3DView(this);
 177         if(imageData){
 178                 if(_pointWidget==NULL)
 179                 {
 180                         _pointWidget = vtkPointWidget::New();
 181                 }
 182 //EED 2017-01-01 Migration VTK7
 183 #if VTK_MAJOR_VERSION <= 5
 184                 _pointWidget->SetInput( imageData );
 185 #else
 186                 _pointWidget->SetInputData( imageData );
 187 #endif
 188                 _myCallback->SetVtkPointWidget(_pointWidget);
 189                 _pointWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 190                 _pointWidget->AllOff();
 191                 _pointWidget->PlaceWidget();
 192                 _pointWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
 193                 if(_planeWidget==NULL)
 194                 {
 195                         _planeWidget = vtkPlaneWidget::New();
 196                 }
 197 //EED 2017-01-01 Migration VTK7
 198 #if VTK_MAJOR_VERSION <= 5
 199                 _planeWidget->SetInput( imageData );
 200 #else
 201                 _planeWidget->SetInputData( imageData );
 202 #endif
 203
 204                 _myCallback->SetVtkPlaneWidget(_planeWidget);
 205                 _planeWidget->NormalToXAxisOn();
 206                 _planeWidget->SetResolution(50);
 207                 _planeWidget->SetRepresentationToOutline();
 208                 int dim[3];
 209                 imageData->GetDimensions(dim);
 210                 int px=(dim[0]/2);
 211                 int py=(dim[1]/2);
 212                 int pz=(dim[2]/2);
 213                 int dd=20;
 214                 _planeWidget->PlaceWidget( px-dd , px+dd , py-dd , py+dd , pz-dd , pz+dd );
 215                 if(_vtkplane==NULL)
 216                 {
 217                         _vtkplane               = vtkPolyData::New();
 218                         _probe                  = vtkProbeFilter::New();
 219                         _contourMapper  = vtkPolyDataMapper::New();
 220
 221 //EED 2017-01-01 Migration VTK7
 222 #if VTK_MAJOR_VERSION <= 5
 223                         _probe->SetInput(_vtkplane);
 224                         _contourMapper->SetInput( _probe->GetPolyDataOutput() );
 225 #else
 226                         _myCallback->SetVtkPolyData(_vtkplane);
 227                         _myCallback->SetVtkProbeFilter(_probe);
 228                         _planeWidget->GetPolyData(_vtkplane);
 229                         _probe->SetInputData(_vtkplane);
 230                         _probe->SetSourceData(imageData);
 231 //                      _probe->Update();
 232                         _contourMapper->SetInputData( _probe->GetPolyDataOutput() );
 233 #endif
 234
 235                         _contourPlaneActor = vtkActor::New();
 236                         _contourPlaneActor->SetMapper(_contourMapper);
 237                         _contourPlaneActor->VisibilityOff();
 238                         _planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 239                         _planeWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
 240                         _planeWidget->Off();
 241                         _wxvtk3Dbaseview->GetRenderer()->AddActor( _contourPlaneActor );
 242                 }
 243                 _planeWidget->GetPolyData(_vtkplane);
 244
 245 //EED 2017-01-01 Migration VTK7
 246 #if VTK_MAJOR_VERSION <= 5
 247                 _probe->SetSource( imageData );
 248 #else
 249                 _probe->SetSourceData( imageData );
 250 #endif
 251                 _contourMapper->SetScalarRange( imageData->GetScalarRange() );
 252         ConfigureFreePlanes();
 253         }
 254
 255 //      _pointWidget->SetPosition( x,y,z );
 256
 257 //EED01 29Mars2009
 258 // MACOS probleme vtk-window out of wx-window
 259 //      _pointWidget->On();
 260
 261 // PlaneWidget
 262
 263
 264 // An initial camera view is created.  The Dolly() method moves
 265         // the camera towards the FocalPoint, thereby enlarging the image.
 266 //      _wxvtk3Dbaseview->GetRenderer()->SetActiveCamera(_wxvtk3Dbaseview->GetCamera());
 267
 268         // Note that when camera movement occurs (as it does in the Dolly()
 269         // method), the clipping planes often need adjusting. Clipping planes
 270         // consist of two planes: near and far along the view direction. The
 271         // near plane clips out objects in front of the plane; the far plane
 272         // clips out objects behind the plane. This way only what is drawn
 273         // between the planes is actually rendered.
 274
 275 //EED 27/05/2013
 276 //      _wxvtk3Dbaseview->GetRenderer()->ResetCameraClippingRange();
 277
 278 //EED 27/05/2013
 279 //      this->ResetCamera();
 280         _wxvtk3Dbaseview->ResetCamera();
 281
 282         // Set a background color for the renderer and set the size of the
 283         // render window (expressed in pixels).
 284         _wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.36 , 0.36 , 0.36 );
 285 //EED 23oct2010 _wxvtk3Dbaseview->GetRenWin()->SetSize(400, 400);
 286 }
 287
 288 //-------------------------------------------------------------------
 289 void wxVtkMPR3DView::ConfigureFreePlanes()
 290 {
 291         // The shared picker enables us to use 3 planes at one time
 292         // and gets the picking order right
 293         vtkCellPicker* picker = vtkCellPicker::New();
 294         picker->SetTolerance(0.005);
 295         // The 3 image plane widgets
 296         _planeWidgetX = GetPlaneWidget('x', 1, 0, 0, picker);
 297         _planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker);
 298         _planeWidgetZ = GetPlaneWidget('z', 0, 0, 1, picker);
 299         _planeWidgetY->SetLookupTable(_planeWidgetX->GetLookupTable());
 300         _planeWidgetZ->SetLookupTable(_planeWidgetX->GetLookupTable());
 301         picker->UnRegister(NULL);
 302
 303         // ColorWindow ColorLevel Callback
 304         vtkWindowLevelImagePlaneWidgetCallback *wlipwc  = vtkWindowLevelImagePlaneWidgetCallback::New();
 305         wlipwc->_wxvtkmpr3Dview                                                 = this;
 306         _planeWidgetX->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 307         _planeWidgetY->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 308         _planeWidgetZ->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 309
 310         // Cursor Callback
 311         _planeWidgetX->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 312         _planeWidgetY->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 313         _planeWidgetZ->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 314 }
 315
 316
 317 //-------------------------------------------------------------------
 318 void wxVtkMPR3DView::SetImage()
 319 {
 320         vtkImageData *imageData                                 = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 321         // Orthogonal planes B&W
 322 //EED 2017-01-01 Migration VTK7
 323 #if VTK_MAJOR_VERSION <= 5
 324         _planeWidgetX->SetInput( imageData );
 325         _planeWidgetY->SetInput( imageData );
 326         _planeWidgetZ->SetInput( imageData );
 327         // -- Plane widget
 328         _probe->SetSource( imageData );
 329 #else
 330         _planeWidgetX->SetInputData( imageData );
 331         _planeWidgetY->SetInputData( imageData );
 332         _planeWidgetZ->SetInputData( imageData );
 333         // -- Plane widget
 334         _probe->SetSourceData( imageData );
 335 #endif
 336
 337         _vtkmpr3Ddataviewer->SetImage();
 338 }
 339
 340 //-------------------------------------------------------------------
 341 vtkImagePlaneWidget* wxVtkMPR3DView::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker)
 342 {
 343         vtkProperty* prop1                                      = 0;
 344         vtkImagePlaneWidget* planeWidget        = 0;
 345         planeWidget                                                     = vtkImagePlaneWidget::New();
 346         planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 347         planeWidget->EnabledOff();
 348         vtkImageData *image                                     = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 349         int xMin = 0, xMax = 0, yMin = 0, yMax = 0, zMin = 0, zMax = 0;
 350         double xSpacing = 0, ySpacing = 0, zSpacing = 0;
 351         if(image)
 352         {
 353
 354 //EED 2017-01-01 Migration VTK7
 355 #if VTK_MAJOR_VERSION <= 5
 356                 planeWidget->SetInput( image );
 357 #else
 358                 planeWidget->SetInputData( image );
 359 #endif
 360                 image->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
 361                 image->GetSpacing(xSpacing, ySpacing, zSpacing);
 362         }
 363
 364         if (activationkey=='x')
 365         {
 366                 planeWidget->SetPlaneOrientationToXAxes();
 367                 planeWidget->SetSlicePosition((xMax+xMin)/2.*xSpacing);
 368         }
 369
 370         if (activationkey=='y')
 371         {
 372                 planeWidget->SetPlaneOrientationToYAxes();
 373                 planeWidget->SetSlicePosition((yMax+yMin)/2.*ySpacing);
 374         }
 375
 376         if (activationkey=='z')
 377         {
 378                 planeWidget->SetPlaneOrientationToZAxes();
 379                 planeWidget->SetSlicePosition((zMax+zMin)/2.*zSpacing);
 380         }
 381
 382         planeWidget->DisplayTextOn();
 383         planeWidget->SetPicker(picker);
 384         planeWidget->SetKeyPressActivationValue(activationkey);
 385         prop1 = planeWidget->GetPlaneProperty();
 386         prop1->SetColor(r, g, b);
 387         return planeWidget;
 388 }
 389
 390 //-------------------------------------------------------------------
 391 void wxVtkMPR3DView::SetFreePlanesOrtho()
 392 {
 393         _planeWidgetX->SetPlaneOrientationToXAxes();
 394         _planeWidgetY->SetPlaneOrientationToYAxes();
 395         _planeWidgetZ->SetPlaneOrientationToZAxes();
 396         vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
 397         isbv->SetParent_refresh_waiting();
 398         isbv->EvaluateToRefresh();
 399 }
 400
 401 //-------------------------------------------------------------------
 402 void wxVtkMPR3DView::FreePlaneVisible(int plane, bool ok)
 403 {
 404         if(_planeWidgetX)
 405         {
 406         if (ok==true)
 407         {
 408 //EED 21/07/2013
 409 //                      1. Make backup of actors
 410 //                      2. Putting out all the actors
 411 //                      3. Add vtkImagageActor in first place
 412 //                      4. Putting again the originals actors
 413 //                      5. Refresh
 414
 415                         // step 1, step 2
 416                 vtkActorCollection *tmpCollection    = vtkActorCollection::New();
 417                         vtkActorCollection *collectionActors =  _wxvtk3Dbaseview->GetRenderer()->GetActors();
 418                 int i,size=collectionActors->GetNumberOfItems();
 419                         for (i=0; i<size; i++)
 420                         {
 421                                 tmpCollection->AddItem ( (vtkActor*)collectionActors->GetItemAsObject(0)  );
 422                                 _wxvtk3Dbaseview->GetRenderer()->RemoveActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
 423                         } // for
 424
 425                         // step 3
 426                         if ((plane==1) && (ok==true)) _planeWidgetX->EnabledOn();
 427                         if ((plane==2) && (ok==true)) _planeWidgetY->EnabledOn();
 428                         if ((plane==3) && (ok==true)) _planeWidgetZ->EnabledOn();
 429
 430                         // step 4
 431                          for (i=0; i<size; i++)
 432                         {
 433                                 _wxvtk3Dbaseview->GetRenderer()->AddActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
 434                         } // for
 435                         tmpCollection->Delete();
 436
 437                         // step 5 Refresh
 438                         _wxvtk3Dbaseview->GetRenderer()->GetRenderWindow()->Render();
 439         } // if ok
 440                 if ((plane==1) && (ok==false)) _planeWidgetX->EnabledOff();
 441                 if ((plane==2) && (ok==false)) _planeWidgetY->EnabledOff();
 442                 if ((plane==3) && (ok==false)) _planeWidgetZ->EnabledOff();
 443         } // if  _planeWidgetX
 444 }
 445
 446 //-------------------------------------------------------------------
 447 void wxVtkMPR3DView::FreePlaneInteraction(bool ok)
 448 {
 449         if(_planeWidgetX)
 450         {
 451                 if (ok==true)
 452                 {
 453                         _planeWidgetX->InteractionOn();
 454                         _planeWidgetY->InteractionOn();
 455                         _planeWidgetZ->InteractionOn();
 456                 } else {
 457                         _planeWidgetX->InteractionOff();
 458                         _planeWidgetY->InteractionOff();
 459                         _planeWidgetZ->InteractionOff();
 460                 }
 461         } // if _planeWidgetX
 462 }
 463
 464
 465
 466 //-------------------------------------------------------------------
 467 void wxVtkMPR3DView::setColorTransferFunction(vtkColorTransferFunction* colortable)
 468 {
 469         GetVtkMPR3DDataViewer ()->setColorTransferFunction(colortable);
 470 }
 471
 472 //-------------------------------------------------------------------
 473 void wxVtkMPR3DView::Refresh()   // virtual
 474 {
 475         _vtkmpr3Ddataviewer->Refresh();
 476         if (_wxvtkmpr3DviewCntrlPanel!=NULL)
 477         {
 478                 _wxvtkmpr3DviewCntrlPanel->Refresh();
 479         }
 480 }
 481
 482 //-------------------------------------------------------------------
 483 void wxVtkMPR3DView::RefreshView()   // virtual
 484 {
 485         double spc[3];
 486         vtkImageData *image = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 487
 488         if(image)
 489         {
 490                 image->GetSpacing(spc);
 491         }else{
 492                 spc[0] = 0;
 493                 spc[1] = 0;
 494                 spc[2] = 0;
 495         }
 496
 497 //      GetVtkMPR3DDataViewer()->Configure();
 498
 499
 500         //double nx=1,ny=0,nz=0; // JPRx
 501         Refresh();
 502         double x=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetX();
 503         double y=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetY();
 504         double z=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetZ();
 505
 506         SetColorWindow( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorWindow() );
 507         SetColorLevel( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorLevel() );
 508
 509         double normal[3];
 510         // FreePlaneX
 511         _planeWidgetX->GetNormal(normal);
 512         if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
 513         {
 514                 _planeWidgetX->SetSlicePosition( x*spc[0] );
 515 //EED 2016/02/19
 516 //              _planeWidgetX->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 517                 _planeWidgetX->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 518         }
 519         // FreePlaneY
 520         _planeWidgetY->GetNormal(normal);
 521         if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
 522         {
 523                 _planeWidgetY->SetSlicePosition( y*spc[1] );
 524 //EED 2016/02/19
 525                 _planeWidgetY->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 526                 _planeWidgetY->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 527         }
 528         // FreePlaneZ
 529         _planeWidgetZ->GetNormal(normal);
 530         if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
 531         {
 532                 _planeWidgetZ->SetSlicePosition( z*spc[2] );
 533 //EED 2016/02/19
 534                 _planeWidgetZ->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 535                 _planeWidgetZ->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
 536         }
 537
 538         x=x*spc[0];
 539         y=y*spc[1];
 540         z=z*spc[2];
 541
 542         if(_pointWidget)
 543         {
 544 //EED 2017-01-01 Migration VTK7
 545 #if VTK_MAJOR_VERSION <= 5
 546                 _pointWidget->SetInput( image );
 547 #else
 548                 _pointWidget->SetInputData( image );
 549 #endif
 550                 _pointWidget->PlaceWidget();
 551
 552                 _pointWidget->SetPosition( x,y,z );
 553                 double in[3];
 554                 double normal[3];
 555                 vtkTransform *transform = GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetTransformOrientation();
 556                 in[0]=1;                in[1]=0;                        in[2]=0;
 557                 transform->TransformPoint(in,normal);
 558                 _planeWidget->SetNormal( normal[0],normal[1],normal[2] );
 559                 _planeWidget->SetCenter( x,y,z );
 560                 _planeWidget->UpdatePlacement();
 561                 _planeWidget->GetPolyData(_vtkplane);
 562 //EED 2017-01-01 Migration VTK7
 563 #if VTK_MAJOR_VERSION <= 5
 564                 // ..
 565 #else
 566                 _probe->Update();
 567 #endif
 568         }
 569
 570         if (_wxvtkmpr3DviewCntrlPanel!=NULL)
 571         {
 572                 _wxvtkmpr3DviewCntrlPanel->UpdateControlPanel(false);
 573         } // if _wxvtkmpr3DviewCntrlPanel
 574 }
 575
 576 //-------------------------------------------------------------------
 577 wxPanel* wxVtkMPR3DView::CreateControlPanel(wxWindow *parent, bool align)
 578 {
 579         _wxvtkmpr3DviewCntrlPanel = new wxVtkMPR3DViewCntrlPanel(parent,this, align);
 580         return _wxvtkmpr3DviewCntrlPanel;
 581 }
 582 //-------------------------------------------------------------------
 583 vtkMPR3DDataViewer* wxVtkMPR3DView::GetVtkMPR3DDataViewer()
 584 {
 585         return _vtkmpr3Ddataviewer;
 586 }
 587 //-------------------------------------------------------------------
 588 wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() throw (char*)
 589 {
 590         if(_wxvtk3Dbaseview==NULL)
 591         {
 592                 throw "wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() _wxvtk3Dbaseview=NULL";
 593         }
 594         return _wxvtk3Dbaseview;
 595 }
 596 //-------------------------------------------------------------------
 597
 598 //-------------------------------------------------------------------
 599 void wxVtkMPR3DView::VisibleImageActor(int idPosition, bool visible){
 600         /*if(visible == false)
 601                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == false \n");
 602         else
 603         {
 604                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == true\n");
 605                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> idPosition == %d\n", idPosition);
 606         }
 607         */
 608         if (visible!=_vtkmpr3Ddataviewer->GetVisiblePosition(idPosition)){
 609                 if (visible==false){
 610                         _wxvtk3Dbaseview->GetRenderer()->RemoveActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
 611                 } else {
 612 //EED 21/07/2013
 613 //                      1. Make backup of actors
 614 //                      2. Putting out all the actors
 615 //                      3. Add vtkImagageActor in first place
 616 //                      4. Putting again the originals actors
 617 //                      5. Refresh
 618
 619 //EED 4 nov 2015  (For Estelle)
 620 // This works in ViewerNV with nTymeView=6
 621 // Because the others modules (5) generate also more actors with an OutLineFilter -> Actor
 622 //showOutlineActor(false);
 623
 624                         // step 1, step 2
 625                         vtkActorCollection *tmpCollection    = vtkActorCollection::New();
 626                         vtkActorCollection *collectionActors =  _wxvtk3Dbaseview->GetRenderer()->GetActors();
 627                         int i,size=collectionActors->GetNumberOfItems();
 628                         for (i=0; i<size; i++)
 629                         {
 630                                 tmpCollection->AddItem ( (vtkActor*)collectionActors->GetItemAsObject(0)  );
 631                                 _wxvtk3Dbaseview->GetRenderer()->RemoveActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
 632                         } // for
 633
 634                         // step 3
 635                         _wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
 636
 637                         // step 4
 638                          for (i=0; i<size; i++)
 639                         {
 640                                 _wxvtk3Dbaseview->GetRenderer()->AddActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
 641                         } // for
 642                         tmpCollection->Delete();
 643
 644                         // step 5 Refresh
 645                         _wxvtk3Dbaseview->GetRenderer()->GetRenderWindow()->Render();
 646
 647                 } // if visible
 648                 _vtkmpr3Ddataviewer->SetVisiblePosition(idPosition,visible);
 649         } // if visible
 650 }
 651
 652
 653 //-------------------------------------------------------------------
 654 void wxVtkMPR3DView::VisiblePointWidget( bool visible )
 655 {
 656         if(_pointWidget){
 657                 if (visible==true)
 658                 {
 659                         _pointWidget->On();
 660                 } else {
 661                         _pointWidget->Off();
 662                 } // visible
 663         } // _pointWidget
 664 }
 665
 666 //-------------------------------------------------------------------
 667 void wxVtkMPR3DView::VisiblePlaneWidget( bool visible )
 668 {
 669         if(_planeWidget){
 670                 if (visible==true)
 671                 {
 672                         _planeWidget->On();
 673                         _contourPlaneActor->VisibilityOn();
 674
 675                 } else {
 676                         _planeWidget->Off();
 677                         _contourPlaneActor->VisibilityOff();
 678                 } // if visible
 679                 _wxvtk3Dbaseview->GetRenderer()->GetRenderWindow()->Render();
 680         } // _planeWidget
 681 }
 682
 683 //CPR: Method added 30 Nov 2009
 684 void wxVtkMPR3DView::showOutlineActor(bool value)
 685 {
 686         vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
 687         if(value == true)
 688         {
 689                 _wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
 690         } else {
 691                 _wxvtk3Dbaseview->GetRenderer()->RemoveActor( _outlineActor );
 692         } // value
 693 }
 694
 695 //-------------------------------------------------------------------
 696 void wxVtkMPR3DView::SetVtkMPR3DDataViewer(vtkMPR3DDataViewer *vtkmpr3Ddataviewer)
 697 {
 698         _vtkmpr3Ddataviewer = vtkmpr3Ddataviewer;
 699 }
 700
 701 //-------------------------------------------------------------------
 702 void wxVtkMPR3DView::InitOrientationPointWidget()
 703 {
 704         this->_planeWidget->SetNormal(1,0,0);
 705 }
 706
 707 //-------------------------------------------------------------------
 708 void wxVtkMPR3DView::SetColorWindow(double colorWindow)
 709 {
 710         double colorLevel = _planeWidgetX->GetLevel();
 711         _planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
 712 }
 713
 714 //-------------------------------------------------------------------
 715 void wxVtkMPR3DView::SetColorLevel(double colorLevel)
 716 {
 717         double colorWindow = _planeWidgetX->GetWindow();
 718         _planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
 719 }
 720
 721 //-------------------------------------------------------------------
 722 // EED 25 Janvier 2007 testLoic
 723 void wxVtkMPR3DView::TestLoic1()
 724 {
 725
 726         int sizeIma     = 128;
 727   double dimIma = sizeIma*2;
 728
 729         vtkImageData *vtkimagedata =  this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 730
 731         double spc[3];
 732         vtkimagedata->GetSpacing(spc);
 733
 734 //EED 2017-01-01 Migration VTK7
 735 #if VTK_MAJOR_VERSION <= 5
 736         vtkimagedata->Update();
 737 #else
 738         // ..
 739 #endif
 740
 741
 742          double p[3], n[3];
 743      this->_planeWidget->GetOrigin(p);
 744      this->_planeWidget->GetNormal(n);
 745
 746
 747
 748     vtkPlaneSource* pSource = vtkPlaneSource::New( );
 749     pSource->SetResolution( sizeIma - 1, sizeIma - 1 );
 750     pSource->SetOrigin( p[ 0 ], p[ 1 ], p[ 2 ] );
 751     pSource->SetPoint1( p[ 0 ] + dimIma - 1.0, p[ 1 ], p[ 2 ] );
 752     pSource->SetPoint2( p[ 0 ], p[ 1 ]+ dimIma - 1.0 , p[ 2 ] );
 753     pSource->Update( );
 754     pSource->SetCenter( p[ 0 ], p[ 1 ], p[ 2 ] );
 755     pSource->SetNormal( n[ 0 ], n[ 1 ], n[ 2 ] );
 756     pSource->Update( );
 757
 758     vtkProbeFilter* slices = vtkProbeFilter::New();
 759
 760 //EED 2017-01-01 Migration VTK7
 761 #if VTK_MAJOR_VERSION <= 5
 762     slices->SetInput( ( vtkDataSet* )pSource->GetOutput( ) );
 763     slices->SetSource( vtkimagedata );
 764     slices->Update( );
 765     pSource->Delete( );
 766 #else
 767     slices->SetInputData( ( vtkDataSet* )pSource->GetOutput( ) );
 768     slices->SetSourceData( vtkimagedata );
 769 #endif
 770
 771         vtkStructuredPoints   *stPoints = vtkStructuredPoints::New();
 772         stPoints -> GetPointData( )->SetScalars(  slices->GetOutput()->GetPointData()->GetScalars()  );
 773         stPoints -> SetDimensions( sizeIma, sizeIma, 1 );
 774
 775 //EED 2017-01-01 Migration VTK7
 776 #if VTK_MAJOR_VERSION <= 5
 777         stPoints -> SetScalarType( vtkimagedata->GetScalarType() );
 778         stPoints -> SetScalarTypeToShort();
 779         stPoints -> Update();
 780 #else
 781         vtkInformation* info=stPoints->GetInformation();
 782         vtkDataObject::SetPointDataActiveScalarInfo(info, VTK_SHORT, 1);
 783 #endif
 784
 785
 786
 787         vtkImageChangeInformation  *change = vtkImageChangeInformation ::New();
 788 //EED 2017-01-01 Migration VTK7
 789 #if VTK_MAJOR_VERSION <= 5
 790         change    -> SetInput( stPoints );
 791 #else
 792         change    -> SetInputData( stPoints );
 793 #endif
 794         change    -> Update();    //important
 795
 796         double _range[2];
 797         vtkimagedata->GetScalarRange(_range);
 798    vtkWindowLevelLookupTable *_bwlookup = vtkWindowLevelLookupTable::New( );
 799    _bwlookup->SetHueRange( 0 , 1 );
 800    _bwlookup->SetNumberOfColors( (int)(_range[1] - _range[0] + 1) );
 801    _bwlookup->SetTableRange( _range[0] , _range[1] );
 802    _bwlookup->SetSaturationRange( 0 , 0 );
 803    _bwlookup->SetValueRange( 0 , 1 );
 804    _bwlookup->SetAlphaRange( 1 , 1 );
 805    _bwlookup->Build( );
 806
 807    vtkLookupTable * _collookup = vtkLookupTable::New( );
 808    _collookup->SetNumberOfColors( 256 );
 809    _collookup->SetTableRange( 0 , 255 );
 810    _collookup->Build( );
 811    _collookup->SetTableValue( 0  , 1 , 0 , 0 , 1 );
 812    _collookup->SetTableValue(128 , 0 , 0 , 1 , 1 );
 813    _collookup->SetTableValue(255 , 0 , 1 , 0 , 1 );
 814
 815
 816         vtkMetaImageWriter *writer = vtkMetaImageWriter::New( );
 817 //EED 2017-01-01 Migration VTK7
 818 #if VTK_MAJOR_VERSION <= 5
 819         writer->SetInput( stPoints );
 820 #else
 821         writer->SetInputData( stPoints );
 822 #endif
 823         writer->SetFileName( "C:/Users/Images/temp_EED/image.mhd" );
 824         writer->SetFileDimensionality( 2 );
 825         writer->Write( );
 826
 827
 828
 829         vtkDataSetMapper *_3DSliceMapper = vtkDataSetMapper::New( );
 830
 831 //EED 2017-01-01 Migration VTK7
 832 #if VTK_MAJOR_VERSION <= 5
 833         _3DSliceMapper->SetInput(change->GetOutput( ) );
 834 #else
 835         _3DSliceMapper->SetInputData(change->GetOutput( ) );
 836 #endif
 837
 838         _3DSliceMapper->SetLookupTable( _bwlookup );
 839         _3DSliceMapper->SetScalarRange( _range );
 840         _3DSliceMapper->ImmediateModeRenderingOn( );
 841
 842         vtkActor *_3DSliceActor = vtkActor::New( );
 843         _3DSliceActor->SetMapper( _3DSliceMapper );
 844
 845  // The usual rendering stuff.
 846   vtkCamera *camera = vtkCamera::New();
 847       camera->SetPosition(1,1,1);
 848       camera->SetFocalPoint(0,0,0);
 849
 850   vtkRenderer *renderer = vtkRenderer::New();
 851   vtkRenderWindow *renWin = vtkRenderWindow::New();
 852     renWin->AddRenderer(renderer);
 853
 854   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
 855     iren->SetRenderWindow(renWin);
 856
 857   renderer->AddActor( _3DSliceActor );
 858       renderer->SetActiveCamera(camera);
 859       renderer->ResetCamera();
 860       renderer->SetBackground(1,1,1);
 861
 862   renWin->SetSize(300,300);
 863
 864   // interact with data
 865   renWin->Render();
 866   iren->Start();
 867
 868
 869
 870
 871
 872 }
 873
 874 //-------------------------------------------------------------------
 875 // EED 25 Janvier 2007 testLoic
 876 void wxVtkMPR3DView::TestLoic2()
 877 {
 878 /*
 879  vtkVectorNorm *vecMagnitude = vtkVectorNorm::New();
 880                 vecMagnitude->SetInput(VtkMainGrid);
 881                 vecMagnitude->NormalizeOff();
 882                 vecMagnitude->Update();
 883 */
 884
 885         vtkPlane *slicePlane = vtkPlane::New();
 886                                 this->_planeWidget->GetPlane( slicePlane );
 887
 888         vtkCutter* sliceCutter = vtkCutter::New();
 889                 vtkImageData *vtkimagedata =  this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 890 //EED 2017-01-01 Migration VTK7
 891 #if VTK_MAJOR_VERSION <= 5
 892                 sliceCutter->SetInput( vtkimagedata );
 893 #else
 894                 sliceCutter->SetInputData( vtkimagedata );
 895 #endif
 896                 sliceCutter->SetCutFunction( slicePlane );
 897
 898 //EED
 899 //        vtkLookupTable *lut = BuildHueWeightBaseMap();
 900
 901         vtkPolyDataMapper *slice = vtkPolyDataMapper::New();
 902 //EED 2017-01-01 Migration VTK7
 903 #if VTK_MAJOR_VERSION <= 5
 904                 slice->SetInput( sliceCutter->GetOutput() );
 905 #else
 906                 slice->SetInputData( sliceCutter->GetOutput() );
 907 #endif
 908         double range[2];
 909 // EED
 910 //                slice->Update();
 911                                   vtkimagedata->GetScalarRange(range);
 912 //                vecMagnitude->GetOutput()->GetScalarRange( range );
 913 //                range[1] *= 0.7; // reduce the upper range by 30%
 914 //                slice->SetScalarRange( range );
 915 //                slice->SetLookupTable( lut );
 916
 917 //       vtkActor *sliceActor = vtkActor::New();
 918 //                sliceActor->SetMapper( slice );
 919
 920         vtkPolyDataMapper *contourMapper = vtkPolyDataMapper::New();
 921
 922 //EED 2017-01-01 Migration VTK7
 923 #if VTK_MAJOR_VERSION <= 5
 924                 contourMapper->SetInput( sliceCutter->GetOutput() );
 925 #else
 926                 contourMapper->SetInputData( sliceCutter->GetOutput() );
 927 #endif
 928
 929                 contourMapper->SetScalarRange( range );
 930   //              contourMapper->SetLookupTable( lut );
 931
 932         vtkActor *contourActor = vtkActor::New();
 933                 contourActor->SetMapper( contourMapper );
 934
 935  // The usual rendering stuff.
 936   vtkCamera *camera = vtkCamera::New();
 937       camera->SetPosition(1,1,1);
 938       camera->SetFocalPoint(0,0,0);
 939
 940   vtkRenderer *renderer = vtkRenderer::New();
 941   vtkRenderWindow *renWin = vtkRenderWindow::New();
 942     renWin->AddRenderer(renderer);
 943
 944   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
 945     iren->SetRenderWindow(renWin);
 946
 947   renderer->AddActor( contourActor );
 948 //  renderer->AddActor(cubeActor);
 949       renderer->SetActiveCamera(camera);
 950       renderer->ResetCamera();
 951       renderer->SetBackground(1,1,1);
 952
 953   renWin->SetSize(300,300);
 954
 955   // interact with data
 956   renWin->Render();
 957   iren->Start();
 958 }
 959
 960 //-------------------------------------------------------------------
 961 void wxVtkMPR3DView::SetBackGroundType(int type)
 962 {
 963         if (_wxvtk3Dbaseview!=NULL)
 964         {
 965                 if (type==0)
 966                 {
 967                         _wxvtk3Dbaseview->GetRenderer()->GradientBackgroundOff();
 968                         _wxvtk3Dbaseview->GetRenderer()->SetBackground( 0 , 0 , 0 );
 969                 } // if type
 970                 if (type==1)
 971                 {
 972                         _wxvtk3Dbaseview->GetRenderer()->GradientBackgroundOn();
 973                         _wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.33 , 0.33 , 0.33 );
 974                         _wxvtk3Dbaseview->GetRenderer()->SetBackground2( 0.66 , 0.66 , 0.66 );
 975                 } // if type
 976                 _wxvtk3Dbaseview->Refresh();
 977         }// if _wxvtk3Dbaseview
 978
 979 }
 980
 981