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

   1 #include "wxVtkMPR3DView.h"
   2
   3 #include "vtkDataSetMapper.h"
   4 #include "vtkPlaneSource.h"
   5 #include "vtkProbeFilter.h"
   6 #include "vtkStructuredPoints.h"
   7 #include "vtkMetaImageWriter.h"
   8 #include "vtkImageChangeInformation.h"
   9 #include "vtkPointData.h"
  10 #include "vtkWindowLevelLookupTable.h"
  11 #include "vtkPlane.h"
  12 #include "vtkCutter.h"
  13 #include "vtkTransform.h"
  14 #include "vtkmyPWCallback_3DPointWidget.h"
  15 #include "vtkProperty.h"
  16 #include "vtkCellPicker.h"
  17 #include "vtkInteractorStyleBaseView.h"
  18 #include <vtkImageReslice.h>
  19
  20 #ifdef _DEBUG
  21 #define new DEBUG_NEW
  22 #endif
  23
  24
  25
  26 //----------------------------------------------------------------------------
  27 //----------------------------------------------------------------------------
  28 //----------------------------------------------------------------------------
  29 class vtkWindowLevelImagePlaneWidgetCallback : public vtkCommand
  30 {
  31          public:
  32            static vtkWindowLevelImagePlaneWidgetCallback *New()
  33            {
  34                    return new vtkWindowLevelImagePlaneWidgetCallback;
  35            }
  36
  37            void Execute( vtkObject *caller, unsigned long ev, void *callData )
  38                 {
  39                         bool needRefresh=false;
  40                         double slicePosition;
  41                         vtkImagePlaneWidget* ipw =dynamic_cast< vtkImagePlaneWidget* >( caller );
  42                     if ( ev == vtkCommand::WindowLevelEvent)
  43                         {
  44                                 double info[2];
  45                                 ipw->GetWindowLevel(info);
  46                                 _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorWindow(info[0]);
  47                                 _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorLevel(info[1]);
  48                                 needRefresh=true;
  49                         } // if ev
  50
  51                    if ( ev == vtkCommand::InteractionEvent)
  52                    {
  53                            needRefresh=true;
  54                            if (ipw->GetCursorDataStatus())
  55                            {
  56                                    double xyzv[4];
  57                                    ipw->GetCursorData(xyzv);
  58                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX(xyzv[0]);
  59                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY(xyzv[1]);
  60                                    _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ(xyzv[2]);
  61                            } else {
  62
  63
  64                                    vtkImageData *image = _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
  65                                    double spc[3];
  66                                    if(image)
  67                                    {
  68                                            image->GetSpacing(spc);
  69                                    }
  70
  71                                    double normal[3];
  72                                    slicePosition = ipw->GetSlicePosition();
  73                                    ipw->GetNormal(normal);
  74
  75                                    // FreePlaneX
  76                                    if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
  77                                    {
  78                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX( slicePosition/spc[0] );
  79                                    }
  80                                    // FreePlaneY
  81                                    if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
  82                                    {
  83                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY( slicePosition/spc[1] );
  84                                    }
  85                                    // FreePlaneZ
  86                                    if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
  87                                    {
  88                                            _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ( slicePosition/spc[2] );
  89                                    }
  90                            }// ipw->GetCursorDataStatus
  91                    }// ev
  92
  93                    if (needRefresh)
  94                    {
  95                            vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)_wxvtkmpr3Dview->GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
  96                            isbv->SetParent_refresh_waiting();
  97                            isbv->EvaluateToRefresh();
  98                    }
  99
 100            } // Execute
 101
 102                 vtkWindowLevelImagePlaneWidgetCallback(){}
 103                 wxVtkMPR3DView *_wxvtkmpr3Dview;
 104 };
 105
 106
 107
 108 //-------------------------------------------------------------------
 109 //-------------------------------------------------------------------
 110 //-------------------------------------------------------------------
 111
 112
 113
 114 wxVtkMPR3DView::wxVtkMPR3DView( wxVtk3DBaseView *wxvtk3Dbaseview )
 115 {
 116         _wxvtk3Dbaseview                        =       wxvtk3Dbaseview;
 117         _vtkmpr3Ddataviewer                     =       NULL;
 118         _wxvtkmpr3DviewCntrlPanel       =       NULL;
 119         //_myCallback                           =       NULL;
 120         _pointWidget                            =       NULL;
 121         _planeWidget                            =       NULL;
 122         _vtkplane                                       =       NULL;
 123         _probe                                          =       NULL;
 124         _contourMapper                          =       NULL;
 125
 126         _planeWidgetX                           =       NULL;
 127         _planeWidgetY                           =       NULL;
 128         _planeWidgetZ                           =       NULL;
 129         //EED Fev 3 2010
 130         //Free Planes
 131 }
 132
 133
 134 //-------------------------------------------------------------------
 135 wxVtkMPR3DView::~wxVtkMPR3DView()
 136 {
 137         VisiblePointWidget( false );
 138         VisiblePlaneWidget( false );
 139
 140         if(_pointWidget!=NULL){
 141                 _pointWidget->Delete();
 142         }
 143         if(_planeWidget!=NULL){
 144                 _planeWidget->Delete();
 145         }
 146
 147         if(_vtkplane!=NULL){
 148                 _vtkplane->Delete();
 149                 _probe->Delete();
 150                 _contourMapper->Delete();
 151                 _contourPlaneActor->Delete();
 152         }
 153
 154 }
 155
 156 //-------------------------------------------------------------------
 157 void wxVtkMPR3DView::RemoveActor(vtkActor* actor)
 158 {
 159         _wxvtk3Dbaseview->GetRenderer()->RemoveActor(actor);
 160 }
 161
 162 //-------------------------------------------------------------------
 163 void wxVtkMPR3DView::ResetCamera(int* ext, double* origin,double* spc)
 164 {
 165         if(ext == NULL){
 166                 _wxvtk3Dbaseview->GetRenderer()->ResetCamera ();
 167                 _wxvtk3Dbaseview->GetCamera()->Dolly(2.0);
 168
 169 //EED 21 mars 2012  FLIP probleme  ..PLOP..
 170                 _wxvtk3Dbaseview->GetCamera()->Roll(180);
 171                 _wxvtk3Dbaseview->GetCamera()->Yaw(180+30);
 172                 _wxvtk3Dbaseview->GetCamera()->Pitch(-22);
 173
 174         }else{
 175
 176                 /*double x = (spc[0])*(origin[0]+(((double)ext[1]-(double)ext[0])/2.0));
 177                 double y = (spc[1])*(origin[1]+(double)ext[3]);
 178                 double z = (spc[2])*(origin[2]+(((double)ext[5]-(double)ext[4])/2.0));*/
 179                 /*double x0=(spc[0])*((double)ext[0]+origin[0]);
 180                 double x1=(spc[0])*((double)ext[1]+origin[0]);
 181                 double y0=(spc[1])*((double)ext[2]+origin[1]);
 182                 double y1=(spc[1])*((double)ext[3]+origin[1]);
 183                 double z0=(spc[2])*((double)ext[4]+origin[2]);
 184                 double z1=(spc[2])*((double)ext[5]+origin[2]);*/
 185
 186                 double x0=(spc[0])*((double)ext[0]);
 187                 double x1=(spc[0])*((double)ext[1]);
 188                 double y0=(spc[1])*((double)ext[2]);
 189                 double y1=(spc[1])*((double)ext[3]);
 190                 double z0=(spc[2])*((double)ext[4]);
 191                 double z1=(spc[2])*((double)ext[5]);
 192
 193                 _wxvtk3Dbaseview->GetRenderer()->ResetCamera(x0,x1,y0,y1,z0,z1);
 194                 //_wxvtk3Dbaseview->GetCamera()->SetPosition(x,y,z);
 195                 _wxvtk3Dbaseview->GetCamera()->Dolly(1.5);
 196
 197         }
 198
 199 }
 200
 201 //-------------------------------------------------------------------
 202 void wxVtkMPR3DView::Configure()
 203 {
 204         vtkImageData *imageData = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 205
 206         _wxvtk3Dbaseview->Configure();
 207         _wxvtkmpr3DviewCntrlPanel->UpdateControlPanel();
 208
 209 //EED 27 Mai 2009
 210 //      _wxvtk3Dbaseview->GetRenderer()->Clear();
 211
 212         // Actors are added to the renderer.
 213         vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
 214         _wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
 215 //      _wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(0)  );    // _saggital
 216 //      _wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(1)  );    // _axial
 217 //      _wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(2)  );    // _coronal
 218
 219
 220         // An initial camera view is created.  The Dolly() method moves
 221         // the camera towards the FocalPoint, thereby enlarging the image.
 222         _wxvtk3Dbaseview->GetRenderer()->SetActiveCamera(_wxvtk3Dbaseview->GetCamera());
 223         this->ResetCamera();
 224
 225
 226         // Set a background color for the renderer and set the size of the
 227         // render window (expressed in pixels).
 228         _wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.36 , 0.36 , 0.36 );
 229 //EED 23oct2010 _wxvtk3Dbaseview->GetRenWin()->SetSize(400, 400);
 230
 231         // Note that when camera movement occurs (as it does in the Dolly()
 232         // method), the clipping planes often need adjusting. Clipping planes
 233         // consist of two planes: near and far along the view direction. The
 234         // near plane clips out objects in front of the plane; the far plane
 235         // clips out objects behind the plane. This way only what is drawn
 236         // between the planes is actually rendered.
 237         _wxvtk3Dbaseview->GetRenderer()->ResetCameraClippingRange();
 238
 239
 240         // vtkPointWidget
 241         //if(_myCallback!=NULL){
 242                 //_myCallback->Delete();
 243
 244         //}
 245         vtkmyPWCallback_3DPointWidget *_myCallback = vtkmyPWCallback_3DPointWidget::New();
 246         _myCallback->SetWxVtkMPR3DView(this);
 247
 248
 249         if(imageData){
 250                 if(_pointWidget==NULL){
 251                         _pointWidget = vtkPointWidget::New();
 252                 }
 253                 _pointWidget->SetInput( imageData );
 254                 _myCallback->SetVtkPointWidget(_pointWidget);
 255                 _pointWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 256                 _pointWidget->AllOff();
 257                 _pointWidget->PlaceWidget();
 258
 259                 _pointWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
 260
 261                 if(_planeWidget==NULL)
 262                 {
 263                         _planeWidget = vtkPlaneWidget::New();
 264                 }
 265                 _planeWidget->SetInput( imageData );
 266
 267                 _myCallback->SetVtkPlaneWidget(_planeWidget);
 268
 269                 _planeWidget->NormalToXAxisOn();
 270                 _planeWidget->SetResolution(50);
 271                 _planeWidget->SetRepresentationToOutline();
 272
 273                 int dim[3];
 274
 275                 imageData->GetDimensions(dim);
 276                 int px=(dim[0]/2);
 277                 int py=(dim[1]/2);
 278                 int pz=(dim[2]/2);
 279                 int dd=20;
 280                 _planeWidget->PlaceWidget( px-dd , px+dd , py-dd , py+dd , pz-dd , pz+dd );
 281
 282                 if(_vtkplane==NULL){
 283                         _vtkplane = vtkPolyData::New();
 284
 285                         _probe = vtkProbeFilter::New();
 286                         _probe->SetInput(_vtkplane);
 287
 288                         _contourMapper = vtkPolyDataMapper::New();
 289
 290                         _contourMapper->SetInput( _probe->GetPolyDataOutput() );
 291
 292                         _contourPlaneActor = vtkActor::New();
 293                         _contourPlaneActor->SetMapper(_contourMapper);
 294
 295                         _contourPlaneActor->VisibilityOff();
 296
 297                         _planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 298                         _planeWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
 299
 300                         _planeWidget->Off();
 301                         _wxvtk3Dbaseview->GetRenderer()->AddActor( _contourPlaneActor );
 302
 303                 }
 304
 305                 _planeWidget->GetPolyData(_vtkplane);
 306                 _probe->SetSource( imageData );
 307                 _contourMapper->SetScalarRange( imageData->GetScalarRange() );
 308
 309         ConfigureFreePlanes();
 310         }
 311
 312 //      _pointWidget->SetPosition( x,y,z );
 313
 314 //EED01 29Mars2009
 315 // MACOS probleme vtk-window out of wx-window
 316 //      _pointWidget->On();
 317
 318 // PlaneWidget
 319
 320 }
 321
 322 //-------------------------------------------------------------------
 323 void wxVtkMPR3DView::ConfigureFreePlanes()
 324 {
 325
 326         // The shared picker enables us to use 3 planes at one time
 327         // and gets the picking order right
 328         vtkCellPicker* picker = vtkCellPicker::New();
 329         picker->SetTolerance(0.005);
 330
 331         // The 3 image plane widgets
 332         _planeWidgetX = GetPlaneWidget('x', 1, 0, 0, picker);
 333         _planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker);
 334         _planeWidgetZ = GetPlaneWidget('z', 0, 0, 1, picker);
 335
 336         _planeWidgetY->SetLookupTable(_planeWidgetX->GetLookupTable());
 337         _planeWidgetZ->SetLookupTable(_planeWidgetX->GetLookupTable());
 338
 339         picker->UnRegister(NULL);
 340
 341         // ColorWindow ColorLevel Callback
 342         vtkWindowLevelImagePlaneWidgetCallback *wlipwc  = vtkWindowLevelImagePlaneWidgetCallback::New();
 343         wlipwc->_wxvtkmpr3Dview                                                 = this;
 344         _planeWidgetX->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 345         _planeWidgetY->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 346         _planeWidgetZ->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 347
 348         // Cursor Callback
 349         _planeWidgetX->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 350         _planeWidgetY->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 351         _planeWidgetZ->AddObserver(vtkCommand::InteractionEvent, wlipwc );
 352 }
 353
 354 //-------------------------------------------------------------------
 355 vtkImagePlaneWidget* wxVtkMPR3DView::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker)
 356 {
 357         vtkProperty* prop1                                      = 0;
 358         vtkImagePlaneWidget* planeWidget        = 0;
 359         planeWidget                                                     = vtkImagePlaneWidget::New();
 360         planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
 361         planeWidget->EnabledOff();
 362         vtkImageData *image                                     = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 363         int xMin = 0, xMax = 0, yMin = 0, yMax = 0, zMin = 0, zMax = 0;
 364         double xSpacing = 0, ySpacing = 0, zSpacing = 0;
 365         if(image)
 366         {
 367                 planeWidget->SetInput( image );
 368                 image->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
 369                 image->GetSpacing(xSpacing, ySpacing, zSpacing);
 370         }
 371
 372         if (activationkey=='x')
 373         {
 374                 planeWidget->SetPlaneOrientationToXAxes();
 375                 planeWidget->SetSlicePosition((xMax+xMin)/2.*xSpacing);
 376         }
 377
 378         if (activationkey=='y')
 379         {
 380                 planeWidget->SetPlaneOrientationToYAxes();
 381                 planeWidget->SetSlicePosition((yMax+yMin)/2.*ySpacing);
 382         }
 383
 384         if (activationkey=='z')
 385         {
 386                 planeWidget->SetPlaneOrientationToZAxes();
 387                 planeWidget->SetSlicePosition((zMax+zMin)/2.*zSpacing);
 388         }
 389
 390         planeWidget->DisplayTextOn();
 391         planeWidget->SetPicker(picker);
 392         planeWidget->SetKeyPressActivationValue(activationkey);
 393         prop1 = planeWidget->GetPlaneProperty();
 394         prop1->SetColor(r, g, b);
 395         return planeWidget;
 396 }
 397
 398 //-------------------------------------------------------------------
 399 void wxVtkMPR3DView::SetFreePlanesOrtho()
 400 {
 401         _planeWidgetX->SetPlaneOrientationToXAxes();
 402         _planeWidgetY->SetPlaneOrientationToYAxes();
 403         _planeWidgetZ->SetPlaneOrientationToZAxes();
 404         vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
 405         isbv->SetParent_refresh_waiting();
 406         isbv->EvaluateToRefresh();
 407 }
 408
 409 //-------------------------------------------------------------------
 410 void wxVtkMPR3DView::FreePlaneVisible(int plane, bool ok)
 411 {
 412         if(_planeWidgetX)
 413         {
 414                 if ((plane==1) && (ok==true)) _planeWidgetX->EnabledOn();
 415                 if ((plane==2) && (ok==true)) _planeWidgetY->EnabledOn();
 416                 if ((plane==3) && (ok==true)) _planeWidgetZ->EnabledOn();
 417                 if ((plane==1) && (ok==false)) _planeWidgetX->EnabledOff();
 418                 if ((plane==2) && (ok==false)) _planeWidgetY->EnabledOff();
 419                 if ((plane==3) && (ok==false)) _planeWidgetZ->EnabledOff();
 420         }
 421 }
 422
 423 //-------------------------------------------------------------------
 424 void wxVtkMPR3DView::FreePlaneInteraction(bool ok)
 425 {
 426         if(_planeWidgetX)
 427         {
 428                 if (ok==true)
 429                 {
 430                         _planeWidgetX->InteractionOn();
 431                         _planeWidgetY->InteractionOn();
 432                         _planeWidgetZ->InteractionOn();
 433                 } else {
 434                         _planeWidgetX->InteractionOff();
 435                         _planeWidgetY->InteractionOff();
 436                         _planeWidgetZ->InteractionOff();
 437                 }
 438         } // if _planeWidgetX
 439 }
 440
 441
 442
 443 //-------------------------------------------------------------------
 444 void wxVtkMPR3DView::setColorTransferFunction(vtkColorTransferFunction* colortable)
 445 {
 446         GetVtkMPR3DDataViewer ()->setColorTransferFunction(colortable);
 447 }
 448
 449 //-------------------------------------------------------------------
 450 void wxVtkMPR3DView::Refresh()   // virtual
 451 {
 452         _vtkmpr3Ddataviewer     ->      Refresh();
 453         if (_wxvtkmpr3DviewCntrlPanel!=NULL)
 454         {
 455                 _wxvtkmpr3DviewCntrlPanel->Refresh();
 456         }
 457
 458 }
 459
 460 //-------------------------------------------------------------------
 461 void wxVtkMPR3DView::RefreshView()   // virtual
 462 {
 463         double spc[3];
 464         vtkImageData *image = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 465
 466         if(image){
 467                 image->GetSpacing(spc);
 468         }else{
 469                 spc[0] = 0;
 470                 spc[1] = 0;
 471                 spc[2] = 0;
 472         }
 473         //this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData()->GetSpacing(spc);
 474
 475         //double nx=1,ny=0,nz=0; // JPRx
 476         Refresh();
 477         double x=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetX();
 478         double y=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetY();
 479         double z=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetZ();
 480
 481         SetColorWindow( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorWindow() );
 482         SetColorLevel( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorLevel() );
 483
 484
 485         double normal[3];
 486         // FreePlaneX
 487         _planeWidgetX->GetNormal(normal);
 488         if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
 489         {
 490                 _planeWidgetX->SetSlicePosition( x*spc[0] );
 491         }
 492         // FreePlaneY
 493         _planeWidgetY->GetNormal(normal);
 494         if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
 495         {
 496                 _planeWidgetY->SetSlicePosition( y*spc[1] );
 497         }
 498         // FreePlaneZ
 499         _planeWidgetZ->GetNormal(normal);
 500         if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
 501         {
 502                 _planeWidgetZ->SetSlicePosition( z*spc[2] );
 503         }
 504
 505         x=x*spc[0];
 506         y=y*spc[1];
 507         z=z*spc[2];
 508
 509         if(_pointWidget){
 510                 _pointWidget->SetPosition( x,y,z );
 511                 double in[3];
 512                 double normal[3];
 513                 vtkTransform *transform = GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetTransformOrientation();
 514                 in[0]=1;                in[1]=0;                        in[2]=0;
 515                 transform->TransformPoint(in,normal);
 516                 _planeWidget->SetNormal( normal[0],normal[1],normal[2] );
 517                 _planeWidget->SetCenter( x,y,z );
 518                 _planeWidget->UpdatePlacement();
 519                 _planeWidget->GetPolyData(_vtkplane);
 520         }
 521 }
 522
 523 //-------------------------------------------------------------------
 524 wxPanel* wxVtkMPR3DView::CreateControlPanel(wxWindow *parent, bool align)
 525 {
 526         _wxvtkmpr3DviewCntrlPanel = new wxVtkMPR3DViewCntrlPanel(parent,this, align);
 527         return _wxvtkmpr3DviewCntrlPanel;
 528 }
 529 //-------------------------------------------------------------------
 530 vtkMPR3DDataViewer* wxVtkMPR3DView::GetVtkMPR3DDataViewer()
 531 {
 532         return _vtkmpr3Ddataviewer;
 533 }
 534 //-------------------------------------------------------------------
 535 wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() throw (char*)
 536 {
 537         if(_wxvtk3Dbaseview==NULL){
 538                 throw "wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() _wxvtk3Dbaseview=NULL";
 539         }
 540         return _wxvtk3Dbaseview;
 541 }
 542 //-------------------------------------------------------------------
 543
 544 //-------------------------------------------------------------------
 545 void wxVtkMPR3DView::VisibleImageActor(int idPosition, bool visible){
 546         /*if(visible == false)
 547                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == false \n");
 548         else
 549         {
 550                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == true\n");
 551                 printf("CPR: wxVtkMPR3DView::VisibleImageActor-> idPosition == %d\n", idPosition);
 552         }
 553         */
 554         if (visible!=_vtkmpr3Ddataviewer->GetVisiblePosition(idPosition)){
 555                 if (visible==false){
 556                         _wxvtk3Dbaseview->GetRenderer()->RemoveActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
 557                 } else {
 558                         _wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
 559                 }
 560                 _vtkmpr3Ddataviewer->SetVisiblePosition(idPosition,visible);
 561         }
 562 }
 563
 564
 565 //-------------------------------------------------------------------
 566 void wxVtkMPR3DView::VisiblePointWidget( bool visible )
 567 {
 568         if(_pointWidget){
 569                 if (visible==true)
 570                 {
 571                         _pointWidget->On();
 572                 } else {
 573                         _pointWidget->Off();
 574                 }
 575         }
 576 }
 577
 578 //-------------------------------------------------------------------
 579 void wxVtkMPR3DView::VisiblePlaneWidget( bool visible )
 580 {
 581         if(_planeWidget){
 582                 if (visible==true)
 583                 {
 584                         _planeWidget->On();
 585                         _contourPlaneActor->VisibilityOn();
 586
 587                 } else {
 588                         _planeWidget->Off();
 589                         _contourPlaneActor->VisibilityOff();
 590                 }
 591         }
 592 }
 593
 594 //CPR: Method added 30 Nov 2009
 595 void wxVtkMPR3DView::showOutlineActor(bool value)
 596 {
 597         vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
 598         if(value == true)
 599         {
 600                 _wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
 601         }
 602         else
 603         {
 604                 _wxvtk3Dbaseview->GetRenderer()->RemoveActor( _outlineActor );
 605         }
 606 }
 607
 608 //-------------------------------------------------------------------
 609 void wxVtkMPR3DView::SetVtkMPR3DDataViewer(vtkMPR3DDataViewer *vtkmpr3Ddataviewer)
 610 {
 611         _vtkmpr3Ddataviewer = vtkmpr3Ddataviewer;
 612 }
 613
 614 //-------------------------------------------------------------------
 615 void wxVtkMPR3DView::InitOrientationPointWidget()
 616 {
 617         this->_planeWidget->SetNormal(1,0,0);
 618 }
 619
 620 //-------------------------------------------------------------------
 621 void wxVtkMPR3DView::SetColorWindow(double colorWindow)
 622 {
 623         double colorLevel = _planeWidgetX->GetLevel();
 624         _planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
 625 }
 626
 627 //-------------------------------------------------------------------
 628 void wxVtkMPR3DView::SetColorLevel(double colorLevel)
 629 {
 630         double colorWindow = _planeWidgetX->GetWindow();
 631         _planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
 632 }
 633
 634 //-------------------------------------------------------------------
 635 // EED 25 Janvier 2007 testLoic
 636 void wxVtkMPR3DView::TestLoic1()
 637 {
 638
 639         int sizeIma     = 128;
 640   double dimIma = sizeIma*2;
 641
 642         vtkImageData *vtkimagedata =  this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 643
 644         double spc[3];
 645         vtkimagedata->GetSpacing(spc);
 646
 647         vtkimagedata->Update();
 648
 649
 650          double p[3], n[3];
 651      this->_planeWidget->GetOrigin(p);
 652      this->_planeWidget->GetNormal(n);
 653
 654
 655
 656     vtkPlaneSource* pSource = vtkPlaneSource::New( );
 657     pSource->SetResolution( sizeIma - 1, sizeIma - 1 );
 658     pSource->SetOrigin( p[ 0 ], p[ 1 ], p[ 2 ] );
 659     pSource->SetPoint1( p[ 0 ] + dimIma - 1.0, p[ 1 ], p[ 2 ] );
 660     pSource->SetPoint2( p[ 0 ], p[ 1 ]+ dimIma - 1.0 , p[ 2 ] );
 661     pSource->Update( );
 662     pSource->SetCenter( p[ 0 ], p[ 1 ], p[ 2 ] );
 663     pSource->SetNormal( n[ 0 ], n[ 1 ], n[ 2 ] );
 664     pSource->Update( );
 665
 666     vtkProbeFilter* slices = vtkProbeFilter::New();
 667     slices->SetInput( ( vtkDataSet* )pSource->GetOutput( ) );
 668     slices->SetSource( vtkimagedata );
 669     slices->Update( );
 670     pSource->Delete( );
 671
 672         vtkStructuredPoints   *stPoints = vtkStructuredPoints::New();
 673         stPoints -> GetPointData( )->SetScalars(  slices->GetOutput()->GetPointData()->GetScalars()  );
 674         stPoints -> SetDimensions( sizeIma, sizeIma, 1 );
 675         stPoints -> SetScalarType( vtkimagedata->GetScalarType() );
 676         stPoints -> SetScalarTypeToShort();
 677         stPoints -> Update();
 678
 679         vtkImageChangeInformation  *change = vtkImageChangeInformation ::New();
 680         change    -> SetInput( stPoints );
 681         change    -> Update();    //important
 682
 683         double _range[2];
 684         vtkimagedata->GetScalarRange(_range);
 685    vtkWindowLevelLookupTable *_bwlookup = vtkWindowLevelLookupTable::New( );
 686    _bwlookup->SetHueRange( 0 , 1 );
 687    _bwlookup->SetNumberOfColors( (int)(_range[1] - _range[0] + 1) );
 688    _bwlookup->SetTableRange( _range[0] , _range[1] );
 689    _bwlookup->SetSaturationRange( 0 , 0 );
 690    _bwlookup->SetValueRange( 0 , 1 );
 691    _bwlookup->SetAlphaRange( 1 , 1 );
 692    _bwlookup->Build( );
 693
 694    vtkLookupTable * _collookup = vtkLookupTable::New( );
 695    _collookup->SetNumberOfColors( 256 );
 696    _collookup->SetTableRange( 0 , 255 );
 697    _collookup->Build( );
 698    _collookup->SetTableValue( 0  , 1 , 0 , 0 , 1 );
 699    _collookup->SetTableValue(128 , 0 , 0 , 1 , 1 );
 700    _collookup->SetTableValue(255 , 0 , 1 , 0 , 1 );
 701
 702
 703         vtkMetaImageWriter *writer = vtkMetaImageWriter::New( );
 704         writer->SetInput( stPoints );
 705         writer->SetFileName( "C:/Users/Images/temp_EED/image.mhd" );
 706         writer->SetFileDimensionality( 2 );
 707         writer->Write( );
 708
 709
 710
 711         vtkDataSetMapper *_3DSliceMapper = vtkDataSetMapper::New( );
 712         _3DSliceMapper->SetInput(change->GetOutput( ) );
 713         _3DSliceMapper->SetLookupTable( _bwlookup );
 714         _3DSliceMapper->SetScalarRange( _range );
 715         _3DSliceMapper->ImmediateModeRenderingOn( );
 716
 717         vtkActor *_3DSliceActor = vtkActor::New( );
 718         _3DSliceActor->SetMapper( _3DSliceMapper );
 719
 720  // The usual rendering stuff.
 721   vtkCamera *camera = vtkCamera::New();
 722       camera->SetPosition(1,1,1);
 723       camera->SetFocalPoint(0,0,0);
 724
 725   vtkRenderer *renderer = vtkRenderer::New();
 726   vtkRenderWindow *renWin = vtkRenderWindow::New();
 727     renWin->AddRenderer(renderer);
 728
 729   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
 730     iren->SetRenderWindow(renWin);
 731
 732   renderer->AddActor( _3DSliceActor );
 733       renderer->SetActiveCamera(camera);
 734       renderer->ResetCamera();
 735       renderer->SetBackground(1,1,1);
 736
 737   renWin->SetSize(300,300);
 738
 739   // interact with data
 740   renWin->Render();
 741   iren->Start();
 742
 743
 744
 745
 746
 747 }
 748
 749 //-------------------------------------------------------------------
 750 // EED 25 Janvier 2007 testLoic
 751 void wxVtkMPR3DView::TestLoic2()
 752 {
 753 /*
 754  vtkVectorNorm *vecMagnitude = vtkVectorNorm::New();
 755                 vecMagnitude->SetInput(VtkMainGrid);
 756                 vecMagnitude->NormalizeOff();
 757                 vecMagnitude->Update();
 758 */
 759
 760         vtkPlane *slicePlane = vtkPlane::New();
 761                                 this->_planeWidget->GetPlane( slicePlane );
 762
 763
 764         vtkCutter* sliceCutter = vtkCutter::New();
 765                 vtkImageData *vtkimagedata =  this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
 766                 sliceCutter->SetInput( vtkimagedata );
 767                 sliceCutter->SetCutFunction( slicePlane );
 768
 769 //EED
 770 //        vtkLookupTable *lut = BuildHueWeightBaseMap();
 771
 772
 773         vtkPolyDataMapper *slice = vtkPolyDataMapper::New();
 774                 slice->SetInput( sliceCutter->GetOutput() );
 775         double range[2];
 776 // EED
 777 //                slice->Update();
 778                                   vtkimagedata->GetScalarRange(range);
 779 //                vecMagnitude->GetOutput()->GetScalarRange( range );
 780 //                range[1] *= 0.7; // reduce the upper range by 30%
 781 //                slice->SetScalarRange( range );
 782 //                slice->SetLookupTable( lut );
 783
 784 //       vtkActor *sliceActor = vtkActor::New();
 785 //                sliceActor->SetMapper( slice );
 786
 787         vtkPolyDataMapper *contourMapper = vtkPolyDataMapper::New();
 788                 contourMapper->SetInput( sliceCutter->GetOutput() );
 789                 contourMapper->SetScalarRange( range );
 790   //              contourMapper->SetLookupTable( lut );
 791
 792
 793         vtkActor *contourActor = vtkActor::New();
 794                 contourActor->SetMapper( contourMapper );
 795
 796
 797
 798  // The usual rendering stuff.
 799   vtkCamera *camera = vtkCamera::New();
 800       camera->SetPosition(1,1,1);
 801       camera->SetFocalPoint(0,0,0);
 802
 803   vtkRenderer *renderer = vtkRenderer::New();
 804   vtkRenderWindow *renWin = vtkRenderWindow::New();
 805     renWin->AddRenderer(renderer);
 806
 807   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
 808     iren->SetRenderWindow(renWin);
 809
 810   renderer->AddActor( contourActor );
 811 //  renderer->AddActor(cubeActor);
 812       renderer->SetActiveCamera(camera);
 813       renderer->ResetCamera();
 814       renderer->SetBackground(1,1,1);
 815
 816   renWin->SetSize(300,300);
 817
 818   // interact with data
 819   renWin->Render();
 820   iren->Start();
 821 }
 822
 823
 824