int ext[6];
bbGetInputIn()->GetExtent(ext);
- int sizeX=ext[1]-ext[0];
- int sizeY=ext[3]-ext[2];
- int sizeZ=ext[5]-ext[4];
+ int sizeX=ext[1]-ext[0]+1;
+ int sizeY=ext[3]-ext[2]+1;
+ int sizeZ=ext[5]-ext[4]+1;
double spc[3];
bbGetInputIn()->GetSpacing(spc);
invSpc[2] = 1/spc[2];
mimage = vtkImageData::New();
- mimage->SetDimensions(sizeX,sizeY,sizeZ);
- mimage->SetOrigin(0,0,0);
- mimage->SetExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
- mimage->SetWholeExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
+ mimage->SetSpacing(bbGetInputIn()->GetSpacing());
+ mimage->SetDimensions(bbGetInputIn()->GetDimensions());
+ mimage->SetExtent(bbGetInputIn()->GetExtent());
+ mimage->SetOrigin(bbGetInputIn()->GetOrigin());
+ //mimage->SetDimensions(sizeX,sizeY,sizeZ);
+ //mimage->SetOrigin(0,0,0);
+ //mimage->SetExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
+ //mimage->SetWholeExtent( 0 , sizeX-1 , 0 , sizeY-1 , 0, sizeZ-1 );
mimage->SetScalarTypeToUnsignedShort();
- mimage->SetSpacing(spc);
+ //mimage->SetSpacing(spc);
mimage->AllocateScalars();
+
int i,j,k;
int sizeLstPointR = bbGetInputlstPointR().size();
int iAxe,sizeAxe=bbGetInputlstPointX().size();
p++;
}
- for (iAxe=0 ; iAxe<sizeAxe; iAxe++)
+ for (iAxe=0 ; iAxe<sizeAxe-1; iAxe++)
{
if (sizeLstPointR<iAxe)
{