TreeExtraction_MH_JFC Class Reference

#include <wxEmptyPanelWidget_2.h>

Collaboration diagram for TreeExtraction_MH_JFC:

[legend]

List of all members.

Public Member Functions
virtual void	SetInput (marExperiment *input)
void	SetBranchLevel (int branchLevel)
marExperiment *	GetInput ()
vtkImageData *	GetVolume ()
void	Update ()
void	SetCleanTreeLevel (int cleanLevel)
Static Public Member Functions
static TreeExtraction_MH_JFC *	New ()
Protected Member Functions
	TreeExtraction_MH_JFC ()
	~TreeExtraction_MH_JFC ()
void	Execute (int count, int sens=1, double *vit=NULL, int cleanLevel=0)
Private Member Functions
void	Reset ()
void	InitImagaDataResult ()
void	BruleVolume (vtkImageData *imagedata, int x, int y, int z, int r, int value)
void	BruleVolume (vtkImageData *imagedata, vtkPolyData *inputAxisRadio, int value)
void	FindVit (marExperiment *newExperiment, int actualAxis, double *pointBif, double *vit)
Private Attributes
double	_cleanLevel
int	_branchLevel
carotidaBifurcacion *	_prgov
vtkImageData *	_imagedataResult
marExperiment *	_marExperiment

---

Detailed Description

Definition at line 24 of file wxEmptyPanelWidget_2.h.

---

Constructor & Destructor Documentation

TreeExtraction_MH_JFC::TreeExtraction_MH_JFC

(

)

[protected]

Definition at line 55 of file wxEmptyPanelWidget_2.cxx.

References _branchLevel, _cleanLevel, and _imagedataResult.

                                            {
        _branchLevel            = 2;
        _imagedataResult        = NULL;
        _cleanLevel                     = 400;
}

TreeExtraction_MH_JFC::~TreeExtraction_MH_JFC

(

)

[protected]

Definition at line 61 of file wxEmptyPanelWidget_2.cxx.

References Reset().

                                             {
        Reset();
}

Here is the call graph for this function:

---

Member Function Documentation

void TreeExtraction_MH_JFC::BruleVolume	(	vtkImageData *	imagedata,
		vtkPolyData *	inputAxisRadio,
		int	value
	)			[private]

Definition at line 113 of file wxEmptyPanelWidget_2.cxx.

References BruleVolume().

{
        int x,y,z,r;
        double point[3],radio;
        double fractpart, intpart;
        int i,maxPoints=inputAxisRadio->GetPoints()->GetNumberOfPoints();
        vtkDoubleArray *allRadios = (vtkDoubleArray*)(inputAxisRadio->GetPointData()->GetScalars("radio"));
        for (i=0; i<maxPoints; i++){    
                inputAxisRadio->GetPoints()->GetPoint(i,point);
                radio=allRadios->GetValue(i);
                x = (int)(point[0]);
                y = (int)(point[1]);    
                z = (int)(point[2]);
                fractpart = modf (radio , &intpart);
                if (fractpart>=0.5){
                 intpart++;
                }
                r = (int)( intpart );
                BruleVolume(imagedata,x,y,z,r,value);
        }
}

Here is the call graph for this function:

void TreeExtraction_MH_JFC::BruleVolume	(	vtkImageData *	imagedata,
		int	x,
		int	y,
		int	z,
		int	r,
		int	value
	)			[private]

Definition at line 76 of file wxEmptyPanelWidget_2.cxx.

Referenced by BruleVolume(), and Execute().

