X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FPrintFile.cxx;h=f1868e72b3c550ae0e176673b980a1abd263f5d9;hb=fa8fa610de8d935491343df2d8a543ff6fdb6e69;hp=b6fd6ce4fbaee82eb4db9f933cecce632bb8dbed;hpb=0a63df7d5417e936948123cb00a80d58addcae66;p=gdcm.git

diff --git a/Example/PrintFile.cxx b/Example/PrintFile.cxx
index b6fd6ce4..f1868e72 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/07/05 12:57:36 $
-  Version:   $Revision: 1.43 $
+  Date:      $Date: 2005/10/01 19:39:14 $
+  Version:   $Revision: 1.61 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -16,15 +16,18 @@
                                                                                 
 =========================================================================*/
 #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 "gdcmGlobal.h"
+#include "gdcmDictSet.h"
 #include "gdcmArgMgr.h"
-
+#include "gdcmOrientation.h"
 #include <iostream>
 
 
@@ -114,22 +117,27 @@ void ShowLutData(gdcm::File *f)
       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] [noshadowseq][noshadow][noseq][debug] ",
-   "        level = 0,1,2 : depending on the amount of details user wants to see",
-   "        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)",
+   " \n PrintFile : \n                                                        ",
+   " Display the header of a ACR-NEMA/PAPYRUS/DICOM File                      ",
+   " usage: PrintFile {filein=inputFileName|dirin=inputDirectoryName}[level=n]",
+   "                       [forceload=listOfElementsToForceLoad]              ",
+   "                       [dict= privateDirectory]                           ",
+   "                       [ { [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 ",
+   "      privateDirectory : source file full path name of Shadow Group elems ",
+   "      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   
@@ -142,29 +150,54 @@ int main(int argc, char *argv[])
       return 0;
    }
 
-   char *fileName = am->ArgMgrWantString("filein",usage);
+   char *fileName = am->ArgMgrGetString("filein",(char *)0);
+   char *dirName  = am->ArgMgrGetString("dirin",(char *)0);
 
-   int loadMode = 0x00000000;
+   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;
+ }
+
+   if (am->ArgMgrDefined("debug"))
+      gdcm::Debug::DebugOn();
+ 
+   int loadMode = gdcm::LD_ALL;
    if ( am->ArgMgrDefined("noshadowseq") )
-      loadMode |= NO_SHADOWSEQ;
+      loadMode |= gdcm::LD_NOSHADOWSEQ;
    else 
    {
    if ( am->ArgMgrDefined("noshadow") )
-         loadMode |= NO_SHADOW;
+         loadMode |= gdcm::LD_NOSHADOW;
       if ( am->ArgMgrDefined("noseq") )
-         loadMode |= NO_SEQ;
+         loadMode |= gdcm::LD_NOSEQ;
    }
 
-   int level = am->ArgMgrGetInt("level", 2);
+   int level = am->ArgMgrGetInt("level", 1);
 
-   if (am->ArgMgrDefined("debug"))
-      gdcm::Debug::DebugOn();
+   int forceLoadNb;
+   uint16_t *elemsToForceLoad 
+                           = am->ArgMgrGetXInt16Enum("forceload", &forceLoadNb);
 
    bool showlut = ( 0 != am->ArgMgrDefined("SHOWLUT") );
- 
+
+   bool ddict = am->ArgMgrDefined("dict") ? true : false;
+   char *dict = 0;
+
+   if (ddict)
+   {
+     dict = am->ArgMgrGetString("dict",(char *)0);
+   }
+
    /* if unused Param we give up */
    if ( am->ArgMgrPrintUnusedLabels() )
