]> Creatis software - bbtk.git/blobdiff - packages/vtk/src/bbvtkImagePlanes.cxx
Try to display planes (from Image Planes)
[bbtk.git] / packages / vtk / src / bbvtkImagePlanes.cxx
index c80aad8bfa36459e54991a59d7e2affac85febbb..79ab3850b81e61815cc83328d6940080d02b13cb 100644 (file)
@@ -1,20 +1,33 @@
-/*=========================================================================
-                                                                                
+/*=========================================================================                                                                               
   Program:   bbtk
   Module:    $RCSfile: bbvtkImagePlanes.cxx,v $
   Language:  C++
-  Date:      $Date: 2008/10/08 14:46:44 $
-  Version:   $Revision: 1.7 $
-                                                                                
-  Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
-  l'Image). All rights reserved. See Doc/License.txt or
-  http://www.creatis.insa-lyon.fr/Public/bbtk/License.html for details.
-                                                                                
-     This software is distributed WITHOUT ANY WARRANTY; without even
-     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
-     PURPOSE.  See the above copyright notices for more information.
-                                                                                
+  Date:      $Date: 2008/12/03 09:34:39 $
+  Version:   $Revision: 1.14 $
 =========================================================================*/
+
+/* ---------------------------------------------------------------------
+
+* Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale)
+* Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux
+*
+*  This software is governed by the CeCILL-B license under French law and 
+*  abiding by the rules of distribution of free software. You can  use, 
+*  modify and/ or redistribute the software under the terms of the CeCILL-B 
+*  license as circulated by CEA, CNRS and INRIA at the following URL 
+*  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html 
+*  or in the file LICENSE.txt.
+*
+*  As a counterpart to the access to the source code and  rights to copy,
+*  modify and redistribute granted by the license, users are provided only
+*  with a limited warranty  and the software's author,  the holder of the
+*  economic rights,  and the successive licensors  have only  limited
+*  liability. 
+*
+*  The fact that you are presently reading this means that you have had
+*  knowledge of the CeCILL-B license and that you accept its terms.
+* ------------------------------------------------------------------------ */                                                                         
+
 /**
  *  \file 
  *  \brief 
@@ -25,6 +38,7 @@
 #include "bbvtkPackage.h"
 #include "vtkCellPicker.h"
 #include "vtkProperty.h"
+#include "vtkMetaImageWriter.h"
 
 #include "bbstdCast.h"
 namespace bbstd
@@ -50,34 +64,38 @@ namespace bbvtk
 namespace bbvtk
 {
 
-
    BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,ImagePlanes)
    BBTK_BLACK_BOX_IMPLEMENTATION(ImagePlanes,bbtk::AtomicBlackBox);
 
    void ImagePlanes::bbUserConstructor() 
    { 
-     Init();
+     planeWidgetX = planeWidgetY = planeWidgetZ = 0;
+     imageX       = imageY       = imageZ       = 0;
      bbSetInputIn(0);
      std::vector<double> vect;
      vect.push_back(512);
      vect.push_back(256);
-     bbSetInputWindowLevel (vect);
+     bbSetInputWindowLevel (vect);  
+
+     bbSetModifiedStatus();
    }
    
-   void ImagePlanes::bbUserCopyConstructor() 
+   void ImagePlanes::bbUserCopyConstructor(bbtk::BlackBox::Pointer
    { 
-     Init();
+     planeWidgetX = planeWidgetY = planeWidgetZ = 0;
+     imageX       = imageY       = imageZ       = 0;
    }
 
    void ImagePlanes::Init() 
-   { 
-     
+   {  
+     if (planeWidgetX != 0) return;
+       
      // The shared picker enables us to use 3 planes at one time
      // and gets the picking order right
      vtkCellPicker* picker = vtkCellPicker::New();
      picker->SetTolerance(0.005);
   
-     // The 3 image plane widgets are used to probe the dataset.
+     // The 3 image plane widgets 
      planeWidgetX = vtkImagePlaneWidget::New();
      planeWidgetX->DisplayTextOn();
      planeWidgetX->SetPicker(picker);
@@ -93,9 +111,6 @@ namespace bbvtk
      prop2->SetColor(1, 1, 0);
      planeWidgetY->SetLookupTable(planeWidgetX->GetLookupTable());
 
-     // for the z-slice, turn off texture interpolation:
-     // interpolation is now nearest neighbour, to demonstrate
-     // cross-hair cursor snapping to pixel centers
      planeWidgetZ = vtkImagePlaneWidget::New();
      planeWidgetZ->DisplayTextOn();
      planeWidgetZ->SetPicker(picker);
@@ -109,20 +124,22 @@ namespace bbvtk
      bbSetOutputPlaneZ(planeWidgetZ);
 
      picker->UnRegister(NULL);
-     //        planeWidgetX->Register(NULL);
-     //        planeWidgetY->Register(NULL);
-     //        planeWidgetZ->Register(NULL);
+     
+     myCallbackPlane *_myCallback = myCallbackPlane::New();
+     _myCallback->SetCurrentBlackBox(this);
+     planeWidgetX->AddObserver(vtkCommand::InteractionEvent,_myCallback);
+     planeWidgetY->AddObserver(vtkCommand::InteractionEvent,_myCallback);
+     planeWidgetZ->AddObserver(vtkCommand::InteractionEvent,_myCallback);  
    }
 
 //---------------------------------------------------------------------
   void ImagePlanes::bbUserDestructor()
   {
-    //planeWidgetX->UnRegister(NULL);//Delete();
-    //planeWidgetY->UnRegister(NULL);//Delete();
-    //planeWidgetZ->UnRegister(NULL);//Delete();
-    planeWidgetX->Delete();
-    planeWidgetY->Delete();
-    planeWidgetZ->Delete();
+
+    if (planeWidgetX) planeWidgetX->UnRegister(NULL);
+    if (planeWidgetY) planeWidgetY->UnRegister(NULL);
+    if (planeWidgetZ) planeWidgetZ->UnRegister(NULL);
+
   }
   
 //---------------------------------------------------------------------  
@@ -130,24 +147,61 @@ namespace bbvtk
    {
      if (bbGetInputIn()!=0)
        {
+        Init();
+
         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);
+        
         planeWidgetX->SetInput(bbGetInputIn());
-        planeWidgetX->SetPlaneOrientationToXAxes();
-        planeWidgetX->SetSlicePosition((xMax+xMin)/2.); 
+        planeWidgetX->SetPlaneOrientationToXAxes();     
+        planeWidgetX->SetSlicePosition((xMax+xMin)/2.*xSpacing);
+        planeWidgetX->GetResliceOutput()->Update();
+        bbSetOutput("ImageX", planeWidgetX->GetResliceOutput());
+        
+        
+        vtkMetaImageWriter *writer = vtkMetaImageWriter::New();
+        
+        std::string newFilenameX("newFilenameX.mhd");
+        writer->SetInput(planeWidgetX->GetResliceOutput());
+        writer->SetFileName(newFilenameX.c_str());
+        writer->Update();
+        
+        //planeWidgetX->GetResliceOutput()->Print(std::cout);
 
         planeWidgetY->SetInput(bbGetInputIn());
         planeWidgetY->SetPlaneOrientationToYAxes();
-        planeWidgetY->SetSlicePosition((yMax+yMin)/2.); 
-
+        planeWidgetY->SetSlicePosition((yMax+yMin)/2.*ySpacing);
+        planeWidgetY->GetResliceOutput()->Update(); 
+        bbSetOutput("ImageY", planeWidgetY->GetResliceOutput());
+        //planeWidgetY->GetResliceOutput()->Print(std::cout);
+        
+        std::string newFilenameY("newFilenameY.mhd");
+        writer->SetInput(planeWidgetY->GetResliceOutput());
+        writer->SetFileName(newFilenameY.c_str());
+        writer->Update();
+                
         planeWidgetZ->SetInput(bbGetInputIn());
         planeWidgetZ->SetPlaneOrientationToZAxes();
-        planeWidgetZ->SetSlicePosition((zMax+zMin)/2.); 
-
+        planeWidgetZ->SetSlicePosition((zMax+zMin)/2.*zSpacing);
+        planeWidgetZ->GetResliceOutput()->Update(); 
+        bbSetOutput("ImageZ", planeWidgetZ->GetResliceOutput());
+       // planeWidgetZ->GetResliceOutput()->Print(std::cout);
+       
        // planeWidgetZ->SetWindowLevel(512,256);
-
+       
+        std::string newFilenameZ("newFilenameZ.mhd");
+        writer->SetInput(planeWidgetZ->GetResliceOutput());
+        writer->SetFileName(newFilenameZ.c_str());
+        writer->Update();
+       
          planeWidgetZ->SetWindowLevel(bbGetInputWindowLevel()[0],bbGetInputWindowLevel()[1]);
+       
+       // writer->vtkMetaImageWriter::Delete();
+        
        }
    }
 }//namespace bbtk