X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=vtk%2FvtkgdcmViewer.cxx;h=47fc46c69a6dca5bccdfae4106d52a0ed5e288f1;hb=9cfc0e09d292826980e3d63728d46b5c6b185bf4;hp=e0a17943fa76e999a7a0ac8f4875d7c034f5eb9e;hpb=c11dce1faf0053eb984aed57f01e6b7d29cbd1ab;p=gdcm.git

diff --git a/vtk/vtkgdcmViewer.cxx b/vtk/vtkgdcmViewer.cxx
index e0a17943..47fc46c6 100644
--- a/vtk/vtkgdcmViewer.cxx
+++ b/vtk/vtkgdcmViewer.cxx
@@ -1,7 +1,36 @@
-#include <iostream>
-
+/*=========================================================================
+                                                                                
+  Program:   gdcm
+  Module:    $RCSfile: vtkgdcmViewer.cxx,v $
+  Language:  C++
+  Date:      $Date: 2005/06/29 16:15:16 $
+  Version:   $Revision: 1.26 $
+                                                                                
+  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/Gdcm/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.
+                                                                                
+=========================================================================*/
+// This example illustrates how the vtkGdcmReader vtk class can be
+// used in order to:
+//  * produce a simple (vtk based) Dicom image STACK VIEWER.
+//  * dump the stack considered as a volume in a vtkStructuredPoints
+//    vtk file: the vtk gdcm wrappers can be seen as a simple way to convert
+//    a stack of Dicom images into a native vtk volume.
+//
+// Usage:
+//  * the filenames of the Dicom images constituting the stack should be
+//    given as command line arguments,
+//  * you can navigate through the stack by hitting any character key,
+//  * the produced vtk file is named "foo.vtk" (in the invocation directory).
+// 
+//----------------------------------------------------------------------------
 #include <vtkRenderWindowInteractor.h>
-#include <vtkImageViewer2.h>
+#include <vtkImageViewer.h>
 #include <vtkStructuredPoints.h>
 #include <vtkStructuredPointsWriter.h>
 #include <vtkCommand.h>
@@ -10,7 +39,7 @@
 #include <vtkLookupTable.h>
 
 #include "vtkGdcmReader.h"
-
+#include "gdcmDocument.h"  // for NO_SHADOWSEQ
 #ifndef vtkFloatingPointType
 #define vtkFloatingPointType float
 #endif
