X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkImagePlanes.cxx;h=92daa652fad7dca3429fe42ab4b561694f1e2eb4;hb=b3225bf9915c36f8e79d9a328b525cde6eed9247;hp=7f0008e41c91c2e9f9f0aefa9c09ba89850f0ae1;hpb=8c690fedb5ea897b6a69584e1e237cb5dbad97a5;p=bbtk.git diff --git a/packages/vtk/src/bbvtkImagePlanes.cxx b/packages/vtk/src/bbvtkImagePlanes.cxx index 7f0008e..92daa65 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/05/28 08:12:16 $ - Version: $Revision: 1.25 $ + Date: $Date: 2010/03/07 13:18:41 $ + Version: $Revision: 1.32 $ =========================================================================*/ /* --------------------------------------------------------------------- @@ -38,6 +38,8 @@ #include "bbvtkPackage.h" #include "vtkCellPicker.h" #include "vtkProperty.h" +#include "vtkPolyData.h" + #include "vtkMetaImageWriter.h" #include "vtkPNGWriter.h" @@ -58,6 +60,8 @@ #include "vtkPlaneWidget.h" +#include + #include "bbstdRelay.h" #include "vtkObjectFactory.h" @@ -133,7 +137,7 @@ namespace bbvtk }; //================================================================ - //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.25 $"); + //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.32 $"); //================================================================ @@ -156,6 +160,20 @@ namespace bbvtk vect.push_back(0); bbSetInputWindowLevel (vect); mVtkCallback = 0; + + + std::vector vectpoints; + + bbSetOutputPlane3Pts(0); + bbSetOutputImage3Pts(0); + bbSetInputPointsX(vectpoints); + bbSetInputPointsY(vectpoints); + bbSetInputPointsZ(vectpoints); + + _imageReslicer = NULL; + image=NULL; + _transform =NULL; + _matrix =NULL; } @@ -172,35 +190,25 @@ 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 = GetPlaneWidget('3', 0, 1, 1, picker); + planeWidget3Pts->SetLookupTable(planeWidgetX->GetLookupTable()); + bbSetOutputPlaneX(planeWidgetX); bbSetOutputPlaneY(planeWidgetY); bbSetOutputPlaneZ(planeWidgetZ); + bbSetOutputPlane3Pts(planeWidget3Pts); bbSetOutputImageX(planeWidgetX->GetResliceOutput()); bbSetOutputImageY(planeWidgetY->GetResliceOutput()); bbSetOutputImageZ(planeWidgetZ->GetResliceOutput()); + bbSetInputInteractor(0); + //bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput()); picker->UnRegister(NULL); @@ -239,59 +247,180 @@ namespace bbvtk } //--------------------------------------------------------------------- - void ImagePlanes::Process() - { - if (bbGetInputIn()!=0) - { - - if ( 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); - - // Initial values : center of the volume (in real world, not in pixels!) - double xSpacing, ySpacing, zSpacing; - bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing); - - bbGetOutputPlaneX()->SetInput(bbGetInputIn()); - bbGetOutputPlaneX()->SetPlaneOrientationToXAxes(); - bbGetOutputPlaneX()->SetSlicePosition((xMax+xMin)/2.*xSpacing); - - bbGetOutputPlaneY()->SetInput(bbGetInputIn()); - bbGetOutputPlaneY()->SetPlaneOrientationToYAxes(); - bbGetOutputPlaneY()->SetSlicePosition((yMax+yMin)/2.*ySpacing); - - bbGetOutputPlaneZ()->SetInput(bbGetInputIn()); - bbGetOutputPlaneZ()->SetPlaneOrientationToZAxes(); - bbGetOutputPlaneZ()->SetSlicePosition((zMax+zMin)/2.*zSpacing); - - if (bbGetInputWindowLevel()[0]!=0) - { - bbGetOutputPlaneZ()->SetWindowLevel(bbGetInputWindowLevel()[0], - bbGetInputWindowLevel()[1]); - } - else - { - double *range = image->GetScalarRange(); - bbGetOutputPlaneZ()->SetWindowLevel(range[1] - range[0], - 0.5*(range[1]+range[0])); - } - } - - // UPDATE DES SORTIES - bbGetOutputPlaneX()->GetResliceOutput()->Update(); - bbGetOutputPlaneY()->GetResliceOutput()->Update(); - bbGetOutputPlaneZ()->GetResliceOutput()->Update(); - } - } - +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(); + + + // Initial values : center of the volume (in real world, not in pixels!) + double xSpacing, ySpacing, zSpacing; + bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing); + + bbGetOutputPlaneX()->SetInput(bbGetInputIn()); + bbGetOutputPlaneX()->SetPlaneOrientationToXAxes(); + bbGetOutputPlaneX()->SetSlicePosition((xMax+xMin)/2.*xSpacing); + + + // bbGetOutputPlaneX()->SetOrigin( 58*xSpacing , 80*ySpacing , 82*zSpacing ); + // bbGetOutputPlaneX()->SetPoint1( 0*xSpacing, 146*ySpacing, 186*zSpacing); + // bbGetOutputPlaneX()->SetPoint2( 126*xSpacing, 146*ySpacing, 0*zSpacing); + + + bbGetOutputPlaneY()->SetInput(bbGetInputIn()); + bbGetOutputPlaneY()->SetPlaneOrientationToYAxes(); + bbGetOutputPlaneY()->SetSlicePosition((yMax+yMin)/2.*ySpacing); + + bbGetOutputPlaneZ()->SetInput(bbGetInputIn()); + bbGetOutputPlaneZ()->SetPlaneOrientationToZAxes(); + bbGetOutputPlaneZ()->SetSlicePosition((zMax+zMin)/2.*zSpacing); + + if (bbGetInputWindowLevel()[0]!=0) + { + bbGetOutputPlaneZ()->SetWindowLevel(bbGetInputWindowLevel()[0], + bbGetInputWindowLevel()[1]); + } + else + { + double *range = image->GetScalarRange(); + bbGetOutputPlaneZ()->SetWindowLevel(range[1] - range[0], + 0.5*(range[1]+range[0])); + } + + updateInteractor(); + + } + // UPDATE DES SORTIES + bbGetOutputPlaneX()->GetResliceOutput()->Update(); + bbGetOutputPlaneY()->GetResliceOutput()->Update(); + bbGetOutputPlaneZ()->GetResliceOutput()->Update(); + + std::vector pointsx = bbGetInputPointsX(); + std::vector pointsy = bbGetInputPointsY(); + std::vector pointsz = bbGetInputPointsZ(); + + //std::cout<=3){ + + + //Get the corresponding three points out of the vectors + double origin[3]; + origin[0] = pointsx[0]; + origin[1] = pointsy[0]; + origin[2] = pointsz[0]; + + double point1[3]; + point1[0] = pointsx[1]; + point1[1] = pointsy[1]; + point1[2] = pointsz[1]; + double point2[3]; + point2[0]= pointsx[2]; + 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); + + double *newx = getCrossProduct(vect2, crossp); + + int ext[6],factor=0; + bbGetInputIn()->GetExtent(ext); + + factor = ext[0]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() { + } //----------------------------------------------------------------- @@ -304,9 +433,85 @@ namespace bbvtk { } + 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; + vectCross = new double[3]; + vectCross[0] = vect0[1]*vect1[2]-(vect0[2]*vect1[1]); + vectCross[1] = -(vect0[0]*vect1[2]-(vect0[2]*vect1[0])); + vectCross[2] = vect0[0]*vect1[1]-(vect0[1]*vect1[0]); + + return vectCross; +} +/** +** Returns the magnitud of the given vector +**/ +double ImagePlanes::getMagnitud(double* vect){ + + double mag; + + mag = sqrt(pow(vect[0],2) + pow(vect[1],2) + pow(vect[2],2)); + + std::cout<<"mag "<