vtk/vtkgdcmViewer.cxx

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