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 _boxWidgetS1->HandlesOn ();
43 vtkStripper* striper = vtkStripper::New();
44 striper->SetInput( _cleanFilter->GetOutput() );
45 //striper->SetInput( _cubesFilter->GetOutput() );
48 _boxWidgetS1->SetInput(striper->GetOutput());
49 _boxWidgetS1->PlaceWidget();
51 _boxWidgetS1->HandlesOn ();
54 _tissuePlanes = vtkPlanes::New();
56 int x1,x2,y1,y2,z1,z2;
57 this->_imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
58 _tissuePlanes->SetBounds (x1,x2,y1,y2,z1,z2);
62 _boxWidgetS1->GetPlanes( _tissuePlanes );
64 _tissueClipper = vtkClipPolyData::New();
65 _tissueClipper->SetInput( striper->GetOutput() );
66 _tissueClipper->SetClipFunction( _tissuePlanes );
67 _tissueClipper->InsideOutOn( );
68 _dataMapper->SetInput( _tissueClipper->GetOutput() );
69 observer->SetPlanes( _tissuePlanes );
70 observer->SetActor( dataActor );
71 _boxWidgetS1->AddObserver( vtkCommand::InteractionEvent , observer );
76 _dataMapper->SetInput(_cleanFilter->GetOutput());
79 dataActor->SetMapper(_dataMapper);
80 this->_prop3D = dataActor;
85 this->changeIsoValue(this->_maxgreylevel);
91 wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd()
93 _cubesFilter->Delete();
94 _cleanFilter->Delete();
95 _dataMapper->Delete();
98 void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface()
100 _cubesFilter->Update();
101 _cleanFilter->Update();
102 _dataMapper->Update();
105 ** changes the isovalue in a prop3D
107 void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){
110 _cubesFilter->SetValue(0,value);
111 _cubesFilter->Update();
112 _cleanFilter->Update();
113 _dataMapper->Update();
118 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
119 return _maxgreylevel;
123 ** Get's the max grey level of the image
125 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){
127 int ext[6], i, j, k,max=0;
130 for(i = ext[0]; i < ext[1];i++){
131 for(j = ext[2]; j < ext[3];j++){
132 for(k = ext[4]; k < ext[5];k++){
133 unsigned short* ptr = (unsigned short*)img->GetScalarPointer(i,j,k);
134 int temp = (int)*ptr;
147 ** Sets the VTK image data
150 void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
151 _imagedata = imagedata;