X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkImagePlanes.cxx;h=426189b031cf5d428edce3cc1b4cb5235d4d22d3;hb=2825c6ab65c1333b6236d54bb157b315097e6fb5;hp=92daa652fad7dca3429fe42ab4b561694f1e2eb4;hpb=b3225bf9915c36f8e79d9a328b525cde6eed9247;p=bbtk.git diff --git a/packages/vtk/src/bbvtkImagePlanes.cxx b/packages/vtk/src/bbvtkImagePlanes.cxx index 92daa65..426189b 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: 2010/03/07 13:18:41 $ - Version: $Revision: 1.32 $ + Date: $Date: 2012/07/09 13:02:42 $ + Version: $Revision: 1.37 $ =========================================================================*/ /* --------------------------------------------------------------------- @@ -65,7 +65,10 @@ #include "bbstdRelay.h" #include "vtkObjectFactory.h" +#include "vtkImageFlip.h" +#include "vtkImageReslice.h" +#include "vtkImageChangeInformation.h" namespace bbstd { @@ -92,7 +95,6 @@ namespace bbtk namespace bbvtk { - //==================================================================== // Add the specialized adaptors to the package typedef vtkImagePlaneWidget* I; @@ -137,15 +139,13 @@ namespace bbvtk }; //================================================================ - //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.32 $"); + //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.37 $"); //================================================================ BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,ImagePlanes) BBTK_BLACK_BOX_IMPLEMENTATION(ImagePlanes,bbtk::AtomicBlackBox); - - void ImagePlanes::bbUserSetDefaultValues() { bbSetOutputPlaneX(0); @@ -161,26 +161,25 @@ namespace bbvtk bbSetInputWindowLevel (vect); mVtkCallback = 0; - - std::vector vectpoints; + std::vector vectpoints; bbSetOutputPlane3Pts(0); - bbSetOutputImage3Pts(0); - bbSetInputPointsX(vectpoints); - bbSetInputPointsY(vectpoints); - bbSetInputPointsZ(vectpoints); - - _imageReslicer = NULL; - image=NULL; - _transform =NULL; - _matrix =NULL; + bbSetOutputImage3Pts(0); + bbSetInputPointsX(vectpoints); + bbSetInputPointsY(vectpoints); + bbSetInputPointsZ(vectpoints); + + _imageReslicer = NULL; + image = NULL; + _transform = NULL; + _matrix = NULL; + + bbSetInputInterpolation(1); } - void ImagePlanes::bbUserInitializeProcessing() { - /// CREATION DES WIDGETS if (bbGetOutputPlaneX() != 0) return; @@ -191,48 +190,75 @@ namespace bbvtk // The 3 image plane widgets vtkImagePlaneWidget* planeWidgetX = GetPlaneWidget('x', 1, 0, 0, picker); - vtkImagePlaneWidget* planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker); + vtkImagePlaneWidget* planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker); planeWidgetY->SetLookupTable(planeWidgetX->GetLookupTable()); vtkImagePlaneWidget* planeWidgetZ = GetPlaneWidget('z', 0, 0, 1, picker); planeWidgetZ->SetLookupTable(planeWidgetX->GetLookupTable()); - vtkImagePlaneWidget* planeWidget3Pts = GetPlaneWidget('3', 0, 1, 1, picker); + vtkImagePlaneWidget* planeWidget3Pts = GetPlaneWidget('3', 0, 1, 1, picker); planeWidget3Pts->SetLookupTable(planeWidgetX->GetLookupTable()); + // EED MPR view orientation correction.. + vtkImageFlip *flipYFilter = vtkImageFlip::New(); + flipYFilter->SetFilteredAxis(1); // flip y axis + flipYFilter->SetInput( planeWidgetX->GetResliceOutput() ); + flipYFilter->Update(); + + vtkImageChangeInformation *image = vtkImageChangeInformation::New(); + image->SetInput( planeWidgetY->GetResliceOutput() ); + image->SetOutputSpacing( 1,1,1 ); + image->CenterImageOn(); + image->Update(); + _imageTransform = vtkTransform::New(); + vtkImageReslice *slicer =vtkImageReslice::New(); + slicer->SetInput( image->GetOutput() ); + slicer->SetInformationInput( image->GetOutput() ); + slicer->SetResliceTransform( _imageTransform ); + slicer->SetOutputOrigin(0 , 0 , 0 ); + slicer->SetInterpolationModeToNearestNeighbor(); + slicer->Update(); + vtkImageChangeInformation *imageResult = vtkImageChangeInformation::New(); + imageResult->SetInput( slicer->GetOutput() ); + double spc[3]; + planeWidgetY->GetResliceOutput()->GetSpacing(spc); + imageResult->SetOutputSpacing( spc[1], spc[0], spc[2] ); + imageResult->SetOutputOrigin( 0,0,0 ); + + bbSetOutputPlaneX(planeWidgetX); bbSetOutputPlaneY(planeWidgetY); bbSetOutputPlaneZ(planeWidgetZ); - bbSetOutputPlane3Pts(planeWidget3Pts); - bbSetOutputImageX(planeWidgetX->GetResliceOutput()); - bbSetOutputImageY(planeWidgetY->GetResliceOutput()); + bbSetOutputPlane3Pts(planeWidget3Pts); + bbSetOutputImageX( flipYFilter->GetOutput() ); // EED MPR view orientation correction.. + bbSetOutputImageY( imageResult->GetOutput() ); // EED MPR view orientation correction.. bbSetOutputImageZ(planeWidgetZ->GetResliceOutput()); - bbSetInputInteractor(0); - //bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput()); + bbSetInputInteractor(0); + //bbSetOutputImage3Pts(planeWidget3Pts->GetResliceOutput()); - picker->UnRegister(NULL); + if(picker != 0) + picker->UnRegister(NULL); mVtkCallback = VtkCallbackType::New(); mVtkCallback->SetBlackBox(this); planeWidgetX->AddObserver(vtkCommand::InteractionEvent,mVtkCallback); planeWidgetY->AddObserver(vtkCommand::InteractionEvent,mVtkCallback); - planeWidgetZ->AddObserver(vtkCommand::InteractionEvent,mVtkCallback); + planeWidgetZ->AddObserver(vtkCommand::InteractionEvent,mVtkCallback); + } //--------------------------------------------------------------------- void ImagePlanes::bbUserFinalizeProcessing() { - + if (bbGetOutputPlaneX()) { /* - bbGetOutputPlaneX()->RemoveObserver(mVtkCallback); bbGetOutputPlaneY()->RemoveObserver(mVtkCallback); bbGetOutputPlaneZ()->RemoveObserver(mVtkCallback); - - + bbGetOutputPlaneX()->Delete(); bbGetOutputPlaneY()->Delete(); bbGetOutputPlaneZ()->Delete(); @@ -241,26 +267,35 @@ namespace bbvtk //bbGetOutputPlaneX()->SetInput(NULL); //bbGetOutputPlaneY()->SetInput(NULL); //bbGetOutputPlaneZ()->SetInput(NULL); - + } - } //--------------------------------------------------------------------- -void ImagePlanes::Process() -{ + void ImagePlanes::Process() + { + if (bbGetInputIn()!=0) { - - int xMin, xMax, yMin, yMax, zMin, zMax; - bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax); - + int dim[3]; + int ext[6]; + bbGetOutputPlaneX()->GetResliceOutput()->GetWholeExtent(ext); + dim[0] = ext[1]-ext[0]+1; + dim[1] = ext[3]-ext[2]+1; + dim[2] = ext[5]-ext[4]+1; + _imageTransform->Identity(); + _imageTransform->PostMultiply(); + _imageTransform->Translate( (int)(-(dim[0]/2)) , (int)(-(dim[1]/2)) ,0); + _imageTransform->RotateZ(90); + + 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!) + int xMin, xMax, yMin, yMax, zMin, zMax; + bbGetInputIn()->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax); double xSpacing, ySpacing, zSpacing; bbGetInputIn()->GetSpacing(xSpacing, ySpacing, zSpacing); @@ -268,12 +303,10 @@ void ImagePlanes::Process() 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); @@ -286,30 +319,41 @@ void ImagePlanes::Process() { bbGetOutputPlaneZ()->SetWindowLevel(bbGetInputWindowLevel()[0], bbGetInputWindowLevel()[1]); - } - else - { + } else { double *range = image->GetScalarRange(); - bbGetOutputPlaneZ()->SetWindowLevel(range[1] - range[0], - 0.5*(range[1]+range[0])); - } - + bbGetOutputPlaneZ()->SetWindowLevel(range[1]-range[0], + 0.5*(range[1]+range[0])); + } // windowlevel updateInteractor(); - - } + + dim[0] = xMax-xMin+1; + dim[1] = yMax-yMin+1; + dim[2] = zMax-zMin+1; + _imageTransform->Identity(); + _imageTransform->PostMultiply(); + _imageTransform->Translate( (int)(-(dim[0]/2)*(1/xSpacing)) , (int)(-(dim[2]/2)*(1/zSpacing)) ,0); + _imageTransform->RotateZ(90); + + } // image + // UPDATE DES SORTIES + bbGetOutputPlaneX()->SetResliceInterpolate( bbGetInputInterpolation() ); + bbGetOutputPlaneY()->SetResliceInterpolate( bbGetInputInterpolation() ); + bbGetOutputPlaneZ()->SetResliceInterpolate( bbGetInputInterpolation() ); + + bbGetOutputPlaneX()->GetResliceOutput()->Update(); bbGetOutputPlaneY()->GetResliceOutput()->Update(); bbGetOutputPlaneZ()->GetResliceOutput()->Update(); - + std::vector pointsx = bbGetInputPointsX(); std::vector pointsy = bbGetInputPointsY(); std::vector pointsz = bbGetInputPointsZ(); //std::cout<=3){ - + if (pointsx.size()==pointsy.size() && pointsx.size()==pointsz.size()&&pointsx.size()>=3) + { //Get the corresponding three points out of the vectors double origin[3]; @@ -332,7 +376,7 @@ void ImagePlanes::Process() double* crossp = getCrossProduct(vect1, vect2); double *newx = getCrossProduct(vect2, crossp); - + int ext[6],factor=0; bbGetInputIn()->GetExtent(ext); @@ -356,10 +400,10 @@ void ImagePlanes::Process() if (_imageReslicer==NULL){ _imageReslicer = vtkImageReslice::New(); _imageReslicer->SetOutputDimensionality(2); - _imageReslicer->SetInterpolationModeToLinear(); _transform = vtkTransform::New(); _matrix = vtkMatrix4x4::New(); } + _imageReslicer->SetInterpolationMode( bbGetInputInterpolation() ); _imageReslicer->SetInput( bbGetInputIn() ); _imageReslicer->SetInformationInput(bbGetInputIn()); //fill out the information with the created vectors and using the spacing of the image @@ -370,11 +414,8 @@ void ImagePlanes::Process() _imageReslicer->GetOutput()->Update(); _imageReslicer->GetOutput()->UpdateInformation(); - bbSetOutputImage3Pts(_imageReslicer->GetOutput()); - - _matrix->Identity(); _matrix->SetElement(0,0,newx[0]*xSpacing); @@ -390,19 +431,15 @@ void ImagePlanes::Process() _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); - - - - } - } -} + bbSetOutputTransform3Pts((vtkLinearTransform*)_transform); + } // pointsx pointsy pointsz + } // bbGetInputIn + } -void ImagePlanes::updateInteractor(){ + void ImagePlanes::updateInteractor(){ vtkRenderWindowInteractor* interactor = bbGetInputInteractor(); @@ -416,25 +453,25 @@ void ImagePlanes::updateInteractor(){ bbGetOutputPlane3Pts()->SetInteractor(interactor); bbGetOutputPlane3Pts()->EnabledOn(); } -} + } //----------------------------------------------------------------- - void vtkImageDataPointerRelay::bbUserSetDefaultValues() + void vtkImageDataPointerRelay::bbUserSetDefaultValues() { } //----------------------------------------------------------------- - void vtkImageDataPointerRelay::bbUserInitializeProcessing() + void vtkImageDataPointerRelay::bbUserInitializeProcessing() { } //----------------------------------------------------------------- - void vtkImageDataPointerRelay::bbUserFinalizeProcessing() + void vtkImageDataPointerRelay::bbUserFinalizeProcessing() { } - vtkImagePlaneWidget* ImagePlanes::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker) - { + vtkImagePlaneWidget* ImagePlanes::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker) + { vtkProperty* prop1 = 0; vtkImagePlaneWidget* planeWidget = 0; @@ -446,9 +483,9 @@ void ImagePlanes::updateInteractor(){ prop1->SetColor(r, g, b); return planeWidget; - } + } -double* ImagePlanes::getCrossProduct(double* vect0,double* vect1){ + double* ImagePlanes::getCrossProduct(double* vect0,double* vect1){ double* vectCross; vectCross = new double[3]; vectCross[0] = vect0[1]*vect1[2]-(vect0[2]*vect1[1]); @@ -456,32 +493,27 @@ double* ImagePlanes::getCrossProduct(double* vect0,double* vect1){ 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 ImagePlanes::getMagnitud(double* vect){ double mag; - mag = sqrt(pow(vect[0],2) + pow(vect[1],2) + pow(vect[2],2)); - - std::cout<<"mag "<