X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkImagePlanes.cxx;h=b62cd077c4d8223d1deb2509503f0a4024e31532;hb=b297410bce4fc9141ade21cbeca71061369b4932;hp=fca5f6a3c9f8531ad744d8a583f658c69a871f73;hpb=8a545acd568e25e30b6255191f69320265fae99d;p=bbtk.git diff --git a/packages/vtk/src/bbvtkImagePlanes.cxx b/packages/vtk/src/bbvtkImagePlanes.cxx index fca5f6a..b62cd07 100644 --- 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/06/23 15:58:56 $ - Version: $Revision: 1.27 $ + Date: $Date: 2010/01/21 16:03:19 $ + Version: $Revision: 1.31 $ =========================================================================*/ /* --------------------------------------------------------------------- @@ -38,6 +38,8 @@ #include "bbvtkPackage.h" #include "vtkCellPicker.h" #include "vtkProperty.h" +#include "vtkPolyData.h" + #include "vtkMetaImageWriter.h" #include "vtkPNGWriter.h" @@ -135,7 +137,7 @@ namespace bbvtk }; //================================================================ - //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.27 $"); + //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.31 $"); //================================================================ @@ -170,6 +172,8 @@ namespace bbvtk _imageReslicer = NULL; image=NULL; + _transform =NULL; + _matrix =NULL; } @@ -186,36 +190,14 @@ namespace bbvtk 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); @@ -225,6 +207,7 @@ namespace bbvtk bbSetOutputImageX(planeWidgetX->GetResliceOutput()); bbSetOutputImageY(planeWidgetY->GetResliceOutput()); bbSetOutputImageZ(planeWidgetZ->GetResliceOutput()); + bbSetInputInteractor(0); //bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput()); picker->UnRegister(NULL); @@ -269,12 +252,13 @@ void ImagePlanes::Process() 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; @@ -310,6 +294,8 @@ void ImagePlanes::Process() 0.5*(range[1]+range[0])); } + updateInteractor(); + } // UPDATE DES SORTIES bbGetOutputPlaneX()->GetResliceOutput()->Update(); @@ -320,36 +306,12 @@ void ImagePlanes::Process() std::vector pointsy = bbGetInputPointsY(); std::vector pointsz = bbGetInputPointsZ(); - std::cout<=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]; @@ -362,41 +324,103 @@ void ImagePlanes::Process() 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 "<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]GetOutput()->Update; + //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() { + } //----------------------------------------------------------------- @@ -409,7 +433,20 @@ void ImagePlanes::Process() { } - + 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; @@ -465,9 +502,9 @@ double* ImagePlanes::makeVector(double podouble0[3], double podouble1[3]){ 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;