{
        int ext[6];
        int i, j, k, ii,jj,kk;  
        imagedata->GetExtent(ext);
        unsigned short *ptr;
        int x1=x-r;
        int x2=x+r;
        int y1=y-r;
        int y2=y+r;
        int z1=z-r;
        int z2=z+r;
        double tmpsqrt;
        for(i=x1;i<=x2;i++){
                for(j=y1;j<=y2;j++){
                        for(k=z1;k<=z2;k++){
                                ptr=(unsigned short *)  imagedata->GetScalarPointer(i,j,k);
                                ii=i-x;
                                jj=j-y;
                                kk=k-z;
                                tmpsqrt=ii*ii + jj*jj + kk*kk;
                                if ((ptr!=NULL) && ( sqrt(tmpsqrt)<=r))
                                {
                                        if (value==0){
                                                *ptr=value;
                                        } else {
                                                if (*ptr==0)  {
                                                        *ptr=value;
                                                }
                                        }
                                } 
                        } // for z
                }  // for y
        } // fox x
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::Execute	(	int	count,
		int	sens = 1,
		double *	vit = NULL,
		int	cleanLevel = 0
	)			[protected]

Definition at line 173 of file wxEmptyPanelWidget_2.cxx.

References _cleanLevel, _imagedataResult, _marExperiment, _prgov, BruleVolume(), marExperiment::extractVascularTree(), FindVit(), marExperiment::getAxis(), marAxis::GetAxisData(), marExperiment::getDynData(), marExperiment::getNumberOfAxes(), marDynData::getVolume(), New(), carotidaBifurcacion::SetHumbral(), carotidaBifurcacion::SetInput(), carotidaBifurcacion::SetInput2(), carotidaBifurcacion::SetMaxPropMasa(), carotidaBifurcacion::SetMaxPropRadio(), carotidaBifurcacion::SetMinPropMasa(), and marExperiment::setStartPoint().

Referenced by Update().

{
        //------------Arbre Extraction --------------
        // Call Maracas axis extraction
        count--;
/*
time_t start,end;
double dif;
time (&start);
*/
                int idActualAxis_Before = _marExperiment->getNumberOfAxes();
                _marExperiment->extractVascularTree( sens , vit, cleanLevel );
/*
time (&end);
dif = difftime (end,start);
FILE *ff;
ff=fopen("c:/temp/MaracasTIME.txt","a+");
fprintf(ff,"\n TreeExtraction_MH_JFC::Execute %d \n",count);
fprintf(ff,"evt %.2lf \n",dif);
fclose(ff);
*/
                int idActualAxis = _marExperiment->getNumberOfAxes();
                // Bifurcation Detection JFC
                marAxis *maraxis = _marExperiment->getAxis( );
                if ((maraxis!=NULL) && (idActualAxis_Before!=idActualAxis)) {
                        vtkImageData *imagedata                 = _marExperiment->getDynData( )->getVolume( )->castVtk();
                        vtkPolyData  *inputAxisRadio    = maraxis->GetAxisData();

                        if (count>0){
/*
time (&start);
*/
                                _prgov = carotidaBifurcacion::New();
                                _prgov->SetMaxPropRadio(5);
                                _prgov->SetHumbral(0.45);
                                _prgov->SetMaxPropMasa(10);
                                _prgov->SetMinPropMasa(0.01);
                                _prgov->SetInput( inputAxisRadio );
                                _prgov->SetInput2( imagedata );
                                _prgov->Update();
/*
time (&end);
dif = difftime (end,start);
ff=fopen("c:/temp/MaracasTIME.txt","a+");
fprintf(ff,"\tcB  %.2lf \n",dif);
fclose(ff);
*/
                        }

                        // Brule Volume
                        BruleVolume(imagedata            , inputAxisRadio , 0                    );
                        BruleVolume(_imagedataResult , inputAxisRadio , idActualAxis );
                        inputAxisRadio->Delete();


                        if (count>0){
                                // LOOP for branch extraction
                                vtkPolyData *bifurcations =  _prgov->GetOutput();
                                int i,maxPoints = bifurcations->GetPoints()->GetNumberOfPoints();
                                int x,y,z;
                                double point[3], vit[3];
                                for (i=0; i<maxPoints; i++){
                                        bifurcations->GetPoints()->GetPoint(i,point);
                                        x = (int)(point[0]);
                                        y = (int)(point[1]);
                                        z = (int)(point[2]);
                                        FindVit( _marExperiment , idActualAxis , point , vit );
                                        _marExperiment->setStartPoint( x , y , z );
                                        //Branch extraction;
                                        Execute(count, 0 , vit , (int)_cleanLevel );
                                }
                                bifurcations->Delete();
//                              _prgov->Delete();
//                              _prgov = NULL;
                        } // of count>0

                } // maraxi=NULL
}

Here is the call graph for this function:

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::FindVit	(	marExperiment *	newExperiment,
		int	actualAxis,
		double *	pointBif,
		double *	vit
	)			[private]

Definition at line 138 of file wxEmptyPanelWidget_2.cxx.

References marExperiment::getAxis(), marAxis::GetAxisData(), marExperiment::getNumberOfAxes(), min, and marExperiment::setAxis().

Referenced by Execute().

{
        double dist,min=-1;
        double point[3];
        double dx,dy,dz;
        int jPoint,maxPoints;
        int iAxis,maxAxis;

        maxAxis=newExperiment->getNumberOfAxes();
        for (iAxis=0;iAxis<maxAxis;iAxis++){
                if (iAxis!=actualAxis){
                        newExperiment->setAxis(iAxis);
                        marAxis *maraxis = newExperiment->getAxis( );
                        if (maraxis!=NULL) {
                                vtkPolyData  *axisVtk   = maraxis->GetAxisData();
                                maxPoints=axisVtk->GetPoints()->GetNumberOfPoints();
                                for (jPoint=0; jPoint<maxPoints; jPoint++){     
                                        axisVtk->GetPoints()->GetPoint(jPoint,point);
                                        dx = pointBif[0]-point[0];
                                        dy = pointBif[1]-point[1];
                                        dz = pointBif[2]-point[2];
                                        dist = sqrt( dx*dx + dy*dy + dz*dz );
                                        if ((min==-1) || (dist<min)) {
                                                min=dist;
                                                vit[0] = dx;
                                                vit[1] = dy;
                                                vit[2] = dz;
                                        } //if 
                                } // for
                        } // if maraxis!=NULL
                } // i!=actualAxis
        } //for
}

Here is the call graph for this function:

Here is the caller graph for this function:

marExperiment * TreeExtraction_MH_JFC::GetInput

(

)

Definition at line 291 of file wxEmptyPanelWidget_2.cxx.

References _marExperiment.

                                              {
        return _marExperiment;
}

vtkImageData * TreeExtraction_MH_JFC::GetVolume

(

)

Definition at line 295 of file wxEmptyPanelWidget_2.cxx.

References _imagedataResult.

Referenced by wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

                                              {
        return _imagedataResult;
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::InitImagaDataResult

(

)

[private]

Definition at line 253 of file wxEmptyPanelWidget_2.cxx.

References _imagedataResult, _marExperiment, marExperiment::getDynData(), marDynData::getVolume(), New(), and Reset().

Referenced by Update().

                                               {
        Reset();
        vtkImageData *imagedata = _marExperiment->getDynData( )->getVolume( )->castVtk();
        imagedata->Update();
        _imagedataResult = vtkImageData::New();
        int dim[3];
        imagedata->GetDimensions(dim);
        _imagedataResult->SetScalarType(imagedata->GetScalarType());

        _imagedataResult->SetDimensions(dim[0],dim[1],dim[2]); 
        _imagedataResult->SetExtent( 0 , dim[0]-1 , 0 , dim[1]-1 , 0 , dim[2]-1 ); 
        _imagedataResult->SetWholeExtent(0 , dim[0]-1 , 0 , dim[1]-1 , 0 , dim[2]-1 );
        _imagedataResult->SetSpacing( imagedata->GetSpacing()  ); 
        _imagedataResult->SetOrigin( 0,0,0  ); 
        _imagedataResult->Update();
        _imagedataResult->UpdateData();
        unsigned short *ptr;
        int i,j,k;
        for ( i=0 ; i<dim[0] ; i++ ){
                for ( j=0 ; j<dim[1] ; j++ ){
                        for ( k=0 ; k<dim[2] ; k++ ){
                                ptr = (unsigned short *)        _imagedataResult->GetScalarPointer(i,j,k);
                                *ptr = 0;
                        }
                }
        }
}

Here is the call graph for this function:

Here is the caller graph for this function:

static TreeExtraction_MH_JFC* TreeExtraction_MH_JFC::New

(

)

[static]

Referenced by Execute(), InitImagaDataResult(), and wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::Reset

(

)

[private]

Definition at line 65 of file wxEmptyPanelWidget_2.cxx.

References _imagedataResult.

Referenced by InitImagaDataResult(), and ~TreeExtraction_MH_JFC().

                                 {
        if (_imagedataResult!=NULL){
                _imagedataResult->Delete();
        }
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::SetBranchLevel

(

int

branchLevel

)

Definition at line 299 of file wxEmptyPanelWidget_2.cxx.

References _branchLevel.

Referenced by wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

                                                         {
        _branchLevel = branchLevel;
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::SetCleanTreeLevel

(

int

cleanLevel

)

Definition at line 303 of file wxEmptyPanelWidget_2.cxx.

References _cleanLevel.

Referenced by wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

                                                             {
        _cleanLevel = (double)cleanLevel;
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::SetInput

(

marExperiment *

input

)

[virtual]

Definition at line 71 of file wxEmptyPanelWidget_2.cxx.

References _marExperiment.

Referenced by wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

                                                        {
        _marExperiment = input;
}

Here is the caller graph for this function:

void TreeExtraction_MH_JFC::Update

(

)

Definition at line 281 of file wxEmptyPanelWidget_2.cxx.

References _branchLevel, Execute(), and InitImagaDataResult().

Referenced by wxEmptyPanelWidget_2::OnBtnExtractTree_MH_JFC().

                                  {
/*
FILE *ff;
ff=fopen("c:/temp/MaracasTIME.txt","w");
fclose(ff);
*/
        InitImagaDataResult();
        Execute(_branchLevel);
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

Member Data Documentation

int TreeExtraction_MH_JFC::_branchLevel [private]

Definition at line 42 of file wxEmptyPanelWidget_2.h.

Referenced by SetBranchLevel(), TreeExtraction_MH_JFC(), and Update().

double TreeExtraction_MH_JFC::_cleanLevel [private]

Definition at line 41 of file wxEmptyPanelWidget_2.h.

Referenced by Execute(), SetCleanTreeLevel(), and TreeExtraction_MH_JFC().

vtkImageData* TreeExtraction_MH_JFC::_imagedataResult [private]

Definition at line 44 of file wxEmptyPanelWidget_2.h.

Referenced by Execute(), GetVolume(), InitImagaDataResult(), Reset(), and TreeExtraction_MH_JFC().

marExperiment* TreeExtraction_MH_JFC::_marExperiment [private]

Definition at line 45 of file wxEmptyPanelWidget_2.h.

Referenced by Execute(), GetInput(), InitImagaDataResult(), and SetInput().

carotidaBifurcacion* TreeExtraction_MH_JFC::_prgov [private]

Definition at line 43 of file wxEmptyPanelWidget_2.h.

Referenced by Execute().

---

The documentation for this class was generated from the following files:

• /home/davila/Creatis/All2/creatools_source/creaMaracasVisu/lib/maracasVisuLib/src/interface/wxWindows/widgets/include/wxEmptyPanelWidget_2.h
• /home/davila/Creatis/All2/creatools_source/creaMaracasVisu/lib/maracasVisuLib/src/interface/wxWindows/widgets/include/wxEmptyPanelWidget_2.cxx