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 $
=========================================================================*/
/* ---------------------------------------------------------------------
};
//================================================================
- //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.29 $");
+ //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.30 $");
//================================================================
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];
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);
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;
(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);
}
_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();
_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);