#include <iostream>
#include "gdcm.h"
#include "gdcmHeaderEntry.h"
#include "gdcmDicomDir.h"

#include <fstream>
#ifdef GDCM_NO_ANSI_STRING_STREAM
#  include <strstream>
#  define  ostringstream ostrstream
# else
#  include <sstream>
#endif
#include <string>
#include <iostream>

#define  ostringstream ostrstream
//extern "C"


int testDicomDir(int argc, char* argv[])
{  
   gdcmDicomDir *e1;
   ListTag::iterator deb , fin;
   
   ListDicomDirPatient::iterator  itPatient;
   ListDicomDirStudy::iterator    itStudy;
   ListDicomDirSerie::iterator    itSerie;
   ListDicomDirImage::iterator    itImage;
   TSKey v;
    
   std::string file; 
   if (argc > 1) 
      file = argv[1];    
   else {
      file += GDCM_DATA_ROOT;
      file += "/DICOMDIR";
    }

   e1 = new gdcmDicomDir(file.c_str());
   if (argc > 2) {
      int level = atoi(argv[2]);   
      e1->SetPrintLevel(level);
   }

   if(e1->GetDicomDirPatients().begin() == e1->GetDicomDirPatients().end() )
   {
      std::cout<<"Empty list"<<std::endl;
      return(1);
   }

// Simple examples of structure exploitation 
// DON'T REMOVE neither the folowing code, nor the commented out lines

/*
  cout << std::endl << std::endl  
	     << " = Liste des PATIENT ==========================================" 
	     << std::endl<< std::endl;	


   itPatient = e1->GetDicomDirPatients().begin();	
   while ( itPatient != e1->GetDicomDirPatients().end() ) {  // on degouline la liste de PATIENT
      std::cout << (*itPatient)->GetEntryByNumber(0x0010, 0x0010) << std::endl; // Patient's Name   
      itPatient ++;    
   }
*/

/*        
   cout << std::endl << std::endl  
        << " = Liste des PATIENT/STUDY ==========================================" 
        << std::endl<< std::endl;	

   itPatient = e1->GetDicomDirPatients().begin();	
   while ( itPatient != e1->GetDicomDirPatients().end() ) {  // on degouline la liste de PATIENT
      std::cout << (*itPatient)->GetEntryByNumber(0x0010, 0x0010) << std::endl; // Patient's Name 
      itStudy = ((*itPatient)->GetDicomDirStudies()).begin();	
      while (itStudy != (*itPatient)->GetDicomDirStudies().end() ) { // on degouline les STUDY de ce patient	
         std::cout << "--- "<< (*itStudy)->GetEntryByNumber(0x0008, 0x1030) << std::endl; // Study Description
         ++itStudy;		
      }
      itPatient ++;    
   }   
 */
 
   std::cout << std::endl << std::endl  
	     << " = Liste des PATIENT/STUDY/SERIE ==========================================" 
	     << std::endl<< std::endl;	
 
   itPatient = e1->GetDicomDirPatients().begin();	
   while ( itPatient != e1->GetDicomDirPatients().end() ) {  // on degouline la liste de PATIENT
       // Patient's Name, Patient ID 
      std::cout << "Pat.Name:[" << (*itPatient)->GetEntryByNumber(0x0010, 0x0010) <<"]";
      std::cout << " Pat.ID:[";
      std::cout << (*itPatient)->GetEntryByNumber(0x0010, 0x0020) << "]" << std::endl;
      itStudy = ((*itPatient)->GetDicomDirStudies()).begin();	
      while (itStudy != (*itPatient)->GetDicomDirStudies().end() ) { // on degouline les STUDY de ce patient	
         std::cout << "--- Stud.descr:["    << (*itStudy)->GetEntryByNumber(0x0008, 0x1030) << "]";// Study Description 
         std::cout << " Stud.ID:["<< (*itStudy)->GetEntryByNumber(0x0020, 0x0010);                 // Study ID
	 std::cout << "]" << std::endl;
         itSerie = ((*itStudy)->GetDicomDirSeries()).begin();
         while (itSerie != (*itStudy)->GetDicomDirSeries().end() ) { // on degouline les SERIES de cette study	   
            std::cout << "--- --- Ser.Descr:["<< (*itSerie)->GetEntryByNumber(0x0008, 0x103e)<< "]";// Serie Description
            std::cout << " Ser.nb:[" <<   (*itSerie)->GetEntryByNumber(0x0020, 0x0011);            // Serie number
            std::cout << "] Mod.:["    <<   (*itSerie)->GetEntryByNumber(0x0008, 0x0060) << "]";   // Modality
	    std::cout << std::endl;	    
            ++itSerie;   
         }
         ++itStudy;		
      }
      itPatient ++;    
   } 
   
    
 /*
   cout << std::endl << std::endl  
	     << " = Liste des PATIENT/STUDY/SERIE/IMAGE ===================================" 
	     << std::endl<< std::endl;
 
   itPatient = e1->GetDicomDirPatients().begin();	
   while ( itPatient != e1->GetDicomDirPatients().end() ) {  // on degouline la liste de PATIENT
      std::cout << (*itPatient)->GetEntryByNumber(0x0010, 0x0010) << std::endl; // Patient's Name
      itStudy = ((*itPatient)->GetDicomDirStudies()).begin();	
      while (itStudy != (*itPatient)->GetDicomDirStudies().end() ) { // on degouline les STUDY de ce patient	
         std::cout << "--- "<< (*itStudy)->GetEntryByNumber(0x0008, 0x1030) << std::endl; // Study Description
         itSerie = ((*itStudy)->GetDicomDirSeries()).begin();
         while (itSerie != (*itStudy)->GetDicomDirSeries().end() ) { // on degouline les SERIES de cette study	   
            std::cout << "--- --- "<< (*itSerie)->GetEntryByNumber(0x0008, 0x103e) << std::endl; // Serie Description
            itImage = ((*itSerie)->GetDicomDirImages()).begin();
            while (itImage != (*itSerie)->GetDicomDirImages().end() ) { // on degouline les SERIES de cette study	   
               std::cout << "--- --- --- "<< (*itImage)->GetEntryByNumber(0x0004, 0x1500) << std::endl; // File name
               ++itImage;   
            }
            ++itSerie;   
         }
         ++itStudy;		
      }  
      itPatient ++;    
   }   
 */  

   std::cout << std::endl << std::endl  
        << " = Contenu Complet du DICOMDIR ==========================================" 
        << std::endl<< std::endl;	
   e1->Print();
   
   std::cout<<std::flush;
   delete e1;

   return(0);
}