]> Creatis software - bbtk.git/blob - packages/vtk/src/bbvtkInversCrop.cxx
Merge branch 'vtk8itk5wx3-mingw64' of ssh://git.creatis.insa-lyon.fr/bbtk into vtk8it...
[bbtk.git] / packages / vtk / src / bbvtkInversCrop.cxx
1 //===== 
2 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
3 //===== 
4 #include "bbvtkInversCrop.h"
5 #include "bbvtkPackage.h"
6 namespace bbvtk
7 {
8
9 BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,InversCrop)
10 BBTK_BLACK_BOX_IMPLEMENTATION(InversCrop,bbtk::AtomicBlackBox);
11 //===== 
12 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
13 //===== 
14 void InversCrop::Process()
15 {
16
17 // THE MAIN PROCESSING METHOD BODY
18 //   Here we simply set the input 'In' value to the output 'Out'
19 //   And print out the output value
20 // INPUT/OUTPUT ACCESSORS ARE OF THE FORM :
21 //    void bbSet{Input|Output}NAME(const TYPE&)
22 //    const TYPE& bbGet{Input|Output}NAME() const 
23 //    Where :
24 //    * NAME is the name of the input/output
25 //      (the one provided in the attribute 'name' of the tag 'input')
26 //    * TYPE is the C++ type of the input/output
27 //      (the one provided in the attribute 'type' of the tag 'input')
28
29 printf("\n\n EED InversCrop::Process Start \n" );
30
31         if (bbGetInputActive()==true)
32         {
33             if ( (bbGetInputImageFix()!=NULL) && (bbGetInputImageMove()!=NULL) )
34             {
35              if ( bbGetInputImageFix()->GetScalarType()==bbGetInputImageMove()->GetScalarType() ) 
36              {
37                 // Creating Image
38                 int dim[3];
39                 int ext[6];
40         
41 //EED 2017-01-01 Migration VTK7
42 #if (VTK_MAJOR_VERSION <= 5) 
43         bbGetInputImageFix()->GetWholeExtent(ext);
44 #endif
45 #if (VTK_MAJOR_VERSION >= 6) 
46         bbGetInputImageFix()->GetExtent(ext);
47 #endif
48                 dim[0]= ext[1]-ext[0]+1;
49                 dim[1]= ext[3]-ext[2]+1;
50                 dim[2] = 1;
51                 if (bbGetInputImageFix()->GetDataDimension()==3) 
52                 {
53                         dim[2] = ext[5]-ext[4]+1;
54                 } 
55                 
56 printf("EED InversCrop::Process 1   %d\n", bbGetInputImageFix()->GetDataDimension()    );
57 printf("EED InversCrop::Process 2   %d %d  %d \n",dim[0],dim[1],dim[2]  );
58
59                 if (bbGetInputType()==0)
60                 {
61                         _imageoutput = vtkImageData::New();
62                         _imageoutput->Initialize();
63                         _imageoutput->SetSpacing( bbGetInputImageFix()->GetSpacing() );
64                         _imageoutput->SetDimensions(  dim[0], dim[1], dim[2] );                 
65 //EED 2017-01-01 Migration VTK7
66 #if (VTK_MAJOR_VERSION <= 5) 
67                         _imageoutput->SetScalarType( bbGetInputImageFix()->GetScalarType() );
68                         _imageoutput->AllocateScalars();
69 #endif
70 #if (VTK_MAJOR_VERSION >= 6) 
71                         _imageoutput->AllocateScalars(bbGetInputImageFix()->GetScalarType() , 1);
72 #endif
73                         
74                         // Duplicating Fix Image
75                         long sizeimage = dim[0]*dim[1]*dim[2]*bbGetInputImageFix()->GetScalarSize();    
76                         memcpy( _imageoutput->GetScalarPointer() , bbGetInputImageFix()->GetScalarPointer() , sizeimage);
77                 }
78                 if (bbGetInputType()==1)
79                 {
80                         _imageoutput=bbGetInputImageFix();
81                 }
82
83                 // Copy the Move Image
84                 int j,k; 
85                 int px,py,pz;
86
87 //EED 2017-01-01 Migration VTK7
88 #if (VTK_MAJOR_VERSION <= 5) 
89                 bbGetInputImageMove()->GetWholeExtent(ext);
90 #endif
91 #if (VTK_MAJOR_VERSION >= 6) 
92                 bbGetInputImageMove()->GetExtent(ext);
93 #endif
94                 int dimMoveX    = ext[1]-ext[0]+1;
95                 int dimMoveY    = ext[3]-ext[2]+1;
96                 int dimMoveZ    = 1;
97                 if (bbGetInputImageMove()->GetDataDimension()==3) 
98                 {
99                         dimMoveZ = ext[5]-ext[4]+1;
100                 } 
101
102                 
103 printf("EED InversCrop::Process 3   %d %d  %d \n",dimMoveX,dimMoveY,dimMoveZ  );
104
105                 
106                 
107                 int spxM                = 0;  // start px MoveImage
108                 int sizeXM              = 0;  // sizeX MoveImage
109 printf("EED InversCrop::Process 3.1   sizeOrigin %d \n", bbGetInputOrigin().size() );
110 printf("EED InversCrop::Process 3.1   Origin %d  %d %d\n", bbGetInputOrigin()[0],bbGetInputOrigin()[1],bbGetInputOrigin()[2] );
111
112                 px                              = bbGetInputOrigin()[0];
113                 spxM                    = 0;
114                 if (px<0)
115                 { 
116                    spxM = px*(-1);
117                    px   = 0;
118                 }
119                 sizeXM  = dimMoveX-spxM;
120                 if (px+sizeXM>=dim[0]) { sizeXM=dim[0]-px; }
121                 sizeXM  = sizeXM * bbGetInputImageFix()->GetScalarSize();
122                 for (k=0; k<dimMoveZ; k++)
123                 {
124                    for (j=0; j<dimMoveY; j++)
125                    {
126                         py = j+bbGetInputOrigin()[1];
127                                 if (bbGetInputImageMove()->GetDataDimension()==3) 
128                                 {
129                                         pz = k+bbGetInputOrigin()[2];
130                                 } else {
131                                         pz = k;
132                                 }// Dimension 3
133                                 if ( (py<dim[1]) && (pz<dim[2]) &&
134                                  (py>=0)    && (pz>=0)  &&
135                                          (sizeXM>0) ) 
136                                 {
137                                         memcpy( _imageoutput->GetScalarPointer(px,py,pz) , bbGetInputImageMove()->GetScalarPointer(spxM,j,k) , sizeXM );
138                                 }
139
140                    } // for j
141                 } // for k
142 printf("EED InversCrop::Process 4   %d %d  %d \n",px,py,pz  );
143
144                 _imageoutput->Modified();
145         } else { // If Image Fixe Move the same GetScalarType
146                    printf ("ERROR: InversCrop  both ImageFixe and ImageMove need the same format.\n");
147                    printf ("       type ImageFix:%d   type ImageMove:%d\n", bbGetInputImageFix()->GetScalarType(), bbGetInputImageMove()->GetScalarType() );
148
149               }  
150             } else { // If Image Fixe Move != NULL
151                         printf ("ERROR: InversCrop  need ImageFixe and ImageMove to run.\n");
152            } 
153          bbSetOutputOut(_imageoutput);
154
155         } // if Active
156 printf("EED InversCrop::Process End \n" );
157
158 }
159
160 //===== 
161 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
162 //===== 
163 void InversCrop::bbUserSetDefaultValues()
164 {
165 //  SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX 
166 //    Here we initialize the input 'In' to 0
167    bbSetInputActive(true);
168    bbSetInputType(0);
169    bbSetInputImageFix(NULL);
170    bbSetInputImageMove(NULL);
171    std::vector<int> origin;
172    origin.push_back(0);
173    origin.push_back(0);
174    origin.push_back(0);
175    bbSetInputOrigin(origin);
176    _imageoutput=NULL;
177 }
178 //===== 
179 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
180 //===== 
181 void InversCrop::bbUserInitializeProcessing()
182 {
183
184 //  THE INITIALIZATION METHOD BODY :
185 //    Here does nothing 
186 //    but this is where you should allocate the internal/output pointers 
187 //    if any 
188
189   
190 }
191 //===== 
192 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
193 //===== 
194 void InversCrop::bbUserFinalizeProcessing()
195 {
196
197 //  THE FINALIZATION METHOD BODY :
198 //    Here does nothing 
199 //    but this is where you should desallocate the internal/output pointers 
200 //    if any
201   
202 }
203 }
204 // EO namespace bbvtk
205
206