2 #include "wxMaracasSurfaceRenderingManagerDataMhd.h"
4 #include "vtkStripper.h"
6 wxMaracasSurfaceRenderingManagerDataMhd::wxMaracasSurfaceRenderingManagerDataMhd(vtkImageData* imagedata, std::string dataname, vtkRenderWindowInteractor* interactor)
7 : wxMaracasSurfaceRenderingManagerData(NULL, dataname)
10 this->setVtkImageData(imagedata);
13 _maxgreylevel = getMaxLevel(imagedata);
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() );
24 _dataMapper = vtkPolyDataMapper::New( );
25 _dataMapper->ScalarVisibilityOff( );
26 _dataMapper->ImmediateModeRenderingOn();
27 vtkActor* dataActor = vtkActor::New();
32 _boxWidgetS1 = vtkBoxWidget::New();
33 _boxWidgetS1->SetInteractor( interactor );
34 _boxWidgetS1->SetPlaceFactor(1.25);
36 _boxWidgetS1->SetInput( this->_imagedata );
37 _boxWidgetS1->PlaceWidget();
38 boxSurfaceObserver* observer = boxSurfaceObserver::New();
40 vtkStripper* striper = vtkStripper::New();
41 striper->SetInput( _cleanFilter->GetOutput() );
42 //striper->SetInput( _cubesFilter->GetOutput() );
45 _boxWidgetS1->SetInput(striper->GetOutput());
46 _boxWidgetS1->PlaceWidget();
48 _boxWidgetS1->HandlesOn ();
51 _tissuePlanes = vtkPlanes::New();
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);
59 _boxWidgetS1->GetPlanes( _tissuePlanes );
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 );
71 _boxWidgetS1->HandlesOn ();
74 _dataMapper->SetInput(_cleanFilter->GetOutput());
77 dataActor->SetMapper(_dataMapper);
78 this->_prop3D = dataActor;
83 this->changeIsoValue(this->_maxgreylevel);
89 wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd()
91 _cubesFilter->Delete();
92 _cleanFilter->Delete();
93 _dataMapper->Delete();
96 void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface()
98 _cubesFilter->Update();
99 _cleanFilter->Update();
100 _dataMapper->Update();
103 ** changes the isovalue in a prop3D
105 void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){
108 _cubesFilter->SetValue(0,value);
109 _cubesFilter->Update();
110 _cleanFilter->Update();
111 _dataMapper->Update();
116 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
117 return _maxgreylevel;
121 ** Get's the max grey level of the image
123 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){
125 int ext[6], i, j, k,max=0;
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;
145 ** Sets the VTK image data
148 void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
149 _imagedata = imagedata;