-#include "CheckBoard.h"
+/*
+# ---------------------------------------------------------------------
+#
+# Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
+# pour la Santé)
+# Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
+# Previous Authors : Laurent Guigues, Jean-Pierre Roux
+# CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
+#
+# This software is governed by the CeCILL-B license under French law and
+# abiding by the rules of distribution of free software. You can use,
+# modify and/ or redistribute the software under the terms of the CeCILL-B
+# license as circulated by CEA, CNRS and INRIA at the following URL
+# http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+# or in the file LICENSE.txt.
+#
+# As a counterpart to the access to the source code and rights to copy,
+# modify and redistribute granted by the license, users are provided only
+# with a limited warranty and the software's author, the holder of the
+# economic rights, and the successive licensors have only limited
+# liability.
+#
+# The fact that you are presently reading this means that you have had
+# knowledge of the CeCILL-B license and that you accept its terms.
+# ------------------------------------------------------------------------ */
+#include "CheckBoard.h"
/*
* Constructor
*/
//------------------------------------------------------------
-CheckBoard::CheckBoard()
+CheckBoard::CheckBoard() :
+ image1(NULL),
+ image2(NULL),
+ squaresX(0),
+ squaresY(0),
+ newImage(NULL),
+ processed(false),
+ resample(NULL),
+ result(NULL)
{
+/*
image1=NULL;
image2=NULL;
squaresX=0;
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(); }
}
//------------------------------------------------------------
if(image1!=NULL && image2!=NULL)
{
-
// IMAGE 1
// Information from image1
image1->GetExtent(extImg1);
type = image1->GetScalarType();
-
+
// The Z dimension is 1 for a 2D image
dimImg1[0] = extImg1[1] - extImg1[0]+1;
dimImg1[1] = extImg1[3] - extImg1[2]+1;
long numPixelsImg1 = dimImg1[0]*dimImg1[1];
long numPixelsImg2 = dimImg2[0]*dimImg2[1];
-
+
double factorX = 1;
double factorY = 1;
- vtkImageResample *resample = vtkImageResample::New();
- vtkImageData *result;
+ resample = vtkImageResample::New();
+ result = vtkImageData::New();
if(numPixelsImg1<numPixelsImg2)
{
resample->SetInput(image1);
result = resample->GetOutput();
result->Update();
createImage(result,image2,dimImg2[0],dimImg2[1]);
-
} //if
else if (numPixelsImg1>numPixelsImg2)
{
initialize(dimImg1, spcImg1);
result = resample->GetOutput();
-
+
result->Update();
createImage(image1,result,dimImg1[0],dimImg1[1]);
result = resample->GetOutput();
result->Update();
-
+
createImage(image1,result,dimImg1[0],dimImg1[1]);
} //else
-
+
resample->Delete();
processed=true;
}
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);
}
-
//std::cout << "The image has been checkboardized!" << std::endl;
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;
+ int counterX=0, counterY; // JPR : uninitialized
+ 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
{
squaresY=rows;
}
-
git://git.creatis.insa-lyon.fr / creaRigidRegistration.git / blobdiff