Add some verbosity to show the troubles with SerieHelper

[gdcm.git] / Example / exSerieHelper.cxx

diff --git a/Example/exSerieHelper.cxx b/Example/exSerieHelper.cxx
index 37ffee32fc30c36ac8fe8c4fa7a91298b59ed1b5..22296f670742d44a6400c9040a1ea542b68fbbdf 100644 (file)
--- a/Example/exSerieHelper.cxx
+++ b/Example/exSerieHelper.cxx
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: exSerieHelper.cxx,v $
   Language:  C++
-  Date:      $Date: 2007/05/30 15:13:25 $
-  Version:   $Revision: 1.9 $
+  Date:      $Date: 2007/10/30 11:45:01 $
+  Version:   $Revision: 1.15 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -17,6 +17,7 @@
 =========================================================================*/
 #include "gdcmSerieHelper.h"
 #include "gdcmFile.h"
+#include "gdcmDirList.h" // for FileList
 #include "gdcmDebug.h"
 #include <iostream>
 
@@ -37,6 +38,10 @@ int main(int argc, char *argv[])
 
   
    std::cout << "Dir Name :[" << dirName << "]" << std::endl;
+   //   
+   // Sometimes using only SerieHelper is not enought !
+   // See also exXcoherentFileSet
+   //
 
    s = GDCM_NAME_SPACE::SerieHelper::New();
    s->SetLoadMode(GDCM_NAME_SPACE::LD_ALL);     // Load everything for each File
@@ -52,45 +57,67 @@ int main(int argc, char *argv[])
    std::cout << " ---------------------------------------- Finish printing (1)"
              << std::endl;
 
+
+   GDCM_NAME_SPACE::FileList::const_iterator it;
+   GDCM_NAME_SPACE::FileList *l;
+   std::cout << std::endl << " ---------------------------------------- Recap"
+             << std::endl;  
+   l = s->GetFirstSingleSerieUIDFileSet();
+   while (l)
+   { 
+      it = l->begin();
+      std::cout << "SerieUID [" <<  (*it)->GetEntryString(0x0020,0x000e) <<"]   Serie Description ["
+                << (*it)->GetEntryString(0x0008,0x103e) << "] "  
+                << " : " << l->size() << " files" << std::endl;
+      l = s->GetNextSingleSerieUIDFileSet();
+   } 
+    std::cout << " ----------------------------------------End Recap"
+             << std::endl << std::endl;
+
    int nbFiles;
+   double zspacing = 0.;
    // For all the Single SerieUID Files Sets of the GDCM_NAME_SPACE::Serie
-   GDCM_NAME_SPACE::FileList *l = s->GetFirstSingleSerieUIDFileSet();
+   l = s->GetFirstSingleSerieUIDFileSet();
    while (l)
    { 
       nbFiles = l->size() ;
-      if ( l->size() > 3 ) // Why not ? Just an example, for testing
+      if ( l->size() > 5 ) // Why not ? Just an example, for testing
       {
-         std::cout << "Sort list : " << nbFiles << " long" << std::endl;
-         s->OrderFileList(l);  // sort the list
+         std::cout << "Sort list : " << nbFiles << " long" << std::endl; 
+ 
+         //---------------------------------------------------------
+         s->OrderFileList(l);  // sort the list (and compute ZSpacing !)
+         //---------------------------------------------------------
+         std::cout << "List after sorting : " << nbFiles << " long" << std::endl;
+           
+          zspacing = s->GetZSpacing();
+         // Just to show : GetZSpacing from a GDCM_NAME_SPACE::SerieHelper is right  
+         std::cout << "GetZSpacing() of sorted SingleSerieUIDFileSet "
+                   << "from GDCM_NAME_SPACE::SerieHelper: " << zspacing << std::endl;
+         std::cout << " ('-1' means all the files have the same position)" << std::endl;
+         int fileCount = 0;
+         for (std::vector<GDCM_NAME_SPACE::File* >::iterator it2 =  l->begin();
+                                            it2 != l->end();
+                                          ++it2)
+         {
+          // Just to show : GetZSpacing from a GDCM_NAME_SPACE::File may be different        
+             std::cout << (*it2)->GetFileName() << " -->  GetZSpacing() from GDCM_NAME_SPACE::File : " 
+                       << (*it2)->GetZSpacing() << std::endl; 
+           fileCount++;      
+         }
+         std::cout << "Iterate trough vector, nb of files : " << fileCount << std::endl;  
+
+         //break; // we only deal with the first one ... Why not ?
       }
       l = s->GetNextSingleSerieUIDFileSet();
    } 
-   std::cout << " -------------------------------------------- Finish sorting"
+   std::cout << std::endl
+             << " ------------------Prints all the Single SerieUID File Sets (sorted or not) -----"
              << std::endl;
    s->Print(); // Prints all the Single SerieUID File Sets (sorted or not)
    std::cout << " -------------------------------------------- Finish printing"
              << std::endl;
-  
-  // Just to show : GetZSpacing from a gdcm::SerieHelper is right  
-    std::cout << "GetZSpacing() from gdcm::SerieHelper: "    << s->GetZSpacing()    << std::endl;
-
-   // Only for the first Coherent File List 
-   // ( Why not ? Just an example, for testing )
-   // Display all the file names
-
-   std::string fileName; 
-   l = s->GetFirstSingleSerieUIDFileSet();
-   for (std::vector<GDCM_NAME_SPACE::File* >::iterator it =  l->begin();
-                                            it != l->end();
-                                          ++it)
-   {
-      fileName = (*it)->GetFileName();
-      std::cout << fileName << std::endl;
- // Just to show : GetZSpacing from a gdcm::File is wrong!        
-      std::cout << "       GetZSpacing() from gdcm::File : " << (*it)->GetZSpacing() << std::endl;      
-   } 
      
-
    s->Delete();
 
    return 0;