lib/maracasVisuLib/src/interface/wxWindows/widgets/vtkPlane2DView.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 "vtkPlane2DView.h"
  27
  28 #include "vtkCellArray.h"
  29
  30 #include "vtkInteractorStyleBaseView2D.h"
  31 #include "vtkPointData.h"
  32
  33
  34 vtkPlane2DView::vtkPlane2DView( wxWindow *parent)
  35 : wxVtk2DBaseView(parent)
  36 {
  37
  38         _backX                                  =       -99999;
  39         _backY                                  =       -99999;
  40         _backZ                                  =       -99999;
  41
  42         _backOrient[0]                  =       -99999;
  43         _backOrient[1]                  =       -99999;
  44         _backOrient[2]                  =       -99999;
  45         _backOrient[3]                  =       -99999;
  46
  47         _active                                 =   true;
  48         _mip_visualization              =       true;
  49         _mip_width                              =       2;
  50
  51     _transform1                         =       vtkTransform::New();
  52     _transform2                         =       vtkTransform::New();
  53         _transform1->Identity();
  54         _transform2->Identity();
  55
  56         _sizeIma                                =       200;
  57
  58         _pSource                                =       NULL;
  59         _3Dslices                               =       NULL;
  60         _stPoints                               =       NULL;
  61         _change                                 =       NULL;
  62         _imageResult                    =       NULL;
  63
  64         // line horizontal
  65     _pts                                        =       NULL;
  66     _lineActor                          =       NULL;
  67         _lineMapper                             =       NULL;
  68         _pd                                             =       NULL;
  69
  70         _interactorstyleplane2D =       NULL;
  71 }
  72 //-------------------------------------------------------------------
  73 vtkPlane2DView::~vtkPlane2DView()
  74 {
  75         ResetPlane();
  76         // Horizontal Line
  77     if (_pts            != NULL)        {       _pts                    -> Delete(); }
  78     if (_lineActor      != NULL)        {       _lineActor              -> Delete(); }
  79     if (_lineMapper     != NULL)        {       _lineMapper             -> Delete(); }
  80     if (_pd                     != NULL)        {       _pd                             -> Delete(); }
  81         _transform1             -> Delete();
  82         _transform2             -> Delete();
  83 }
  84
  85 //-------------------------------------------------------------------
  86 void vtkPlane2DView::ResetBack()
  87 {
  88                 _backX=-1;
  89                 _backY=-1;
  90                 _backZ=-1;
  91                 _backOrient[0]=-1;
  92                 _backOrient[1]=-1;
  93                 _backOrient[2]=-1;
  94                 _backOrient[3]=-1;
  95 }
  96
  97 //-------------------------------------------------------------------
  98 void vtkPlane2DView::ResetPlane()
  99 {
 100         if (_pSource    !=NULL) { _pSource      -> Delete();    }
 101         if (_3Dslices   !=NULL) { _3Dslices     -> Delete();    }
 102         if (_stPoints   !=NULL) { _stPoints     -> Delete();    }
 103         if (_change             !=NULL) { _change       -> Delete();    }
 104 }
 105
 106 //-------------------------------------------------------------------
 107 vtkMPRBaseData *vtkPlane2DView::GetVtkmprbasedata()
 108 {
 109         return (vtkMPRBaseData*)GetVtkBaseData();
 110 }
 111
 112 //-------------------------------------------------------------------
 113 void vtkPlane2DView::SetPSource(int sizeIma){
 114         int dimIma      = sizeIma;
 115
 116         double x = GetVtkmprbasedata()->GetX();
 117         double y = GetVtkmprbasedata()->GetY();
 118         double z = GetVtkmprbasedata()->GetZ();
 119
 120         double spc[3];
 121         vtkImageData *imagedata = GetVtkmprbasedata()->GetImageData();
 122         imagedata->GetSpacing(spc);
 123         x=x*spc[0];
 124         y=y*spc[1];
 125         z=z*spc[2];
 126
 127
 128         vtkMPRBaseData  *mprbasedata = (vtkMPRBaseData*)this->GetVtkBaseData();
 129         vtkTransform    *transform       = mprbasedata->GetTransformOrientation();
 130
 131         double in[3];  // temp
 132         double pA[3];
 133         double pB[3];
 134         in[0]=1;                in[1]=0;                        in[2]=0;
 135         transform->TransformPoint(in,_n);
 136
 137         in[0]=0;        in[1]=dimIma-1;         in[2] = 0;
 138         transform->TransformPoint(in,pA);
 139
 140         in[0]=0;        in[1]=0;                in[2]=dimIma-1;
 141         transform->TransformPoint(in,pB);
 142
 143     _pSource -> SetPoint1( pA   );
 144     _pSource -> SetPoint2( pB   );
 145
 146         _pSource -> SetOrigin( 0        , 0             , 0             );
 147     _pSource -> SetResolution( sizeIma-1 , sizeIma -1 );
 148         _pSource -> Update();
 149         _pSource -> SetCenter( x, y, z );
 150     _pSource -> SetNormal( _n );
 151     _pSource -> Update( );
 152
 153 // EED Borrame
 154 //      transform->Delete();
 155
 156 }
 157
 158 //-------------------------------------------------------------------
 159 void vtkPlane2DView::ExtractPlane()
 160 {
 161         double x = GetVtkmprbasedata()->GetX();
 162         double y = GetVtkmprbasedata()->GetY();
 163         double z = GetVtkmprbasedata()->GetZ();
 164
 165         double spc[3];
 166         vtkImageData *imagedata = GetVtkmprbasedata()->GetImageData();
 167         imagedata->GetSpacing(spc);
 168         x=x*spc[0];
 169         y=y*spc[1];
 170         z=z*spc[2];
 171
 172         vtkTransform *transform = GetVtkmprbasedata()->GetTransformOrientation();
 173         double orientation[4];
 174         transform->GetOrientationWXYZ(orientation);
 175
 176         bool okOrientation=true;
 177         if ((orientation[0]!=_backOrient[0]) || (orientation[1]!=_backOrient[1]) ||
 178                 (orientation[2]!=_backOrient[2]) || (orientation[3]!=_backOrient[3]))
 179         {
 180                 okOrientation=false;
 181         }
 182
 183         bool okPosicion=true;
 184         if ( (x!=_backX) || (y!=_backY) || (z!=_backZ) )
 185         {
 186                 okPosicion=false;
 187         }
 188
 189         if ((okPosicion==false) || (okOrientation==false) ) {
 190                 if (_active==true){
 191
 192                         bool ok = false;
 193
 194                         if ( _mip_visualization==true )
 195                         {
 196                                 if (_interactorstyleplane2D!=NULL)
 197                                 {
 198                                         if (_interactorstyleplane2D->GetStateRotate()==false )
 199                                         {
 200                                                 ok=true;
 201                                         }
 202                                 }
 203                         }
 204 //                      ok=true;
 205
 206                         if (ok==true)
 207                         {
 208                                 Extract_MIP_PlaneVTK();
 209                         } else {
 210                                 Extract_One_PlaneVTK();
 211                         } // ok
 212
 213                 } // active
 214                 _backX=x;
 215                 _backY=y;
 216                 _backZ=z;
 217                 _backOrient[0]=orientation[0];
 218                 _backOrient[1]=orientation[1];
 219                 _backOrient[2]=orientation[2];
 220                 _backOrient[3]=orientation[3];
 221         } //okPosition okOrientation
 222 }
 223
 224 //-------------------------------------------------------------------
 225 void vtkPlane2DView::Extract_One_PlaneVTK()
 226 {
 227         vtkImageData *imagedata = GetVtkmprbasedata()->GetImageData();
 228         SetPSource(_sizeIma);
 229
 230 //EED 2017-01-01 Migration VTK7
 231 #if VTK_MAJOR_VERSION <= 5
 232         _3Dslices -> SetInput( ( vtkDataSet* )_pSource->GetOutput( ) );
 233         _3Dslices -> SetSource( imagedata );
 234 #else
 235         _3Dslices -> SetInputData( ( vtkDataSet* )_pSource->GetOutput( ) );
 236         _3Dslices -> SetSourceData( imagedata );
 237 #endif
 238
 239         _3Dslices -> Update( );
 240         _stPoints -> GetPointData( )->SetScalars(  _3Dslices->GetOutput()->GetPointData()->GetScalars()  );
 241         _stPoints -> SetDimensions( _sizeIma, _sizeIma, 1 );
 242
 243 //EED 2017-01-01 Migration VTK7
 244 #if VTK_MAJOR_VERSION <= 5
 245         _stPoints -> SetScalarType( imagedata->GetScalarType() );
 246         _stPoints -> SetScalarTypeToShort();
 247         _stPoints -> Update();
 248 //      _change   -> SetInput( _stPoints );
 249 //      _change   -> Update();    //important
 250         _imageViewer2XYZ->GetVtkImageViewer2()->SetInput ( _stPoints );
 251
 252 #else
 253         _stPoints -> SetScalarType( imagedata->GetScalarType() );
 254         _stPoints -> SetScalarTypeToShort();
 255         _imageViewer2XYZ->GetVtkImageViewer2()->SetInputData( _stPoints );
 256
 257 #endif
 258
 259
 260 //      _imageViewer2XYZ->GetVtkImageViewer2()->SetInput ( _change->GetOutput() );
 261         //      vtkImageActor *imageActor = _imageViewer2XYZ->GetVtkImageViewer2()->GetImageActor();
 262 }
 263
 264 //-------------------------------------------------------------------
 265 void vtkPlane2DView::Extract_MIP_PlaneVTK( /*double heightDefinition*/ )
 266 {
 267         int mipWidth;
 268         double sp;
 269         int sizeWidth = (_mip_width*2)+1 ;
 270         int deltaPixel;
 271         int iWidth,itmp,tmpSizeWith;
 272         double spc[3];
 273         vtkImageData *imagedata = GetVtkmprbasedata()->GetImageData();
 274         SetPSource(_sizeIma);
 275         imagedata->GetSpacing(spc);
 276         bool heightDefinition=false;
 277         if (_mip_width<3)
 278         {
 279                 heightDefinition=true;
 280         }
 281         if (heightDefinition==true)
 282         {
 283                 mipWidth        =       _mip_width;
 284                 sp                      =       spc[0];
 285                 deltaPixel      =       1;
 286         } else {
 287                 mipWidth        =       2;
 288                 tmpSizeWith =   (mipWidth*2) + 1;
 289                 sp                      =       (spc[0]*sizeWidth)/tmpSizeWith;
 290                 sizeWidth       =       tmpSizeWith;
 291                 deltaPixel      =       4;
 292         }
 293         std::vector< vtkProbeFilter* > slicesLST;
 294         _pSource->Push( -mipWidth * sp );
 295         _pSource->Update();
 296         for ( iWidth=0 ; iWidth<sizeWidth ; iWidth++ )
 297         {
 298                 vtkProbeFilter *slice = vtkProbeFilter::New();
 299 //EED 2017-01-01 Migration VTK7
 300 #if VTK_MAJOR_VERSION <= 5
 301                 slice -> SetInput( ( vtkDataSet* )_pSource->GetOutput( ) );
 302                 slice -> SetSource( imagedata );
 303 #else
 304                 slice -> SetInputData( ( vtkDataSet* )_pSource->GetOutput( ) );
 305                 slice -> SetSourceData( imagedata );
 306 #endif
 307                 slice -> Update( );
 308                 slicesLST.push_back( slice );
 309                 _pSource->Push( sp );
 310         }
 311         if (_imageResult ==NULL)
 312         {
 313                 _imageResult = vtkImageData::New();
 314                 _imageResult -> SetDimensions(_sizeIma,_sizeIma,1);
 315                 _imageResult -> SetSpacing(1,1,1);
 316                 _imageResult -> SetExtent(0,_sizeIma-1,0,_sizeIma-1,0,0);
 317 //EED 2017-01-01 Migration VTK7
 318 #if VTK_MAJOR_VERSION <= 5
 319                 _imageResult -> SetScalarType( imagedata->GetScalarType() );
 320                 _imageResult -> SetWholeExtent(0,_sizeIma-1,0,_sizeIma-1,0,0);
 321                 _imageResult -> AllocateScalars();
 322                 _imageResult -> Update();
 323 #else
 324                 _imageResult -> AllocateScalars(imagedata->GetScalarType(),1);
 325 #endif
 326
 327         }
 328         unsigned short *pTemp;
 329         unsigned short *pResult;
 330         pResult = (unsigned short*)_imageResult->GetScalarPointer( 0 , 0 , 0 );
 331         int iPixels , sizePixels = _sizeIma*_sizeIma;
 332         for(iPixels=0 ; iPixels<sizePixels ; iPixels=iPixels+deltaPixel)
 333         {
 334                 pTemp = (unsigned short*)slicesLST[0]->GetOutput()->GetPointData()->GetScalars()->GetVoidPointer(0);
 335                 pResult = (unsigned short*)_imageResult->GetScalarPointer( 0 , 0 , 0 );
 336                 pResult[iPixels] = pTemp[iPixels];
 337                 for (iWidth=1;iWidth<sizeWidth;iWidth++)
 338                 {
 339                         pTemp = (unsigned short*)slicesLST[iWidth]->GetOutput()->GetPointData()->GetScalars()->GetVoidPointer(0);
 340
 341                         if (pResult[iPixels]< pTemp[iPixels])
 342                         {
 343                                 pResult[iPixels] = pTemp[iPixels];
 344                         }
 345                 }
 346                 if (deltaPixel!=1)
 347                 {
 348                         for (itmp=1;itmp<deltaPixel;itmp++)
 349                         {
 350                                 pResult[iPixels+itmp] = pResult[iPixels];
 351                         }
 352                 }
 353         }
 354         for (iWidth=0;iWidth<sizeWidth;iWidth++)
 355         {
 356                 slicesLST[iWidth]->Delete();
 357         }
 358         _imageResult->Modified();
 359 //EED 2017-01-01 Migration VTK7
 360 #if VTK_MAJOR_VERSION <= 5
 361         _imageViewer2XYZ->GetVtkImageViewer2()->SetInput( _imageResult );
 362 #else
 363         _imageViewer2XYZ->GetVtkImageViewer2()->SetInputData( _imageResult );
 364 #endif
 365 }
 366
 367 //-------------------------------------------------------------------
 368 void vtkPlane2DView::Configure( )
 369 {
 370         wxVtk2DBaseView::Configure(false);
 371
 372         HorizontalLine();
 373 // Borrame
 374 //      CircleLine();
 375
 376         _pSource         = vtkPlaneSource::New( );
 377     _3Dslices    = vtkProbeFilter::New( ) ;
 378         _stPoints        = vtkStructuredPoints::New( );
 379         _change          = vtkImageChangeInformation::New();
 380
 381         wxVTKRenderWindowInteractor *iren               = GetWxVTKRenderWindowInteractor();
 382         ExtractPlane();
 383 //      SetActive(false);
 384         _imageViewer2XYZ -> GetVtkImageViewer2()        -> SetupInteractor ( iren );
 385
 386         //////////////
 387         //RaC 04-2010 Look for the description of new changes in vtkInteractorStyleBaseView2D constructor.
 388         //
 389         // Previous version
 390         // SetInteractorStyleImage( vtkInteractorStyleBaseView2D::New() );
 391         //
 392         // Actual version
 393         vtkInteractorStyleBaseView2D *style2D = vtkInteractorStyleBaseView2D::New();
 394         manualInteractorWindowLevel *_manualinteractorwindowlevel= new manualInteractorWindowLevel();
 395         style2D->SetInteractorWindowLevel( _manualinteractorwindowlevel );
 396         vtkInteractorScrollZ *_vtkInteractorScrollZ = new vtkInteractorScrollZ();
 397         style2D->SetInteractorScrollZ(_vtkInteractorScrollZ);
 398         SetInteractorStyleImage( style2D );
 399
 400         // RaC
 401         //////////////
 402
 403         _interactorstyleplane2D = new vtkInteractorStylePlane2D();
 404         ((vtkInteractorStyleBaseView*)GetInteractorStyleBaseView())->AddInteractorStyleMaracas( _interactorstyleplane2D );
 405
 406         vtkCamera *camera =_imageViewer2XYZ->GetVtkImageViewer2()->GetRenderer()->GetActiveCamera();
 407
 408         camera->SetViewUp               (       0                       ,       1                               ,       0       );
 409         camera->SetPosition             ((0+_sizeIma)/2 , (0+_sizeIma)/2        , 10000 );
 410
 411         camera->SetFocalPoint   ((0+_sizeIma)/2 , (0+_sizeIma)/2        ,       0       );
 412         camera->SetClippingRange(0.01, 100000);
 413         camera->ComputeViewPlaneNormal();
 414         camera->SetParallelScale( _sizeIma/3.0 );
 415
 416         // text information over the graphic window
 417         _vtkInfoTextImage                                               = new vtkInfoTextImage();
 418         _vtkInfoTextImageInteractorPlane2D              = new vtkInfoTextImageInteractorPlane2D();
 419         _vtkInfoTextImage->SetWxVtk2DBaseView(this);
 420         _vtkInfoTextImage->SetMarImageData(  GetVtkmprbasedata()->GetMarImageData() );
 421         _vtkInfoTextImageInteractorPlane2D->SetModelVtkInfoTextImage(_vtkInfoTextImage);
 422         _vtkInfoTextImage->Configure();
 423         ((vtkInteractorStyleBaseView*)this->GetInteractorStyleBaseView())->AddInteractorStyleMaracas(_vtkInfoTextImageInteractorPlane2D);
 424 }
 425
 426 //-------------------------------------------------------------------
 427 void vtkPlane2DView::HorizontalLine()
 428 {
 429 // Axe Horizontal
 430         _pts = vtkPoints::New();
 431         _pts->SetNumberOfPoints(2);
 432         _pts->SetPoint(0, -1000 , -1000 , -1000 );
 433         _pts->SetPoint(1,  1000 ,  1000 ,  1000 );
 434         vtkCellArray *lines = vtkCellArray::New();
 435         lines->InsertNextCell(2);
 436         lines->InsertCellPoint(0);
 437         lines->InsertCellPoint(1);
 438         _pd = vtkPolyData::New();
 439     _pd->SetPoints( _pts );
 440     _pd->SetLines( lines );
 441         lines->Delete();  //do not delete lines ??
 442         _lineActor                                      =       vtkActor::New();
 443     _lineMapper                                 =       vtkPolyDataMapper::New();
 444 //EED 2017-01-01 Migration VTK7
 445 #if VTK_MAJOR_VERSION <= 5
 446         _lineMapper->SetInput(_pd);
 447 #else
 448         _lineMapper->SetInputData(_pd);
 449 #endif
 450         _lineMapper->ImmediateModeRenderingOn();
 451         _lineActor->SetMapper(_lineMapper);
 452         _lineActor->GetProperty()->BackfaceCullingOn();
 453         _lineActor->GetProperty()->SetDiffuseColor(0,0,1);
 454         _lineActor->GetProperty()->SetLineWidth(2);
 455         _lineActor->GetProperty()->SetOpacity(0);
 456     _imageViewer2XYZ->GetVtkImageViewer2()->GetRenderer()->AddActor( _lineActor );
 457 }
 458
 459 //-------------------------------------------------------------------
 460 void vtkPlane2DView::RotationStart()
 461 {
 462         vtkMPRBaseData  *mprbasedata = (vtkMPRBaseData*)this->GetVtkBaseData();
 463         vtkTransform    *transform       = mprbasedata->GetTransformOrientation();
 464         _transform1->SetMatrix( transform->GetMatrix() );
 465 }
 466
 467 //-------------------------------------------------------------------
 468 void vtkPlane2DView::RotationDrag(double vx, double vy, bool ok_v, bool ok_ang)
 469 {
 470         if (ok_ang==false)
 471         {
 472                 _ang =sqrt( vx*vx + vy*vy ) / 1.5;
 473         }
 474
 475         if (ok_v==false){
 476                 _vxb=-vy;
 477                 _vyb=vx;
 478         }
 479         _transform2->Identity();
 480         _transform2->RotateWXYZ(_ang,0,_vxb,_vyb);
 481
 482         vtkMPRBaseData  *mprbasedata = (vtkMPRBaseData*)this->GetVtkBaseData();
 483         vtkTransform *transform          =      vtkTransform::New();
 484         transform->Identity();
 485         transform->Concatenate(_transform1);
 486         transform->Concatenate(_transform2);
 487         mprbasedata->InitTransformOrientation(transform);
 488         transform->Delete();
 489
 490         // Refresh Horizontal Line
 491         _pts->SetPoint( 0 , (_sizeIma/2) - _vxb*2       , (_sizeIma/2) - _vyb*2 , 1 );
 492         _pts->SetPoint( 1 , (_sizeIma/2) + _vxb*2       , (_sizeIma/2) + _vyb*2 , 1 );
 493 //      RefreshCircleLine();
 494 }
 495
 496 //-------------------------------------------------------------------
 497 void vtkPlane2DView::Refresh(  )
 498 {
 499         ExtractPlane();
 500         UpdateColorWindowLevel();
 501         wxVtkBaseView::Refresh();
 502 }
 503 //-------------------------------------------------------------------
 504 void vtkPlane2DView::SetImgSize( int imgSize )
 505 {
 506         _sizeIma = imgSize;
 507 }
 508 //-------------------------------------------------------------------
 509 int vtkPlane2DView::GetImgSize()
 510 {
 511         return _sizeIma;
 512 }
 513 //-------------------------------------------------------------------
 514 int     vtkPlane2DView::GetActualSlice()  // virtual
 515 {
 516         _cx = GetVtkmprbasedata()->GetX();
 517         _cy = GetVtkmprbasedata()->GetY();
 518         _cz = GetVtkmprbasedata()->GetZ();
 519         return 0;
 520 }
 521
 522 //-------------------------------------------------------------------
 523
 524 bool vtkPlane2DView::GetMipVisualization()
 525 {
 526         return _mip_visualization;
 527 }
 528
 529 //-------------------------------------------------------------------
 530 int vtkPlane2DView::GetMipWidth()
 531 {
 532         return _mip_width;
 533 }
 534
 535 //-------------------------------------------------------------------
 536 void vtkPlane2DView::SetActualSlice(int slice)  // Virtual
 537 {
 538         double dir=(double)slice/3.0;
 539         GetVtkmprbasedata()->SetX( _cx + (_n[0]*dir) );
 540         GetVtkmprbasedata()->SetY( _cy + (_n[1]*dir) );
 541         GetVtkmprbasedata()->SetZ( _cz + (_n[2]*dir) );
 542 }
 543 //---------------------------------------------------------------------------
 544 vtkInteractorStylePlane2D       *vtkPlane2DView::GetInteractorstyleplane2D()
 545 {
 546         return this->_interactorstyleplane2D;
 547 }
 548 //---------------------------------------------------------------------------
 549 void vtkPlane2DView::SetActive(bool active)
 550 {
 551         _active = active;
 552 }
 553
 554 //---------------------------------------------------------------------------
 555
 556 void vtkPlane2DView::SetMipVisualization(bool ok)
 557 {
 558         _mip_visualization=ok;
 559 }
 560
 561 //---------------------------------------------------------------------------
 562
 563 void vtkPlane2DView::SetMipWidth(int value)
 564 {
 565         _mip_width=value;
 566 }
 567
 568 //---------------------------------------------------------------------------
 569 void vtkPlane2DView::SetVisibleLine(bool ok)
 570 {
 571         double opacity;
 572         if (ok==true)
 573         {
 574                 opacity=1;
 575         } else {
 576                 opacity=0;
 577         }
 578         _lineActor->GetProperty()->SetOpacity(opacity);
 579 }
 580
 581 // ----------------------------------------------------------------------------
 582 void vtkPlane2DView::TransfromCoordViewWorld2(double &X, double &Y, double &Z)
 583 {
 584         double spc[3];
 585         GetVtkmprbasedata()->GetImageData()->GetSpacing(spc);
 586         double xx = X;
 587         double yy = Y;
 588         double zz = 0;
 589         TransFromCoordScreenToWorld(xx,yy,zz);
 590         vtkTransform *transf1 = vtkTransform::New();
 591         transf1->Identity();
 592         vtkTransform *transf2 = GetVtkmprbasedata()->GetTransformOrientation();
 593         transf1->Concatenate(transf2->GetMatrix());
 594         double in[4], out[4];
 595         double center = GetImgSize() / 2;
 596         in[0] = 0;
 597         in[1] = xx - center;
 598         in[2] = yy - center;
 599         in[3] = 0;
 600         transf1->MultiplyPoint(in,out);
 601         transf1->Delete();
 602         X = out[0] + GetVtkmprbasedata()->GetX() ;
 603         Y = out[1] + GetVtkmprbasedata()->GetY() ;
 604         Z = out[2] + GetVtkmprbasedata()->GetZ() ;
 605 }
 606