]> Creatis software - gdcm.git/blobdiff - src/gdcmHeader.cxx
Just a few modif on Dicom Dict, according to 2003 Releaase.
[gdcm.git] / src / gdcmHeader.cxx
index 10f5700dac5c3445f5e847d8d1f4f801d823b978..4f2d2affb872a154c12ec4bd148249eb939ec6dd 100644 (file)
@@ -13,7 +13,7 @@
 #include "gdcmUtil.h"
 #include "gdcmHeader.h"
 
-#include <iddcmjpeg.h>
+#include "iddcmjpeg.h"
 
 // Refer to gdcmHeader::CheckSwap()
 #define HEADER_LENGTH_TO_READ       256
@@ -43,7 +43,7 @@ gdcmHeader::gdcmHeader(const char *InFilename, bool exception_on_error)
   else
     dbg.Error(!fp, "gdcmHeader::gdcmHeader cannot open file", InFilename);
   ParseHeader();
-  AddAndDefaultElements();
+  LoadElements();
 }
 
 
@@ -121,14 +121,7 @@ void gdcmHeader::CheckSwap()
    
    entCur = deb + 128;
    if(memcmp(entCur, "DICM", (size_t)4) == 0) {
-      filetype = TrueDicom;
       dbg.Verbose(1, "gdcmHeader::CheckSwap:", "looks like DICOM Version3");
-   } else {
-      filetype = Unknown;
-      dbg.Verbose(1, "gdcmHeader::CheckSwap:", "not a DICOM Version3 file");
-   }
-
-   if(filetype == TrueDicom) {
       // Next, determine the value representation (VR). Let's skip to the
       // first element (0002, 0000) and check there if we find "UL" 
       // - or "OB" if the 1st one is (0002,0001) -,
@@ -145,8 +138,12 @@ void gdcmHeader::CheckSwap()
       // * the 4 bytes of the first tag (0002, 0000),or (0002, 0001)
       // i.e. a total of  136 bytes.
       entCur = deb + 136;
+      // FIXME
+      // Use gdcmHeader::dicom_vr to test all the possibilities
+      // instead of just checking for UL, OB and UI !?
       if(  (memcmp(entCur, "UL", (size_t)2) == 0) ||
-          (memcmp(entCur, "OB", (size_t)2) == 0) )
+         (memcmp(entCur, "OB", (size_t)2) == 0) ||
+         (memcmp(entCur, "UI", (size_t)2) == 0) )
        {
          filetype = ExplicitVR;
          dbg.Verbose(1, "gdcmHeader::CheckSwap:",
@@ -172,11 +169,12 @@ void gdcmHeader::CheckSwap()
       rewind(fp);
       fseek (fp, 132L, SEEK_SET);
       return;
-   } // End of TrueDicom
+   } // End of DicomV3
 
    // Alas, this is not a DicomV3 file and whatever happens there is no file
    // preamble. We can reset the file position indicator to where the data
    // is (i.e. the beginning of the file).
+    dbg.Verbose(1, "gdcmHeader::CheckSwap:", "not a DICOM Version3 file");
    rewind(fp);
 
    // Our next best chance would be to be considering a 'clean' ACR/NEMA file.
@@ -216,6 +214,7 @@ void gdcmHeader::CheckSwap()
    // It is time for despaired wild guesses. So, let's assume this file
    // happens to be 'dirty' ACR/NEMA, i.e. the length of the group is
    // not present. Then the only info we have is the net2host one.
+   filetype = Unknown;
    if (! net2host )
       sw = 0;
    else
