X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FPrintFile.cxx;h=1801f17e80e4294c0f5e0a4e8a76727742496263;hb=2e1f797ea8181180ebba6a48f0a4463bf43acedc;hp=b2eb2dec424c8eb2a6be387e06e9ad9a4d480c5a;hpb=3910d4f8047161daeb3a0992d00e641f1d11a668;p=gdcm.git

diff --git a/Example/PrintFile.cxx b/Example/PrintFile.cxx
index b2eb2dec..1801f17e 100644
--- a/Example/PrintFile.cxx
+++ b/Example/PrintFile.cxx
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: PrintFile.cxx,v $
   Language:  C++
-  Date:      $Date: 2005/06/29 16:00:13 $
-  Version:   $Revision: 1.42 $
+  Date:      $Date: 2005/08/29 12:25:37 $
+  Version:   $Revision: 1.48 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -16,27 +16,29 @@
                                                                                 
 =========================================================================*/
 #include "gdcmFile.h"
+#include "gdcmDocument.h"
 #include "gdcmSeqEntry.h"
 #include "gdcmSQItem.h"
 #include "gdcmBinEntry.h"
 
 #include "gdcmFileHelper.h"
 #include "gdcmDebug.h"
+#include "gdcmDirList.h"
 
 #include "gdcmArgMgr.h"
 
 #include <iostream>
 
 
-void ShowLutData(gdcm::File *e1);
+void ShowLutData(gdcm::File *f);
 