-   { 
+   {
       am->ArgMgrUsage(usage);
       delete am;
       return 0;
@@ -173,155 +206,316 @@ 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 *f = new gdcm::File();
-   f->SetLoadMode(loadMode);
 
-   bool res = f->Load( fileName );
-   if ( !res )
+ 
+   if (ddict)
    {
-      delete f;
-      return 0;
+      gdcm::Global::GetDicts()->GetDefaultPubDict()->AddDict(dict);   
    }
 
-   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= [" 
+   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)
+
+      gdcm::File *f = new gdcm::File();
+      f->SetLoadMode(loadMode);
+      f->SetFileName( fileName );
+
+      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;
+                << "] "<< std::endl;
 
-   int numberOfScalarComponents=f->GetNumberOfScalarComponents();
-   std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents 
-             <<std::endl
-             << " LUT = " << (f->HasLUT() ? "TRUE" : "FALSE")
-             << std::endl;
+      int numberOfScalarComponents=f->GetNumberOfScalarComponents();
+      std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents 
+                <<std::endl
+                << " LUT = " << (f->HasLUT() ? "TRUE" : "FALSE")
+                << std::endl;
 
-   if ( f->GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) 
-   {
-      std::cout << "Transfer Syntax not loaded. " << std::endl
-                << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE"
+      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;
-   }
+         return 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 )
-   {
-      uint8_t* lutrgba = fh->GetLutRGBA();
-      if ( lutrgba == 0 )
-      {
-         std::cout << "Lut RGBA (Palette Color) not built " << std::endl;
+      std::string transferSyntaxName = f->GetTransferSyntaxName();
+      std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" 
+                << std::endl;
+      std::cout << " SwapCode= " << f->GetSwapCode() << std::endl;
+      std::cout << " ------" << 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;
+
+      // Lets's get and print some usefull fields about 'Orientation'
+      // ------------------------------------------------------------
+
+      std::string strPatientPosition = 
+                                      f->GetEntryValue(0x0018,0x5100);
+      if ( strPatientPosition != gdcm::GDCM_UNFOUND 
+        && strPatientPosition != "" )  
+            std::cout << "PatientPosition (0x0010,0x5100)= [" 
+                      << strPatientPosition << "]" << std::endl;
+ 
+     std::string strPatientOrientation = 
+                                      f->GetEntryValue(0x0020,0x0020);
+      if ( strPatientOrientation != gdcm::GDCM_UNFOUND
+        && strPatientOrientation != "")  
+         std::cout << "PatientOrientation (0x0020,0x0020)= [" 
+                   << strPatientOrientation << "]" << std::endl;
+
+      std::string strImageOrientationPatient = 
+                                      f->GetEntryValue(0x0020,0x0037);  
+      if ( strImageOrientationPatient != gdcm::GDCM_UNFOUND
+        && strImageOrientationPatient != "" )  
+         std::cout << "ImageOrientationPatient (0x0020,0x0037)= [" 
+                   << strImageOrientationPatient << "]" << std::endl;
+
+      std::string strImageOrientationRET = 
+                                      f->GetEntryValue(0x0020,0x0035);
+      if ( strImageOrientationRET != gdcm::GDCM_UNFOUND
+        && strImageOrientationRET != "" )  
+         std::cout << "ImageOrientationRET (0x0020,0x0035)= [" 
+                   << strImageOrientationRET << "]" << std::endl;
+
+      // Let's compute 'user friendly' results about 'Orientation'
+      // ---------------------------------------------------------
  
-        // Nothing is written yet to get LUT Data user friendly
-        // The following is to be moved into a PixelRedaConvert method
+      gdcm::Orientation o;
+
+      if ( strImageOrientationPatient != gdcm::GDCM_UNFOUND ||
+           strImageOrientationRET     != gdcm::GDCM_UNFOUND )
+      {
   
-         gdcm::SeqEntry *modLutSeq = f->GetSeqEntry(0x0028,0x3000);
-         if ( modLutSeq !=0 )
+         gdcm::OrientationType orient = o.GetOrientationType( f );
+ 
+         std::cout << "TypeOrientation = " << orient << " (-> " 
+                   << o.GetOrientationTypeString(orient) << " )" << std::endl;
+      }
+
+      std::string ori = o.GetOrientation ( f );
+      if (ori != "\\" )
+         std::cout << "Orientation [" << ori << "]" << std::endl;
+
+
+      // Display the LUT as an int array (for debugging purpose)
+      if ( f->HasLUT() && showlut )
+      {
+         uint8_t* lutrgba = fh->GetLutRGBA();
+         if ( lutrgba == 0 )
          {
-            gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
-            if ( !sqi )
+            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
+  
+            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 )
+                  {
+                      //gdcmWarningMacro( "Wrong LUT descriptor" );
+                      std::cout << "Wrong LUT descriptor" << std::endl;
+                  }                                                  
+                  gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006);
+                  if ( b != 0 )
                   {
-                     std::cout << "---------------------------------------"
+                     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 ( fh->GetLutItemSize() == 8 )
          {
-            for (int i=0;i<fh->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<fh->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->gdcm::Document::IsReadable())
+      // Try downcast to please MSVC
+     if ( !((gdcm::Document *)f)->IsReadable() )
+     {
+         std::cout <<std::endl<<fileName<<" is NOT 'gdcm parsable'"<<std::endl;
+      }
+     
+      if (f->IsReadable())
+         std::cout <<std::endl<<fileName<<" is Readable"<<std::endl;
+      else if ( f->GetSeqEntry(0x0041,0x1010) )
+      {
+         std::cout <<std::endl<<fileName<<" looks like a 'PAPYRUS image' file"
+                   <<std::endl;
+      }
+      else if ( f->GetSeqEntry(0x0004,0x1220) )
+      {
+         std::cout <<std::endl<<fileName<<" looks like a 'DICOMDIR file'"
+                   <<std::endl;
+      }
+      else 
+      {
+         std::cout <<std::endl<<fileName<<" doesn't look like an image file "
+             <<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(f);
+      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 )
+      {
+         std::cout << std::endl<<" Start processing :[" << it->c_str() << "]"
+                   << std::endl;
+         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();
+         gdcm::Orientation o;
+         std::string ori = o.GetOrientation ( f );
+         if (ori != gdcm::GDCM_UNFOUND )
+            std::cout << "- Orientation [" << ori << "]" << std::endl;
+
+         gdcm::OrientationType d = o.GetOrientationType( f );
+ 
+         std::cout << "TypeOrientation = " << d << " (-> " 
+                   << o.GetOrientationTypeString(d) << std::endl;
+        
+         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 (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;   
 }