vtk/vtkgdcmViewer.cxx

   1 #include <vtkRenderWindowInteractor.h>
   2 #include <vtkImageViewer2.h>
   3 #include <vtkStructuredPoints.h>
   4 #include <vtkStructuredPointsWriter.h>
   5 #include <vtkCommand.h>
   6 #include <vtkRenderer.h>
   7
   8 #include "vtkGdcmReader.h"
   9
  10 //----------------------------------------------------------------------------
  11 // Callback for the interaction
  12 class vtkgdcmObserver : public vtkCommand
  13 {
  14   public:
  15   virtual char const *GetClassName() const { return "vtkgdcmObserver";}
  16   static vtkgdcmObserver *New()
  17     { return new vtkgdcmObserver; }
  18   vtkgdcmObserver()
  19     {
  20       this->ImageViewer = NULL;
  21     }
  22   virtual void Execute(vtkObject *wdg, unsigned long event, void* calldata)
  23     {
  24       if ( this->ImageViewer )
  25         {
  26           if ( event == vtkCommand::CharEvent )
  27           {
  28             int max = ImageViewer->GetWholeZMax();
  29             int slice = (ImageViewer->GetZSlice() + 1 ) % ++max;
  30             ImageViewer->SetZSlice( slice );
  31             ImageViewer->GetRenderer()->ResetCameraClippingRange();
  32             ImageViewer->Render();
  33           }
  34         }
  35     }
  36     vtkImageViewer2 *ImageViewer;
  37 };
  38
  39
  40 int main(int argc, char *argv[])
  41 {
  42   if( argc < 2 )
  43     return 0;
  44
  45   vtkGdcmReader *reader = vtkGdcmReader::New();
  46
  47   if( argc == 2 )
  48     reader->SetFileName( argv[1] );
  49   else
  50     for(int i=1; i< argc; i++)
  51       reader->AddFileName( argv[i] );
  52
  53   reader->DebugOn();
  54   reader->Update();
  55
  56   //print debug info:
  57   reader->GetOutput()->Print( std::cout );
  58
  59   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
  60
  61   vtkImageViewer2 *viewer = vtkImageViewer2::New();
  62   viewer->SetInput ( reader->GetOutput() );
  63   viewer->SetupInteractor (iren);
  64
  65 //  float *range = reader->GetOutput()->GetScalarRange();
  66 //  viewer->SetColorWindow (range[1] - range[0]);
  67 //  viewer->SetColorLevel (0.5 * (range[1] + range[0]));
  68
  69   // Here is where we setup the observer,
  70   vtkgdcmObserver *obs = vtkgdcmObserver::New();
  71   obs->ImageViewer = viewer;
  72   iren->AddObserver(vtkCommand::CharEvent,obs);
  73   obs->Delete();
  74
  75   iren->Initialize();
  76   iren->Start();
  77
  78   //if you wish you can export dicom to a vtk file
  79   vtkStructuredPointsWriter *writer = vtkStructuredPointsWriter::New();
  80   writer->SetInput( reader->GetOutput());
  81   writer->SetFileName( "foo.vtk" );
  82   writer->SetFileTypeToBinary();
  83   writer->Write();
  84
  85   reader->Delete();
  86   iren->Delete();
  87   viewer->Delete();
  88   writer->Delete();
  89
  90   return 0;
  91 }