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();
31 _boxWidgetS1 = vtkBoxWidget::New();
32 _boxWidgetS1->SetInteractor( interactor );
33 _boxWidgetS1->SetPlaceFactor(1.25);
35 _boxWidgetS1->SetInput( this->_imagedata );
36 _boxWidgetS1->PlaceWidget();
37 boxSurfaceObserver* observer = boxSurfaceObserver::New();
39 _boxWidgetS1->HandlesOn ();
42 vtkStripper* striper = vtkStripper::New();
43 striper->SetInput( _cleanFilter->GetOutput() );
44 //striper->SetInput( _cubesFilter->GetOutput() );
47 _boxWidgetS1->SetInput(striper->GetOutput());
48 _boxWidgetS1->PlaceWidget();
50 _tissuePlanes = vtkPlanes::New();
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);
58 _boxWidgetS1->GetPlanes( _tissuePlanes );
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 );
72 _dataMapper->SetInput(_cleanFilter->GetOutput());
75 dataActor->SetMapper(_dataMapper);
76 this->_prop3D = dataActor;
81 this->changeIsoValue(this->_maxgreylevel);
86 wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd(){
87 _cubesFilter->Delete();
88 _cleanFilter->Delete();
89 _dataMapper->Delete();
92 void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface(){
93 _cubesFilter->Update();
94 _cleanFilter->Update();
95 _dataMapper->Update();
98 ** changes the isovalue in a prop3D
100 void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){
103 _cubesFilter->SetValue(0,value);
104 _cubesFilter->Update();
105 _cleanFilter->Update();
106 _dataMapper->Update();
111 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
112 return _maxgreylevel;
116 ** Get's the max grey level of the image
118 int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){
120 int ext[6], i, j, k,max=0;
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;
140 ** Sets the VTK image data
143 void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
144 _imagedata = imagedata;