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