2075 creaRigidRegistration Feature New Normal Registration 3D 1P

[creaRigidRegistration.git] / PackRecalage / src / bbPackRecalageReSlicerBox.cxx

diff --git a/PackRecalage/src/bbPackRecalageReSlicerBox.cxx b/PackRecalage/src/bbPackRecalageReSlicerBox.cxx
index 29953b0481c3882fb7742b3943776ff83d6a86d5..984f4549b516ecd42adb04886e9f2f7deecaf3da 100644 (file)
--- 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
 }