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

diff --git a/PackRecalage/src/bbPackRecalageReSlicerBox.cxx b/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
index 29953b0..984f454 100644
--- a/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
+++ b/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
@@ -30,6 +30,50 @@
 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()
@@ -45,45 +89,50 @@ void ReSlicerBox::Process()
 			image->CenterImageOn();
 		}
 
-		//slicer =vtkImageReslice::New();
-		slicer->SetInput( image->GetOutput() );
-		slicer->SetInformationInput( image->GetOutput() );
-		slicer->SetResliceTransform(bbGetInputTransform());
- 		slicer->SetOutputOrigin( -(bbGetInputOrigin()[0]) , -(bbGetInputOrigin()[1]) , -(bbGetInputOrigin()[2]) );
-
+                double tmpbackInfo[19];
+		GetBackInfo(tmpbackInfo,bbGetInputTransform()->GetMatrix(),bbGetInputOrigin() );
 
-		if(bbGetInputInterpolate())
+		if ( CompareBackInfo(backInfoA,tmpbackInfo)==false ) 
 		{
-			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 ); 
+			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() );
-
-		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
-	}
-	else
-	{
-		bbSetOutputOut( NULL );
+			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 );
 	}  
 }
 
@@ -92,12 +141,18 @@ void ReSlicerBox::bbUserSetDefaultValues()
 { 
 	std::vector<int> empty;
 	bbSetInputOrigin(empty); 
-    bbSetInputIn(NULL); 
-    bbSetInputTransform(NULL);
-    bbSetOutputOut(NULL);
+    	bbSetInputIn(NULL); 
+    	bbSetInputTransform(NULL);
+    	bbSetOutputOut(NULL);
+
 	bbSetInputCentered(false);
 	bbSetInputInterpolate(false);
   
+	sizeBackInfo=19;
+	for (int i=0; i<sizeBackInfo; i++)
+	{
+		backInfoA[i]=0;
+	} // for
 }