]> Creatis software - gdcm.git/blobdiff - src/gdcmDataEntry.cxx
upgrades for 4DSplitter
[gdcm.git] / src / gdcmDataEntry.cxx
index 4f8937e9faf603d99db15640e62da7abaa3b076d..78f186c59b3f9f45ac7c13143ba78e09a5917367 100644 (file)
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: gdcmDataEntry.cxx,v $
   Language:  C++
-  Date:      $Date: 2005/10/25 14:27:37 $
-  Version:   $Revision: 1.11 $
+  Date:      $Date: 2011/03/29 07:36:00 $
+  Version:   $Revision: 1.57 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
 
 #include <fstream>
 
-namespace gdcm 
+#if defined(__BORLANDC__)
+ #include <mem.h>    // for memcpy
+ #include <stdlib.h> // for atof
+ #include <ctype.h>  // for isdigit
+#endif
+#include <string.h> // memcpy
+#include <stdlib.h> // atof
+
+// Could be defined like MAX_SIZE_LOAD_ELEMENT_VALUE
+#define GDCM_MAX_LENGTH_TO_CONVERT_TO_HEXA 8
+
+namespace GDCM_NAME_SPACE 
 {
 //-----------------------------------------------------------------------------
 #define MAX_SIZE_PRINT_ELEMENT_VALUE 0x7fffffff
@@ -34,17 +45,22 @@ uint32_t DataEntry::MaxSizePrintEntry = MAX_SIZE_PRINT_ELEMENT_VALUE;
 //-----------------------------------------------------------------------------
 // Constructor / Destructor
 /**
- * \brief   Constructor for a given DictEntry
- * @param   e Pointer to existing dictionary entry
+ * \brief   Constructor for a given DataEntry
+ * @param   group group number of the Data Entry to be created
+ * @param   elem element number of the Data Entry to be created
+ * @param   vr Value Representation of the Data Entry to be created 
  */
-DataEntry::DataEntry(DictEntry *e) 
-            : DocEntry(e)
+DataEntry::DataEntry(uint16_t group,uint16_t elem,
+                                     VRKey const &vr) 
+            : DocEntry(group,elem,vr)
 {
    State = STATE_LOADED;
    Flag = FLAG_NONE;
 
-   BinArea = 0;
-   SelfArea = true;
+   StrArea     = 0;
+   StrHexaArea = 0;   
+   BinArea     = 0;
+   SelfArea    = true;
 }
 
 /**
@@ -52,10 +68,12 @@ DataEntry::DataEntry(DictEntry *e)
  * @param   e Pointer to existing Doc entry
  */
 DataEntry::DataEntry(DocEntry *e)
-            : DocEntry(e->GetDictEntry())
+            //: DocEntry(e->GetDictEntry())
+            : DocEntry(e->GetGroup(),e->GetElement(), e->GetVR()  )
 {
-   Flag = 0;
+   Flag = FLAG_NONE;
    BinArea = 0;
+
    SelfArea = true;
 
    Copy(e);
@@ -69,13 +87,15 @@ DataEntry::~DataEntry ()
    DeleteBinArea();
 }
 
-//-----------------------------------------------------------------------------
-// Print
-
 //-----------------------------------------------------------------------------
 // Public
 /**
- * \brief Sets the value (non string) of the current Dicom Header Entry
+ * \brief Sets the value (non string) of the current DataEntry
+ * @param area area
+ * @param self self=true : The area : *belongs" to the DataEntry 
+ *                                  : will be delete with the DataEntry
+ *             self=false  The area *is not* deleted with the DataEntry
+ *              
  */
 void DataEntry::SetBinArea( uint8_t *area, bool self )  
 { 
@@ -87,7 +107,9 @@ void DataEntry::SetBinArea( uint8_t *area, bool self )
    State = STATE_LOADED;
 }
 /**
- * \brief Inserts the value (non string) into the current Dicom Header Entry
+ * \brief Inserts the value (non string) into the current DataEntry
+ * @param area area
+ * @param length length 
  */
 void DataEntry::CopyBinArea( uint8_t *area, uint32_t length )
 {
@@ -101,12 +123,110 @@ void DataEntry::CopyBinArea( uint8_t *area, uint32_t length )
       NewBinArea();
       memcpy(BinArea,area,length);
       if( length!=lgh )
-         BinArea[length]=0;
+         BinArea[length]=0; // padd with zero
 
       State = STATE_LOADED;
    }
 }
 
