]> Creatis software - bbtk.git/blob - packages/vtk/src/bbvtkTransform.cxx
Clean code
[bbtk.git] / packages / vtk / src / bbvtkTransform.cxx
1 /*
2  # ---------------------------------------------------------------------
3  #
4  # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
5  #                        pour la SantÈ)
6  # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
7  # Previous Authors : Laurent Guigues, Jean-Pierre Roux
8  # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9  #
10  #  This software is governed by the CeCILL-B license under French law and
11  #  abiding by the rules of distribution of free software. You can  use,
12  #  modify and/ or redistribute the software under the terms of the CeCILL-B
13  #  license as circulated by CEA, CNRS and INRIA at the following URL
14  #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
15  #  or in the file LICENSE.txt.
16  #
17  #  As a counterpart to the access to the source code and  rights to copy,
18  #  modify and redistribute granted by the license, users are provided only
19  #  with a limited warranty  and the software's author,  the holder of the
20  #  economic rights,  and the successive licensors  have only  limited
21  #  liability.
22  #
23  #  The fact that you are presently reading this means that you have had
24  #  knowledge of the CeCILL-B license and that you accept its terms.
25  # ------------------------------------------------------------------------ */
26
27
28 #include "bbvtkTransform.h"
29 #include "bbvtkPackage.h"
30 #include "vtkMath.h"
31 namespace bbvtk
32 {
33   BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,Transform);
34   BBTK_BLACK_BOX_IMPLEMENTATION(Transform,bbtk::AtomicBlackBox);
35
36         // --------------------------------------------------------------       
37
38         void Transform::bbUserSetDefaultValues()
39         {
40         bbSetInputIn(NULL);
41                 bbSetInputInverse(false);
42
43                 vecScale.push_back(1);  // scale x
44                 vecScale.push_back(1);  // scale y
45                 vecScale.push_back(1);  // scale z
46                 bbSetInputScale(vecScale);
47
48                 vecRotateWXYZ.push_back(0);  //angle 
49                 vecRotateWXYZ.push_back(1);  //vx
50                 vecRotateWXYZ.push_back(0);  //vy
51                 vecRotateWXYZ.push_back(0);  //vz
52                 bbSetInputRotateWXYZ(vecRotateWXYZ);
53
54                 vecTranslate.push_back(0);  //tx
55                 vecTranslate.push_back(0);  //ty
56                 vecTranslate.push_back(0);  //tz
57                 bbSetInputTranslate(vecTranslate);
58
59                 vecSpacing.push_back(1);  //spacing x
60                 vecSpacing.push_back(1);  //spacing y
61                 vecSpacing.push_back(1);  //spacing z
62                 bbSetInputSpacing(vecSpacing);
63         
64                 result = NULL;
65         }
66
67         // --------------------------------------------------------------       
68         void Transform::bbUserInitializeProcessing()
69         {
70                 bbUserFinalizeProcessing();
71                 result = vtkTransform::New();
72                 result->Update();
73                 resultInverse = vtkTransform::New();
74                 resultInverse->Update();
75         }
76
77         // --------------------------------------------------------------       
78         void Transform::bbUserFinalizeProcessing()
79         {
80                 if (result!=NULL)
81                 {
82                         result->Delete();
83                         result=NULL;
84                 }
85         }
86
87
88         // --------------------------------------------------------------       
89   void Transform::Process()
90   {
91           bbUserInitializeProcessing();
92           if (bbGetInputIn()!=NULL)
93       {
94                   result->Concatenate( bbGetInputIn()->GetMatrix() );
95       } // if In
96           if ((bbGetInputTranslate().size()>=3) && (bbGetInputSpacing().size()>=3))
97       {
98                   double tx = bbGetInputTranslate()[0] * bbGetInputSpacing()[0];
99                   double ty = bbGetInputTranslate()[1] * bbGetInputSpacing()[1];
100                   double tz = bbGetInputTranslate()[2] * bbGetInputSpacing()[2];
101                   result->Translate(tx,ty,tz);
102       }  // if translate spacing  size >=3
103           if (bbGetInputScale().size()>=3)
104       {
105                   result->Scale(bbGetInputScale()[0], bbGetInputScale()[1], bbGetInputScale()[2]);
106       }  // if scale size >=3
107           if (bbGetInputRotateWXYZ().size()>=4)
108       {
109                   result->RotateWXYZ(bbGetInputRotateWXYZ()[0],bbGetInputRotateWXYZ()[1], bbGetInputRotateWXYZ()[2], bbGetInputRotateWXYZ()[3]);
110       } // if rotation size >=4
111  
112       if (bbGetInputRotateToNormal().size()>=3)
113       {
114           double v1[3],v2[3],v3[3];
115           v1[0]             = 0;
116           v1[1]             = 0;
117           v1[2]             = 1;
118           v2[0]             = bbGetInputRotateToNormal()[0];
119           v2[1]             = bbGetInputRotateToNormal()[1];
120           v2[2]             = bbGetInputRotateToNormal()[2];
121           
122           double mag= sqrt( v2[0]*v2[0] + v2[1]*v2[1]+ v2[2]*v2[2] ) ;
123                   v2[0] = v2[0] / mag;
124                   v2[1] = v2[1] / mag;
125                   v2[2] = v2[2] / mag;
126           vtkMath *vtkmath  = vtkMath::New();
127           double ang        = acos ( vtkmath->Dot(v1,v2)  );
128           ang               = vtkmath->DegreesFromRadians( ang );
129           vtkmath->Cross( v1,v2,v3 );
130           vtkmath->Delete();
131           result->RotateWXYZ( ang  , v3[0], v3[1], v3[2] );
132       } // if rotation size >=4
133
134           if (bbGetInputInverse()==false)
135           {
136                 bbSetOutputOut(result);
137           } else {
138                 vtkMatrix4x4 *matrix;
139                 matrix=vtkMatrix4x4::New();
140                 result->GetInverse(matrix);
141                 resultInverse->SetMatrix( matrix );
142                 bbSetOutputOut(resultInverse);
143           }  // if Inverse
144   }
145   
146 }// EO namespace bbvtk