]> Creatis software - gdcm.git/blobdiff - src/gdcmDocument.cxx
* Test/PrintAllDocument.cxx: looping on files is now effective. It used to
[gdcm.git] / src / gdcmDocument.cxx
index fa1d4a7d8d494cd55aa2c6aaba2f0a204cd94227..e0d531a7e382d6cb1268fb6741a96759cf59d497 100644 (file)
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: gdcmDocument.cxx,v $
   Language:  C++
-  Date:      $Date: 2004/06/21 12:38:29 $
-  Version:   $Revision: 1.21 $
+  Date:      $Date: 2004/06/25 12:56:39 $
+  Version:   $Revision: 1.33 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -66,7 +66,7 @@
 #define UI1_2_840_10008_1_2_4_91 "1.2.840.10008.1.2.4.91"
 // RLE Lossless
 #define UI1_2_840_10008_1_2_5    "1.2.840.10008.1.2.5"
-
+// UI1_1_2_840_10008_1_2_5
 #define str2num(str, typeNum) *((typeNum *)(str))
 
 //-----------------------------------------------------------------------------
@@ -74,9 +74,9 @@
 const unsigned int gdcmDocument::HEADER_LENGTH_TO_READ = 256;
 
 // Refer to gdcmDocument::SetMaxSizeLoadEntry()
-const unsigned int gdcmDocument::MAX_SIZE_LOAD_ELEMENT_VALUE = 4096;
+const unsigned int gdcmDocument::MAX_SIZE_LOAD_ELEMENT_VALUE = 0x7fffffff;// 4096;
 
-const unsigned int gdcmDocument::MAX_SIZE_PRINT_ELEMENT_VALUE = 64;
+const unsigned int gdcmDocument::MAX_SIZE_PRINT_ELEMENT_VALUE = 0x7fffffff;//64;
 
 //-----------------------------------------------------------------------------
 // Constructor / Destructor
