X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=PackRecalage%2Fsrc%2FbbPackRecalageReSlicerBox.cxx;h=984f4549b516ecd42adb04886e9f2f7deecaf3da;hb=5b6e06ea55193ecce4fe2c59b1f3cb698835020b;hp=9aa94249fbea1fc5ca53cba907d9228f21fcea55;hpb=683ff9386e374db4f94ae5d4c32859cc808c3ad4;p=creaRigidRegistration.git

diff --git a/PackRecalage/src/bbPackRecalageReSlicerBox.cxx b/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
index 9aa9424..984f454 100644
--- a/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
+++ b/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
@@ -1,53 +1,161 @@
+/*
+# ---------------------------------------------------------------------
+#
+# Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image 
+#                        pour la Santé)
+# Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
+#
+#  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 "bbPackRecalageReSlicerBox.h"
 #include "bbPackRecalagePackage.h"
 
 namespace bbPackRecalage
 {
 
+
+//-----------------------------------------------------------------------------
+void ReSlicerBox::GetBackInfo(double *backInfo, vtkMatrix4x4 *vtkmatrix, std::vector<int> point )
+{
+
+		backInfo[0]=vtkmatrix->GetElement(0,0);
+		backInfo[1]=vtkmatrix->GetElement(0,1);
+		backInfo[2]=vtkmatrix->GetElement(0,2);
+		backInfo[3]=vtkmatrix->GetElement(0,3);
+		backInfo[4]=vtkmatrix->GetElement(1,0);
+		backInfo[5]=vtkmatrix->GetElement(1,1);
+		backInfo[6]=vtkmatrix->GetElement(1,2);
+		backInfo[7]=vtkmatrix->GetElement(1,3);
+		backInfo[8]=vtkmatrix->GetElement(2,0);
+		backInfo[9]=vtkmatrix->GetElement(2,1);
+		backInfo[10]=vtkmatrix->GetElement(2,2);
+		backInfo[11]=vtkmatrix->GetElement(2,3);
+		backInfo[12]=vtkmatrix->GetElement(3,0);
+		backInfo[13]=vtkmatrix->GetElement(3,1);
+		backInfo[14]=vtkmatrix->GetElement(3,2);
+		backInfo[15]=vtkmatrix->GetElement(3,3);
+
+		backInfo[16]=-point[0];
+		backInfo[17]=-point[1];;
+		backInfo[18]=-point[2];;
+}
+
+//-----------------------------------------------------------------------------
+bool ReSlicerBox::CompareBackInfo(double* backInfoA, double * backInfoB)
+{
+	bool ok=true;
+	int i;
+	for (i=0; i<sizeBackInfo; i++ )
+	{
+           if (backInfoA[i]!=backInfoB[i])
+	   {
+		ok=false;
+	   } // if
+	} //for
+	return ok;
+}
+
+
+//-----------------------------------------------------------------------------
 BBTK_ADD_BLACK_BOX_TO_PACKAGE(PackRecalage,ReSlicerBox)
 BBTK_BLACK_BOX_IMPLEMENTATION(ReSlicerBox,bbtk::AtomicBlackBox);
 void ReSlicerBox::Process()
 {
- 
-	if(bbGetInputOrigin()[0] != NULL)
+	if(!bbGetInputOrigin().empty()  && bbGetInputTransform() != NULL)
 	{
-		image = vtkImageChangeInformation::New();
+		//image = vtkImageChangeInformation::New();
 		image->SetInput( bbGetInputIn() );
-		image->SetOutputSpacing( 1,1,1 ); 
-	
-		slicer =vtkImageReslice::New();
-		slicer->SetInput( image->GetOutput() );
-		slicer->SetResliceTransform(bbGetInputTransform());
- 		slicer->SetOutputOrigin( 0 , 0 , 0 );
-		slicer->SetOutputOrigin( -(bbGetInputOrigin()[0]) , -(bbGetInputOrigin()[1]) , -(bbGetInputOrigin()[2]) );
-		slicer->Update();
-	
-		imageResult = vtkImageChangeInformation::New();
-		imageResult->SetInput( slicer->GetOutput() );
-		double spc[3];
-		bbGetInputIn()->GetSpacing(spc);
-		imageResult->SetOutputSpacing( spc ); 
-		imageResult->SetOutputOrigin( 0,0,0 ); 
+		image->SetOutputSpacing( 1,1,1 );
+
+		if(bbGetInputCentered())
+		{
+			image->CenterImageOn();
+		}
+
+                double tmpbackInfo[19];
+		GetBackInfo(tmpbackInfo,bbGetInputTransform()->GetMatrix(),bbGetInputOrigin() );
+
+		if ( CompareBackInfo(backInfoA,tmpbackInfo)==false ) 
+		{
+			GetBackInfo(backInfoA,bbGetInputTransform()->GetMatrix(),bbGetInputOrigin() );
+
+			//slicer =vtkImageReslice::New();
+			slicer->SetInput( image->GetOutput() );
+			slicer->SetInformationInput( image->GetOutput() );
+			slicer->SetResliceTransform(bbGetInputTransform());
+ 			slicer->SetOutputOrigin( -(bbGetInputOrigin()[0]) , -(bbGetInputOrigin()[1]) , -(bbGetInputOrigin()[2]) );
+
+			if(bbGetInputInterpolate())
+			{
+				slicer->InterpolateOn();
+			}		
+			slicer->Update();
+
+			//imageResult = vtkImageChangeInformation::New();
+			imageResult->SetInput( slicer->GetOutput() );
+			double spc[3];
+			bbGetInputIn()->GetSpacing(spc);
+			imageResult->SetOutputSpacing( spc ); 
+			imageResult->SetOutputOrigin( 0,0,0 ); 
 	
-		bbSetOutputOut( imageResult->GetOutput() );		
-	}
-	else
-	{
-		bbSetOutputOut( bbGetInputIn() );
-	}
-	  
+			bbSetOutputOut( imageResult->GetOutput() );
+
+			if (bbGetInputTransform()!=NULL)
+			{ 
+		  	 	bbGetInputTransform()->Update();
+		   		vtkMatrix4x4 *m = bbGetInputTransform()->GetMatrix();
+		   		if (m!=NULL)
+		   		{
+		     			printf("EED ReSlicerBox::Process Translation %d %d %d \n",-(bbGetInputOrigin()[0]),-(bbGetInputOrigin()[1]),-(bbGetInputOrigin()[2]));
+		     			printf("EED ReSlicerBox::Process Matrix %f %f %f %f\n", m->GetElement(0,0),m->GetElement(0,1),m->GetElement(0,2),m->GetElement(0,3));
+		     			printf("EED ReSlicerBox::Process Matrix %f %f %f %f\n", m->GetElement(1,0),m->GetElement(1,1),m->GetElement(1,2),m->GetElement(1,3));
+		     			printf("EED ReSlicerBox::Process Matrix %f %f %f %f\n", m->GetElement(2,0),m->GetElement(2,1),m->GetElement(2,2),m->GetElement(2,3));
+		     			printf("EED ReSlicerBox::Process Matrix %f %f %f %f\n", m->GetElement(3,0),m->GetElement(3,1),m->GetElement(3,2),m->GetElement(3,3));
+	     	   		} // if m 
+			} // if Transform
+		} // Compare BackInfo
+	} else {
+			bbSetOutputOut( NULL );
+	}  
 }
+
+
 void ReSlicerBox::bbUserSetDefaultValues()
-{
- 
-	std::vector<int> nuevo(3,0);
-	nuevo.push_back(0);
-	bbSetInputOrigin(nuevo); 
-    bbSetInputIn(NULL); 
-    bbSetInputTransform(vtkIdentityTransform::New());
-    bbSetOutputOut(NULL);
+{ 
+	std::vector<int> empty;
+	bbSetInputOrigin(empty); 
+    	bbSetInputIn(NULL); 
+    	bbSetInputTransform(NULL);
+    	bbSetOutputOut(NULL);
+
+	bbSetInputCentered(false);
+	bbSetInputInterpolate(false);
   
+	sizeBackInfo=19;
+	for (int i=0; i<sizeBackInfo; i++)
+	{
+		backInfoA[i]=0;
+	} // for
 }
+
+
 void ReSlicerBox::bbUserInitializeProcessing()
 {
  
@@ -59,16 +167,13 @@ void ReSlicerBox::bbUserInitializeProcessing()
     slicer =vtkImageReslice::New();
 	imageResult = vtkImageChangeInformation::New();
 }
+	
 void ReSlicerBox::bbUserFinalizeProcessing()
-{
- 
-//  THE FINALIZATION METHOD BODY : 
-//    Here does nothing  
-//    but this is where you should desallocate the internal/output pointers  
-//    if any 
+{ 
 	image->Delete();
 	slicer->Delete();
     imageResult->Delete();
 }
+
 }
-// EO namespace bbPackRecalage
\ No newline at end of file
+// EO namespace bbPackRecalage