]> Creatis software - creaRigidRegistration.git/blobdiff - lib/Transparency.cxx
git://git.creatis.insa-lyon.fr / creaRigidRegistration.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Added template methods
[creaRigidRegistration.git] / lib / Transparency.cxx
diff --git a/lib/Transparency.cxx b/lib/Transparency.cxx
index 0cf822d90a6b36df7954c8ac5c243472748af1f3..1459e2da459d5d6dfb3cd9c4877af42e0bb81108 100644 (file)
--- a/lib/Transparency.cxx
+++ b/lib/Transparency.cxx
@@ -138,20 +138,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (char)*dataImagePointer1*_percent + (char)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_SIGNED_CHAR)
        {
@@ -160,20 +147,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (signed char)*dataImagePointer1*_percent + (signed char)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_UNSIGNED_CHAR)
        {
@@ -182,20 +156,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (unsigned char)*dataImagePointer1*_percent + (unsigned char)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_SHORT)
        {
@@ -204,20 +165,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (short)*dataImagePointer1*_percent + (short)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_UNSIGNED_SHORT)
        {
@@ -226,20 +174,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (unsigned short)*dataImagePointer1*_percent + (unsigned short)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_INT)
        {
@@ -248,20 +183,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (int)*dataImagePointer1*_percent + (int)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_UNSIGNED_INT)
        {
@@ -270,20 +192,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (unsigned int)*dataImagePointer1*_percent + (unsigned int)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_LONG)
        {
@@ -292,20 +201,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (long)*dataImagePointer1*_percent + (long)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_UNSIGNED_LONG)
        {
@@ -314,20 +210,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (unsigned long)*dataImagePointer1*_percent + (unsigned long)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);              
        }
        else if(_type == VTK_FLOAT)
        {
@@ -336,20 +219,7 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                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);
-               for(i = 0; i < sizeX; i++)
-               {
-                       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);
-                               
-                               *dataImageResultPointer = (float)*dataImagePointer1*_percent + (float)*dataImagePointer2*(1-_percent);
-                       }
-               }
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
        }
        else if(_type == VTK_DOUBLE)
        {
@@ -357,24 +227,32 @@ void Transparency::createImage(vtkImageData *img1, vtkImageData *img2, int sizeX
                double* dataImagePointer1 = NULL;
                double* dataImagePointer2 = NULL;
                double* dataImageResultPointer = NULL;
+
+               createImageByType(dataImagePointer1, dataImagePointer2, dataImageResultPointer, img1, img2, sizeX, sizeY);
+       }
+
+       _newImage->Update();
+}
+
+template <class T>
+void Transparency::createImageByType(T* dataImagePointer1, T* dataImagePointer2, T* dataImageResultPointer, vtkImageData *img1, vtkImageData *img2, int sizeX, int sizeY)
+{
+       int i, j;
        
-               dataImagePointer1 = (double*) img1->GetScalarPointer(0,0,0);
-               dataImagePointer2 = (double*) img2->GetScalarPointer(0,0,0);
-               dataImageResultPointer = (double*) _newImage->GetScalarPointer(0,0,0);
-               for(i = 0; i < sizeX; i++)
+       dataImagePointer1 = (T*) img1->GetScalarPointer(0,0,0);
+       dataImagePointer2 = (T*) img2->GetScalarPointer(0,0,0);
+       dataImageResultPointer = (T*) _newImage->GetScalarPointer(0,0,0);
+       for(i = 0; i < sizeX; i++)
+       {
+               for(j = 0; j < sizeY; j++)
                {
-                       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);
-                               
-                               *dataImageResultPointer = (double)*dataImagePointer1*_percent + (double)*dataImagePointer2*(1-_percent);                                
-                       }
+                       dataImagePointer1 = (T*)img1->GetScalarPointer(i,j,0);
+                       dataImagePointer2 = (T*)img2->GetScalarPointer(i,j,0);
+                       dataImageResultPointer = (T*)_newImage->GetScalarPointer(i,j,0);
+                       
+                       *dataImageResultPointer = (T)*dataImagePointer1*_percent + (T)*dataImagePointer2*(1-_percent);                          
                }
        }
-
-       _newImage->Update();
 }
 
 /*
@@ -407,4 +285,5 @@ void Transparency::setImage2(vtkImageData *image)
 void Transparency::setPercent(int percent)
 {
     _percent=(double)percent/100.0;    
-}
\ No newline at end of file
+}
+
Manual rigid registration