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