Testing/TestDicomDir.cxx

   1 /*=========================================================================
   2
   3   Program:   gdcm
   4   Module:    $RCSfile: TestDicomDir.cxx,v $
   5   Language:  C++
   6   Date:      $Date: 2006/05/31 16:51:41 $
   7   Version:   $Revision: 1.43 $
   8
   9   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
  10   l'Image). All rights reserved. See Doc/License.txt or
  11   http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
  12
  13      This software is distributed WITHOUT ANY WARRANTY; without even
  14      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  15      PURPOSE.  See the above copyright notices for more information.
  16
  17 =========================================================================*/
  18 #include "gdcmDocEntry.h"
  19 #include "gdcmDataEntry.h"
  20 #include "gdcmDicomDir.h"
  21 #include "gdcmDicomDirPatient.h"
  22 #include "gdcmDicomDirStudy.h"
  23 #include "gdcmDicomDirSerie.h"
  24 #include "gdcmDicomDirImage.h"
  25 #include "gdcmTS.h"
  26
  27 #include <iostream>
  28 #include <fstream>
  29
  30 #include <stdlib.h>
  31
  32 // check *all* the dicom elements (gdcm::DocEntry)
  33 // of this gdcm::DicomDirObject
  34 int CompareSQItem(gdcm::SQItem *pa1, gdcm::SQItem *pa2 )
  35 {
  36    gdcm::DocEntry *e1;
  37    gdcm::DocEntry *e2;
  38
  39    e2 = pa2->GetFirstEntry();
  40    while (!e2)
  41    {
  42       // locate the corresponding element in 'source' file
  43       e1 = pa1->GetDocEntry( e2->GetGroup(),e2->GetElement() );
  44
  45       // an element doesn't exist in origin file
  46       if (!e1)
  47       {
  48        std::cout << "DicomDir element " << std::hex
  49                  << e2->GetGroup() << "," <<e2->GetElement() << std::endl;
  50        return 1;
  51       }
  52       // skip SeqEntries (I don't want to deal with 'recursion pbs' here)
  53       if ( !dynamic_cast<gdcm::DataEntry *>(e1) ||
  54            !dynamic_cast<gdcm::DataEntry *>(e2) )
  55          continue;
  56
  57       // a value is read as GDCM_UNFOUND
  58       if ( ((gdcm::DataEntry *)e1)->GetString() == gdcm::GDCM_UNFOUND )
  59       {
  60          std::cout << "for gdcm source DicomDir : element (" << std::hex
  61                    << e1->GetGroup() << "," <<e1->GetElement()
  62                    << ") has values [" << gdcm::GDCM_UNFOUND << "]"
  63                    << std::endl;
  64          return 1;
  65       }
  66
  67       // values differ in source file and destination file
  68       if ( ((gdcm::DataEntry *)e1)->GetString() !=
  69            ((gdcm::DataEntry *)e2)->GetString() )
  70       {
  71
  72          // serious trouble : values differ in source and destination file
  73          std::cout << "for gdcm DicomDir element (" << std::hex
  74                    << e2->GetGroup() << "," <<e2->GetElement()
  75                    << ") values differ ["
  76                    << ((gdcm::DataEntry *)e1)->GetString() << "] != ["
  77                    << ((gdcm::DataEntry *)e2)->GetString() << "]"
  78                    << std::endl;
  79           return 1;
  80       }
  81    }
  82    return 0;
  83 }
  84
  85 int TestDicomDir(int argc, char *argv[])
  86 {
  87    gdcm::DicomDir *dicomdir;
  88
  89    gdcm::DicomDirPatient *pa1;
  90    gdcm::DicomDirStudy   *st1;
  91    gdcm::DicomDirSerie   *se1;
  92    gdcm::DicomDirImage   *im1;
  93
  94    gdcm::TSKey v;
  95
  96    std::string file;
  97    if (argc > 1)
  98       file = argv[1];
  99    else
 100    {
 101       file += GDCM_DATA_ROOT;
 102       file += "/DICOMDIR";
 103    }
 104
 105    dicomdir = gdcm::DicomDir::New();
 106    dicomdir->SetFileName(file);
 107    dicomdir->Load();
 108    if (argc > 2)
 109    {
 110       int level = atoi(argv[2]);
 111       dicomdir->SetPrintLevel(level);
 112    }
 113
 114    // Test if the DicomDir is readable
 115    if( !dicomdir->IsReadable() )
 116    {
 117       std::cout<<"          DicomDir '"<<file
 118                <<"' is not readable"<<std::endl
 119                <<"          ...Failed"<<std::endl;
 120
 121       dicomdir->Delete();
 122       return 1;
 123    }
 124    else
 125    {
 126       std::cout<<"          DicomDir '"<<file
 127                <<"' is readable"<<std::endl;
 128    }
 129
 130    // Test if the DicomDir contains any patient
 131    if( !dicomdir->GetFirstPatient() )
 132    {
 133       std::cout<<"          DicomDir '"<<file
 134                <<" has no patient"<<std::endl
 135                <<"          ...Failed"<<std::endl;
 136
 137       dicomdir->Delete();
 138       return 1;
 139    }
 140
 141    // step by step structure full exploitation
 142    std::cout << std::endl << std::endl
 143              << " = PATIENT/STUDY/SERIE/IMAGE List ============================"
 144              << std::endl<< std::endl;
 145
 146    pa1 = dicomdir->GetFirstPatient();
 147    while ( pa1 )
 148    {  // we process all the PATIENT of this DICOMDIR
 149       std::cout << pa1->GetEntryString(0x0010, 0x0010) << std::endl; // Patient's Name
 150
 151       st1 = pa1->GetFirstStudy();
 152       while ( st1 )
 153       { // we process all the STUDY of this patient
 154          std::cout << "--- "<< st1->GetEntryString(0x0008, 0x1030) // Study Description
 155          << std::endl;
 156          std::cout << " Stud.ID:[" << st1->GetEntryString(0x0020, 0x0010) // Study ID
 157          << "]";
 158
 159          se1 = st1->GetFirstSerie();
 160          while ( se1 )
 161          { // we process all the SERIES of this study
 162             std::cout << "--- --- "<< se1->GetEntryString(0x0008, 0x103e) << std::endl;      // Serie Description
 163             std::cout << " Ser.nb:["         <<  se1->GetEntryString(0x0020, 0x0011);        // Series number
 164             std::cout << "] Mod.:["          <<  se1->GetEntryString(0x0008, 0x0060) << "]"; // Modality
 165
 166             im1 = se1->GetFirstImage();
 167             while ( im1 ) { // we process all the IMAGE of this serie
 168                std::cout << "--- --- --- "<< im1->GetEntryString(0x0004, 0x1500) << std::endl; // File name
 169                im1 = se1->GetNextImage();
 170             }
 171             se1 = st1->GetNextSerie();
 172          }
 173          st1 = pa1->GetNextStudy();
 174       }
 175       pa1 = dicomdir->GetNextPatient();
 176    }
 177
 178    std::cout << std::endl << std::endl
 179              << " = DICOMDIR full content ===================================="
 180              << std::endl<< std::endl;
 181   // dicomdir->Print();
 182
 183    // ------------------------- second stage ---------------------------
 184
 185    // Write on disc what we read
 186    dicomdir->Write("NewDICOMDIR");
 187
 188    std::cout << std::endl << std::endl
 189              << "NewDICOMDIR written on disc ================================="
 190              << std::endl<< std::endl;
 191   // Read what we wrote
 192    gdcm::DicomDir *d2 = gdcm::DicomDir::New();
 193    d2->SetFileName("NewDICOMDIR");
 194    d2->Load();
 195    if (!d2->IsReadable())
 196    {
 197       std::cout << std::endl << std::endl
 198                 << "Read NewDicomDir from disc failed ========================"
 199                 << std::endl<< std::endl;
 200       return 1;
 201    }
 202    std::cout << std::endl << std::endl
 203              << "NewDICOMDIR successfully read from disc ================================="
 204              << std::endl<< std::endl;
 205
 206    gdcm::DicomDirPatient *pa2;
 207    gdcm::DicomDirStudy   *st2;
 208    gdcm::DicomDirSerie   *se2;
 209    gdcm::DicomDirImage   *im2;
 210
 211    pa1 = dicomdir->GetFirstPatient();
 212    pa2 = d2->GetFirstPatient();
 213
 214    if (!d2)
 215    {
 216       std::cout << "NewDICOMDIR contains no Patient ?!?" << std::endl;
 217       dicomdir->Delete();
 218       d2->Delete();
 219       return 1;
 220    }
 221
 222    while ( pa1 && pa2 )
 223    {  // we process all the PATIENT of this DICOMDIR
 224
 225       if ( CompareSQItem(pa2,pa1) == 1 )
 226       {
 227          dicomdir->Delete();
 228          d2->Delete();
 229          return 1;
 230       }
 231
 232       // just to allow human reader to be sure ...
 233       std::cout << pa2->GetEntryString(0x0010, 0x0010)
 234                 << std::endl; // Patient's Name
 235
 236       st1 = pa1->GetFirstStudy();
 237       st2 = pa2->GetFirstStudy();
 238
 239       while ( st1 && st2 )
 240       {
 241          if ( CompareSQItem(st2,st1) == 1 )
 242          {
 243             dicomdir->Delete();
 244             d2->Delete();
 245             return 1;
 246          }
 247
 248          // just to allow human reader to be sure ...
 249          std::cout << "--- "<< st2->GetEntryString(0x0008, 0x1030);
 250          // << std::endl;    // Study Description
 251          std::cout << " Stud.ID:["
 252                    << st2->GetEntryString(0x0020, 0x0010)
 253                    << "]" << std::endl; // Study ID
 254
 255          se1 = st1->GetFirstSerie();
 256          se2 = st2->GetFirstSerie();
 257
 258          while ( se1 && se2 )
 259          { // we process all the SERIE of this study
 260             if ( CompareSQItem(se2,se1) == 1 )
 261               return 1;
 262
 263             std::cout << "--- --- " << se2->GetEntryString(0x0008, 0x103e);      // Serie Description
 264             std::cout << " Ser.nb:["<< se2->GetEntryString(0x0020, 0x0011);        // Series number
 265             std::cout << "] Mod.:[" << se2->GetEntryString(0x0008, 0x0060) << "]" << std::endl; // Modality
 266             im1 = se1->GetFirstImage();
 267             im2 = se2->GetFirstImage();
 268
 269             while ( im1 && im2 ) // we process all the IMAGE of this serie
 270             {
 271                if ( CompareSQItem(im2,im1) == 1 )
 272                {
 273                   dicomdir->Delete();
 274                   d2->Delete();
 275                   return 1;
 276                }
 277
 278                im1 = se1->GetNextImage();
 279                im2 = se2->GetNextImage();
 280             }
 281             se1 = st1->GetNextSerie();
 282             se2 = st2->GetNextSerie();
 283          }
 284          st1 = pa1->GetNextStudy();
 285          st2 = pa2->GetNextStudy();
 286       }
 287       pa1 = dicomdir->GetNextPatient();
 288       pa2 = dicomdir->GetNextPatient();
 289    }
 290
 291    std::cout << std::flush;
 292    dicomdir->Delete();
 293    d2->Delete();
 294
 295    return 0;
 296 }