-void ShowLutData(gdcm::File *e1)
+void ShowLutData(gdcm::File *f)
 {
      // Nothing is written yet to get LUT Data user friendly
      // The following is to be moved into a PixelReadConvert method
      // Let here, waiting for a clever idea on the way to do it.
   
-      gdcm::SeqEntry *modLutSeq = e1->GetSeqEntry(0x0028,0x3000);
+      gdcm::SeqEntry *modLutSeq = f->GetSeqEntry(0x0028,0x3000);
       if ( modLutSeq !=0 )
       {
          gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
@@ -71,7 +73,7 @@ void ShowLutData(gdcm::File *e1)
             gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006); 
             if ( b != 0 )
             { 
-               int BitsAllocated = e1->GetBitsAllocated();
+               int BitsAllocated = f->GetBitsAllocated();
                if ( BitsAllocated <= 8 )
                { 
                   int mult;
@@ -114,57 +116,77 @@ void ShowLutData(gdcm::File *e1)
       else
          std::cout << "No LUT Data SeqEntry (0x0028,0x3000) found " 
                    << std::endl;
-
    }
 
 int main(int argc, char *argv[])
 {
 
    START_USAGE(usage)
-   " \n PrintFile : \n",
-   " Display the header of a ACR-NEMA/PAPYRUS/DICOM File",
-   " usage: PrintFile filein=fileName [level=n] [noshadow] [noseq] [debug] ",
-   "        level = 0,1,2 : depending on the amount of details user wants to see",
-   "        noshadow : user doesn't want to load Private groups (odd number)",
-   "        noseq    : user doesn't want to load Sequences ",
-   "        debug    : user wants to run the program in 'debug mode' ",
-   "        showlut  : user wants to display the Palette Color (as an int array)",
+   " \n PrintFile : \n                                                        ",
+   " Display the header of a ACR-NEMA/PAPYRUS/DICOM File                      ",
+   " usage: PrintFile {filein=inputFileName|dirin=inputDirectoryName}[level=n]",
+   "                       [forceload=listOfElementsToForceLoad]              ",
+   "                       [ { [noshadowseq] | [noshadow][noseq] } ] [debug]  ",
+   "      level = 0,1,2 : depending on the amount of details user wants to see",
+   "      listOfElementsToForceLOad : group-elem,g2-e2,... (in hexa, no space)",
+   "                                of Elements to load whatever their length ",
+   "        noshadowseq: user doesn't want to load Private Sequences          ",
+   "        noshadow   : user doesn't want to load Private groups (odd number)",
+   "        noseq      : user doesn't want to load Sequences                  ",
+   "        debug      : user wants to run the program in 'debug mode'        ",
+   "        showlut :user wants to display the Palette Color (as an int array)",
    FINISH_USAGE
 
    // Initialize Arguments Manager   
    gdcm::ArgMgr *am= new gdcm::ArgMgr(argc, argv);
   
-   if (argc == 1)
+   if (argc == 1 || am->ArgMgrDefined("usage") )
    {
       am->ArgMgrUsage(usage); // Display 'usage'
       delete am;
       return 0;
    }
 
-   char *fileName = am->ArgMgrWantString("filein",usage);
+   char *fileName = am->ArgMgrGetString("filein",(char *)0);
+   char *dirName  = am->ArgMgrGetString("dirin",(char *)0);
+
+   if ( (fileName == 0 && dirName == 0)
+        ||
+        (fileName != 0 && dirName != 0) )
+   {
+       std::cout <<std::endl
+                 << "Either 'filein=' or 'dirin=' must be present;" 
+                 << std::endl << "Not both" << std::endl;
+       am->ArgMgrUsage(usage); // Display 'usage'  
+       delete am;
+       return 0;
+ }
 
-   int loadMode;
+   if (am->ArgMgrDefined("debug"))
+      gdcm::Debug::DebugOn();
+ 
+   int loadMode = 0x00000000;
    if ( am->ArgMgrDefined("noshadowseq") )
-      loadMode = NO_SHADOWSEQ;
-   else if ( am->ArgMgrDefined("noshadow") && am->ArgMgrDefined("noseq") )
-      loadMode = NO_SEQ | NO_SHADOW;  
-   else if ( am->ArgMgrDefined("noshadow") )
-      loadMode = NO_SHADOW;
-   else if ( am->ArgMgrDefined("noseq") )
-      loadMode = NO_SEQ;
-   else
-      loadMode = 0;
+      loadMode |= NO_SHADOWSEQ;
+   else 
+   {
+   if ( am->ArgMgrDefined("noshadow") )
+         loadMode |= NO_SHADOW;
+      if ( am->ArgMgrDefined("noseq") )
+         loadMode |= NO_SEQ;
+   }
 
    int level = am->ArgMgrGetInt("level", 2);
 
-   if (am->ArgMgrDefined("debug"))
-      gdcm::Debug::DebugOn();
+   int forceLoadNb;
+   uint16_t *elemsToForceLoad 
+                        = am->ArgMgrGetXInt16Enum("forceload", &forceLoadNb);
 
    bool showlut = ( 0 != am->ArgMgrDefined("SHOWLUT") );
- 
+
    /* if unused Param we give up */
    if ( am->ArgMgrPrintUnusedLabels() )
-   { 
+   {
       am->ArgMgrUsage(usage);
       delete am;
       return 0;
@@ -173,150 +195,225 @@ int main(int argc, char *argv[])
    delete am;  // we don't need Argument Manager any longer
 
    // ----------- End Arguments Manager ---------
- 
-   // gdcm::File::IsReadable() is no usable here, because we deal with
-   // any kind of gdcm-Parsable *document* 
-   // not only gdcm::File (as opposed to gdcm::DicomDir)
 
-   gdcm::File *e1 = new gdcm::File();
-   e1->SetLoadMode(loadMode);
+   if ( fileName != 0 ) // ====== Deal with a single file ======
+   { 
+      // gdcm::File::IsReadable() is no usable here, because we deal with
+      // any kind of gdcm-Parsable *document* 
+      // not only gdcm::File (as opposed to gdcm::DicomDir)
 
-   bool res = e1->Load( fileName );
-   if ( !res )
-   {
-      delete e1;
-      return 0;
-   }
+      gdcm::File *f = new gdcm::File();
+      f->SetLoadMode(loadMode);
+      f->SetFileName( fileName );
 
-   gdcm::FileHelper *f1 = new gdcm::FileHelper(e1);
-   f1->SetPrintLevel( level );
-
-   f1->Print();   
-
-   std::cout << "\n\n" << std::endl; 
-
-   std::cout <<std::endl;
-   std::cout <<" dataSize    " << f1->GetImageDataSize()    << std::endl;
-   std::cout <<" dataSizeRaw " << f1->GetImageDataRawSize() << std::endl;
-
-   int nX,nY,nZ,sPP,planarConfig;
-   std::string pixelType;
-   nX=e1->GetXSize();
-   nY=e1->GetYSize();
-   nZ=e1->GetZSize();
-   std::cout << " DIMX=" << nX << " DIMY=" << nY << " DIMZ=" << nZ << std::endl;
-
-   pixelType    = e1->GetPixelType();
-   sPP          = e1->GetSamplesPerPixel();
-   planarConfig = e1->GetPlanarConfiguration();
-   
-   std::cout << " pixelType= ["            << pixelType 
-             << "] SamplesPerPixel= ["     << sPP
-             << "] PlanarConfiguration= [" << planarConfig 
-             << "] "<< std::endl 
-             << " PhotometricInterpretation= [" 
-                                << e1->GetEntryValue(0x0028,0x0004)
-             << "] "<< std::endl;
-
-   int numberOfScalarComponents=e1->GetNumberOfScalarComponents();
-   std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents 
-             <<std::endl
-             << " LUT = " << (e1->HasLUT() ? "TRUE" : "FALSE")
-             << std::endl;
-
-   if ( e1->GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) 
-   {
-      std::cout << "Transfer Syntax not loaded. " << std::endl
-                << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE"
+      for (int ri=0; ri<forceLoadNb; ri++)
+      {
+         f->AddForceLoadElement((uint32_t)elemsToForceLoad[2*ri], 
+                                (uint32_t)elemsToForceLoad[2*ri+1] ); 
+      }
+
+      bool res = f->Load();
+
+      if ( !res )
+      {
+         std::cout << "Cannot process file [" << fileName << "]" << std::endl;
+         std::cout << "Either it doesn't exist, or it's read protected " 
+                   << std::endl;
+         std::cout << "or it's not a Dicom File, or its 'header' is bugged" 
+                   << std::endl;
+         std::cout << "use 'PrintFile filein=... debug' to try to guess the pb"
+                   << std::endl;
+         delete f;
+         return 0;
+      }
+
+      gdcm::FileHelper *fh = new gdcm::FileHelper(f);
+      fh->SetPrintLevel( level );
+
+      fh->Print();
+
+      std::cout << "\n\n" << std::endl; 
+
+      std::cout <<std::endl;
+      std::cout <<" dataSize    " << fh->GetImageDataSize()    << std::endl;
+      std::cout <<" dataSizeRaw " << fh->GetImageDataRawSize() << std::endl;
+
+      int nX,nY,nZ,sPP,planarConfig;
+      std::string pixelType;
+      nX=f->GetXSize();
+      nY=f->GetYSize();
+      nZ=f->GetZSize();
+      std::cout << " DIMX=" << nX << " DIMY=" << nY << " DIMZ=" << nZ << std::endl;
+
+      pixelType    = f->GetPixelType();
+      sPP          = f->GetSamplesPerPixel();
+      planarConfig = f->GetPlanarConfiguration();
+
+      std::cout << " pixelType= ["            << pixelType 
+                << "] SamplesPerPixel= ["     << sPP
+                << "] PlanarConfiguration= [" << planarConfig 
+                << "] "<< std::endl 
+                << " PhotometricInterpretation= [" 
+                                << f->GetEntryValue(0x0028,0x0004)
+                << "] "<< std::endl;
+
+      int numberOfScalarComponents=f->GetNumberOfScalarComponents();
+      std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents 
+                <<std::endl
+                << " LUT = " << (f->HasLUT() ? "TRUE" : "FALSE")
                 << std::endl;
-      return 0;
-   }
+
+      if ( f->GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) 
+      {
+         std::cout << "Transfer Syntax not loaded. " << std::endl
+                   << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE"
+                << std::endl;
+         return 0;
+      }
   
-   std::string transferSyntaxName = e1->GetTransferSyntaxName();
-   std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" 
-             << std::endl;
-   std::cout << " SwapCode= " << e1->GetSwapCode() << std::endl;
- 
-   // Display the LUT as an int array (for debugging purpose)
-   if ( e1->HasLUT() && showlut )
-   {
-      uint8_t* lutrgba = f1->GetLutRGBA();
-      if ( lutrgba == 0 )
+      std::string transferSyntaxName = f->GetTransferSyntaxName();
+      std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" 
+                << std::endl;
+      std::cout << " SwapCode= " << f->GetSwapCode() << std::endl;
+
+      //std::cout << "\n\n" << std::endl; 
+      //std::cout << "X spacing " << f->GetXSpacing() << std::endl;
+      //std::cout << "Y spacing " << f->GetYSpacing() << std::endl;
+      //std::cout << "Z spacing " << f->GetZSpacing() << std::endl;
+
+      // Display the LUT as an int array (for debugging purpose)
+      if ( f->HasLUT() && showlut )
       {
-         std::cout << "Lut RGBA (Palette Color) not built " << std::endl;
+         uint8_t* lutrgba = fh->GetLutRGBA();
+         if ( lutrgba == 0 )
+         {
+            std::cout << "Lut RGBA (Palette Color) not built " << std::endl;
  
-        // Nothing is written yet to get LUT Data user friendly
-        // The following is to be moved into a PixelRedaConvert method
+           // Nothing is written yet to get LUT Data user friendly
+           // The following is to be moved into a PixelRedaConvert method
   
-         gdcm::SeqEntry *modLutSeq = e1->GetSeqEntry(0x0028,0x3000);
-         if ( modLutSeq !=0 )
-         {
-            gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
-            if ( !sqi )
+            gdcm::SeqEntry *modLutSeq = f->GetSeqEntry(0x0028,0x3000);
+            if ( modLutSeq !=0 )
             {
-               std::string lutDescriptor = sqi->GetEntryValue(0x0028,0x3002);
-               int length;   // LUT length in Bytes
-               int deb;      // Subscript of the first Lut Value
-               int nbits;    // Lut item size (in Bits)
-               int nbRead;    // nb of items in LUT descriptor (must be = 3)
-
-               nbRead = sscanf( lutDescriptor.c_str(),
-                                 "%d\\%d\\%d",
-                                  &length, &deb, &nbits );
-               if ( nbRead != 3 )
-               {
-                   //gdcmWarningMacro( "Wrong LUT descriptor" );
-                   std::cout << "Wrong LUT descriptor" << std::endl;
-               }                                                  
-               gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006);
-               if ( b != 0 )
+               gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
+               if ( !sqi )
                {
-                  if ( b->GetLength() != 0 )
+               std::string lutDescriptor = sqi->GetEntryValue(0x0028,0x3002);
+                  int length;   // LUT length in Bytes
+                  int deb;      // Subscript of the first Lut Value
+                  int nbits;    // Lut item size (in Bits)
+                  int nbRead;    // nb of items in LUT descriptor (must be = 3)
+
+                  nbRead = sscanf( lutDescriptor.c_str(),
+                                    "%d\\%d\\%d",
+                                     &length, &deb, &nbits );
+                  if ( nbRead != 3 )
                   {
-                     std::cout << "---------------------------------------"
+                      //gdcmWarningMacro( "Wrong LUT descriptor" );
+                      std::cout << "Wrong LUT descriptor" << std::endl;
+                  }                                                  
+                  gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006);
+                  if ( b != 0 )
+                  {
+                     if ( b->GetLength() != 0 )
+                     {
+                        std::cout << "---------------------------------------"
                                << " We should never reach this point      "
                                << std::endl;
-                     //LoadEntryBinArea(b);    //LUT Data (CTX dependent)
-                  }   
-              }
-           }      
-         }
+                        //LoadEntryBinArea(b);    //LUT Data (CTX dependent)
+                     }   
+                 }
+              }      
+            }
+            else
+               std::cout << "No LUT Data (0x0028,0x3000) found " << std::endl;
+        }
          else
-             std::cout << "No LUT Data (0x0028,0x3000) found " << std::endl;
-     }
-      else
-      {
-         if ( f1->GetLutItemSize() == 8 )
          {
-            for (int i=0;i<f1->GetLutItemNumber();i++)
-               std::cout << i << " : \t"
+            if ( fh->GetLutItemSize() == 8 )
+            {
+               for (int i=0;i<fh->GetLutItemNumber();i++)
+                  std::cout << i << " : \t"
                          << (int)(lutrgba[i*4])   << " "
                          << (int)(lutrgba[i*4+1]) << " "
                          << (int)(lutrgba[i*4+2]) << std::endl;
-         }
-         else // LutItemSize assumed to be = 16
-         {
-            uint16_t* lutrgba16 = (uint16_t*)lutrgba;
-            for (int i=0;i<f1->GetLutItemNumber();i++)
-               std::cout << i << " : \t"
+            }
+            else // LutItemSize assumed to be = 16
+            {
+               uint16_t* lutrgba16 = (uint16_t*)lutrgba;
+               for (int i=0;i<fh->GetLutItemNumber();i++)
+                  std::cout << i << " : \t"
                          << (int)(lutrgba16[i*4])   << " "
                          << (int)(lutrgba16[i*4+1]) << " "
                          << (int)(lutrgba16[i*4+2]) << std::endl;
+            }
          }
       }
+      else if (showlut)
+      {
+         std::cout << "Try LUT Data "<< std::endl;
+         ShowLutData(f);
+      }
+     
+      if (f->IsReadable())
+         std::cout <<std::endl<<fileName<<" is Readable"<<std::endl;
+      else
+         std::cout <<std::endl<<fileName<<" is NOT Readable"<<std::endl;
+      std::cout<<std::flush;
+      delete f;
+      delete fh;
+      return 0;
    }
-   else if (showlut)
+   else  // ====== Deal with a Directory ======
    {
-      std::cout << "Try LUT Data "<< std::endl;
-      ShowLutData(e1);
+      std::cout << "dirName [" << dirName << "]" << std::endl;
+      gdcm::DirList dirList(dirName,1); // gets recursively the file list
+      gdcm::DirListType fileList = dirList.GetFilenames();
+      gdcm::File *f;
+      bool res;
+      for( gdcm::DirListType::iterator it  = fileList.begin();
+                                 it != fileList.end();
+                                 ++it )
+      {
+         f = new gdcm::File();
+         f->SetLoadMode(loadMode);
+         f->SetFileName( it->c_str() );
+
+         for (int ri=0; ri<forceLoadNb; ri++)
+         {
+            printf("%04x,%04x\n",elemsToForceLoad[2*ri], elemsToForceLoad[2*ri+1]);
+            f->AddForceLoadElement((uint32_t)elemsToForceLoad[2*ri], 
+                                   (uint32_t)elemsToForceLoad[2*ri+1]); 
+         }
+         res = f->Load();
+
+         if ( !res )
+         {
+            std::cout << "Cannot process file [" << it->c_str() << "]" 
+                      << std::endl;
+            std::cout << "Either it doesn't exist, or it's read protected " 
+                      << std::endl;
+            std::cout << "or it's not a Dicom File, or its 'header' is bugged" 
+                      << std::endl;
+            std::cout << "use 'PrintFile filein=... debug' to try to guess the pb"
+                   << std::endl;
+            delete f;
+            continue;
+         }
+
+         gdcm::FileHelper *fh = new gdcm::FileHelper(f);
+         fh->SetPrintLevel( level );
+
+         fh->Print();
+
+         if (f->IsReadable())
+            std::cout <<std::endl<<it->c_str()<<" is Readable"<<std::endl;
+         else
+            std::cout <<std::endl<<it->c_str()<<" is NOT Readable"<<std::endl;
+         std::cout << "\n\n" << std::endl;
+         delete f;
+         delete fh;
+      }
+      std::cout<<std::flush;
    }
-     
-   if (e1->IsReadable())
-      std::cout <<std::endl<<fileName<<" is Readable"<<std::endl;
-   else
-      std::cout <<std::endl<<fileName<<" is NOT Readable"<<std::endl;
-   std::cout<<std::flush;
-   delete e1;
-   delete f1;
-   return 0;   
 }