@@ -126,6 +126,25 @@ gdcmDocument::gdcmDocument(std::string const & inFilename,
    
    long l=ParseDES( this, beg, lgt, false); // le Load sera fait a la volee
    (void)l; //is l used anywhere ?
+
+   rewind(fp); 
+   
+   // Load 'non string' values
+      
+   std::string PhotometricInterpretation = GetEntryByNumber(0x0028,0x0004);   
+   if( PhotometricInterpretation == "PALETTE COLOR " ) {
+      LoadEntryVoidArea(0x0028,0x1200);  // gray LUT   
+      LoadEntryVoidArea(0x0028,0x1201);  // R    LUT
+      LoadEntryVoidArea(0x0028,0x1202);  // G    LUT
+      LoadEntryVoidArea(0x0028,0x1203);  // B    LUT
+      
+      LoadEntryVoidArea(0x0028,0x1221);  // Segmented Red   Palette Color LUT Data
+      LoadEntryVoidArea(0x0028,0x1222);  // Segmented Green Palette Color LUT Data
+      LoadEntryVoidArea(0x0028,0x1223);  // Segmented Blue  Palette Color LUT Data
+   } 
+   //FIXME later : how to use it?
+   LoadEntryVoidArea(0x0028,0x3006);  //LUT Data (CTX dependent) 
+
    CloseFile(); 
   
    // --------------------------------------------------------------
@@ -175,8 +194,8 @@ gdcmDocument::~gdcmDocument (void) {
 
    // Recursive clean up of sequences
    for (TagDocEntryHT::iterator it = tagHT.begin(); it != tagHT.end(); ++it )
-   {
-      delete (it->second);
+   { 
+         delete it->second;
    }
    tagHT.clear();
 }
@@ -206,14 +225,14 @@ void gdcmDocument::PrintShaDict(std::ostream & os) {
  * \brief   Get the public dictionary used
  */
 gdcmDict *gdcmDocument::GetPubDict(void) {
-   return(RefPubDict);
+   return RefPubDict;
 }
 
 /**
  * \brief   Get the shadow dictionary used
  */
 gdcmDict *gdcmDocument::GetShaDict(void) {
-   return(RefShaDict);
+   return RefShaDict;
 }
 
 /**
@@ -222,7 +241,7 @@ gdcmDict *gdcmDocument::GetShaDict(void) {
  */
 bool gdcmDocument::SetShaDict(gdcmDict *dict){
    RefShaDict=dict;
-   return(!RefShaDict);
+   return !RefShaDict;
 }
 
 /**
@@ -231,7 +250,7 @@ bool gdcmDocument::SetShaDict(gdcmDict *dict){
  */
 bool gdcmDocument::SetShaDict(DictKey dictName){
    RefShaDict=gdcmGlobal::GetDicts()->GetDict(dictName);
-   return(!RefShaDict);
+   return !RefShaDict;
 }
 
 /**
@@ -245,17 +264,18 @@ bool gdcmDocument::SetShaDict(DictKey dictName){
 bool gdcmDocument::IsReadable(void) { 
    if(Filetype==gdcmUnknown) {
       dbg.Verbose(0, "gdcmDocument::IsReadable: wrong filetype");
-      return(false);
+      return false;
    }
    if(!tagHT.empty()<=0) { 
       dbg.Verbose(0, "gdcmDocument::IsReadable: no tags in internal"
                      " hash table.");
-      return(false);
+      return false;
    }
 
    return(true);
 }
 
+
 /**
  * \brief   Internal function that checks whether the Transfer Syntax given
  *          as argument is the one present in the current document.
@@ -389,11 +409,12 @@ bool gdcmDocument::IsRLELossLessTransferSyntax(void)
  * @return  True when RLE Lossless found. False in all
  *          other cases.
  */
 bool gdcmDocument::IsJPEGLossless(void)
 {
    return (   IsGivenTransferSyntax(UI1_2_840_10008_1_2_4_55)
            || IsGivenTransferSyntax(UI1_2_840_10008_1_2_4_57)
-           || IsGivenTransferSyntax(UI1_2_840_10008_1_2_4_90) );
+           || IsGivenTransferSyntax(UI1_2_840_10008_1_2_4_70) );
 }
                                                                                 
 /**
@@ -490,34 +511,28 @@ bool gdcmDocument::CloseFile(void) {
 
 /**
  * \brief Writes in a file all the Header Entries (Dicom Elements) 
- *        of the Chained List
  * @param fp file pointer on an already open file
- * @param type Type of the File to be written 
+ * @param filetype Type of the File to be written 
  *          (ACR-NEMA, ExplicitVR, ImplicitVR)
  * \return Always true.
  */
-bool gdcmDocument::Write(FILE *fp, FileType type) {
+bool gdcmDocument::WriteF(FileType filetype) {
 /// \todo
-///  ==============
-///      The stuff will have to be rewritten using the SeQuence based 
-///       tree-like stucture instead  of the chained list .
-///      (so we shall remove the Group*HT from the gdcmDocument)
+/// ==============
+///      The stuff is rewritten using the SeQuence based 
+///       tree-like stucture (cf : Print )
 ///      To be checked
 /// =============
 
    /// \todo move the following lines (and a lot of others, to be written)
    /// to a future function CheckAndCorrectHeader
    
-   /// \todo
-   /// Question :
-   /// Comment pourrait-on savoir si le DcmHeader vient d'un fichier
-   /// DicomV3 ou non (FileType est un champ de gdcmDocument ...)
    /// WARNING : Si on veut ecrire du DICOM V3 a partir d'un DcmHeader ACR-NEMA
-   /// no way 
+   /// no way (check : FileType est un champ de gdcmDocument ...)
    /// a moins de se livrer a un tres complique ajout des champs manquants.
    /// faire un CheckAndCorrectHeader (?)  
 
-   if (type == gdcmImplicitVR) 
+   if (filetype == gdcmImplicitVR) 
    {
       std::string implicitVRTransfertSyntax = UI1_2_840_10008_1_2;
       ReplaceOrCreateByNumber(implicitVRTransfertSyntax,0x0002, 0x0010);
@@ -530,7 +545,7 @@ bool gdcmDocument::Write(FILE *fp, FileType type) {
       SetEntryLengthByNumber(18, 0x0002, 0x0010);
    } 
 
-   if (type == gdcmExplicitVR) 
+   if (filetype == gdcmExplicitVR)
    {
       std::string explicitVRTransfertSyntax = UI1_2_840_10008_1_2_1;
       ReplaceOrCreateByNumber(explicitVRTransfertSyntax,0x0002, 0x0010);
@@ -545,18 +560,20 @@ bool gdcmDocument::Write(FILE *fp, FileType type) {
 
 /**
  * \todo rewrite later, if really usefull
+ *               ('Group Length' element is optional in DICOM)
  *
- *       --> Warning : un-updated odd groups lengths can causes pb 
+ *       --> Warning : un-updated odd groups lengthes can causes pb
  *       -->           (xmedcon breaks)
  *       --> to be re- written with future org.
  *
- * if ( (type == ImplicitVR) || (type == ExplicitVR) )
- *    UpdateGroupLength(false,type);
- * if ( type == ACR)
+ * if ( (filetype == ImplicitVR) || (filetype == ExplicitVR) )
+ *    UpdateGroupLength(false,filetype);
+ * if ( filetype == ACR)
  *    UpdateGroupLength(true,ACR);
  */
+   gdcmElementSet::Write(fp,filetype);
 
-   WriteEntries(fp,type);
    return true;
 }
 
@@ -639,14 +656,14 @@ gdcmBinEntry * gdcmDocument::ReplaceOrCreateByNumber(
                                          guint16 Elem)
 {
    gdcmDocEntry* a;
-   gdcmBinEntry* b;    
+   gdcmBinEntry* b;
    a = GetDocEntryByNumber( Group, Elem);
    if (a == NULL) {
       a =NewBinEntryByNumber(Group, Elem);
       if (a == NULL) 
          return NULL;
 
-               b = new gdcmBinEntry(a);                        
+      b = new gdcmBinEntry(a);
       AddEntry(b);
    }   
    SetEntryByNumber(voidArea, lgth, Group, Elem);
@@ -716,7 +733,7 @@ std::string gdcmDocument::GetEntryByName(std::string tagName) {
  * @return  Corresponding element value representation when it exists,
  *          and the string GDCM_UNFOUND ("gdcm::Unfound") otherwise.
  */
-std::string gdcmDocument::GetEntryVRByName(std::string tagName) {
+std::string gdcmDocument::GetEntryVRByName(TagName tagName) {
    gdcmDictEntry *dictEntry = RefPubDict->GetDictEntryByName(tagName); 
    if( dictEntry == NULL)
       return GDCM_UNFOUND;
@@ -848,7 +865,7 @@ bool gdcmDocument::SetEntryByNumber(void *content,
    TagKey key = gdcmDictEntry::TranslateToKey(group, element);
    if ( ! tagHT.count(key))
       return false;
-               
+
 /* Hope Binray field length is never wrong    
    if(lgth%2) // Non even length are padded with a space (020H).
    {  
@@ -950,13 +967,37 @@ void *gdcmDocument::LoadEntryVoidArea(guint16 Group, guint16 Elem)
       delete[] a;
       return NULL;
    }
-
+   return a;  
+}
+/**
+ * \brief         Loads (from disk) the element content 
+ *                when a string is not suitable
+ * @param Element  Entry whose voidArea is going to be loaded
+ */
+void *gdcmDocument::LoadEntryVoidArea(gdcmBinEntry *Element) 
+{
+   size_t o =(size_t)Element->GetOffset();
+   fseek(fp, o, SEEK_SET);
+   size_t l = Element->GetLength();
+   char* a = new char[l];
+   if(!a) {
+      dbg.Verbose(0, "gdcmDocument::LoadEntryVoidArea cannot allocate a");
+      return NULL;
+   }
+   Element->SetVoidArea((void *)a);
+   /// \todo check the result 
+   size_t l2 = fread(a, 1, l ,fp);
+   if(l != l2) 
+   {
+      delete[] a;
+      return NULL;
+   }
    return a;  
 }
 
 /**
  * \brief   Sets a 'non string' value to a given Dicom Element
- * @param   area
+ * @param   area area containing the 'non string' value
  * @param   group     Group number of the searched Dicom Element 
  * @param   element Element number of the searched Dicom Element 
  * @return  
@@ -1048,7 +1089,7 @@ void gdcmDocument::UpdateShaEntries(void) {
  */
 gdcmDocEntry* gdcmDocument::GetDocEntryByNumber(guint16 group, guint16 element) 
 {
-   TagKey key = gdcmDictEntry::TranslateToKey(group, element);   
+   TagKey key = gdcmDictEntry::TranslateToKey(group, element);
    if ( ! tagHT.count(key))
       return NULL;
    return tagHT.find(key)->second;
@@ -1086,233 +1127,6 @@ void gdcmDocument::LoadDocEntrySafe(gdcmDocEntry * entry) {
    fseek(fp, PositionOnEntry, SEEK_SET);
 }
 
-
-/**
- * \brief Writes in a file (according to the requested format)
- *        the group, the element, the value representation and the length
- *        of a single gdcmDocEntry passed as argument.
- * @param tag  pointer on the gdcmDocEntry to be written
- * @param _fp  already open file pointer
- * @param type type of the File to be written 
- */
-void gdcmDocument::WriteEntryTagVRLength(gdcmDocEntry *tag,
-                                       FILE *_fp,
-                                       FileType type)
-{
-   guint16 group  = tag->GetGroup();
-   VRKey   vr     = tag->GetVR();
-   guint16 el     = tag->GetElement();
-   guint32 lgr    = tag->GetReadLength();
-
-   if ( (group == 0xfffe) && (el == 0x0000) ) 
-     // Fix in order to make some MR PHILIPS images e-film readable
-     // see gdcmData/gdcm-MR-PHILIPS-16-Multi-Seq.dcm:
-     // we just *always* ignore spurious fffe|0000 tag !   
-      return; 
-
-   fwrite ( &group,(size_t)2 ,(size_t)1 ,_fp);  //group
-   fwrite ( &el,(size_t)2 ,(size_t)1 ,_fp);     //element
-      
-   if ( type == gdcmExplicitVR ) {
-
-      // Special case of delimiters:
-      if (group == 0xfffe) {
-         // Delimiters have NO Value Representation and have NO length.
-         // Hence we skip writing the VR and length and we pad by writing
-         // 0xffffffff
-
-         int ff=0xffffffff;
-         fwrite (&ff,(size_t)4 ,(size_t)1 ,_fp);
-         return;
-      }
-
-      guint16 z=0;
-      guint16 shortLgr = lgr;
-      if (vr == "unkn") {     // Unknown was 'written'
-         // deal with Little Endian            
-         fwrite ( &shortLgr,(size_t)2 ,(size_t)1 ,_fp);
-         fwrite ( &z,  (size_t)2 ,(size_t)1 ,_fp);
-      } else {
-         fwrite (vr.c_str(),(size_t)2 ,(size_t)1 ,_fp);                     
-         if ( (vr == "OB") || (vr == "OW") || (vr == "SQ") )
-         {
-            fwrite ( &z,  (size_t)2 ,(size_t)1 ,_fp);
-            fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp);
-         } else {
-            fwrite ( &shortLgr,(size_t)2 ,(size_t)1 ,_fp);
-         }
-      }
-   } 
-   else // IMPLICIT VR 
-   { 
-      fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp);
-   }
-}
-      
-/**
- * \brief Writes in a file (according to the requested format)
- *        the value of a single gdcmDocEntry passed as argument.
- * @param tag  Pointer on the gdcmDocEntry to be written
- * @param _fp  Already open file pointer
- * @param type type of the File to be written
- */
-// \todo TODO : to be re -written recursively !
-void gdcmDocument::WriteEntryValue(gdcmDocEntry *Entry, FILE *_fp,FileType type)
-{
-   (void)type;
-   guint16 group  = Entry->GetGroup();
-   VRKey   vr     = Entry->GetVR();
-   guint32 lgr    = Entry->GetReadLength();
-
-   if (vr == "SQ")
-      // SeQuences have no value:
-      return;
-   if (group == 0xfffe)
-      // Delimiters have no associated value:
-      return;
-               
-// if (gdcmBinEntry* BinEntry = dynamic_cast< gdcmBinEntry* >(Entry) ) {
-      void *voidArea;
-      gdcmBinEntry *BinEntry= (gdcmBinEntry *)Entry;;
-      voidArea = BinEntry->GetVoidArea();
-      if (voidArea != NULL) 
-      { // there is a 'non string' LUT, overlay, etc
-         fwrite ( voidArea,(size_t)lgr ,(size_t)1 ,_fp); // Elem value
-         return;            
-      }
-// } 
-
-   if (vr == "US" || vr == "SS") 
-   {
-      // some 'Short integer' fields may be mulivaluated
-      // each single value is separated from the next one by '\'
-      // we split the string and write each value as a short int
-      std::vector<std::string> tokens;
-      tokens.erase(tokens.begin(),tokens.end()); // clean any previous value
-      Tokenize (((gdcmValEntry *)Entry)->GetValue(), tokens, "\\");
-      for (unsigned int i=0; i<tokens.size();i++) 
-      {
-         guint16 val_uint16 = atoi(tokens[i].c_str());
-         void *ptr = &val_uint16;
-         fwrite ( ptr,(size_t)2 ,(size_t)1 ,_fp);
-      }
-      tokens.clear();
-      return;
-   }
-   if (vr == "UL" || vr == "SL") 
-   {
-      // Some 'Integer' fields may be multivaluated (multiple instances 
-      // of integers). But each single integer value is separated from the
-      // next one by '\' (backslash character). Hence we split the string
-      // along the '\' and write each value as an int:
-      std::vector<std::string> tokens;
-      tokens.erase(tokens.begin(),tokens.end()); // clean any previous value
-      Tokenize (((gdcmValEntry *)Entry)->GetValue(), tokens, "\\");
-      for (unsigned int i=0; i<tokens.size();i++) 
-      {
-         guint32 val_uint32 = atoi(tokens[i].c_str());
-         void *ptr = &val_uint32;
-         fwrite ( ptr,(size_t)4 ,(size_t)1 ,_fp);
-      }
-      tokens.clear();
-      return;
-   }           
-   fwrite (((gdcmValEntry *)Entry)->GetValue().c_str(),
-           (size_t)lgr ,(size_t)1, _fp); // Elem value
-}
-
-/**
- * \brief Writes in a file (according to the requested format)
- *        a single gdcmDocEntry passed as argument.
- * \sa    WriteEntryValue, WriteEntryTagVRLength.
- * @param tag  Pointer on the gdcmDocEntry to be written
- * @param _fp  Already open file pointer
- * @param type type of the File to be written
- */
-
-bool gdcmDocument::WriteEntry(gdcmDocEntry *tag, FILE *_fp,FileType type)
-{
-   guint32 length = tag->GetLength();
-
-   if (gdcmValEntry* ValEntry = dynamic_cast< gdcmValEntry* >(tag) )
-   {
-      // The value of a tag MUST (see the DICOM norm) be an odd number of
-      // bytes. When this is not the case, pad with an additional byte:
-      if(length%2==1) {
-         ValEntry->SetValue(ValEntry->GetValue()+"\0");
-         ValEntry->SetLength(ValEntry->GetReadLength()+1);
-      }
-      WriteEntryTagVRLength(ValEntry, _fp, type);
-      WriteEntryValue(ValEntry, _fp, type);
-      return true;
-   }
-
-   if (gdcmBinEntry* BinEntry = dynamic_cast< gdcmBinEntry* >(tag) )
-   {
-    (void)BinEntry; //not used
-      /// \todo FIXME : when voidArea belong to gdcmBinEntry only, fix
-      /// voidArea length
-      //
-      // The value of a tag MUST (see the DICOM norm) be an odd number of
-      // bytes. When this is not the case, pad with an additional byte:
-/*
-      if(length%2==1) { 
-         tag->SetValue(tag->GetValue()+"\0");
-         tag->SetLength(tag->GetReadLength()+1);
-      }
-*/
-      WriteEntryTagVRLength(tag, _fp, type);
-      WriteEntryValue(tag, _fp, type);
-      return true;
-   }
-   
-  return false; //default behavior ?
-}
-
-/**
- * \brief   writes on disc according to the requested format
- *          (ACR-NEMA, ExplicitVR, ImplicitVR) the image
- *          using the Chained List
- * \warning does NOT add the missing elements in the header :
- *           it's up to the user doing it !
- *           (function CheckHeaderCoherence to be written)
- * \warning DON'T try, right now, to write a DICOM image
- *           from an ACR Header (meta elements will be missing!)
- * @param   type type of the File to be written 
- *          (ACR-NEMA, ExplicitVR, ImplicitVR)
- * @param   _fp already open file pointer
- */
-
-bool gdcmDocument::WriteEntries(FILE *_fp,FileType type)
-{ 
-   /// \todo FIXME : explore recursively the whole structure...
-   /// \todo (?) check write failures (after *each* fwrite)
-   dbg.Verbose(0, "gdcmDocument::WriteEntries: entering.");
-   for (TagDocEntryHT::iterator it = tagHT.begin(); it != tagHT.end(); ++it )
-   {
-      gdcmDocEntry * entry = it->second;
-
-      if ( type == gdcmACR ){ 
-         if (entry->GetGroup() < 0x0008)
-            // Ignore pure DICOM V3 groups
-            continue;
-         if (entry->GetGroup() %2)
-            // Ignore the "shadow" groups
-            continue;
-         if (entry->GetVR() == "SQ" ) // ignore Sequences
-            continue;    
-      } 
-      if (! WriteEntry(entry, _fp, type) ) {
-         dbg.Verbose(0, "gdcmDocument::WriteEntries: write failure.");
-         return false;
-      }
-   }
-   return true;
-}   
-
-
 /**
  * \brief   Swaps back the bytes of 4-byte long integer accordingly to
  *          processor order.
@@ -1339,8 +1153,8 @@ guint32 gdcmDocument::SwapLong(guint32 a) {
          dbg.Error(" gdcmDocument::SwapLong : unset swap code");
          a=0;
    }
-   return(a);
-}
+   return a;
+} 
 
 /**
  * \brief   Unswaps back the bytes of 4-byte long integer accordingly to
@@ -1348,7 +1162,7 @@ guint32 gdcmDocument::SwapLong(guint32 a) {
  * @return  The properly unswaped 32 bits integer.
  */
 guint32 gdcmDocument::UnswapLong(guint32 a) {
-   return (SwapLong(a));
+   return SwapLong(a);
 }
 
 /**
@@ -1358,7 +1172,7 @@ guint32 gdcmDocument::UnswapLong(guint32 a) {
 guint16 gdcmDocument::SwapShort(guint16 a) {
    if ( (sw==4321)  || (sw==2143) )
       a =(((a<<8) & 0x0ff00) | ((a>>8)&0x00ff));
-   return (a);
+   return a;
 }
 
 /**
@@ -1366,7 +1180,7 @@ guint16 gdcmDocument::SwapShort(guint16 a) {
  * @return  The properly unswaped 16 bits integer.
  */
 guint16 gdcmDocument::UnswapShort(guint16 a) {
-   return (SwapShort(a));
+   return SwapShort(a);
 }
 
 //-----------------------------------------------------------------------------
@@ -1374,7 +1188,7 @@ guint16 gdcmDocument::UnswapShort(guint16 a) {
 
 /**
  * \brief   Parses a DocEntrySet (Zero-level DocEntries or SQ Item DocEntries)
- * @return  false if file is not ACR-NEMA / PAPYRUS / DICOM 
+ * @return  length of the parsed set. 
  */ 
 
 long gdcmDocument::ParseDES(gdcmDocEntrySet *set, long offset, long l_max, bool delim_mode) {
@@ -1409,7 +1223,7 @@ long gdcmDocument::ParseDES(gdcmDocEntrySet *set, long offset, long l_max, bool
             NewValEntry->SetDepthLevel(depth);
             set->AddEntry(NewValEntry);
             LoadDocEntry(NewValEntry);
-            if (/*!delim_mode && */NewValEntry->isItemDelimitor())
+            if (NewValEntry->isItemDelimitor())
                break;
             if ( !delim_mode && ftell(fp)-offset >= l_max)
             {
@@ -1435,10 +1249,17 @@ long gdcmDocument::ParseDES(gdcmDocEntrySet *set, long offset, long l_max, bool
          if (NewDocEntry->GetGroup()   == 0x7fe0 && 
              NewDocEntry->GetElement() == 0x0010 )
          {
-             if (NewDocEntry->GetLength()==0xffffffff)
+             if (NewDocEntry->GetReadLength()==0xffffffff)
+             {
                 // Broken US.3405.1.dcm
                 Parse7FE0(); // to skip the pixels 
                              // (multipart JPEG/RLE are trouble makers)
+             }
+             else
+             {
+                SkipToNextDocEntry(NewDocEntry);
+                l = NewDocEntry->GetFullLength(); 
+             }
          }
          else
          {
@@ -1470,7 +1291,6 @@ long gdcmDocument::ParseDES(gdcmDocEntrySet *set, long offset, long l_max, bool
                                 l, delim_mode);
             (void)lgt;  //not used...
          }
-         // FIXME : on en fait quoi, de lgt ?
          set->AddEntry(sq);
          if ( !delim_mode && ftell(fp)-offset >= l_max)
          {
@@ -1479,7 +1299,7 @@ long gdcmDocument::ParseDES(gdcmDocEntrySet *set, long offset, long l_max, bool
       }
       delete NewDocEntry;
    }
-   return l; // ?? 
+   return l; // Probably useless 
 }
 
 /**
@@ -1490,7 +1310,7 @@ long gdcmDocument::ParseSQ(gdcmSeqEntry *set,
                            long offset, long l_max, bool delim_mode)
 {
    int SQItemNumber = 0;
-               
+
    gdcmDocEntry *NewDocEntry = (gdcmDocEntry *)0;
    gdcmSQItem *itemSQ;
    bool dlm_mod;
@@ -1500,11 +1320,12 @@ long gdcmDocument::ParseSQ(gdcmSeqEntry *set,
    (void)depth; //not used
 
    while (true) {
+
       NewDocEntry = ReadNextDocEntry();   
+      if (!NewDocEntry)
+         break;
       if(delim_mode) {   
          if (NewDocEntry->isSequenceDelimitor()) {
-            /// \todo add the Sequence Delimitor
-            /// \todo find the trick to put it properly !
             set->SetSequenceDelimitationItem(NewDocEntry);
             break;
           }
@@ -1515,7 +1336,6 @@ long gdcmDocument::ParseSQ(gdcmSeqEntry *set,
 
       itemSQ = new gdcmSQItem(set->GetDepthLevel());
       itemSQ->AddEntry(NewDocEntry);
-      /// \todo  no value, no voidArea. Think of it while printing !
       l= NewDocEntry->GetReadLength();
       
       if (l == 0xffffffff)
@@ -1525,14 +1345,14 @@ long gdcmDocument::ParseSQ(gdcmSeqEntry *set,
    
       lgr=ParseDES(itemSQ, NewDocEntry->GetOffset(), l, dlm_mod);
       
-      set->AddEntry(itemSQ,SQItemNumber);       
+      set->AddEntry(itemSQ,SQItemNumber); 
       SQItemNumber ++;
       if (!delim_mode && (ftell(fp)-offset) >= l_max) {
          break;
       }
    }
    lgth = ftell(fp) - offset;
-   return(lgth);
+   return lgth;
 }
 
 /**
@@ -1555,9 +1375,6 @@ void gdcmDocument::LoadDocEntry(gdcmDocEntry *Entry)
    //          (fffe e0dd) tells us the current SeQuence just ended
    if( group == 0xfffe ) {
       // NO more value field for SQ !
-      //Entry->SetValue("gdcm::Skipped");
-      // appel recursif de Load Value
-      // (meme pb que pour le parsing)
       return;
    }
 
@@ -1570,21 +1387,43 @@ void gdcmDocument::LoadDocEntry(gdcmDocEntry *Entry)
    // The elements whose length is bigger than the specified upper bound
    // are not loaded. Instead we leave a short notice of the offset of
    // the element content and it's length.
+
+   std::ostringstream s;
    if (length > MaxSizeLoadEntry) {
+      if (gdcmBinEntry* BinEntryPtr = dynamic_cast< gdcmBinEntry* >(Entry) )
+      {         
+         s << "gdcm::NotLoaded (BinEntry)";
+         s << " Address:" << (long)Entry->GetOffset();
+         s << " Length:"  << Entry->GetLength();
+         s << " x(" << std::hex << Entry->GetLength() << ")";
+         BinEntryPtr->SetValue(s.str());
+      }
+      // to be sure we are at the end of the value ...
+      fseek(fp,(long)Entry->GetOffset()+(long)Entry->GetLength(),SEEK_SET);      
+      return;
+       // Be carefull : a BinEntry IS_A ValEntry ... 
       if (gdcmValEntry* ValEntryPtr = dynamic_cast< gdcmValEntry* >(Entry) )
       {
-         std::ostringstream s;
-         s << "gdcm::NotLoaded.";
+         s << "gdcm::NotLoaded. (ValEntry)";
          s << " Address:" << (long)Entry->GetOffset();
          s << " Length:"  << Entry->GetLength();
          s << " x(" << std::hex << Entry->GetLength() << ")";
          ValEntryPtr->SetValue(s.str());
       }
       // to be sure we are at the end of the value ...
-      fseek(fp,(long)Entry->GetOffset()+(long)Entry->GetLength(),SEEK_SET);
-      
+      fseek(fp,(long)Entry->GetOffset()+(long)Entry->GetLength(),SEEK_SET);      
       return;
    }
+
+   // When we find a BinEntry not very much can be done :
+   if (gdcmBinEntry* BinEntryPtr = dynamic_cast< gdcmBinEntry* >(Entry) ) {
+
+      LoadEntryVoidArea(BinEntryPtr);
+      s << "gdcm::Loaded (BinEntry)";
+      BinEntryPtr->SetValue(s.str());
+      return;
+   }
     
    // Any compacter code suggested (?)
    if ( IsDocEntryAnInteger(Entry) ) {   
@@ -1634,16 +1473,22 @@ void gdcmDocument::LoadDocEntry(gdcmDocEntry *Entry)
    // We need an additional byte for storing \0 that is not on disk
    std::string NewValue(length,0);
    item_read = fread(&(NewValue[0]), (size_t)length, (size_t)1, fp);
-   if ( item_read != 1 ) {
-      dbg.Verbose(1, "gdcmDocument::LoadElementValue","unread element value");
-      ((gdcmValEntry *)Entry)->SetValue("gdcm::UnRead");
-      return;
+   if (gdcmValEntry* ValEntry = dynamic_cast< gdcmValEntry* >(Entry) ) {  
+      if ( item_read != 1 ) {
+         dbg.Verbose(1, "gdcmDocument::LoadElementValue","unread element value");
+         ValEntry->SetValue("gdcm::UnRead");
+         return;
+      }
+
+      if( (vr == "UI") ) // Because of correspondance with the VR dic
+         ValEntry->SetValue(NewValue.c_str());
+      else
+         ValEntry->SetValue(NewValue);
+   } else {
+   // fusible
+      std::cout << "Should have a ValEntry, here !" << std::endl;
    }
 
-   if( (vr == "UI") ) // Because of correspondance with the VR dic
-      ((gdcmValEntry *)Entry)->SetValue(NewValue.c_str());
-   else
-      ((gdcmValEntry *)Entry)->SetValue(NewValue);
 }
 
 
@@ -1733,14 +1578,12 @@ void gdcmDocument::LoadDocEntry(gdcmDocEntry *Entry)
          Entry->SetDictEntry(NewTag);
       }
        
-      // Heuristic: well some files are really ill-formed.
+      // Heuristic: well, some files are really ill-formed.
       if ( length16 == 0xffff) 
       {
          length16 = 0;
-         //dbg.Verbose(0, "gdcmDocument::FindLength",
-         //            "Erroneous element length fixed.");
-         // Actually, length= 0xffff means that we deal with
-         // Unknown Sequence Length 
+         // Length16= 0xffff means that we deal with
+         // 'Unknown Length' Sequence  
       }
       FixDocEntryFoundLength(Entry, (guint32)length16);
       return;
@@ -1839,7 +1682,7 @@ bool gdcmDocument::CheckDocEntryVR(gdcmDocEntry *Entry, VRKey vr)
                                    "UL","FIXME","Group Length");
          Entry->SetDictEntry(NewEntry);     
       }
-      return(false);
+      return false;
    }
 
    if ( Entry->IsVRUnknown() ) 
@@ -1863,7 +1706,7 @@ bool gdcmDocument::CheckDocEntryVR(gdcmDocEntry *Entry, VRKey vr)
                                  vr,"FIXME",Entry->GetName());
       Entry->SetDictEntry(NewEntry);
    }
-   return(true)
+   return true
 }
 
 /**
@@ -1929,10 +1772,10 @@ std::string gdcmDocument::GetDocEntryValue(gdcmDocEntry *Entry)
 #ifdef GDCM_NO_ANSI_STRING_STREAM
       s << std::ends; // to avoid oddities on Solaris
 #endif //GDCM_NO_ANSI_STRING_STREAM
-      return(s.str());
+      return s.str();
    }
 
-   return(((gdcmValEntry *)Entry)->GetValue());
+   return ((gdcmValEntry *)Entry)->GetValue();
 }
 
 /**
@@ -1982,10 +1825,10 @@ std::string gdcmDocument::GetDocEntryUnvalue(gdcmDocEntry *Entry)
 #ifdef GDCM_NO_ANSI_STRING_STREAM
       s << std::ends; // to avoid oddities on Solaris
 #endif //GDCM_NO_ANSI_STRING_STREAM
-      return(s.str());
+      return s.str();
    }
 
-   return(((gdcmValEntry *)Entry)->GetValue());
+   return ((gdcmValEntry *)Entry)->GetValue();
 }
 
 /**
@@ -2009,16 +1852,6 @@ void gdcmDocument::SkipToNextDocEntry(gdcmDocEntry *entry)
    (void)fseek(fp, (long)(entry->GetReadLength()), SEEK_CUR);
 }
 
-/**
- * \brief   Loads the value for a a given VLEntry 
- * \warning NOT end user intended method !
- * @param   entry 
- */
-void gdcmDocument::LoadVLEntry(gdcmDocEntry *entry) 
-{
-    //SkipBytes(entry->GetLength());
-    LoadDocEntry(entry);
-}
 /**
  * \brief   When the length of an element value is obviously wrong (because
  *          the parser went Jabberwocky) one can hope improving things by
@@ -2529,117 +2362,6 @@ gdcmDocEntry *gdcmDocument::ReadNextDocEntry(void) {
    return NewEntry;
 }
 
-/**
- * \brief   Build a new Element Value from all the low level arguments. 
- *          Check for existence of dictionary entry, and build
- *          a default one when absent.
- * @param   Name    Name of the underlying DictEntry
- */
-gdcmDocEntry *gdcmDocument::NewDocEntryByName(std::string Name) 
-{
-   gdcmDictEntry *NewTag = GetDictEntryByName(Name);
-   if (!NewTag)
-      NewTag = NewVirtualDictEntry(0xffff, 0xffff, "LO", "unkn", Name);
-
-   gdcmDocEntry* NewEntry = new gdcmDocEntry(NewTag);
-   if (!NewEntry) 
-   {
-      dbg.Verbose(1, "gdcmDocument::ObtainDocEntryByName",
-                  "failed to allocate gdcmDocEntry");
-      return (gdcmDocEntry *)0;
-   }
-   return NewEntry;
-}  
-
-/**
- * \brief   Request a new virtual dict entry to the dict set
- * @param   group     group  number of the underlying DictEntry
- * @param   element  element number of the underlying DictEntry
- * @param   vr     VR of the underlying DictEntry
- * @param   fourth owner group
- * @param   name   english name
- */
-gdcmDictEntry *gdcmDocument::NewVirtualDictEntry(guint16 group, guint16 element,
-                                               std::string vr,
-                                               std::string fourth,
-                                               std::string name)
-{
-   return gdcmGlobal::GetDicts()->NewVirtualDictEntry(group,element,vr,fourth,name);
-}
-
-/**
- * \brief   Build a new Element Value from all the low level arguments. 
- *          Check for existence of dictionary entry, and build
- *          a default one when absent.
- * @param   Group group   number of the underlying DictEntry
- * @param   Elem  element number of the underlying DictEntry
- */
-gdcmDocEntry* gdcmDocument::NewDocEntryByNumber(guint16 Group, guint16 Elem) 
-{
-   // Find out if the tag we encountered is in the dictionaries:
-   gdcmDictEntry *DictEntry = GetDictEntryByNumber(Group, Elem);
-   if (!DictEntry)
-      DictEntry = NewVirtualDictEntry(Group, Elem);
-
-   gdcmDocEntry *NewEntry = new gdcmDocEntry(DictEntry);
-   if (!NewEntry) 
-   {
-      dbg.Verbose(1, "gdcmDocument::NewDocEntryByNumber",
-                  "failed to allocate gdcmDocEntry");
-      return (gdcmDocEntry*)0;
-   }
-   return NewEntry;
-}
-
-
-/**
- * \brief   Build a new Element Value from all the low level arguments. 
- *          Check for existence of dictionary entry, and build
- *          a default one when absent.
- * @param   Group group   number of the underlying DictEntry
- * @param   Elem  element number of the underlying DictEntry
- */
-gdcmValEntry *gdcmDocument::NewValEntryByNumber(guint16 Group, guint16 Elem) 
-{
-   // Find out if the tag we encountered is in the dictionaries:
-   gdcmDictEntry *DictEntry = GetDictEntryByNumber(Group, Elem);
-   if (!DictEntry)
-      DictEntry = NewVirtualDictEntry(Group, Elem);
-
-   gdcmValEntry *NewEntry = new gdcmValEntry(DictEntry);
-   if (!NewEntry) 
-   {
-      dbg.Verbose(1, "gdcmDocument::NewValEntryByNumber",
-                  "failed to allocate gdcmValEntry");
-      return NULL;
-   }
-   return NewEntry;
-}
-
-
-/**
- * \brief   Build a new Element Value from all the low level arguments. 
- *          Check for existence of dictionary entry, and build
- *          a default one when absent.
- * @param   Group group   number of the underlying DictEntry
- * @param   Elem  element number of the underlying DictEntry
- */
-gdcmBinEntry *gdcmDocument::NewBinEntryByNumber(guint16 Group, guint16 Elem) 
-{
-   // Find out if the tag we encountered is in the dictionaries:
-   gdcmDictEntry *DictEntry = GetDictEntryByNumber(Group, Elem);
-   if (!DictEntry)
-      DictEntry = NewVirtualDictEntry(Group, Elem);
-
-   gdcmBinEntry *NewEntry = new gdcmBinEntry(DictEntry);
-   if (!NewEntry) 
-   {
-      dbg.Verbose(1, "gdcmDocument::NewBinEntryByNumber",
-                  "failed to allocate gdcmBinEntry");
-      return NULL;
-   }
-   return NewEntry;
-}
 
 /**
  * \brief   Generate a free TagKey i.e. a TagKey that is not present
@@ -2659,68 +2381,6 @@ guint32 gdcmDocument::GenerateFreeTagKeyInGroup(guint16 group)
 }
 
 
-/**
- * \brief   Searches both the public and the shadow dictionary (when they
- *          exist) for the presence of the DictEntry with given name.
- *          The public dictionary has precedence on the shadow one.
- * @param   Name name of the searched DictEntry
- * @return  Corresponding DictEntry when it exists, NULL otherwise.
- */
-gdcmDictEntry *gdcmDocument::GetDictEntryByName(std::string Name) 
-{
-   gdcmDictEntry *found = (gdcmDictEntry *)0;
-   if (!RefPubDict && !RefShaDict) 
-   {
-      dbg.Verbose(0, "gdcmDocument::GetDictEntry",
-                     "we SHOULD have a default dictionary");
-   }
-   if (RefPubDict) 
-   {
-      found = RefPubDict->GetDictEntryByName(Name);
-      if (found)
-         return found;
-   }
-   if (RefShaDict) 
-   {
-      found = RefShaDict->GetDictEntryByName(Name);
-      if (found)
-         return found;
-   }
-   return found;
-}
-
-/**
- * \brief   Searches both the public and the shadow dictionary (when they
- *          exist) for the presence of the DictEntry with given
- *          group and element. The public dictionary has precedence on the
- *          shadow one.
- * @param   group   group number of the searched DictEntry
- * @param   element element number of the searched DictEntry
- * @return  Corresponding DictEntry when it exists, NULL otherwise.
- */
-gdcmDictEntry *gdcmDocument::GetDictEntryByNumber(guint16 group,guint16 element) 
-{
-   gdcmDictEntry *found = (gdcmDictEntry *)0;
-   if (!RefPubDict && !RefShaDict) 
-   {
-      dbg.Verbose(0, "gdcmDocument::GetDictEntry",
-                     "we SHOULD have a default dictionary");
-   }
-   if (RefPubDict) 
-   {
-      found = RefPubDict->GetDictEntryByNumber(group, element);
-      if (found)
-         return found;
-   }
-   if (RefShaDict) 
-   {
-      found = RefShaDict->GetDictEntryByNumber(group, element);
-      if (found)
-         return found;
-   }
-   return found;
-}
-
 /**
  * \brief   Assuming the internal file pointer \ref gdcmDocument::fp 
  *          is placed at the beginning of a tag (TestGroup, TestElement),
@@ -2788,7 +2448,7 @@ guint32 gdcmDocument::ReadItemTagLength(void)
 }
 
 /**
- * \brief   Read the length of an exptected Sequence Delimiter tag i.e.
+ * \brief   Read the length of an expected Sequence Delimiter tag i.e.
  *          (0xfffe, 0xe0dd).
  * \sa      \ref gdcmDocument::ReadTagLength
  * \warning See warning of \ref gdcmDocument::ReadTagLength
@@ -2806,6 +2466,7 @@ guint32 gdcmDocument::ReadSequenceDelimiterTagLength(void)
  *          No other way so 'skip' the Data
  *
  */
+
 void gdcmDocument::Parse7FE0 (void)
 {
    gdcmDocEntry* Element = GetDocEntryByNumber(0x0002, 0x0010);
@@ -2898,6 +2559,8 @@ void gdcmDocument::Parse7FE0 (void)
    }
 }
 
+
+
 /**
  * \brief   Compares two documents, according to \ref gdcmDicomDir rules
  * \warning Does NOT work with ACR-NEMA files
@@ -2946,7 +2609,6 @@ bool gdcmDocument::operator<(gdcmDocument &document)
          }
       }
    }
-
    return false;
 }