+/**
+ * \brief Checks wether the current DataEntry contains number(s)
+ */
+
+bool DataEntry::IsNumerical()
+{
+   const VRKey &vr = GetVR();
+
+   return 
+          vr == "DS" ||
+          vr == "FL" ||
+          vr == "FD" ||
+          vr == "IS" || 
+          vr == "SH" ||
+          vr == "SL" ||
+          vr == "SS" ||
+          vr == "UI" ||
+          vr == "UL" ||
+          vr == "US" ;      
+}
+
+/**
+ * \brief Gets a std::vector of 'double' holding the value(s) of any 'numerical' DataEntry
+ * @param valueVector std::vector double of value(s)
+ * \return false if VR not "a 'numerical" one
+ */
+ bool DataEntry::GetNumerical(std::vector <double> &valueVector)
+ {
+    valueVector.clear();
+    
+    if (!IsNumerical()) // never trust a user !
+       return false;    
+
+    const VRKey &vr = GetVR();
+    int loop;
+    if (vr == "IS" || vr == "DS")
+    {
+       /// \todo rewrite the whole method, in order *not to use* std::string !
+       std::vector<std::string> tokens;
+
+       Util::Tokenize ( GetString().c_str(), tokens, "\\" );
+
+       int nbValues= tokens.size();
+       if (nbValues == 0)
+          return false;
+
+       if (vr == "DS")        
+          for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(atof(tokens[loop].c_str()));
+       else
+          for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(atoi(tokens[loop].c_str()));
+
+       return true;     
+    }    
+
+    uint32_t nbValues = GetValueCount();
+    if (nbValues == 0)
+       return false;
+
+    if( vr == "US") {
+       for (loop=0; loop<nbValues; loop++)
+             valueVector.push_back(((uint16_t *)BinArea)[loop]);  
+       return true;
+    }
+    if( vr == "SS" ) {
+       for (loop=0; loop<nbValues; loop++)
+             valueVector.push_back(((int16_t *)BinArea)[loop]);  
+       return true;
+    }
+    if( vr == "UL") {
+       for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(((uint32_t *)BinArea)[loop]);      
+       return true;
+    }
+    if(vr == "SL" ) {
+       for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(((int32_t *)BinArea)[loop]);      
+       return true;
+    }       
+    if( vr == "FL" ) {
+       for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(((float *)BinArea)[loop]);     
+       return true;
+    }
+    if( vr == "FD" ) {
+       for (loop=0; loop<nbValues; loop++) 
+             valueVector.push_back(((double *)BinArea)[loop]);
+       return true;
+    }     
+    return false;   // ?!? 
+ }
+
+/**
+ * \brief Inserts the elementary (non string) value into the current DataEntry
+ * @param id index of the elementary value to be set
+ * @param val value, passed as a double 
+ */
 void DataEntry::SetValue(const uint32_t &id, const double &val)
 {
    if( !BinArea )
@@ -115,7 +235,7 @@ void DataEntry::SetValue(const uint32_t &id, const double &val)
 
    if( id > GetValueCount() )
    {
-      gdcmErrorMacro("Index (" << id << ")is greater than the data size");
+      gdcmErrorMacro("Index (" << id << ") is greater than the data size");
       return;
    }
 
@@ -150,30 +270,35 @@ void DataEntry::SetValue(const uint32_t &id, const double &val)
    }
 }
 /**
- * \brief returns, as a double (?!?) one of the values 
//      (when entry is multivaluated), identified by its index.
//      Returns 0.0 if index is wrong
- //     FIXME : warn the user there was a problem ! 
+ * \brief returns, as a double one of the values 
*      (when entry is multivaluated), identified by its index.
*      Returns 0.0 if index is wrong
+ * @param id id
  */
 double DataEntry::GetValue(const uint32_t &id) const
 {
    if( !BinArea )
    {
-      gdcmErrorMacro("BinArea not set. Can't get the value");
+      if (GetLength() != 0) // avoid stupid messages
+   /// \todo warn the user there was a problem !
+         gdcmErrorMacro("BinArea not set " << std::hex 
+                     << GetGroup() << " " << GetElement() 
+                     << " Can't get the value");
       return 0.0;
    }
 
    uint32_t count = GetValueCount();
    if( id > count )
    {
-      gdcmErrorMacro("Index (" << id << ")is greater than the data size");
+      gdcmErrorMacro("Index (" << id << ") is greater than the data size");
       return 0.0;
    }
 
-   // FIX the API : user *knows* that entry contains a US
-   //               and he receives a double ?!?
+   // if user *knows* that entry contains a US, 
+   // he just has to cast the double he receives
    
    const VRKey &vr = GetVR();
+
    if( vr == "US" || vr == "SS" )
       return ((uint16_t *)BinArea)[id];
    else if( vr == "UL" || vr == "SL" )
@@ -184,6 +309,7 @@ double DataEntry::GetValue(const uint32_t &id) const
       return ((double *)BinArea)[id];
    else if( Global::GetVR()->IsVROfStringRepresentable(vr) )
    {
+      // this is for VR = "DS", ...
       if( GetLength() )
       {
          // Don't use std::string to accelerate processing
@@ -232,19 +358,21 @@ double DataEntry::GetValue(const uint32_t &id) const
 /**
  * \brief Checks if the multiplicity of the value follows Dictionary VM
  */
-bool DataEntry::IsValueCountValid() const
+bool DataEntry::IsValueCountValid() /*const*/
 {
-  bool valid;
   uint32_t vm;
   const std::string &strVM = GetVM();
   uint32_t vc = GetValueCount();
+  bool valid = vc == 0;
+  if( valid )
+    return true;
   
   // FIXME : what shall we do with VM = "2-n", "3-n", etc
   
   if( strVM == "1-n" )
   {
-    // make sure it is at least one ??? FIXME
-    valid = vc >= 1 || vc == 0;
+    // make sure there is at least one ??? FIXME
+    valid = vc >= 1;
   }
   else
   {
@@ -252,7 +380,7 @@ bool DataEntry::IsValueCountValid() const
     os.str( strVM );
     os >> vm;
     // Two cases:
-    // vm respect the one from the dict
+    // vm respects the one from the dict
     // vm is 0 (we need to check if this element is allowed to be empty) FIXME
 
     // note  (JPR)
@@ -269,7 +397,7 @@ bool DataEntry::IsValueCountValid() const
     // Problem : entry type may depend on the modality and/or the Sequence
     //           it's embedded in !
     //          (Get the information in the 'Conformance Statements' ...)  
-    valid = vc == vm || vc == 0;
+    valid = vc == vm;
   }
   return valid;
 }
@@ -284,14 +412,15 @@ uint32_t DataEntry::GetValueCount( ) const
       return GetLength()/sizeof(uint16_t);
    else if( vr == "UL" || vr == "SL" )
       return GetLength()/sizeof(uint32_t);
-   else if( vr == "FL" )
-      return GetLength()/sizeof(float);
+   else if( vr == "FL" || vr == "OF" )
+      return GetLength()/4 ; // FL has a *4* length! sizeof(float);
    else if( vr == "FD" )
-      return GetLength()/sizeof(double);
+      return GetLength()/8;  // FD has a *8* length! sizeof(double);
    else if( Global::GetVR()->IsVROfStringRepresentable(vr) )
    {
       // Some element in DICOM are allowed to be empty
-      if( !GetLength() ) return 0;
+      if( !GetLength() ) 
+         return 0;
       // Don't use std::string to accelerate processing
       uint32_t count = 1;
       for(uint32_t i=0;i<GetLength();i++)
@@ -301,17 +430,41 @@ uint32_t DataEntry::GetValueCount( ) const
       }
       return count;
    }
-
    return GetLength();
 }
+
 /**
- * \brief Sets the 'value' of an Entry, passed as a std::string
+ * \brief Gets a std::vector of 'double' holding the value(s) of a DS DataEntry
+ * @param valueVector std::vector double of value(s)
+ * \return false if VR not "DS" or DataEntry empty
+ */
+ bool DataEntry::GetDSValue(std::vector <double> &valueVector)
+ {
+    /// \todo rewrite the whole method, in order *not to use* std::string !
+    std::vector<std::string> tokens;
+    
+    if (GetVR() != "DS") // never trust a user !
+       return false;    
+       
+    Util::Tokenize ( GetString().c_str(), tokens, "\\" );
+        
+    int nbValues= tokens.size();
+    if (nbValues == 0)
+       return false;
+               
+    for (int loop=0; loop<nbValues; loop++) 
+       valueVector.push_back(atof(tokens[loop].c_str()));
+    
+    return true;  
+ }
+/**
+ * \brief Sets the 'value' of a DataEntry, passed as a std::string
  * @param value string representation of the value to be set
  */ 
 void DataEntry::SetString(std::string const &value)
 {
    DeleteBinArea();
-
    const VRKey &vr = GetVR();
    if ( vr == "US" || vr == "SS" )
    {
@@ -362,27 +515,18 @@ void DataEntry::SetString(std::string const &value)
       tokens.clear();
    }
    else
-   {
-      if( value.size() > 0 )
-      {
-         // FIXME : should be quicker if we don't create one more std::string
-         //         just to make even the length of a char array ...
-
-         /*
-         std::string finalVal = Util::DicomString( value.c_str() );
-         SetLength(finalVal.size());
-         NewBinArea();
-
-         memcpy(BinArea, &(finalVal[0]), finalVal.size());
-         */
-
-         size_t l =  value.size();    
-         SetLength(l + l%2);
-         NewBinArea();
-         memcpy(BinArea, &(value[0]), l);
-         if (l%2)
-            BinArea[l] = '\0';
-      }
+   {      
+      size_t l =  value.size();    
+      SetLength(l + l%2);
+      NewBinArea();
+      memcpy(BinArea, value.c_str(), l);
+      if (l%2) // padded with blank except for UI
+        {
+        if ( vr == "UI" )
+          BinArea[l] = '\0';
+        else
+          BinArea[l] = ' ';
+        }
    }
    State = STATE_LOADED;
 }
@@ -391,75 +535,222 @@ void DataEntry::SetString(std::string const &value)
  */
 std::string const &DataEntry::GetString() const
 {
-   static std::ostringstream s;
-   const VRKey &vr = GetVR();
-
-   s.str("");
-   StrArea="";
-
-   if( !BinArea )
-      return StrArea;
-      
-   // When short integer(s) are stored, convert the following (n * 2) characters 
-   // as a displayable string, the values being separated by a back-slash
-
-   if( vr == "US" || vr == "SS" )
-   {
-      uint16_t *data=(uint16_t *)BinArea;
-
-      for (unsigned int i=0; i < GetValueCount(); i++)
-      {
-         if( i!=0 )
-            s << '\\';
-         s << data[i];
-      }
-      StrArea=s.str();
+  static std::ostringstream s;
+  const VRKey &vr = GetVR();
+  s.str("");
+  
+  if (!StrArea)
+     StrArea = new std::string();
+  else 
+     *StrArea="";
+
+  if( !BinArea )
+     return *StrArea;
+  // When short integer(s) are stored, convert the following (n * 2) characters
+  // as a displayable string, the values being separated by a back-slash
+  if( vr == "US" )
+  {
+     uint16_t *data=(uint16_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }
+  else if (vr == "SS" )
+  {
+     int16_t *data=(int16_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }      // See above comment on multiple short integers (mutatis mutandis).
+  else if( vr == "UL" )
+  {
+     uint32_t *data=(uint32_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }
+  else if( vr == "SL" )
+  {
+     int32_t *data=(int32_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }     // See above comment on multiple short integers (mutatis mutandis).
+  else if( vr == "FL" )
+  {
+     float *data=(float *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }     // See above comment on multiple short integers (mutatis mutandis).
+  else if( vr == "FD" )
+  {
+     double *data=(double *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        if( i!=0 )
+           s << '\\';
+        s << data[i];
+     }
+     *StrArea=s.str();
+  }
+  else
+  {
+     StrArea->append((const char *)BinArea,GetLength());
+     // to avoid gdcm to propagate oddities in lengthes
+     if ( GetLength()%2)
+        StrArea->append(" ",1);
    }
-   // See above comment on multiple short integers (mutatis mutandis).
-   else if( vr == "UL" || vr == "SL" )
-   {
-      uint32_t *data=(uint32_t *)BinArea;
+  return *StrArea;
+}
 
-      for (unsigned int i=0; i < GetValueCount(); i++)
-      {
+/**
+ * \brief   returns an hexadecimal representation of the DataEntry value
+ */
+std::string const &DataEntry::GetHexaRepresentation() const
+{ 
+  static std::ostringstream s2;
+  const VRKey &vr = GetVR();
+   
+  s2.str("");
+  if (!StrHexaArea)
+     StrHexaArea = new std::string();
+  else 
+     *StrHexaArea="";
+  if( !BinArea )
+     return *StrHexaArea;
+  // When short integer(s) are stored, convert the following (n * 2) characters
+  // as a displayable string, the values being separated by a back-slash
+  
+  s2 << std::hex;
+  
+  if( vr == "US" )
+  {
+     uint16_t *data=(uint16_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        s2  << std::setw( 2 ) << std::setfill( '0' );     
+        if( i!=0 )
+           s2 << '\\';
+        s2 << data[i];
+     }
+     *StrHexaArea=s2.str();
+  }
+  else if (vr == "SS" )
+  {
+     int16_t *data=(int16_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        s2  << std::setw( 4 ) << std::setfill( '0' );
+        if( i!=0 )
+           s2 << '\\';
+        s2  << data[i];
+     }
+     *StrHexaArea=s2.str();
+  }      // See above comment on multiple short integers (mutatis mutandis).
+  else if( vr == "UL" )
+  {
+     uint32_t *data=(uint32_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        s2  << std::setw( 4 ) << std::setfill( '0' );
+        if( i!=0 ) 
+           s2  << '\\';
+        s2 << data[i];
+     }
+     *StrHexaArea=s2.str();
+  }
+  else if( vr == "SL" )
+  {
+     int32_t *data=(int32_t *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        s2  << std::setw( 4 ) << std::setfill( '0' );
          if( i!=0 )
-            s << '\\';
-         s << data[i];
-      }
-      StrArea=s.str();
-   }
-   else if( vr == "FL" )
-   {
-      float *data=(float *)BinArea;
-
-      for (unsigned int i=0; i < GetValueCount(); i++)
-      {
+           s2  << '\\';
+        s2 << data[i];
+     }
+     *StrHexaArea=s2.str();
+  }
+  else if( vr == "FL" )
+  {
+     unsigned char *toto=(unsigned char *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+         s2.str("");
          if( i!=0 )
-            s << '\\';
-         s << data[i];
-      }
-      StrArea=s.str();
-   }
-   else if( vr == "FD" )
-   {
-      double *data=(double *)BinArea;
+           s2 << '\\';
+         unsigned int a4;
+         for(int iif=0; iif<4; iif++)
+         {
+            a4=toto[iif];
+            s2 << a4; 
+         }
+     }
+     *StrHexaArea=s2.str();
+  }
+  else if( vr == "FD" )
+  {
+     //double *data=(double *)BinArea;
+     unsigned char *toto=(unsigned char *)BinArea;
+     for (unsigned int i=0; i < GetValueCount(); i++)
+     {
+        s2.str("");
+        if( i!=0 )
+           s2 << '\\';
+        //s2 << data[i];
+
+         unsigned int a4;
+         for(int iid=0; iid<8; iid++)
+         {
+            a4=toto[iid];
+            s2 << a4; 
+         }
 
-      for (unsigned int i=0; i < GetValueCount(); i++)
-      {
-         if( i!=0 )
-            s << '\\';
-         s << data[i];
-      }
-      StrArea=s.str();
+     }
+     *StrHexaArea=s2.str();
+  }
+  else
+  {
+     unsigned int l = (Length > GDCM_MAX_LENGTH_TO_CONVERT_TO_HEXA) ? GDCM_MAX_LENGTH_TO_CONVERT_TO_HEXA : Length;
+     uint8_t *data=(uint8_t *)BinArea;
+     for (unsigned int i=0; i < l; i++)
+     {
+        if( i!=0 )
+           s2 << '\\';
+        s2 << std::setw( 2 ) << (int)(data[i]);
+     }
+     if (Length > 16)
+        s2 << "\\...";
+     *StrHexaArea=s2.str();
    }
-   else
-      StrArea.append((const char *)BinArea,GetLength());
-
-   return StrArea;
+  return *StrHexaArea;
 }
+
 /**
- * \brief   Copies all the attributes from an other DocEntry 
+ * \brief Copies all the attributes from an other DocEntry 
  * @param doc entry to copy from
+ * @remarks The content BinArea is copied too (StrArea is not)
  */
 void DataEntry::Copy(DocEntry *doc)
 {
@@ -470,100 +761,132 @@ void DataEntry::Copy(DocEntry *doc)
    {
       State = entry->State;
       Flag = entry->Flag;
-      CopyBinArea(entry->BinArea,entry->GetLength());
+      CopyBinArea(entry->BinArea,entry->GetLength());      
    }
 }
+
 /**
- * \brief   Writes the value of a DataEntry
+ * \brief   Writes the 'common part' + the 'value' area of a DataEntry
  * @param fp already open ofstream pointer
  * @param filetype type of the file (ACR, ImplicitVR, ExplicitVR, ...)
  */
-void DataEntry::WriteContent(std::ofstream *fp, FileType filetype)
+void DataEntry::WriteContent(std::ofstream *fp, FileType filetype, 
+                                                      bool insideMetaElements, bool insideSequence)
 { 
-   DocEntry::WriteContent(fp, filetype);
+   // writes the 'common part'
+   DocEntry::WriteContent(fp, filetype, insideMetaElements, insideSequence);
 
    if ( GetGroup() == 0xfffe )
    {
       return; //delimitors have NO value
    }
-
-   uint8_t *binArea8 = BinArea; //safe notation
-   size_t lgr = GetLength();
-   if (BinArea) // the binArea was *actually* loaded
-   {
-
-   //  The same operation should be done if we wanted 
-   //  to write image with Big Endian Transfer Syntax, 
-   //  while working on Little Endian Processor
+   
+   // --> We only deal with Little Endian writting.
    // --> forget Big Endian Transfer Syntax writting!
    //     Next DICOM version will give it up ...
-
-   // --> FIXME 
-   //    The stuff looks nice, but it's probably bugged,
-   //    since troubles occur on big endian processors (SunSparc, Motorola)
-   //    while reading the pixels of a 
-   //    gdcm-written Little-Endian 16 bits per pixel image
-
+   // WARNING - For Implicit VR private element,
+   //           we have *no choice* but considering them as
+   //           something like 'OB' values.
+   //           we rewrite them as we found them on disc.
+   //           Some trouble will occur if element was 
+   //           *actually* OW, if image was produced 
+   //           on Big endian based processor, read and writen 
+   //           on Little endian based processor
+   //           and, later on, somebody needs
+   //           this 'OW' Implicit VR private element (?!?)
+   //           (Same stuff, mutatis mutandis, for Little/Big)
+   // 8/16 bits Pixels problem should be solved automatiquely,
+   // since we ensure the VR (OB vs OW) is conform to Pixel size.
+        
+   uint8_t *data = BinArea; //safe notation
+   size_t l = GetLength(); 
+//   gdcmDebugMacro("in DataEntry::WriteContent " << GetKey() << " AtomicLength: "
+//              << Global::GetVR()->GetAtomicElementLength(this->GetVR() ) // << " BinArea in :" << &BinArea
+//             );
+   if (BinArea) // the binArea was *actually* loaded
+   {
 #if defined(GDCM_WORDS_BIGENDIAN) || defined(GDCM_FORCE_BIGENDIAN_EMULATION)
-
-      /// \todo FIXME : Right now, we only care of Pixels element
-      ///       we should deal with *all* the BinEntries
-      ///       Well, not really since we are not interpreting values read...
-
-      // 8 Bits Pixels *are* OB, 16 Bits Pixels *are* OW
-      // -value forced while Reading process-
-      
-      // -->  WARNING
-      // -->        the following lines *looked* very clever, 
-      // -->        but they don't work on big endian processors.
-      // -->        since I've no access for the moment to a big endian proc :-(
-      // -->        I comment them out, to see the result on the dash board 
-      // -->     
-      
-      // --> Revert to initial code : TestWriteSimple hangs on Darwin :-(     
-      if (GetGroup() == 0x7fe0 && GetVR() == "OW")
-      {  
-         uint16_t *binArea16 = (uint16_t*)binArea8;
-         binary_write (*fp, binArea16, lgr );
-      }
-      else
-      { 
-         // For any other VR, DataEntry is re-written as-is
-         binary_write (*fp, binArea8, lgr );
+      unsigned short vrLgth = 
+                        Global::GetVR()->GetAtomicElementLength(this->GetVR());
+      unsigned int i;
+      switch(vrLgth)
+      {
+         case 1:
+         {
+            binary_write (*fp, data, l );           
+            break;
+         }     
+         case 2:
+         {
+            uint16_t *data16 = (uint16_t *)data;
+            for(i=0;i<l/vrLgth;i++)
+               binary_write( *fp, data16[i]);
+            break;
+         }
+         case 4:
+         {
+            uint32_t *data32 = (uint32_t *)data;
+            for(i=0;i<l/vrLgth;i++)
+               binary_write( *fp, data32[i]);
+            break;
+         }
+         case 8:
+         {
+            double *data64 = (double *)data;
+            for(i=0;i<l/vrLgth;i++)
+               binary_write( *fp, data64[i]);
+            break;
+         }
       }
-
-      // -->  WARNING      
-      // -->         remove the following line, an uncomment the previous ones, 
-      // -->         if it doesn't work better
-      // -->     
-      /*binary_write ( *fp, binArea8, lgr ); // Elem value*/
-      
 #else
-      binary_write ( *fp, binArea8, lgr ); // Elem value
+   binary_write (*fp, data, l );
 #endif //GDCM_WORDS_BIGENDIAN
+
    }
    else
    {
-      // nothing was loaded, but we need to skip space on disc
-      
+      // nothing was loaded, but we need to skip space on disc     
+      if (l != 0)
+      {
       //  --> WARNING : nothing is written; 
       //  --> the initial data (on the the source image) is lost
-      //  --> user is *not* informed !
-      
-      fp->seekp(lgr, std::ios::cur);
+      //  --> user is *not* informed !      
+         gdcmDebugMacro ("Nothing was loaded, but we need to skip space on disc. "
+                      << "Length =" << l << " for " << GetKey() );   
+         fp->seekp(l, std::ios::cur); // At Write time, for unloaded elems
+      }
    }
+   // to avoid gdcm to propagate oddities
+   // (length was already modified)  
+   if (l%2)
+      fp->seekp(1, std::ios::cur);  // At Write time, for non even length elems
+}
+
+/**
+ * \brief   Compute the full length of the elementary DataEntry (not only value
+ *          length) depending on the VR.
+ */
+uint32_t DataEntry::ComputeFullLength()
+{
+   return GetFullLength();
 }
 
 //-----------------------------------------------------------------------------
 // Protected
-void DataEntry::NewBinArea(void)
+
+/// \brief Creates a DataEntry owned BinArea 
+///       (remove previous one if any and relevant StrArea if any)
+void DataEntry::NewBinArea( )
 {
    DeleteBinArea();
    if( GetLength() > 0 )
       BinArea = new uint8_t[GetLength()];
    SelfArea = true;
 }
-
+/// \brief Removes the BinArea, if owned by the DataEntry, 
+///        and the relevant StrArea if any
 void DataEntry::DeleteBinArea(void)
 {
    if (BinArea && SelfArea)
@@ -571,6 +894,16 @@ void DataEntry::DeleteBinArea(void)
       delete[] BinArea;
       BinArea = NULL;
    }
+   if (StrArea)
+   {
+      delete StrArea;
+      StrArea = 0;
+   }
+   if (StrHexaArea)
+   {
+      delete StrHexaArea;
+      StrHexaArea = 0;
+   }   
 }
 
 //-----------------------------------------------------------------------------
@@ -585,7 +918,9 @@ void DataEntry::DeleteBinArea(void)
  */
 void DataEntry::Print(std::ostream &os, std::string const & )
 {
-   os << "D ";
+   //os << "D ";
+   
+   // First, Print the common part (vr [length offset] name).
    DocEntry::Print(os);
 
    uint16_t g = GetGroup();
@@ -602,14 +937,9 @@ void DataEntry::Print(std::ostream &os, std::string const & )
       v = GetString();
       const VRKey &vr = GetVR();
 
-      if( vr == "US" || vr == "SS" )
-         s << " [" << GetString() << "]";
-      else if( vr == "UL" || vr == "SL" )
-         s << " [" << GetString() << "]";
-      else if ( vr == "FL" )
-         s << " [" << GetString() << "]";
-      else if ( vr == "FD" )
-         s << " [" << GetString() << "]";
+      if( vr == "US" || vr == "SS" || vr == "UL" || vr == "SL" 
+       || vr == "FL" || vr == "FD")
+         s << " [" << GetString() << "] =0x(" << GetHexaRepresentation() << ")";
       else
       { 
          if(Global::GetVR()->IsVROfStringRepresentable(vr))
@@ -628,7 +958,7 @@ void DataEntry::Print(std::ostream &os, std::string const & )
             }
             else
             {
-               s << " [gdcm::too long for print (" << cleanString.length() << ") ]";
+               s << " [GDCM_NAME_SPACE::too long for print (" << cleanString.length() << ") ]";
             }
          }
          else
@@ -648,7 +978,7 @@ void DataEntry::Print(std::ostream &os, std::string const & )
             else
             {
                s << " [" << GDCM_BINLOADED << ";"
-               << "length = " << GetLength() << "]";
+               << "length = " << GetLength() << "] =0x(" << GetHexaRepresentation() << ")";      
             }
          }
       }