@@ -343,7 +342,7 @@ void gdcmHeader::FindVR( gdcmElValue *ElVal) {
  * @return  True when ImplicitVRLittleEndian found. False in all other cases.
  */
 bool gdcmHeader::IsImplicitVRLittleEndianTransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -361,7 +360,7 @@ bool gdcmHeader::IsImplicitVRLittleEndianTransferSyntax(void) {
  * @return  True when ExplicitVRLittleEndian found. False in all other cases.
  */
 bool gdcmHeader::IsExplicitVRLittleEndianTransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -379,7 +378,7 @@ bool gdcmHeader::IsExplicitVRLittleEndianTransferSyntax(void) {
  * @return  True when DeflatedExplicitVRLittleEndian found. False in all other cases.
  */
 bool gdcmHeader::IsDeflatedExplicitVRLittleEndianTransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -397,7 +396,7 @@ bool gdcmHeader::IsDeflatedExplicitVRLittleEndianTransferSyntax(void) {
  * @return  True when big endian found. False in all other cases.
  */
 bool gdcmHeader::IsExplicitVRBigEndianTransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -415,7 +414,7 @@ bool gdcmHeader::IsExplicitVRBigEndianTransferSyntax(void) {
  * @return  True when JPEGBaseLineProcess1found. False in all other cases.
  */
 bool gdcmHeader::IsJPEGBaseLineProcess1TransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -428,7 +427,7 @@ bool gdcmHeader::IsJPEGBaseLineProcess1TransferSyntax(void) {
 // faire qq chose d'intelligent a la place de ça
 
 bool gdcmHeader::IsJPEGLossless(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -447,7 +446,7 @@ bool gdcmHeader::IsJPEGLossless(void) {
  * @return  True when JPEGExtendedProcess2-4 found. False in all other cases.
  */
 bool gdcmHeader::IsJPEGExtendedProcess2_4TransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -465,7 +464,7 @@ bool gdcmHeader::IsJPEGExtendedProcess2_4TransferSyntax(void) {
  * @return  True when JPEGExtendedProcess3-5 found. False in all other cases.
  */
 bool gdcmHeader::IsJPEGExtendedProcess3_5TransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -484,7 +483,7 @@ bool gdcmHeader::IsJPEGExtendedProcess3_5TransferSyntax(void) {
  *          other cases.
  */
 bool gdcmHeader::IsJPEGSpectralSelectionProcess6_8TransferSyntax(void) {
-   gdcmElValue* Element = PubElVals.GetElementByNumber(0x0002, 0x0010);
+   gdcmElValue* Element = PubElValSet.GetElementByNumber(0x0002, 0x0010);
    if ( !Element )
       return false;
    LoadElementValueSafe(Element);
@@ -493,6 +492,17 @@ bool gdcmHeader::IsJPEGSpectralSelectionProcess6_8TransferSyntax(void) {
       return true;
    return false;
 }
+/**
+ * \ingroup gdcmHeader
+ * \brief   Predicate for dicom version 3 file.
+ * @return  True when the file is a dicom version 3.
+ */
+bool gdcmHeader::IsDicomV3(void) {
+   if (   (filetype == ExplicitVR)
+       || (filetype == ImplicitVR) )
+      return true;
+   return false;
+}
 
 /**
  * \ingroup gdcmHeader
@@ -644,7 +654,6 @@ void gdcmHeader::FindLength(gdcmElValue * ElVal) {
  * @return  The suggested integer.
  */
 guint32 gdcmHeader::SwapLong(guint32 a) {
-   // FIXME: il pourrait y avoir un pb pour les entiers negatifs ...
    switch (sw) {
    case    0 :
       break;
@@ -713,13 +722,14 @@ void gdcmHeader::LoadElementValue(gdcmElValue * ElVal) {
 
    fseek(fp, (long)ElVal->GetOffset(), SEEK_SET);
    
-   // Sequences not treated yet !
+   // FIXME Sequences not treated yet !
    //
    // Ne faudrait-il pas au contraire trouver immediatement
    // une maniere 'propre' de traiter les sequences (vr = SQ)
    // car commencer par les ignorer risque de conduire a qq chose
    // qui pourrait ne pas etre generalisable
-   //
+   // Well, I'm expecting your code !!!
+    
    if( vr == "SQ" )
       SkipLoad = true;
 
@@ -729,16 +739,7 @@ void gdcmHeader::LoadElementValue(gdcmElValue * ElVal) {
    if( group == 0xfffe )
       SkipLoad = true;
 
-   // The group length doesn't represent data to be loaded in memory, since
-   // each element of the group shall be loaded individualy.
-   if( elem == 0 )
-      //SkipLoad = true;       // modif sauvage JPR
-                               // On charge la longueur du groupe
-                               // quand l'element 0x0000 est présent !
-
    if ( SkipLoad ) {
-           // FIXME the following skip is not necessary
-      SkipElementValue(ElVal);
       ElVal->SetLength(0);
       ElVal->SetValue("gdcm::Skipped");
       return;
@@ -768,26 +769,6 @@ void gdcmHeader::LoadElementValue(gdcmElValue * ElVal) {
        // pour les elements de Value Multiplicity > 1
        // on aura en fait une serie d'entiers
        
-       // code original
-       
-       //if ( IsAnInteger(ElVal) ) {
-       //      guint32 NewInt;
-       //      if( length == 2 ) {
-       //              NewInt = ReadInt16();
-       //      } else if( length == 4 ) {
-       //
-       //      NewInt = ReadInt32();
-       //      } else
-       //              dbg.Error(true, "LoadElementValue: Inconsistency when reading Int.");
-       //      
-       //      //FIXME: make the following an util fonction
-       //      ostringstream s;
-       //      s << NewInt;
-       //      ElVal->SetValue(s.str());
-       //      return;
-       //}
-
-       // modif proposee. JPR
        // on devrait pouvoir faire + compact (?)
                
        if ( IsAnInteger(ElVal) ) {
@@ -803,7 +784,6 @@ void gdcmHeader::LoadElementValue(gdcmElValue * ElVal) {
                                        s << '\\';
                                        NewInt = ReadInt16();
                                        s << NewInt;
-                                       //printf("%s\n", s.str().c_str());
                                }
                        }
                        
@@ -823,8 +803,6 @@ void gdcmHeader::LoadElementValue(gdcmElValue * ElVal) {
                return; 
        }
    
-   
-   
    // FIXME The exact size should be length if we move to strings or whatever
    char* NewValue = (char*)malloc(length+1);
    if( !NewValue) {
@@ -886,6 +864,16 @@ guint32 gdcmHeader::ReadInt32(void) {
    return g;
 }
 
+
+
+
+//
+// TODO : JPR Pour des raisons d'homogeneité de noms, 
+// remplacer les quelques GetxxxByKey 
+// par des GetxxxByNumber, lorsqu'on passe gr, el en param.
+// il peut etre interessant de rajouter des GetxxxByKey, auxquels on passe *vraiment* une TagKey
+//
+
 /**
  * \ingroup gdcmHeader
  * \brief   Build a new Element Value from all the low level arguments. 
@@ -909,6 +897,20 @@ gdcmElValue* gdcmHeader::NewElValueByKey(guint16 Group, guint16 Elem) {
    return NewElVal;
 }
 
+/**
+ * \ingroup gdcmHeader
+ * \brief   TODO
+ * @param   
+ */
+int gdcmHeader::ReplaceOrCreateByNumber(string Value, guint16 Group, guint16 Elem ) {
+
+       gdcmElValue* nvElValue=NewElValueByKey(Group, Elem);
+       PubElValSet.Add(nvElValue);     
+       PubElValSet.SetElValueByNumber(Value, Group, Elem);
+       return(1);
+}   
+
+
 /**
  * \ingroup gdcmHeader
  * \brief   Build a new Element Value from all the low level arguments. 
@@ -937,8 +939,8 @@ gdcmElValue* gdcmHeader::NewElValueByName(string Name) {
  * @return  On succes the newly created ElValue, NULL on failure.      
  */
 gdcmElValue * gdcmHeader::ReadNextElement(void) {
-   guint16 g;
-   guint16 n;
+  
+   guint16 g,n;
    gdcmElValue * NewElVal;
    
    g = ReadInt16();
@@ -977,48 +979,17 @@ bool gdcmHeader::IsAnInteger(gdcmElValue * ElVal) {
    if ( element == 0 )  {  // This is the group length of the group
       if (length == 4)
          return true;
-      else
+      else {
+         printf("Erroneous Group Length element length %d\n",length);
+                    
          dbg.Error("gdcmHeader::IsAnInteger",
-                   "Erroneous Group Length element length.");
+                   "Erroneous Group Length element length.");     
+      }
    }
  
-   /*          
-   // on le traite tt de même (VR peut donner l'info)
-               // faire qq chose + ruse (pas de test si pas de VR)  
-   if ( group % 2 != 0 )
-      // We only have some semantics on documented elements, which are
-      // the even ones.
-      return false; 
-    */
-   
-   /*
-   if ( (length != 4) && ( length != 2) )
-      // Swapping only make sense on integers which are 2 or 4 bytes long.
-               
-               // En fait, pour les entiers de 'Value Multiplicity' supérieur a 1
-               // la longueur n'est pas forcement 2 ou 4 
-               // ET il faudra swapper.
-      return false;
-    */
-   
    if ( (vr == "UL") || (vr == "US") || (vr == "SL") || (vr == "SS") )
       return true;
    
-   
-   // est-ce encore utile?
-   // mieux vaut modifier le source du Dicom Dictionnaty 
-   // et remplacer pour ces 2 cas  RET par US
-   
-   if ( (group == 0x0028) && (element == 0x0005) )
-      // The "Image Dimensions" tag is retained from ACR/NEMA and contains
-      // the number of dimensions of the contained object (1 for Signal,
-      // 2 for Image, 3 for Volume, 4 for Sequence).
-      return true;
-   
-   if ( (group == 0x0028) && (element == 0x0200) )
-      // This tag is retained from ACR/NEMA
-      return true;
-   
    return false;
 }
 
@@ -1033,7 +1004,7 @@ size_t gdcmHeader::GetPixelOffset(void) {
    // is found by indirection through the "Image Location").
    // Inside the group pointed by "Image Location" the searched element
    // is conventionally the element 0x0010 (when the norm is respected).
-   //    When the "Image Location" is absent we default to group 0x7fe0.
+   // When the "Image Location" is absent we default to group 0x7fe0.
    guint16 grPixel;
    guint16 numPixel;
    string ImageLocation = GetPubElValByName("Image Location");
@@ -1043,12 +1014,11 @@ size_t gdcmHeader::GetPixelOffset(void) {
       grPixel = (guint16) atoi( ImageLocation.c_str() );
    }
    if (grPixel != 0x7fe0)
-      // FIXME is this still necessary ?
-      // Now, this looks like an old dirty fix for Philips imager
+      // This is a kludge for old dirty Philips imager.
       numPixel = 0x1010;
    else
       numPixel = 0x0010;
-   gdcmElValue* PixelElement = PubElVals.GetElementByNumber(grPixel, numPixel);
+   gdcmElValue* PixelElement = PubElValSet.GetElementByNumber(grPixel, numPixel);
    if (PixelElement)
       return PixelElement->GetOffset();
    else
@@ -1068,7 +1038,7 @@ size_t gdcmHeader::GetPixelOffset(void) {
 gdcmDictEntry * gdcmHeader::GetDictEntryByKey(guint16 group, guint16 element) {
    gdcmDictEntry * found = (gdcmDictEntry*)0;
    if (!RefPubDict && !RefShaDict) {
-      dbg.Verbose(0, "FIXME in gdcmHeader::GetDictEntry",
+      dbg.Verbose(0, "gdcmHeader::GetDictEntry",
                      "we SHOULD have a default dictionary");
    }
    if (RefPubDict) {
@@ -1095,7 +1065,7 @@ gdcmDictEntry * gdcmHeader::GetDictEntryByKey(guint16 group, guint16 element) {
 gdcmDictEntry * gdcmHeader::GetDictEntryByName(string Name) {
    gdcmDictEntry * found = (gdcmDictEntry*)0;
    if (!RefPubDict && !RefShaDict) {
-      dbg.Verbose(0, "FIXME in gdcmHeader::GetDictEntry",
+      dbg.Verbose(0, "gdcmHeader::GetDictEntry",
                      "we SHOULD have a default dictionary");
    }
    if (RefPubDict) {
@@ -1121,7 +1091,7 @@ gdcmDictEntry * gdcmHeader::GetDictEntryByName(string Name) {
  *          "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetPubElValByNumber(guint16 group, guint16 element) {
-   return PubElVals.GetElValueByNumber(group, element);
+   return PubElValSet.GetElValueByNumber(group, element);
 }
 
 /**
@@ -1139,7 +1109,7 @@ string gdcmHeader::GetPubElValByNumber(guint16 group, guint16 element) {
  *          and the string "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetPubElValRepByNumber(guint16 group, guint16 element) {
-   gdcmElValue* elem =  PubElVals.GetElementByNumber(group, element);
+   gdcmElValue* elem =  PubElValSet.GetElementByNumber(group, element);
    if ( !elem )
       return "gdcm::Unfound";
    return elem->GetVR();
@@ -1154,7 +1124,7 @@ string gdcmHeader::GetPubElValRepByNumber(guint16 group, guint16 element) {
  *          "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetPubElValByName(string TagName) {
-   return PubElVals.GetElValueByName(TagName);
+   return PubElValSet.GetElValueByName(TagName);
 }
 
 /**
@@ -1171,7 +1141,7 @@ string gdcmHeader::GetPubElValByName(string TagName) {
  *          and the string "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetPubElValRepByName(string TagName) {
-   gdcmElValue* elem =  PubElVals.GetElementByName(TagName);
+   gdcmElValue* elem =  PubElValSet.GetElementByName(TagName);
    if ( !elem )
       return "gdcm::Unfound";
    return elem->GetVR();
@@ -1187,7 +1157,7 @@ string gdcmHeader::GetPubElValRepByName(string TagName) {
  *          and the string "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetShaElValByNumber(guint16 group, guint16 element) {
-   return ShaElVals.GetElValueByNumber(group, element);
+   return ShaElValSet.GetElValueByNumber(group, element);
 }
 
 /**
@@ -1205,7 +1175,7 @@ string gdcmHeader::GetShaElValByNumber(guint16 group, guint16 element) {
  *          and the string "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetShaElValRepByNumber(guint16 group, guint16 element) {
-   gdcmElValue* elem =  ShaElVals.GetElementByNumber(group, element);
+   gdcmElValue* elem =  ShaElValSet.GetElementByNumber(group, element);
    if ( !elem )
       return "gdcm::Unfound";
    return elem->GetVR();
@@ -1220,7 +1190,7 @@ string gdcmHeader::GetShaElValRepByNumber(guint16 group, guint16 element) {
  *          "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetShaElValByName(string TagName) {
-   return ShaElVals.GetElValueByName(TagName);
+   return ShaElValSet.GetElValueByName(TagName);
 }
 
 /**
@@ -1237,7 +1207,7 @@ string gdcmHeader::GetShaElValByName(string TagName) {
  *          and the string "gdcm::Unfound" otherwise.
  */
 string gdcmHeader::GetShaElValRepByName(string TagName) {
-   gdcmElValue* elem =  ShaElVals.GetElementByName(TagName);
+   gdcmElValue* elem =  ShaElValSet.GetElementByName(TagName);
    if ( !elem )
       return "gdcm::Unfound";
    return elem->GetVR();
@@ -1321,7 +1291,7 @@ string gdcmHeader::GetElValRepByName(string TagName) {
 
 /**
  * \ingroup gdcmHeader
- * \brief   Accesses an existing gdcmElValue in the PubElVals of this instance
+ * \brief   Accesses an existing gdcmElValue in the PubElValSet of this instance
  *          through it's (group, element) and modifies it's content with
  *          the given value.
  * @param   content new value to substitute with
@@ -1331,60 +1301,61 @@ string gdcmHeader::GetElValRepByName(string TagName) {
 int gdcmHeader::SetPubElValByNumber(string content, guint16 group,
                                     guint16 element)
 {
-   return (  PubElVals.SetElValueByNumber (content, group, element) );
+   return (  PubElValSet.SetElValueByNumber (content, group, element) );
 }
 
 /**
  * \ingroup gdcmHeader
- * \brief   Accesses an existing gdcmElValue in the PubElVals of this instance
+ * \brief   Accesses an existing gdcmElValue in the PubElValSet of this instance
  *          through tag name and modifies it's content with the given value.
  * @param   content new value to substitute with
  * @param   TagName name of the tag to be modified
  */
 int gdcmHeader::SetPubElValByName(string content, string TagName) {
-   return (  PubElVals.SetElValueByName (content, TagName) );
+   return (  PubElValSet.SetElValueByName (content, TagName) );
 }
 
 /**
  * \ingroup gdcmHeader
- * \brief   Accesses an existing gdcmElValue in the PubElVals of this instance
+ * \brief   Accesses an existing gdcmElValue in the PubElValSet of this instance
  *          through it's (group, element) and modifies it's length with
  *          the given value.
- *             NOT FOR BOZOs !
- * @param   contents new length to substitute with
+ * \warning Use with extreme caution.
+ * @param   length new length to substitute with
  * @param   group   group of the ElVal to modify
  * @param   element element of the ElVal to modify
+ * @return  1 on success, 0 otherwise.
  */
-int gdcmHeader::SetPubElValLengthByNumber(guint32 lgr, guint16 group,
-                                    guint16 element)
-{
-       return (  PubElVals.SetElValueLengthByNumber (lgr, group, element) );
+
+int gdcmHeader::SetPubElValLengthByNumber(guint32 length, guint16 group,
+                                    guint16 element) {
+       return (  PubElValSet.SetElValueLengthByNumber (length, group, element) );
 }
 
 /**
  * \ingroup gdcmHeader
- * \brief   Accesses an existing gdcmElValue in the ShaElVals of this instance
+ * \brief   Accesses an existing gdcmElValue in the ShaElValSet of this instance
  *          through it's (group, element) and modifies it's content with
  *          the given value.
  * @param   content new value to substitute with
  * @param   group   group of the ElVal to modify
  * @param   element element of the ElVal to modify
+ * @return  1 on success, 0 otherwise.
  */
 int gdcmHeader::SetShaElValByNumber(string content,
-                                    guint16 group, guint16 element)
-{
-   return (  ShaElVals.SetElValueByNumber (content, group, element) );
+                                    guint16 group, guint16 element) {
+   return (  ShaElValSet.SetElValueByNumber (content, group, element) );
 }
 
 /**
  * \ingroup gdcmHeader
- * \brief   Accesses an existing gdcmElValue in the ShaElVals of this instance
+ * \brief   Accesses an existing gdcmElValue in the ShaElValSet of this instance
  *          through tag name and modifies it's content with the given value.
  * @param   content new value to substitute with
  * @param   TagName name of the tag to be modified
  */
 int gdcmHeader::SetShaElValByName(string content, string TagName) {
-   return (  ShaElVals.SetElValueByName (content, TagName) );
+   return (  ShaElValSet.SetElValueByName (content, TagName) );
 }
 
 /**
@@ -1398,36 +1369,159 @@ void gdcmHeader::ParseHeader(bool exception_on_error) throw(gdcmFormatError) {
    CheckSwap();
    while ( (newElValue = ReadNextElement()) ) {
       SkipElementValue(newElValue);
-      PubElVals.Add(newElValue);
+      PubElValSet.Add(newElValue);
    }
 }
 
 /**
  * \ingroup gdcmHeader
- * \brief   Once the header is parsed add some gdcm convenience/helper elements
- *          in the gdcmElValSet. For example add:
- *          - gdcmImageType which is an entry containing a short for the
- *            type of image and whose value ranges in 
- *               I8   (unsigned 8 bit image)
- *               I16  (unsigned 8 bit image)
- *               IS16 (signed 8 bit image)
- *          - gdcmXsize, gdcmYsize, gdcmZsize whose values are respectively
- *            the ones of the official DICOM fields Rows, Columns and Planes.
+ * \brief   Retrieve the number of columns of image.
+ * @return  The encountered size when found, 0 by default.
  */
-void gdcmHeader::AddAndDefaultElements(void) {
-   gdcmElValue* NewEntry = (gdcmElValue*)0;
+int gdcmHeader::GetXSize(void) {
+   // We cannot check for "Columns" because the "Columns" tag is present
+   // both in IMG (0028,0011) and OLY (6000,0011) sections of the dictionary.
+   string StrSize = GetPubElValByNumber(0x0028,0x0011);
+   if (StrSize == "gdcm::Unfound")
+      return 0;
+   return atoi(StrSize.c_str());
+}
+
+/**
+ * \ingroup gdcmHeader
+ * \brief   Retrieve the number of lines of image.
+ * \warning The defaulted value is 1 as opposed to gdcmHeader::GetXSize()
+ * @return  The encountered size when found, 1 by default.
+ */
+int gdcmHeader::GetYSize(void) {
+   // We cannot check for "Rows" because the "Rows" tag is present
+   // both in IMG (0028,0010) and OLY (6000,0010) sections of the dictionary.
+   string StrSize = GetPubElValByNumber(0x0028,0x0010);
+   if (StrSize != "gdcm::Unfound")
+      return atoi(StrSize.c_str());
+   if ( IsDicomV3() )
+      return 0;
+   else
+      // The Rows (0028,0010) entry is optional for ACR/NEMA. It might
+      // hence be a signal (1d image). So we default to 1:
+      return 1;
+}
+
+/**
+ * \ingroup gdcmHeader
+ * \brief   Retrieve the number of planes of volume or the number
+ *          of frames of a multiframe.
+ * \warning When present we consider the "Number of Frames" as the third
+ *          dimension. When absent we consider the third dimension as
+ *          being the "Planes" tag content.
+ * @return  The encountered size when found, 1 by default.
+ */
+int gdcmHeader::GetZSize(void) {
+   // Both in DicomV3 and ACR/Nema the consider the "Number of Frames"
+   // as the third dimension.
+   string StrSize = GetPubElValByNumber(0x0028,0x0008);
+   if (StrSize != "gdcm::Unfound")
+      return atoi(StrSize.c_str());
+
+   // We then consider the "Planes" entry as the third dimension [we
+   // cannot retrieve by name since "Planes tag is present both in
+   // IMG (0028,0012) and OLY (6000,0012) sections of the dictionary]. 
+   StrSize = GetPubElValByNumber(0x0028,0x0012);
+   if (StrSize != "gdcm::Unfound")
+      return atoi(StrSize.c_str());
+   return 1;
+}
+
+/**
+ * \ingroup gdcmHeader
+ * \brief   Build the Pixel Type of the image.
+ *          Possible values are:
+ *          - U8  unsigned  8 bit,
+ *          - S8    signed  8 bit,
+ *          - U16 unsigned 16 bit,
+ *          - S16   signed 16 bit,
+ *          - U32 unsigned 32 bit,
+ *          - S32   signed 32 bit,
+ * \warning 12 bit images appear as 16 bit.
+ * @return 
+ */
+string gdcmHeader::GetPixelType(void) {
+   string BitsAlloc;
+   BitsAlloc = GetElValByName("Bits Allocated");
+   if (BitsAlloc == "gdcm::Unfound") {
+      dbg.Verbose(0, "gdcmHeader::GetPixelType: unfound Bits Allocated");
+      BitsAlloc = string("16");
+   }
+   if (BitsAlloc == "12")
+      BitsAlloc = string("16");
+
+   string Signed;
+   Signed = GetElValByName("Pixel Representation");
+   if (Signed == "gdcm::Unfound") {
+      dbg.Verbose(0, "gdcmHeader::GetPixelType: unfound Pixel Representation");
+      BitsAlloc = string("0");
+   }
+   if (Signed == "0")
+      Signed = string("U");
+   else
+      Signed = string("S");
 
-   NewEntry = NewElValueByName("gdcmXSize");
-   NewEntry->SetValue(GetElValByName("Rows"));
-   PubElVals.Add(NewEntry);
+   return( BitsAlloc + Signed);
+}
+
+/**
+ * \ingroup gdcmHeader
+ * \brief  This predicate, based on hopefully reasonnable heuristics,
+ *         decides whether or not the current gdcmHeader was properly parsed
+ *         and contains the mandatory information for being considered as
+ *         a well formed and usable image.
+ * @return true when gdcmHeader is the one of a reasonable Dicom file,
+ *         false otherwise. 
+ */
+bool gdcmHeader::IsReadable(void) {
+   if (   GetElValByName("Image Dimensions") != "gdcm::Unfound"
+      && atoi(GetElValByName("Image Dimensions").c_str()) > 4 ) {
+      return false;
+   }
+   if (  GetElValByName("Bits Allocated") == "gdcm::Unfound" )
+      return false;
+   if (  GetElValByName("Bits Stored") == "gdcm::Unfound" )
+      return false;
+   if (  GetElValByName("High Bit") == "gdcm::Unfound" )
+      return false;
+   if (  GetElValByName("Pixel Representation") == "gdcm::Unfound" )
+      return false;
+   return true;
+}
 
-   NewEntry = NewElValueByName("gdcmYSize");
-   NewEntry->SetValue(GetElValByName("Columns"));
-   PubElVals.Add(NewEntry);
+/**
+ * \ingroup gdcmHeader
+ * \brief   Small utility function that creates a new manually crafted
+ *          (as opposed as read from the file) gdcmElValue with user
+ *          specified name and adds it to the public tag hash table.
+ * \note    A fake TagKey is generated so the PubDict can keep it's coherence.
+ * @param   NewTagName The name to be given to this new tag.
+ * @param   VR The Value Representation to be given to this new tag.
+ * @ return The newly hand crafted Element Value.
+ */
+gdcmElValue* gdcmHeader::NewManualElValToPubDict(string NewTagName, string VR) {
+   gdcmElValue* NewElVal = (gdcmElValue*)0;
+   guint32 StuffGroup = 0xffff;   // Group to be stuffed with additional info
+   guint32 FreeElem = 0;
+   gdcmDictEntry* NewEntry = (gdcmDictEntry*)0;
+
+   FreeElem = PubElValSet.GenerateFreeTagKeyInGroup(StuffGroup);
+   if (FreeElem == UINT32_MAX) {
+      dbg.Verbose(1, "gdcmHeader::NewManualElValToPubDict",
+                     "Group 0xffff in Public Dict is full");
+      return (gdcmElValue*)0;
+   }
+   NewEntry = new gdcmDictEntry(StuffGroup, FreeElem,
+                                VR, "GDCM", NewTagName);
+   NewElVal = new gdcmElValue(NewEntry);
+   PubElValSet.Add(NewElVal);
+   return NewElVal;
 
-   NewEntry = NewElValueByName("gdcmZSize");
-   NewEntry->SetValue(GetElValByName("Planes"));
-   PubElVals.Add(NewEntry);
 }
 
 /**
@@ -1437,16 +1531,20 @@ void gdcmHeader::AddAndDefaultElements(void) {
  */
 void gdcmHeader::LoadElements(void) {
    rewind(fp);   
-   TagElValueHT ht = PubElVals.GetTagHt();
+   TagElValueHT ht = PubElValSet.GetTagHt();
    for (TagElValueHT::iterator tag = ht.begin(); tag != ht.end(); ++tag) {
       LoadElementValue(tag->second);
       }
 }
 
 void gdcmHeader::PrintPubElVal(ostream & os) {
-   PubElVals.Print(os);
+   PubElValSet.Print(os);
 }
 
 void gdcmHeader::PrintPubDict(ostream & os) {
    RefPubDict->Print(os);
 }
+
+int gdcmHeader::Write(FILE * fp, FileType type) {
+   return PubElValSet.Write(fp, type);
+}