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             std::cout << "Slice:" << slice << "," << max << std::endl;
  31             ImageViewer->SetZSlice( slice );
  32             ImageViewer->GetRenderer()->ResetCameraClippingRange();
  33             ImageViewer->Render();
  34           }
  35         }
  36     }
  37     vtkImageViewer2 *ImageViewer;
  38 };
  39
  40
  41 int main(int argc, char *argv[])
  42 {
  43
  44   vtkGdcmReader *reader = vtkGdcmReader::New();
  45   reader->SetFileName( argv[1] );
  46   reader->DebugOn();
  47   reader->Update();
  48
  49   //print debug info:
  50   reader->GetOutput()->Print( std::cout );
  51
  52   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
  53
  54   vtkImageViewer2 *viewer = vtkImageViewer2::New();
  55   viewer->SetInput ( reader->GetOutput() );
  56   viewer->SetupInteractor (iren);
  57
  58 //  float *range = reader->GetOutput()->GetScalarRange();
  59 //  viewer->SetColorWindow (range[1] - range[0]);
  60 //  viewer->SetColorLevel (0.5 * (range[1] + range[0]));
  61
  62   // Here is where we setup the observer,
  63   vtkgdcmObserver *obs = vtkgdcmObserver::New();
  64   obs->ImageViewer = viewer;
  65   iren->AddObserver(vtkCommand::CharEvent,obs);
  66   obs->Delete();
  67
  68   iren->Initialize();
  69   iren->Start();
  70
  71   //if you wish you can export dicom to a vtk file
  72   vtkStructuredPointsWriter *writer = vtkStructuredPointsWriter::New();
  73   writer->SetInput( reader->GetOutput());
  74   writer->SetFileName( "foo.vtk" );
  75   writer->SetFileTypeToBinary();
  76   writer->Write();
  77
  78   reader->Delete();
  79   iren->Delete();
  80   viewer->Delete();
  81   writer->Delete();
  82
  83   return 0;
  84 }