@@ -19,94 +48,107 @@
 // Callback for the interaction
 class vtkgdcmObserver : public vtkCommand
 {
-  public:
-  virtual char const *GetClassName() const { return "vtkgdcmObserver";}
-  static vtkgdcmObserver *New() 
-    { return new vtkgdcmObserver; }
-  vtkgdcmObserver()
-    {
+public:
+   virtual char const *GetClassName() const 
+   { 
+      return "vtkgdcmObserver";
+   }
+   static vtkgdcmObserver *New() 
+   { 
+      return new vtkgdcmObserver; 
+   }
+   vtkgdcmObserver()
+   {
       this->ImageViewer = NULL;
-    }
-  virtual void Execute(vtkObject *wdg, unsigned long event, void* calldata)
-    {
+   }
+   virtual void Execute(vtkObject *, unsigned long event, void* )
+   {
       if ( this->ImageViewer )
-        {
-          if ( event == vtkCommand::CharEvent )
-          {
+      {
+         if ( event == vtkCommand::CharEvent )
+         {
             int max = ImageViewer->GetWholeZMax();
             int slice = (ImageViewer->GetZSlice() + 1 ) % ++max;
             ImageViewer->SetZSlice( slice );
             ImageViewer->GetRenderer()->ResetCameraClippingRange();
             ImageViewer->Render();
-          }
-        }
-    }
-    vtkImageViewer2 *ImageViewer;
+         }
+      }
+   }
+   vtkImageViewer *ImageViewer;
 };
 
 
 int main(int argc, char *argv[])
 {
-  if( argc < 2 )
-    return 0;
+   if( argc < 2 )
+      return 0;
   
-  vtkGdcmReader *reader = vtkGdcmReader::New();
+   vtkGdcmReader *reader = vtkGdcmReader::New();
+   reader->AllowLookupTableOff();
 
-  if( argc == 2 )
-    reader->SetFileName( argv[1] );
-  else
-    for(int i=1; i< argc; i++)
-      reader->AddFileName( argv[i] );
+   if( argc == 2 )
+      reader->SetFileName( argv[1] );
+   else
+      for(int i=1; i< argc; i++)
+         reader->AddFileName( argv[i] );
 
-  reader->DebugOn();
-  reader->Update();
-  
-  //print debug info:
-  reader->GetOutput()->Print( std::cout );
-        
-  vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
-
-  vtkImageViewer2 *viewer = vtkImageViewer2::New();
-  if( reader->GetLookupTable() )
-  {
-    //convert to color:
-    vtkImageMapToColors *map = vtkImageMapToColors::New ();
-    map->SetInput (reader->GetOutput());
-    map->SetLookupTable (reader->GetLookupTable());
-    map->SetOutputFormatToRGB();
-    viewer->SetInput ( map->GetOutput() );
-    map->Delete();
-  }
-  else
-  {
-    viewer->SetInput ( reader->GetOutput() );
-  }
-  viewer->SetupInteractor (iren);
-  
-  vtkFloatingPointType *range = reader->GetOutput()->GetScalarRange();
-  viewer->SetColorWindow (range[1] - range[0]);
-  viewer->SetColorLevel (0.5 * (range[1] + range[0]));
-
-  // Here is where we setup the observer, 
-  vtkgdcmObserver *obs = vtkgdcmObserver::New();
-  obs->ImageViewer = viewer;
-  iren->AddObserver(vtkCommand::CharEvent,obs);
-  obs->Delete();
-
-  iren->Initialize();
-  iren->Start();
-
-  //if you wish you can export dicom to a vtk file  
-  vtkStructuredPointsWriter *writer = vtkStructuredPointsWriter::New();
-  writer->SetInput( reader->GetOutput());
-  writer->SetFileName( "foo.vtk" );
-  writer->SetFileTypeToBinary();
-  writer->Write();
-
-  reader->Delete();
-  iren->Delete();
-  viewer->Delete();
-  writer->Delete();
+// TODO : allow user to choose Load Mode
+   reader->SetLoadMode(NO_SHADOWSEQ);  
+   reader->Update();
+
+   //print debug info:
+   reader->GetOutput()->Print( cout );
+
+   vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
+
+   vtkImageViewer *viewer = vtkImageViewer::New();
+
+   if( reader->GetLookupTable() )
+   {
+      //convert to color:
+      vtkImageMapToColors *map = vtkImageMapToColors::New ();
+      map->SetInput (reader->GetOutput());
+      map->SetLookupTable (reader->GetLookupTable());
+      map->SetOutputFormatToRGB();
+      viewer->SetInput ( map->GetOutput() );
+      map->Delete();
+   }
+   else
+   {
+      vtkFloatingPointType *range = reader->GetOutput()->GetScalarRange();
+      viewer->SetColorLevel (0.5 * (range[1] + range[0]));
+      viewer->SetColorWindow (range[1] - range[0]);
+
+      viewer->SetInput ( reader->GetOutput() );
+   }
+   viewer->SetupInteractor (iren);
   
-  return 0;
+   //vtkFloatingPointType *range = reader->GetOutput()->GetScalarRange();
+   //viewer->SetColorWindow (range[1] - range[0]);
+   //viewer->SetColorLevel (0.5 * (range[1] + range[0]));
+
+   // Here is where we setup the observer, 
+   vtkgdcmObserver *obs = vtkgdcmObserver::New();
+   obs->ImageViewer = viewer;
+   iren->AddObserver(vtkCommand::CharEvent,obs);
+   obs->Delete();
+
+   //viewer->Render();
+   iren->Initialize();
+   iren->Start();
+
+   //if you wish you can export dicom to a vtk file  
+   vtkStructuredPointsWriter *writer = vtkStructuredPointsWriter::New();
+   writer->SetInput( reader->GetOutput());
+   writer->SetFileName( "foo.vtk" );
+   writer->SetFileTypeToBinary();
+   //writer->Write();
+
+   reader->Delete();
+   iren->Delete();
+   viewer->Delete();
+   writer->Delete();
+
+   return 0;
 }