1 /*# ---------------------------------------------------------------------
3 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
5 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
6 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
7 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9 # This software is governed by the CeCILL-B license under French law and
10 # abiding by the rules of distribution of free software. You can use,
11 # modify and/ or redistribute the software under the terms of the CeCILL-B
12 # license as circulated by CEA, CNRS and INRIA at the following URL
13 # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
14 # or in the file LICENSE.txt.
16 # As a counterpart to the access to the source code and rights to copy,
17 # modify and redistribute granted by the license, users are provided only
18 # with a limited warranty and the software's author, the holder of the
19 # economic rights, and the successive licensors have only limited
22 # The fact that you are presently reading this means that you have had
23 # knowledge of the CeCILL-B license and that you accept its terms.
24 # ------------------------------------------------------------------------ */
26 #include "bbmaracasvisuAxeVolume.h"
27 #include "bbcreaMaracasVisuPackage.h"
28 namespace bbcreaMaracasVisu
31 BBTK_ADD_BLACK_BOX_TO_PACKAGE(creaMaracasVisu,AxeVolume)
32 BBTK_BLACK_BOX_IMPLEMENTATION(AxeVolume,bbtk::AtomicBlackBox);
33 void AxeVolume::Process()
36 printf("EED AxeVolume::Process start \n");
44 bbGetInputIn()->GetExtent(ext);
45 int sizeX=ext[1]-ext[0]+1;
46 int sizeY=ext[3]-ext[2]+1;
47 int sizeZ=ext[5]-ext[4]+1;
50 bbGetInputIn()->GetSpacing(spc);
56 mimage = vtkImageData::New();
57 mimage->SetSpacing(bbGetInputIn()->GetSpacing());
58 mimage->SetDimensions(bbGetInputIn()->GetDimensions());
59 mimage->SetExtent(bbGetInputIn()->GetExtent());
60 mimage->SetOrigin(bbGetInputIn()->GetOrigin());
61 //mimage->SetDimensions(sizeX,sizeY,sizeZ);
62 //mimage->SetOrigin(0,0,0);
63 //mimage->SetExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
64 //mimage->SetWholeExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
65 mimage->SetScalarTypeToUnsignedShort();
66 //mimage->SetSpacing(spc);
67 mimage->AllocateScalars();
72 int sizeLstPointR = bbGetInputlstPointR().size();
73 int iAxe,sizeAxe = bbGetInputlstPointX().size();
77 int sizeImage = sizeX*sizeY*sizeZ;
83 p = (unsigned short*)mimage->GetScalarPointer (0, 0, 0);
84 for ( i=0 ; i<sizeImage ; i++)
90 for (iAxe=0 ; iAxe<sizeAxe; iAxe++)
92 if (iAxe<sizeLstPointR)
94 r = bbGetInputlstPointR()[ iAxe ]* invSpc[0];
96 if (bbGetInputlstPointR().size()>=1)
98 r = bbGetInputlstPointR()[ bbGetInputlstPointR().size()-1 ] * invSpc[0];
103 px = bbGetInputlstPointX()[iAxe] * invSpc[0];
104 py = bbGetInputlstPointY()[iAxe] * invSpc[1];
105 pz = bbGetInputlstPointZ()[iAxe] * invSpc[2];
114 for ( i=px1 ; i<=px2 ; i++ )
118 for ( j=py1 ; j<py2 ; j++ )
122 for ( k=pz1 ; k<pz2 ; k++ )
124 if ( (i>=0) && (i<sizeX) && (j>=0) && (j<sizeY) &&(k>=0) && (k<sizeZ) ){
125 p = (unsigned short*)mimage->GetScalarPointer (i, j, k);
130 if ( rx + ry + rz <= rr )
140 bbSetOutputOut( mimage );
142 printf("EED AxeVolume::Process end \n");
147 //-----------------------------------------------------------------
148 void AxeVolume::bbUserSetDefaultValues()
153 //-----------------------------------------------------------------
154 void AxeVolume::bbUserInitializeProcessing()
158 //-----------------------------------------------------------------
159 void AxeVolume::bbUserFinalizeProcessing()
163 //-----------------------------------------------------------------
166 // EO namespace bbcreaMaracasVisu