squaresY=0;
newImage = NULL;
processed=false;
+ newImage = NULL;
+ resample = NULL;
+ result = NULL;
}
/*
//if (image1 != NULL ) { image1->Delete(); }
//if (image2 != NULL ) { image2->Delete(); }
if (newImage != NULL ) { newImage->Delete(); }
+ if (result != NULL ) { result->Delete(); }
}
//------------------------------------------------------------
double factorX = 1;
double factorY = 1;
- vtkImageResample *resample = vtkImageResample::New();
- vtkImageData *result;
+ resample = vtkImageResample::New();
+ result = vtkImageData::New();
if(numPixelsImg1<numPixelsImg2)
{
resample->SetInput(image1);
void CheckBoard::initialize(int dimensions[], double spacing[])
{
// Setting the new image
+ if(newImage != NULL){
+ newImage->Delete();
+ }
newImage = vtkImageData::New();
newImage->SetScalarType(type);
newImage->SetSpacing(spacing);
newImage->SetDimensions(dimensions);
newImage->AllocateScalars();
newImage->Update();
+
}
char* dataImagePointer1 = NULL;
char* dataImagePointer2 = NULL;
char* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (char*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (char*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (char*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
-
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (char*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (char*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (char*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (char) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (char) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_SIGNED_CHAR)
{
signed char* dataImagePointer1 = NULL;
signed char* dataImagePointer2 = NULL;
signed char* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (signed char*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (signed char*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (signed char*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (signed char*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (signed char*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (signed char*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (signed char) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (signed char) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_UNSIGNED_CHAR)
{
unsigned char* dataImagePointer1 = NULL;
unsigned char* dataImagePointer2 = NULL;
unsigned char* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (unsigned char*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (unsigned char*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (unsigned char*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (unsigned char*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (unsigned char*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (unsigned char*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned char) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned char) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_SHORT)
{
short* dataImagePointer1 = NULL;
short* dataImagePointer2 = NULL;
short* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (short*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (short*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (short*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (short*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (short*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (short*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (short) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (short) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_UNSIGNED_SHORT)
{
unsigned short* dataImagePointer1 = NULL;
unsigned short* dataImagePointer2 = NULL;
unsigned short* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (unsigned short*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (unsigned short*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (unsigned short*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (unsigned short*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (unsigned short*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (unsigned short*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned short) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned short) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_INT)
int* dataImagePointer1 = NULL;
int* dataImagePointer2 = NULL;
int* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (int*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (int*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (int*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (int*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (int*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (int*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (int) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (int) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_UNSIGNED_INT)
{
unsigned int* dataImagePointer1 = NULL;
unsigned int* dataImagePointer2 = NULL;
unsigned int* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (unsigned int*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (unsigned int*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (unsigned int*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (unsigned int*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (unsigned int*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (unsigned int*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned int) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned int) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_LONG)
{
long* dataImagePointer1 = NULL;
long* dataImagePointer2 = NULL;
long* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (long*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (long*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (long*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (long*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (long*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (long*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (long) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (long) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_UNSIGNED_LONG)
{
unsigned long* dataImagePointer1 = NULL;
unsigned long* dataImagePointer2 = NULL;
unsigned long* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (unsigned long*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (unsigned long*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (unsigned long*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (unsigned long*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (unsigned long*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (unsigned long*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned long) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (unsigned long) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_FLOAT)
{
float* dataImagePointer1 = NULL;
float* dataImagePointer2 = NULL;
float* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (float*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (float*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (float*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (float*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (float*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (float*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (float) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (float) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
-
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
else if(imageType == VTK_DOUBLE)
{
double* dataImagePointer1 = NULL;
double* dataImagePointer2 = NULL;
double* dataImageResultPointer = NULL;
-
- dataImagePointer1 = (double*) img1->GetScalarPointer(0,0,0);
- dataImagePointer2 = (double*) img2->GetScalarPointer(0,0,0);
- dataImageResultPointer = (double*) newImage->GetScalarPointer(0,0,0);
-
- if(squaresX == 0)
- {
- squaresX = 1;
- }
- if(squaresY == 0)
- {
- squaresY = 1;
- }
- int divX = sizeX/squaresX;
- int divY = sizeY/squaresY;
- int i, j, counterX, counterY;
- for(i = 0; i < sizeX; i++)
- {
- counterY = 0;
- for(j = 0; j < sizeY; j++)
- {
- dataImagePointer1 = (double*)img1->GetScalarPointer(i,j,0);
- dataImagePointer2 = (double*)img2->GetScalarPointer(i,j,0);
- dataImageResultPointer = (double*)newImage->GetScalarPointer(i,j,0);
-
- if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (double) *dataImagePointer1;
- }
- else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
- {
- *dataImageResultPointer = (double) *dataImagePointer2;
- }
-
- if(counterX >= 2*divX)
- {
- counterX = 0;
- }
- else if(counterY >= 2*divY)
- {
- counterY = 0;
- }
- counterY++;
- }
- counterX++;
- }
+ createImageByType(dataImagePointer1,dataImagePointer2,dataImageResultPointer,img1,img2,sizeX,sizeY);
}
newImage->Update();
// newImage->Modified();
}
+/**Template used for the different image types
+*/
+template <class T>
+void CheckBoard::createImageByType(T* dataImagePointer1, T* dataImagePointer2, T* dataImageResultPointer, vtkImageData *img1, vtkImageData *img2, int sizeX, int sizeY)
+{
+ dataImagePointer1 = (T*) img1->GetScalarPointer(0,0,0);
+ dataImagePointer2 = (T*) img2->GetScalarPointer(0,0,0);
+ dataImageResultPointer = (T*) newImage->GetScalarPointer(0,0,0);
+
+ if(squaresX == 0)
+ {
+ squaresX = 1;
+ }
+ if(squaresY == 0)
+ {
+ squaresY = 1;
+ }
+
+ int divX = sizeX/squaresX;
+ int divY = sizeY/squaresY;
+ int i, j, counterX, counterY;
+ for(i = 0; i < sizeX; i++)
+ {
+ counterY = 0;
+ for(j = 0; j < sizeY; j++)
+ {
+ dataImagePointer1 = (T*)img1->GetScalarPointer(i,j,0);
+ dataImagePointer2 = (T*)img2->GetScalarPointer(i,j,0);
+ dataImageResultPointer = (T*)newImage->GetScalarPointer(i,j,0);
+ if( (counterX <= divX && counterY <= divY) || ((counterX >= divX && counterX <= 2*divX) && (counterY >= divY && counterY <= 2*divY)) )
+ {
+ *dataImageResultPointer = (T) *dataImagePointer1;
+ }
+ else //if( ((counterX >= divX && counterX <= 2*divX) && counterY <= divY) || (counterX <= divX && (counterY >= divY && counterY <= 2*divY)) )
+ {
+ *dataImageResultPointer = (T) *dataImagePointer2;
+ }
+
+ if(counterX >= 2*divX)
+ {
+ counterX = 0;
+ }
+ else if(counterY >= 2*divY)
+ {
+ counterY = 0;
+ }
+ counterY++;
+ }
+ counterX++;
+ }
+}
/*
* Get new image
}
}
-
//------------------------------------------------------------
void CheckBoard::setInputImage1(vtkImageData *_image)
{
void CheckBoard::setRows(int rows)
{
squaresY=rows;
-}
\ No newline at end of file
+}
+
git://git.creatis.insa-lyon.fr / creaRigidRegistration.git / blobdiff