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

   1
   2 #include "wxMaracasSurfaceRenderingManagerDataMhd.h"
   3
   4 #include "vtkStripper.h"
   5
   6 wxMaracasSurfaceRenderingManagerDataMhd::wxMaracasSurfaceRenderingManagerDataMhd(vtkImageData* imagedata, std::string dataname, vtkRenderWindowInteractor*  interactor)
   7 : wxMaracasSurfaceRenderingManagerData(NULL, dataname)
   8 {
   9
  10         this->setVtkImageData(imagedata);
  11         _dataname = dataname;
  12
  13         _maxgreylevel = getMaxLevel(imagedata);
  14         _prop3D=NULL;
  15
  16         _cubesFilter = vtkMarchingCubes::New();
  17         _cubesFilter->SetInput(this->_imagedata);
  18         _cubesFilter->ComputeGradientsOn ();
  19         _cubesFilter->ComputeScalarsOn ();
  20         _cubesFilter->SetNumberOfContours( 1 );
  21         _cleanFilter = vtkCleanPolyData::New();
  22         _cleanFilter->SetInput ( _cubesFilter->GetOutput() );
  23
  24         _dataMapper = vtkPolyDataMapper::New( );
  25         _dataMapper->ScalarVisibilityOff( );
  26         _dataMapper->ImmediateModeRenderingOn();
  27         vtkActor* dataActor = vtkActor::New();
  28
  29         if(interactor){
  30
  31                 _boxWidgetS1 = vtkBoxWidget::New();
  32                 _boxWidgetS1->SetInteractor( interactor );
  33                 _boxWidgetS1->SetPlaceFactor(1.25);
  34
  35                 _boxWidgetS1->SetInput( this->_imagedata );
  36                 _boxWidgetS1->PlaceWidget();
  37                 boxSurfaceObserver* observer = boxSurfaceObserver::New();
  38
  39                 _boxWidgetS1->HandlesOn ();
  40                 _boxWidgetS1->On();
  41
  42                 vtkStripper* striper = vtkStripper::New();
  43             striper->SetInput( _cleanFilter->GetOutput() );
  44                 //striper->SetInput( _cubesFilter->GetOutput() );
  45
  46                 striper->Update();
  47                 _boxWidgetS1->SetInput(striper->GetOutput());
  48                 _boxWidgetS1->PlaceWidget();
  49
  50                 _tissuePlanes  = vtkPlanes::New();
  51
  52                 int x1,x2,y1,y2,z1,z2;
  53                 this->_imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
  54                 _tissuePlanes->SetBounds  (x1,x2,y1,y2,z1,z2);
  55
  56
  57
  58                 _boxWidgetS1->GetPlanes( _tissuePlanes );
  59
  60                 _tissueClipper = vtkClipPolyData::New();
  61                 _tissueClipper->SetInput( striper->GetOutput() );
  62                 _tissueClipper->SetClipFunction( _tissuePlanes );
  63                 _tissueClipper->InsideOutOn( );
  64                 _dataMapper->SetInput( _tissueClipper->GetOutput() );
  65                 observer->SetPlanes( _tissuePlanes );
  66                 observer->SetActor( dataActor );
  67                 _boxWidgetS1->AddObserver( vtkCommand::InteractionEvent          , observer );
  68
  69
  70
  71         }else{
  72                 _dataMapper->SetInput(_cleanFilter->GetOutput());
  73         }
  74
  75         dataActor->SetMapper(_dataMapper);
  76         this->_prop3D = dataActor;
  77
  78
  79
  80
  81         this->changeIsoValue(this->_maxgreylevel);
  82
  83
  84 }
  85
  86 wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd(){
  87         _cubesFilter->Delete();
  88         _cleanFilter->Delete();
  89         _dataMapper->Delete();
  90 }
  91
  92 void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface(){
  93         _cubesFilter->Update();
  94         _cleanFilter->Update();
  95         _dataMapper->Update();
  96 }
  97 /**
  98 **      changes the isovalue in a prop3D
  99 **/
 100 void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){
 101
 102
 103     _cubesFilter->SetValue(0,value);
 104         _cubesFilter->Update();
 105         _cleanFilter->Update();
 106         _dataMapper->Update();
 107
 108
 109
 110 }
 111 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
 112         return _maxgreylevel;
 113 }
 114
 115 /**
 116         ** Get's the max grey level of the image
 117         **/
 118 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){
 119
 120         int ext[6], i, j, k,max=0;
 121         img->GetExtent(ext);
 122
 123         for(i = ext[0]; i < ext[1];i++){
 124                 for(j = ext[2]; j < ext[3];j++){
 125                         for(k = ext[4]; k < ext[5];k++){
 126                 unsigned short* ptr = (unsigned short*)img->GetScalarPointer(i,j,k);
 127                                 int temp = (int)*ptr;
 128                                 if(temp > max){
 129                     max = temp;
 130                                 }
 131                         }
 132                 }
 133         }
 134         return max;
 135
 136
 137 }
 138
 139 /**
 140 **      Sets the VTK image data
 141 **/
 142
 143 void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
 144         _imagedata = imagedata;
 145 }
 146