X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FPrintFile.cxx;h=4310831a6d90de76c4f8fa5a598fa17109f30771;hb=f3cc252fad08c3e8dbc395f2ca7cabf47bf95302;hp=36bb2ac352a5eff3bd5391f1d071b4c001253822;hpb=978f699b97794e209f9d20e2a33a72404bda54ef;p=gdcm.git

diff --git a/Example/PrintFile.cxx b/Example/PrintFile.cxx
index 36bb2ac3..4310831a 100644
--- a/Example/PrintFile.cxx
+++ b/Example/PrintFile.cxx
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: PrintFile.cxx,v $
   Language:  C++
-  Date:      $Date: 2004/11/16 04:26:17 $
-  Version:   $Revision: 1.10 $
+  Date:      $Date: 2005/08/30 15:13:05 $
+  Version:   $Revision: 1.50 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -15,91 +15,405 @@
      PURPOSE.  See the above copyright notices for more information.
                                                                                 
 =========================================================================*/
+#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>
-#include "gdcm.h"
 
-int main(int argc, char* argv[])
+
+void ShowLutData(gdcm::File *f);
+
+void ShowLutData(gdcm::File *f)
 {
- 
-   gdcm::Header *e1;
-   gdcm::File   *f1;
-   std::string fileName;   
-   if (argc != 2) {
-      std::cout << " usage : PrintDocument fileName" << std::endl;
-   }
+     // 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 = f->GetSeqEntry(0x0028,0x3000);
+      if ( modLutSeq !=0 )
+      {
+         gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
+         if ( sqi != 0 )
+         {
+            std::string lutDescriptor = sqi->GetEntryValue(0x0028,0x3002);
+           if (   /*lutDescriptor   == GDCM_UNFOUND*/ 0 )
+           {
+              //gdcmWarningMacro( "LUT Descriptor is missing" );
+              std::cout << "LUT Descriptor is missing" << std::endl;
+              return;
+            }
+            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)
 
-   if (argc > 1) {
-      fileName=argv[1];
-   } else {
-      fileName += GDCM_DATA_ROOT;
-      fileName += "/test.acr";
+            nbRead = sscanf( lutDescriptor.c_str(),
+                              "%d\\%d\\%d",
+                               &length, &deb, &nbits );
+           std::cout << "length " << length 
+                     << " deb " << deb 
+                     << " nbits " << nbits
+                     << std::endl;
+            if ( nbRead != 3 )
+            {
+                //gdcmWarningMacro( "Wrong LUT descriptor" );
+                std::cout << "Wrong LUT descriptor" << std::endl;
+            }
+            //LUT Data (CTX dependent)    
+            gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006); 
+            if ( b != 0 )
+            { 
+               int BitsAllocated = f->GetBitsAllocated();
+               if ( BitsAllocated <= 8 )
+               { 
+                  int mult;
+                  if ( ( nbits == 16 ) && ( BitsAllocated == 8 ) )
+                  {
+                  // when LUT item size is different than pixel size
+                     mult = 2; // high byte must be = low byte
+                  }
+                  else
+                  {
+                  // See PS 3.3-2003 C.11.1.1.2 p 619
+                     mult = 1;
+                  }
+                  uint8_t *lut = b->GetBinArea();
+                  for( int i=0; i < length; ++i )
+                  {
+                     std::cout << i+deb << " : \t"
+                               << (int) (lut[i*mult + 1]) << std::endl;
+                  }
+               }
+               else
+               {
+                  uint16_t *lut = (uint16_t *)(b->GetBinArea());  
+                  for( int i=0; i < length; ++i )
+                  {
+                     std::cout << i+deb << " : \t"
+                               << (int) (((uint16_t *)lut)[i])
+                               << std::endl;
+                  }             
+               }
+            }  
+            else
+               std::cout << "No LUT Data BinEntry (0x0028,0x3006) found?!? " 
+                         << std::endl;
+         }
+         else
+            std::cout << "No First SQ Item within (0x0028,0x3000) ?!? " 
+                      << std::endl;      
+      }
+      else
+         std::cout << "No LUT Data SeqEntry (0x0028,0x3000) found " 
+                   << std::endl;
    }
-   
-   e1= new gdcm::Header( fileName.c_str() );
-
-   f1 = new gdcm::File(e1);
-
-   e1->SetPrintLevel(2);
-   
-   e1->Print();
-      
-   std::cout << "\n\n" << std::endl; 
-
-   int dataSize = f1->GetImageDataSize();
-   std::cout <<std::endl <<" dataSize " << dataSize << 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 
-             << " SampleserPixel="      << sPP
-             << " PlanarConfiguration=" << planarConfig 
-             << std::endl 
-             << " PhotometricInterpretation=" 
-                                << e1->GetEntryByNumber(0x0028,0x0004)
-             << std::endl;
-
-   int numberOfScalarComponents=e1->GetNumberOfScalarComponents();
-   std::cout << " NumberOfScalarComponents " << numberOfScalarComponents <<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=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 ( e1->GetEntryByNumber(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) {
-     std::cout << "Transfert Syntax not loaded. " << std::endl
-               << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE"
-               << std::endl;
+   if (argc == 1 || am->ArgMgrDefined("usage") )
+   {
+      am->ArgMgrUsage(usage); // Display 'usage'
+      delete am;
+      return 0;
+   }
+
+   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;
+ }
+
+   if (am->ArgMgrDefined("debug"))
+      gdcm::Debug::DebugOn();
+ 
+   int loadMode = gdcm::LD_ALL;
+   if ( am->ArgMgrDefined("noshadowseq") )
+      loadMode |= gdcm::LD_NOSHADOWSEQ;
+   else 
+   {
+   if ( am->ArgMgrDefined("noshadow") )
+         loadMode |= gdcm::LD_NOSHADOW;
+      if ( am->ArgMgrDefined("noseq") )
+         loadMode |= gdcm::LD_NOSEQ;
+   }
+
+   int level = am->ArgMgrGetInt("level", 2);
+
+   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;
-  }
+   } 
+
+   delete am;  // we don't need Argument Manager any longer
+
+   // ----------- End Arguments Manager ---------
+
+   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;
+
+      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"
+                << std::endl;
+         return 0;
+      }
   
-   std::string transferSyntaxName = e1->GetTransfertSyntaxName();
-   std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" << std::endl;
-   
-   if (  transferSyntaxName != "Implicit VR - Little Endian"
-      && transferSyntaxName != "Explicit VR - Little Endian"     
-      && transferSyntaxName != "Deflated Explicit VR - Little Endian"      
-      && transferSyntaxName != "Explicit VR - Big Endian"
-      && transferSyntaxName != "Uncompressed ACR-NEMA"     )
-  {
-      std::cout << std::endl << "==========================================="
-                  << std::endl; 
-         f1->GetPixelConverter()->Print();
-      std::cout << std::endl << "==========================================="
-                  << std::endl; 
-   }      
-   
-   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;
-
-   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;
+ 
+           // 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 )
+            {
+               gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem();
+               if ( !sqi )
+               {
+               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 )
+                  {
+                     if ( b->GetLength() != 0 )
+                     {
+                        std::cout << "---------------------------------------"
+                               << " We should never reach this point      "
+                               << std::endl;
+                        //LoadEntryBinArea(b);    //LUT Data (CTX dependent)
+                     }   
+                 }
+              }      
+            }
+            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"
+                         << (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"
+                         << (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  // ====== Deal with a Directory ======
+   {
+      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;
+   }
 }