wxMaracasSurfaceRenderingManagerDataMhd.cxx

Go to the documentation of this file.

#include "wxMaracasSurfaceRenderingManagerDataMhd.h"

#include "vtkStripper.h"

wxMaracasSurfaceRenderingManagerDataMhd::wxMaracasSurfaceRenderingManagerDataMhd(vtkImageData* imagedata, std::string dataname, vtkRenderWindowInteractor*  interactor)
: wxMaracasSurfaceRenderingManagerData(NULL, dataname)
{

        this->setVtkImageData(imagedata);
        _dataname = dataname;
        
        _maxgreylevel = getMaxLevel(imagedata);
        _prop3D=NULL;                   

        _cubesFilter = vtkMarchingCubes::New(); 
        _cubesFilter->SetInput(this->_imagedata);               
        _cubesFilter->ComputeGradientsOn ();
        _cubesFilter->ComputeScalarsOn ();
        _cubesFilter->SetNumberOfContours( 1 );
        _cleanFilter = vtkCleanPolyData::New();         
        _cleanFilter->SetInput ( _cubesFilter->GetOutput() );

        _dataMapper = vtkPolyDataMapper::New( );
        _dataMapper->ScalarVisibilityOff( );
        _dataMapper->ImmediateModeRenderingOn();
        vtkActor* dataActor = vtkActor::New();

        //if(_boxWidgetS1){
        if(interactor){

                _boxWidgetS1 = vtkBoxWidget::New();
                _boxWidgetS1->SetInteractor( interactor );
                _boxWidgetS1->SetPlaceFactor(1.25); 

                _boxWidgetS1->SetInput( this->_imagedata );
                _boxWidgetS1->PlaceWidget();    
                boxSurfaceObserver* observer = boxSurfaceObserver::New();               

                vtkStripper* striper = vtkStripper::New();
            striper->SetInput( _cleanFilter->GetOutput() );
                //striper->SetInput( _cubesFilter->GetOutput() );
                
                striper->Update();
                _boxWidgetS1->SetInput(striper->GetOutput());
                _boxWidgetS1->PlaceWidget();

                _boxWidgetS1->HandlesOn ();
                _boxWidgetS1->On();

                _tissuePlanes  = vtkPlanes::New();
   
                int x1,x2,y1,y2,z1,z2;
                this->_imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
                _tissuePlanes->SetBounds  (x1,x2,y1,y2,z1,z2);


                
                _boxWidgetS1->GetPlanes( _tissuePlanes );
                
                _tissueClipper = vtkClipPolyData::New();
                _tissueClipper->SetInput( striper->GetOutput() );
                _tissueClipper->SetClipFunction( _tissuePlanes );
                _tissueClipper->InsideOutOn( );
                _dataMapper->SetInput( _tissueClipper->GetOutput() );
                observer->SetPlanes( _tissuePlanes );
                observer->SetActor( dataActor );
                _boxWidgetS1->AddObserver( vtkCommand::InteractionEvent          , observer );


                _boxWidgetS1->HandlesOn ();
                _boxWidgetS1->On();
        }else{          
                _dataMapper->SetInput(_cleanFilter->GetOutput());                               
        }
        
        dataActor->SetMapper(_dataMapper);      
        this->_prop3D = dataActor;

        
        

        this->changeIsoValue(this->_maxgreylevel);      


}


wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd()
{
        _cubesFilter->Delete();
        _cleanFilter->Delete();
        _dataMapper->Delete();
}
        
void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface()
{
        _cubesFilter->Update();    
        _cleanFilter->Update();
        _dataMapper->Update();  
}
void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){     
                
   
    _cubesFilter->SetValue(0,value);            
        _cubesFilter->Update();    
        _cleanFilter->Update();
        _dataMapper->Update();  
        
        
        
}
int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
        return _maxgreylevel;
}

int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){

        int ext[6], i, j, k,max=0;
        img->GetExtent(ext);

        for(i = ext[0]; i < ext[1];i++){
                for(j = ext[2]; j < ext[3];j++){
                        for(k = ext[4]; k < ext[5];k++){
                unsigned short* ptr = (unsigned short*)img->GetScalarPointer(i,j,k);
                                int temp = (int)*ptr;
                                if(temp > max){
                    max = temp;
                                }
                        }
                }
        }
        return max;


}

void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
        _imagedata = imagedata;
}