X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkImagePlanes.cxx;h=c31a38d2583e85fdc3059f9b9f2a1ab97db0af54;hb=731d9ee262e286a24e3ea2776d8e914950840a2c;hp=cd5ba230af6e795378b67065ccd8c90b29aab53f;hpb=17ac9460f5ea5600d1065fbe3ee4ad826d7bb4fc;p=bbtk.git

diff --git a/packages/vtk/src/bbvtkImagePlanes.cxx b/packages/vtk/src/bbvtkImagePlanes.cxx
index cd5ba23..c31a38d 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/15 14:58:01 $
-  Version:   $Revision: 1.24 $
+  Date:      $Date: 2010/08/18 12:36:45 $
+  Version:   $Revision: 1.35 $
 =========================================================================*/
 
 /* ---------------------------------------------------------------------
@@ -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 <vtkImplicitPlaneWidget.h>
+
 #include "bbstdRelay.h"
 
 #include "vtkObjectFactory.h"
@@ -88,7 +92,6 @@ namespace bbtk
 namespace bbvtk
 {
 
-
   //====================================================================
   // Add the specialized adaptors to the package
   typedef vtkImagePlaneWidget* I;
@@ -133,15 +136,13 @@ namespace bbvtk
  };
   //================================================================
 
-  //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.24 $");
+  //vtkCxxRevisionMacro(ImagePlanes::VtkCallbackType, "$Revision: 1.35 $");
 
   //================================================================
 
    BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,ImagePlanes)
    BBTK_BLACK_BOX_IMPLEMENTATION(ImagePlanes,bbtk::AtomicBlackBox);
 
-	
-	
    void ImagePlanes::bbUserSetDefaultValues() 
    { 
      bbSetOutputPlaneX(0);
@@ -156,13 +157,24 @@ namespace bbvtk
      vect.push_back(0);
      bbSetInputWindowLevel (vect);  
      mVtkCallback = 0;
+
+     std::vector<int> vectpoints;
+
+     bbSetOutputPlane3Pts(0);
+     bbSetOutputImage3Pts(0);
+     bbSetInputPointsX(vectpoints);
+     bbSetInputPointsY(vectpoints);
+     bbSetInputPointsZ(vectpoints);
+
+     _imageReslicer = NULL;
+     image          = NULL;
+     _transform     = NULL;
+     _matrix        = NULL;
    }
-   
 
 	
    void ImagePlanes::bbUserInitializeProcessing() 
    {  
-
      /// CREATION DES WIDGETS
      if (bbGetOutputPlaneX() != 0) return;
        
@@ -172,141 +184,304 @@ 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);
+	 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()->SetInput(NULL);
-	//	if (mVtkCallback) bbGetOutputPlaneX()->RemoveObserver(mVtkCallback);
+
+	/*
+	  bbGetOutputPlaneX()->RemoveObserver(mVtkCallback);
+	  bbGetOutputPlaneY()->RemoveObserver(mVtkCallback);
+	  bbGetOutputPlaneZ()->RemoveObserver(mVtkCallback);
+
 	bbGetOutputPlaneX()->Delete();
-      }
-    if (bbGetOutputPlaneY()) 
-      {
-	bbGetOutputPlaneY()->SetInput(NULL);
-	//	if (mVtkCallback) bbGetOutputPlaneY()->RemoveObserver(mVtkCallback);
 	bbGetOutputPlaneY()->Delete();
-      }
-    if (bbGetOutputPlaneZ()) 
-      {
-	bbGetOutputPlaneZ()->SetInput(NULL);
-	//	if (mVtkCallback) bbGetOutputPlaneZ()->RemoveObserver(mVtkCallback);
 	bbGetOutputPlaneZ()->Delete();
-      }
-    if (mVtkCallback) 
-      {
-	//	mVtkCallback->Delete();
+	mVtkCallback->Delete();
+	*/
+	//bbGetOutputPlaneX()->SetInput(NULL);
+	//bbGetOutputPlaneY()->SetInput(NULL);
+	//bbGetOutputPlaneZ()->SetInput(NULL);
+
       }
   }
   
 //---------------------------------------------------------------------  
-   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<int> pointsx = bbGetInputPointsX();
+		std::vector<int> pointsy = bbGetInputPointsY();
+		std::vector<int> pointsz = bbGetInputPointsZ();
+
+		//std::cout<<pointsx.size()<<pointsy.size()<<pointsz.size()<<std::endl;
+
+		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];
+			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]<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()
+  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)
+  {
+		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 "<<mag <<std::endl;
+	return mag;
+  }
+/**
+**	returns the unitary vector of the given vector
+**	u = 1/|vect| . vect
+**/
+  double* ImagePlanes::getNormal(double* vect){
+
+	double* vectnorm;
+	double mag = getMagnitud(vect);
+
+	vectnorm = new double[3];
+
+	if(mag!=0){
+		vectnorm[0] = vect[0]/mag;
+		vectnorm[1] = vect[1]/mag;
+		vectnorm[2] = vect[2]/mag;
+	}else{
+		vectnorm[0] = 0;
+		vectnorm[1] = 0;
+		vectnorm[2] = 0;
+	}
+
+	return vectnorm;
+  }
+
+  double* ImagePlanes::makeVector(double podouble0[3], double podouble1[3]){
+	double *vect;
+	vect = new double[3];
+
+	vect[0]= podouble1[0]-podouble0[0];
+	vect[1]= podouble1[1]-podouble0[1];
+	vect[2]= podouble1[2]-podouble0[2];
+
+	return vect;
+  }
 	
 }//namespace bbtk
 
 #endif // _USE_VTK_
-