*** empty log message ***

[bbtk.git] / packages / vtk / src / bbvtkImagePlanes.cxx

diff --git a/packages/vtk/src/bbvtkImagePlanes.cxx b/packages/vtk/src/bbvtkImagePlanes.cxx
index 24d9f36a6243a5d205e6dfe44e76beb7cfc963ad..9f61f6ea0efc81d02221ef80ab3d25cd10549369 100644 (file)
--- a/packages/vtk/src/bbvtkImagePlanes.cxx
+++ b/packages/vtk/src/bbvtkImagePlanes.cxx
@@ -2,8 +2,8 @@
   Program:   bbtk
   Module:    $RCSfile: bbvtkImagePlanes.cxx,v $
   Language:  C++
-  Date:      $Date: 2009/07/02 07:17:33 $
-  Version:   $Revision: 1.29 $
+  Date:      $Date: 2009/09/01 13:59:48 $
+  Version:   $Revision: 1.30 $
 =========================================================================*/
 
 /* ---------------------------------------------------------------------
@@ -137,7 +137,7 @@ namespace bbvtk
  };
   //================================================================
 
-  //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.29 $");
+  //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.30 $");
 
   //================================================================
 
@@ -330,6 +330,7 @@ void ImagePlanes::Process()
                if(pointsx.size()==pointsy.size() && pointsx.size()==pointsz.size()&&pointsx.size()>=3){
 
 
+                       //Get the corresponding three points out of the vectors
                        double origin[3];
                        origin[0] = pointsx[0];
                        origin[1] = pointsy[0];
@@ -344,6 +345,7 @@ void ImagePlanes::Process()
                        point2[1]= pointsy[2];
                        point2[2]= pointsz[2];  
 
+                       //With the three points we create the corresponding X, Y and Z vectors all orthogonal to each other
                        double* vect1= getNormal(makeVector(origin, point1));
                        double* vect2= getNormal(makeVector(origin, point2));                           
                        double* crossp = getCrossProduct(vect1, vect2);
@@ -356,7 +358,7 @@ void ImagePlanes::Process()
                        factor = ext[0]<ext[3]? ext[3] : ext[0];
                        factor = factor<ext[5]? ext[5] : factor;
 
-
+       //for the plane widgets
                        vtkImagePlaneWidget* plane3pts = (vtkImagePlaneWidget*)bbGetOutputPlane3Pts();
                        plane3pts->SetInput(bbGetInputIn());                    
                        double xSpacing, ySpacing, zSpacing;
@@ -369,7 +371,7 @@ void ImagePlanes::Process()
                                                                        (origin[1]+vect2[1]*factor)*ySpacing,
                                                                        (origin[2]+vect2[2]*factor)*zSpacing);
                        plane3pts->GetResliceOutput()->Update();
-
+//To get the slice of image out of the selected volume
                        if (_imageReslicer==NULL){
                                _imageReslicer = vtkImageReslice::New();                                        
                                _imageReslicer->SetOutputDimensionality(2);
@@ -379,9 +381,10 @@ void ImagePlanes::Process()
                        }
                        _imageReslicer->SetInput( bbGetInputIn() );
                        _imageReslicer->SetInformationInput(bbGetInputIn());    
-                       _imageReslicer->SetResliceAxesDirectionCosines(newx[0],newx[1],newx[2],
-                                                                       vect2[0],vect2[1],vect2[2],
-                                                                       crossp[0],crossp[1],crossp[2]);                 
+                       //fill out the information with the created vectors and using the spacing of the image
+                       _imageReslicer->SetResliceAxesDirectionCosines(newx[0]*xSpacing,newx[1]*xSpacing,newx[2]*xSpacing,
+                                                                       vect2[0]*ySpacing,vect2[1]*ySpacing,vect2[2]*ySpacing,
+                                                                       crossp[0]*zSpacing,crossp[1]*zSpacing,crossp[2]*zSpacing);                      
                        _imageReslicer->SetResliceAxesOrigin(origin[0]*xSpacing,origin[1]*ySpacing,origin[2]*zSpacing);
                        _imageReslicer->GetOutput()->Update();
                        _imageReslicer->GetOutput()->UpdateInformation();
@@ -393,22 +396,23 @@ void ImagePlanes::Process()
                         
                        _matrix->Identity();    
 
-                       _matrix->SetElement(0,0,newx[0]);
-                       _matrix->SetElement(1,0,newx[1]);
-                       _matrix->SetElement(2,0,newx[2]);
-                       _matrix->SetElement(0,1,vect2[0]);
-                       _matrix->SetElement(1,1,vect2[1]);
-                       _matrix->SetElement(2,1,vect2[2]);
-                       _matrix->SetElement(0,2,crossp[0]);
-                       _matrix->SetElement(1,2,crossp[1]);
-                       _matrix->SetElement(2,2,crossp[2]);
-                       _matrix->SetElement(0,3,origin[0]);
-                       _matrix->SetElement(1,3,origin[1]);
-                       _matrix->SetElement(2,3,origin[2]);
+                       _matrix->SetElement(0,0,newx[0]*xSpacing);
+                       _matrix->SetElement(1,0,newx[1]*xSpacing);
+                       _matrix->SetElement(2,0,newx[2]*xSpacing);
+                       _matrix->SetElement(0,1,vect2[0]*ySpacing);
+                       _matrix->SetElement(1,1,vect2[1]*ySpacing);
+                       _matrix->SetElement(2,1,vect2[2]*ySpacing);
+                       _matrix->SetElement(0,2,crossp[0]*zSpacing);
+                       _matrix->SetElement(1,2,crossp[1]*zSpacing);
+                       _matrix->SetElement(2,2,crossp[2]*zSpacing);
+                       _matrix->SetElement(0,3,origin[0]*xSpacing);
+                       _matrix->SetElement(1,3,origin[1]*ySpacing);
+                       _matrix->SetElement(2,3,origin[2]*zSpacing);
 
                        
                        _transform->SetMatrix(_matrix);
 
+                       //set the transformation out to be used by other bbBoxes
                        bbSetOutputTransform3Pts((vtkLinearTransform*)_transform);