Program: bbtk
Module: $RCSfile: bbvtkImagePlanes.cxx,v $
Language: C++
- Date: $Date: 2009/06/24 07:56:02 $
- Version: $Revision: 1.28 $
+ Date: $Date: 2010/03/07 13:18:41 $
+ Version: $Revision: 1.32 $
=========================================================================*/
/* ---------------------------------------------------------------------
#include "bbvtkPackage.h"
#include "vtkCellPicker.h"
#include "vtkProperty.h"
+#include "vtkPolyData.h"
#include "vtkMetaImageWriter.h"
};
//================================================================
- //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.28 $");
+ //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.32 $");
//================================================================
_imageReslicer = NULL;
image=NULL;
+ _transform =NULL;
+ _matrix =NULL;
}
picker->SetTolerance(0.005);
// The 3 image plane widgets
- vtkImagePlaneWidget* planeWidgetX = vtkImagePlaneWidget::New();
- planeWidgetX->DisplayTextOn();
- planeWidgetX->SetPicker(picker);
- planeWidgetX->SetKeyPressActivationValue('x');
- vtkProperty* prop1 = planeWidgetX->GetPlaneProperty();
- prop1->SetColor(1, 0, 0);
-
- vtkImagePlaneWidget* planeWidgetY = vtkImagePlaneWidget::New();
- planeWidgetY->DisplayTextOn();
- planeWidgetY->SetPicker(picker);
- planeWidgetY->SetKeyPressActivationValue('y');
- vtkProperty* prop2 = planeWidgetY->GetPlaneProperty();
- prop2->SetColor(1, 1, 0);
+ vtkImagePlaneWidget* planeWidgetX = GetPlaneWidget('x', 1, 0, 0, picker);
+ vtkImagePlaneWidget* planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker);
planeWidgetY->SetLookupTable(planeWidgetX->GetLookupTable());
- vtkImagePlaneWidget* planeWidgetZ = vtkImagePlaneWidget::New();
- planeWidgetZ->DisplayTextOn();
- planeWidgetZ->SetPicker(picker);
- planeWidgetZ->SetKeyPressActivationValue('z');
- vtkProperty* prop3 = planeWidgetZ->GetPlaneProperty();
- prop3->SetColor(0, 0, 1);
+ vtkImagePlaneWidget* planeWidgetZ = GetPlaneWidget('z', 0, 0, 1, picker);
planeWidgetZ->SetLookupTable(planeWidgetX->GetLookupTable());
- vtkImagePlaneWidget* planeWidget3Pts = vtkImagePlaneWidget::New();
- //vtkPlaneWidget* planeWidget3Pts = vtkPlaneWidget::New();
- planeWidget3Pts->DisplayTextOn();
- planeWidget3Pts->SetPicker(picker);
- planeWidget3Pts->SetKeyPressActivationValue('3');
- vtkProperty* prop3Pts = planeWidget3Pts->GetPlaneProperty();
- prop3Pts->SetColor(0, 1, 1);
+ vtkImagePlaneWidget* planeWidget3Pts = GetPlaneWidget('3', 0, 1, 1, picker);
planeWidget3Pts->SetLookupTable(planeWidgetX->GetLookupTable());
bbSetOutputPlaneX(planeWidgetX);
bbSetOutputImageX(planeWidgetX->GetResliceOutput());
bbSetOutputImageY(planeWidgetY->GetResliceOutput());
bbSetOutputImageZ(planeWidgetZ->GetResliceOutput());
+ bbSetInputInteractor(0);
//bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput());
picker->UnRegister(NULL);
if (bbGetInputIn()!=0)
{
+ int xMin, xMax, yMin, yMax, zMin, zMax;
+ bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
+
if ( image != bbGetInputIn()){//bbGetInputStatus("In") != bbtk::UPTODATE ){
// Input image has changed : reinitialize planes
- image = bbGetInputIn();
-
- int xMin, xMax, yMin, yMax, zMin, zMax;
- bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
+ image = bbGetInputIn();
+
// Initial values : center of the volume (in real world, not in pixels!)
double xSpacing, ySpacing, zSpacing;
0.5*(range[1]+range[0]));
}
+ updateInteractor();
+
}
// UPDATE DES SORTIES
bbGetOutputPlaneX()->GetResliceOutput()->Update();
std::vector<int> pointsy = bbGetInputPointsY();
std::vector<int> pointsz = bbGetInputPointsZ();
- std::cout<<pointsx.size()<<pointsy.size()<<pointsz.size()<<std::endl;
+ //std::cout<<pointsx.size()<<pointsy.size()<<pointsz.size()<<std::endl;
if(pointsx.size()==pointsy.size() && pointsx.size()==pointsz.size()&&pointsx.size()>=3){
- vtkImagePlaneWidget* plane3pts = (vtkImagePlaneWidget*)bbGetOutputPlane3Pts();
- //vtkPlaneWidget* plane3pts = (vtkPlaneWidget*)bbGetOutputPlane3Pts();
-
- plane3pts->SetInput(bbGetInputIn());
-
- //xSpacing = ySpacing = zSpacing = 1;
- double xSpacing, ySpacing, zSpacing;
- bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing);
-
- plane3pts->SetOrigin(pointsx[0]*xSpacing,pointsy[0]*ySpacing,pointsz[0]*zSpacing);
- plane3pts->SetPoint1(pointsx[1]*xSpacing,pointsy[1]*ySpacing,pointsz[1]*zSpacing);
- //plane3pts->SetPoint1((pointsx[1]-pointsx[0])*xSpacing,(pointsy[1]-pointsy[0])*ySpacing,(pointsz[1]-pointsz[0])*zSpacing);
- plane3pts->SetPoint2(pointsx[2]*xSpacing,pointsy[2]*ySpacing,pointsz[2]*zSpacing);
- //plane3pts->SetPoint2((pointsx[2]-pointsx[0])*xSpacing,(pointsy[2]-pointsy[0])*ySpacing,(pointsz[2]-pointsz[0])*zSpacing);
- plane3pts->GetResliceOutput()->Update();
-
- if (_imageReslicer==NULL){
- _imageReslicer = vtkImageReslice::New();
- _imageReslicer->SetOutputDimensionality(2);
- _imageReslicer->SetInterpolationModeToLinear();
- }
-
- _imageReslicer->SetInput( bbGetInputIn() );
- _imageReslicer->SetInformationInput(bbGetInputIn());
-
+ //Get the corresponding three points out of the vectors
double origin[3];
origin[0] = pointsx[0];
origin[1] = pointsy[0];
double point2[3];
point2[0]= pointsx[2];
point2[1]= pointsy[2];
- point2[2]= pointsz[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);
- /*std::cout<<"origin "<<origin[0]<<" "<<origin[1]<<" "<<origin[2]<<" "<<std::endl;
- std::cout<<"point1 "<<point1[0]<<" "<<point1[1]<<" "<<point1[2]<<" "<<std::endl;
- std::cout<<"point2 "<<point2[0]<<" "<<point2[1]<<" "<<point2[2]<<" "<<std::endl;
- std::cout<<"vect1 "<<vect1[0]<<" "<<vect1[1]<<" "<<vect1[2]<<" "<<std::endl;
- std::cout<<"vect2 "<<vect2[0]<<" "<<vect2[1]<<" "<<vect2[2]<<" "<<std::endl;
- std::cout<<"crossp "<<crossp[0]<<" "<<crossp[1]<<" "<<crossp[2]<<" "<<std::endl;*/
-
+ double *newx = getCrossProduct(vect2, crossp);
+
+ int ext[6],factor=0;
+ bbGetInputIn()->GetExtent(ext);
- _imageReslicer->SetResliceAxesDirectionCosines(vect1[0],vect1[1],vect1[2],
- vect2[0],vect2[1],vect2[2],
- crossp[0],crossp[1],crossp[2]);
- //_imageReslicer->SetResliceAxesOrigin(0,0,0);
- _imageReslicer->SetResliceAxesOrigin(origin[0],origin[1],origin[2]);
+ 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;
+ bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing);
+ plane3pts->SetOrigin(pointsx[0]*xSpacing,pointsy[0]*ySpacing,pointsz[0]*zSpacing);
+ plane3pts->SetPoint1((origin[0]+newx[0]*factor)*xSpacing,
+ (origin[1]+newx[1]*factor)*ySpacing,
+ (origin[2]+newx[2]*factor)*zSpacing);
+ plane3pts->SetPoint2((origin[0]+vect2[0]*factor)*xSpacing,
+ (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->SetInterpolationModeToLinear();
+ _transform = vtkTransform::New();
+ _matrix = vtkMatrix4x4::New();
+ }
+ _imageReslicer->SetInput( bbGetInputIn() );
+ _imageReslicer->SetInformationInput(bbGetInputIn());
+ //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();
+
bbSetOutputImage3Pts(_imageReslicer->GetOutput());
+
+
+
+ _matrix->Identity();
+
+ _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);
+
+
+
}
}
}
-
+void ImagePlanes::updateInteractor(){
+
+ vtkRenderWindowInteractor* interactor = bbGetInputInteractor();
+
+ if(interactor){
+ bbGetOutputPlaneX()->SetInteractor(interactor);
+ bbGetOutputPlaneX()->EnabledOn();
+ bbGetOutputPlaneY()->SetInteractor(interactor);
+ bbGetOutputPlaneY()->EnabledOn();
+ bbGetOutputPlaneZ()->SetInteractor(interactor);
+ bbGetOutputPlaneZ()->EnabledOn();
+ bbGetOutputPlane3Pts()->SetInteractor(interactor);
+ bbGetOutputPlane3Pts()->EnabledOn();
+ }
+}
//-----------------------------------------------------------------
void vtkImageDataPointerRelay::bbUserSetDefaultValues()
{
+
}
//-----------------------------------------------------------------
{
}
-
+ vtkImagePlaneWidget* ImagePlanes::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker)
+ {
+ vtkProperty* prop1 = 0;
+ vtkImagePlaneWidget* planeWidget = 0;
+
+ planeWidget = vtkImagePlaneWidget::New();
+ planeWidget->DisplayTextOn();
+ planeWidget->SetPicker(picker);
+ planeWidget->SetKeyPressActivationValue(activationkey);
+ prop1 = planeWidget->GetPlaneProperty();
+ prop1->SetColor(r, g, b);
+
+ return planeWidget;
+ }
double* ImagePlanes::getCrossProduct(double* vect0,double* vect1){
double* vectCross;
double *vect;
vect = new double[3];
- vect[0]= podouble0[0]-podouble1[0];
- vect[1]= podouble0[1]-podouble1[1];
- vect[2]= podouble0[2]-podouble1[2];
+ vect[0]= podouble1[0]-podouble0[0];
+ vect[1]= podouble1[1]-podouble0[1];
+ vect[2]= podouble1[2]-podouble0[2];
return vect;