]> Creatis software - gdcm.git/blobdiff - src/gdcmFile.cxx
Coding Style + Doxygenation
[gdcm.git] / src / gdcmFile.cxx
index 60c86429471f044929d325fd96a35d72c336410d..cd6ec36b48bc77078dffe8723c58d6c785a33137 100644 (file)
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: gdcmFile.cxx,v $
   Language:  C++
-  Date:      $Date: 2004/11/25 10:24:34 $
-  Version:   $Revision: 1.163 $
+  Date:      $Date: 2005/01/06 20:03:27 $
+  Version:   $Revision: 1.183 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
 =========================================================================*/
 
 #include "gdcmFile.h"
+#include "gdcmDocument.h"
 #include "gdcmDebug.h"
 #include "gdcmUtil.h"
 #include "gdcmBinEntry.h"
+#include "gdcmHeader.h"
+#include "gdcmPixelReadConvert.h"
+#include "gdcmPixelWriteConvert.h"
+#include "gdcmDocEntryArchive.h"
+
 #include <fstream>
 
 namespace gdcm 
@@ -28,6 +34,26 @@ typedef std::pair<TagDocEntryHT::iterator,TagDocEntryHT::iterator> IterHT;
 
 //-------------------------------------------------------------------------
 // Constructor / Destructor
+/**
+ * \brief Constructor dedicated to deal with the *pixels* area of a ACR/DICOMV3
+ *        file (Header only deals with the ... header)
+ *        Opens (in read only and when possible) an existing file and checks
+ *        for DICOM compliance. Returns NULL on failure.
+ *        It will be up to the user to load the pixels into memory
+ *        (see GetImageData, GetImageDataRaw)
+ * \note  the in-memory representation of all available tags found in
+ *        the DICOM header is post-poned to first header information access.
+ *        This avoid a double parsing of public part of the header when
+ *        one sets an a posteriori shadow dictionary (efficiency can be
+ *        seen as a side effect).   
+ */
+File::File( )
+{
+   HeaderInternal = new Header( );
+   SelfHeader = true;
+   Initialise();
+}
+
 /**
  * \brief Constructor dedicated to deal with the *pixels* area of a ACR/DICOMV3
  *        file (Header only deals with the ... header)
@@ -70,26 +96,6 @@ File::File(std::string const & filename )
    Initialise();
 }
 
-/**
- * \brief Factorization for various forms of constructors.
- */
-void File::Initialise()
-{
-   WriteMode = WMODE_DECOMPRESSED;
-   WriteType = WTYPE_IMPL_VR;
-
-   PixelConverter = new PixelConvert;
-   Archive = new DocEntryArchive( HeaderInternal );
-
-   if ( HeaderInternal->IsReadable() )
-   {
-      PixelConverter->GrabInformationsFromHeader( HeaderInternal );
-   }
-
-   Pixel_Data = 0;
-   ImageDataSize = 0;
-}
-
 /**
  * \brief canonical destructor
  * \note  If the Header was created by the File constructor,
@@ -97,9 +103,13 @@ void File::Initialise()
  */
 File::~File()
 { 
-   if( PixelConverter )
+   if( PixelReadConverter )
    {
-      delete PixelConverter;
+      delete PixelReadConverter;
+   }
+   if( PixelWriteConverter )
+   {
+      delete PixelWriteConverter;
    }
    if( Archive )
    {
@@ -115,6 +125,14 @@ File::~File()
 
 //-----------------------------------------------------------------------------
 // Print
+void File::Print(std::ostream &os)
+{
+   HeaderInternal->SetPrintLevel(PrintLevel);
+   HeaderInternal->Print(os);
+
+   PixelReadConverter->SetPrintLevel(PrintLevel);
+   PixelReadConverter->Print(os);
+}
 
 //-----------------------------------------------------------------------------
 // Public
@@ -127,7 +145,12 @@ File::~File()
  */
 size_t File::GetImageDataSize()
 {
-   return PixelConverter->GetRGBSize();
+   if ( PixelWriteConverter->GetUserData() )
+   {
+      return PixelWriteConverter->GetUserDataSize();
+   }
+
+   return PixelReadConverter->GetRGBSize();
 }
 
 /**
@@ -139,7 +162,12 @@ size_t File::GetImageDataSize()
  */
 size_t File::GetImageDataRawSize()
 {
-   return PixelConverter->GetDecompressedSize();
+   if ( PixelWriteConverter->GetUserData() )
+   {
+      return PixelWriteConverter->GetUserDataSize();
+   }
+
+   return PixelReadConverter->GetRawSize();
 }
 
 /**
@@ -152,25 +180,44 @@ size_t File::GetImageDataRawSize()
  * @return  Pointer to newly allocated pixel data.
  *          NULL if alloc fails 
  */
-uint8_tFile::GetImageData()
+uint8_t *File::GetImageData()
 {
-   if ( ! GetDecompressed() )
+   if ( PixelWriteConverter->GetUserData() )
+   {
+      return PixelWriteConverter->GetUserData();
+   }
+
+   if ( ! GetRaw() )
    {
       // If the decompression failed nothing can be done.
       return 0;
    }
 
-   if ( HeaderInternal->HasLUT() && PixelConverter->BuildRGBImage() )
+   if ( HeaderInternal->HasLUT() && PixelReadConverter->BuildRGBImage() )
    {
-      return PixelConverter->GetRGB();
+      return PixelReadConverter->GetRGB();
    }
    else
    {
-      // When no LUT or LUT conversion fails, return the decompressed
-      return PixelConverter->GetDecompressed();
+      // When no LUT or LUT conversion fails, return the Raw
+      return PixelReadConverter->GetRaw();
    }
 }
 
+/**
+ * \brief   Allocates necessary memory, 
+ *          Transforms YBR pixels (if any) into RGB pixels
+ *          Transforms 3 planes R, G, B  (if any) into a single RGB Plane
+ *          Copies the pixel data (image[s]/volume[s]) to newly allocated zone. 
+ *          DOES NOT transform Grey plane + 3 Palettes into a RGB Plane
+ * @return  Pointer to newly allocated pixel data.
+ * \        NULL if alloc fails 
+ */
+uint8_t *File::GetImageDataRaw ()
+{
+   return GetRaw();
+}
+
 /**
  * \brief
  *          Read the pixels from disk (uncompress if necessary),
@@ -196,80 +243,43 @@ uint8_t* File::GetImageData()
  * @return  On success, the number of bytes actually copied. Zero on
  *          failure e.g. MaxSize is lower than necessary.
  */
-size_t File::GetImageDataIntoVector (voiddestination, size_t maxSize)
+size_t File::GetImageDataIntoVector (void *destination, size_t maxSize)
 {
-   if ( ! GetDecompressed() )
+   if ( ! GetRaw() )
    {
       // If the decompression failed nothing can be done.
       return 0;
    }
 
-   if ( HeaderInternal->HasLUT() && PixelConverter->BuildRGBImage() )
+   if ( HeaderInternal->HasLUT() && PixelReadConverter->BuildRGBImage() )
    {
-      if ( PixelConverter->GetRGBSize() > maxSize )
+      if ( PixelReadConverter->GetRGBSize() > maxSize )
       {
          dbg.Verbose(0, "File::GetImageDataIntoVector: pixel data bigger"
                         "than caller's expected MaxSize");
          return 0;
       }
       memcpy( destination,
-              (void*)PixelConverter->GetRGB(),
-              PixelConverter->GetRGBSize() );
-      return PixelConverter->GetRGBSize();
+              (void*)PixelReadConverter->GetRGB(),
+              PixelReadConverter->GetRGBSize() );
+      return PixelReadConverter->GetRGBSize();
    }
 
    // Either no LUT conversion necessary or LUT conversion failed
-   if ( PixelConverter->GetDecompressedSize() > maxSize )
+   if ( PixelReadConverter->GetRawSize() > maxSize )
    {
       dbg.Verbose(0, "File::GetImageDataIntoVector: pixel data bigger"
                      "than caller's expected MaxSize");
       return 0;
    }
    memcpy( destination,
-           (void*)PixelConverter->GetDecompressed(),
-           PixelConverter->GetDecompressedSize() );
-   return PixelConverter->GetDecompressedSize();
+           (void*)PixelReadConverter->GetRaw(),
+           PixelReadConverter->GetRawSize() );
+   return PixelReadConverter->GetRawSize();
 }
 
 /**
- * \brief   Allocates necessary memory, 
- *          Transforms YBR pixels (if any) into RGB pixels
- *          Transforms 3 planes R, G, B  (if any) into a single RGB Plane
- *          Copies the pixel data (image[s]/volume[s]) to newly allocated zone. 
- *          DOES NOT transform Grey plane + 3 Palettes into a RGB Plane
- * @return  Pointer to newly allocated pixel data.
- * \        NULL if alloc fails 
- */
-uint8_t* File::GetImageDataRaw ()
-{
-   return GetDecompressed();
-}
-
-uint8_t* File::GetDecompressed()
-{
-   uint8_t* decompressed = PixelConverter->GetDecompressed();
-   if ( ! decompressed )
-   {
-      // The decompressed image migth not be loaded yet:
-      std::ifstream* fp = HeaderInternal->OpenFile();
-      PixelConverter->ReadAndDecompressPixelData( fp );
-      if(fp) 
-         HeaderInternal->CloseFile();
-
-      decompressed = PixelConverter->GetDecompressed();
-      if ( ! decompressed )
-      {
-         dbg.Verbose(0, "File::GetDecompressed: read/decompress of "
-                        "pixel data apparently went wrong.");
-         return 0;
-      }
-   }
-
-   return decompressed;
-}
-
-/**
- * \brief   Points the internal Pixel_Data pointer to the callers inData
+ * \brief   Points the internal pointer to the callers inData
  *          image representation, BUT WITHOUT COPYING THE DATA.
  *          'image' Pixels are presented as C-like 2D arrays : line per line.
  *          'volume'Pixels are presented as C-like 3D arrays : plane per plane 
@@ -281,13 +291,72 @@ uint8_t* File::GetDecompressed()
  *
  * @return boolean
  */
-bool File::SetImageData(uint8_t* inData, size_t expectedSize)
+void File::SetImageData(uint8_t *inData, size_t expectedSize)
 {
-// FIXME : if already allocated, memory leak !
-   Pixel_Data     = inData;
-   ImageDataSize = expectedSize;
-// FIXME : 7fe0, 0010 IS NOT set ...
-   return true;
+   SetUserData(inData,expectedSize);
+}
+
+/**
+ * \brief   Set the image datas defined by the user
+ * \warning When writting the file, this datas are get as default datas to write
+ */
+void File::SetUserData(uint8_t *data, size_t expectedSize)
+{
+   PixelWriteConverter->SetUserData(data,expectedSize);
+}
+
+/**
+ * \brief   Get the image datas defined by the user
+ * \warning When writting the file, this datas are get as default data to write
+ */
+uint8_t *File::GetUserData()
+{
+   return PixelWriteConverter->GetUserData();
+}
+
+/**
+ * \brief   Get the image data size defined by the user
+ * \warning When writting the file, this datas are get as default data to write
+ */
+size_t File::GetUserDataSize()
+{
+   return PixelWriteConverter->GetUserDataSize();
+}
+
+/**
+ * \brief   Get the image datas from the file.
+ *          If a LUT is found, the data are expanded to be RGB
+ */
+uint8_t *File::GetRGBData()
+{
+   return PixelReadConverter->GetRGB();
+}
+
+/**
+ * \brief   Get the image data size from the file.
+ *          If a LUT is found, the data are expanded to be RGB
+ */
+size_t File::GetRGBDataSize()
+{
+   return PixelReadConverter->GetRGBSize();
+}
+
+/**
+ * \brief   Get the image datas from the file.
+ *          If a LUT is found, the datas are not expanded !
+ */
+uint8_t *File::GetRawData()
+{
+   return PixelReadConverter->GetRaw();
+}
+
+/**
+ * \brief   Get the image data size from the file.
+ *          If a LUT is found, the data are not expanded !
+ */
+size_t File::GetRawDataSize()
+{
+   return PixelReadConverter->GetRawSize();
 }
 
 /**
@@ -299,7 +368,7 @@ bool File::SetImageData(uint8_t* inData, size_t expectedSize)
  * @return false if write fails
  */
 
-bool File::WriteRawData(std::string const & fileName)
+bool File::WriteRawData(std::string const &fileName)
 {
   std::ofstream fp1(fileName.c_str(), std::ios::out | std::ios::binary );
    if (!fp1)
@@ -307,7 +376,14 @@ bool File::WriteRawData(std::string const & fileName)
       dbg.Verbose(2, "Fail to open (write) file:", fileName.c_str());
       return false;
    }
-   fp1.write((char*)Pixel_Data, ImageDataSize);
+
+   if(PixelWriteConverter->GetUserData())
+      fp1.write((char*)PixelWriteConverter->GetUserData(), PixelWriteConverter->GetUserDataSize());
+   else if(PixelReadConverter->GetRGB())
+      fp1.write((char*)PixelReadConverter->GetRGB(), PixelReadConverter->GetRGBSize());
+   else if(PixelReadConverter->GetRaw())
+      fp1.write((char*)PixelReadConverter->GetRaw(), PixelReadConverter->GetRawSize());
+
    fp1.close();
 
    return true;
@@ -322,7 +398,7 @@ bool File::WriteRawData(std::string const & fileName)
  * @return false if write fails
  */
 
-bool File::WriteDcmImplVR (std::string const & fileName)
+bool File::WriteDcmImplVR (std::string const &fileName)
 {
    SetWriteTypeToDcmImplVR();
    return Write(fileName);
@@ -331,12 +407,13 @@ bool File::WriteDcmImplVR (std::string const & fileName)
 /**
 * \brief Writes on disk A SINGLE Dicom file, 
  *        using the Explicit Value Representation convention
- *        NO test is performed on  processor "Endiannity". * @param fileName name of the file to be created
+ *        NO test is performed on  processor "Endiannity". 
+ * @param fileName name of the file to be created
  *                 (any already existing file is overwritten)
  * @return false if write fails
  */
 
-bool File::WriteDcmExplVR (std::string const & fileName)
+bool File::WriteDcmExplVR (std::string const &fileName)
 {
    SetWriteTypeToDcmExplVR();
    return Write(fileName);
@@ -356,65 +433,34 @@ bool File::WriteDcmExplVR (std::string const & fileName)
  * @return false if write fails
  */
 
-bool File::WriteAcr (std::string const & fileName)
+bool File::WriteAcr (std::string const &fileName)
 {
    SetWriteTypeToAcr();
    return Write(fileName);
 }
 
-bool File::Write(std::string const& fileName)
-{
-   switch(WriteType)
-   {
-      case WTYPE_IMPL_VR:
-         return WriteBase(fileName,ImplicitVR);
-      case WTYPE_EXPL_VR:
-         return WriteBase(fileName,ExplicitVR);
-      case WTYPE_ACR:
-         return WriteBase(fileName,ACR);
-   }
-   return(false);
-}
-
 /**
- * \brief Access to the underlying \ref PixelConverter RGBA LUT
- */
-uint8_t* File::GetLutRGBA()
-{
-   return PixelConverter->GetLutRGBA();
-}
-
-//-----------------------------------------------------------------------------
-// Protected
-/**
- * \brief NOT a end user inteded function
- *        (used by WriteDcmExplVR, WriteDcmImplVR, WriteAcr, etc)
+ * \brief Writes on disk A SINGLE Dicom file, 
  * @param fileName name of the file to be created
  *                 (any already existing file is overwritten)
- * @param  type file type (ExplicitVR, ImplicitVR, ...)
  * @return false if write fails
  */
-bool File::WriteBase (std::string const & fileName, FileType type)
+bool File::Write(std::string const &fileName)
 {
-   std::ofstream* fp1 = new std::ofstream(fileName.c_str(), 
-                              std::ios::out | std::ios::binary);
-   if (fp1 == NULL)
-   {
-      dbg.Verbose(2, "Failed to open (write) File: " , fileName.c_str());
-      return false;
-   }
-
-   switch(WriteMode)
+   switch(WriteType)
    {
-      case WMODE_NATIVE :
-         SetWriteToNative();
+      case ImplicitVR:
+         SetWriteFileTypeToImplicitVR();
          break;
-      case WMODE_DECOMPRESSED :
-         SetWriteToDecompressed();
+      case ExplicitVR:
+         SetWriteFileTypeToExplicitVR();
          break;
-      case WMODE_RGB :
-         SetWriteToRGB();
+      case ACR:
+      case ACR_LIBIDO:
+         SetWriteFileTypeToACR();
          break;
+      default:
+         SetWriteFileTypeToExplicitVR();
    }
 
    // --------------------------------------------------------------
@@ -425,37 +471,118 @@ bool File::WriteBase (std::string const & fileName, FileType type)
    // just before writting ...
    /// \todo the best trick would be *change* the recognition code
    ///       but pb expected if user deals with, e.g. COMPLEX images
-/*   if ( HeaderInternal->GetFileType() == ACR_LIBIDO)
+   if( WriteType == ACR_LIBIDO )
    {
       SetWriteToLibido();
-   }*/
+   }
+   else
+   {
+      SetWriteToNoLibido();
+   }
    // ----------------- End of Special Patch ----------------
+  
+   switch(WriteMode)
+   {
+      case WMODE_RAW :
+         SetWriteToRaw();
+         break;
+      case WMODE_RGB :
+         SetWriteToRGB();
+         break;
+   }
 
-   bool check=CheckWriteIntegrity();
+   bool check = CheckWriteIntegrity();
    if(check)
    {
-      HeaderInternal->Write(fp1,type);
+      check = HeaderInternal->Write(fileName,WriteType);
    }
 
+   RestoreWrite();
+   RestoreWriteFileType();
+
    // --------------------------------------------------------------
    // Special Patch to allow gdcm to re-write ACR-LibIDO formated images
    // 
    // ...and we restore the Header to be Dicom Compliant again 
    // just after writting
-/*   if ( HeaderInternal->GetFileType() == ACR_LIBIDO )
-   {
-      RestoreWriteFromLibido();
-   }*/
+   RestoreWriteOfLibido();
    // ----------------- End of Special Patch ----------------
 
-   RestoreWrite();
+   return check;
+}
 
-   fp1->close();
-   delete fp1;
+/**
+ * \brief   Accesses an existing DocEntry (i.e. a Dicom Element)
+ *          through it's (group, element) and modifies it's content with
+ *          the given value.
+ * @param   content new value (string) to substitute with
+ * @param   group     group number of the Dicom Element to modify
+ * @param   element element number of the Dicom Element to modify
+ */
+bool File::SetEntryByNumber(std::string const &content,
+                            uint16_t group, uint16_t element)
+{ 
+   return HeaderInternal->SetEntryByNumber(content,group,element);
+}
 
-   return check;
+
+/**
+ * \brief   Accesses an existing DocEntry (i.e. a Dicom Element)
+ *          through it's (group, element) and modifies it's content with
+ *          the given value.
+ * @param   content new value (void*  -> uint8_t*) to substitute with
+ * @param   lgth new value length
+ * @param   group     group number of the Dicom Element to modify
+ * @param   element element number of the Dicom Element to modify
+ */
+bool File::SetEntryByNumber(uint8_t *content, int lgth,
+                            uint16_t group, uint16_t element)
+{
+   return HeaderInternal->SetEntryByNumber(content,lgth,group,element);
+}
+
+/**
+ * \brief   Modifies the value of a given Doc Entry (Dicom Element)
+ *          when it exists. Create it with the given value when unexistant.
+ * @param   content (string) Value to be set
+ * @param   group   Group number of the Entry 
+ * @param   element  Element number of the Entry
+ * \return  pointer to the modified/created Header Entry (NULL when creation
+ *          failed).
+ */ 
+bool File::ReplaceOrCreateByNumber(std::string const &content,
+                                   uint16_t group, uint16_t element)
+{
+   return HeaderInternal->ReplaceOrCreateByNumber(content,group,element) != NULL;
+}
+
+/*
+ * \brief   Modifies the value of a given Header Entry (Dicom Element)
+ *          when it exists. Create it with the given value when unexistant.
+ *          A copy of the binArea is made to be kept in the Document.
+ * @param   binArea (binary) value to be set
+ * @param   group   Group number of the Entry 
+ * @param   element  Element number of the Entry
+ * \return  pointer to the modified/created Header Entry (NULL when creation
+ *          failed).
+ */
+bool File::ReplaceOrCreateByNumber(uint8_t *binArea, int lgth,
+                                   uint16_t group, uint16_t element)
+{
+   return HeaderInternal->ReplaceOrCreateByNumber(binArea,lgth,group,element) != NULL;
 }
 
+/**
+ * \brief Access to the underlying \ref PixelReadConverter RGBA LUT
+ */
+uint8_t* File::GetLutRGBA()
+{
+   return PixelReadConverter->GetLutRGBA();
+}
+
+//-----------------------------------------------------------------------------
+// Protected
+
 /**
  * \brief Check the write integrity
  *
@@ -466,23 +593,40 @@ bool File::WriteBase (std::string const & fileName, FileType type)
  */
 bool File::CheckWriteIntegrity()
 {
-   if(Pixel_Data)
+   if(PixelWriteConverter->GetUserData())
    {
+      int numberBitsAllocated = HeaderInternal->GetBitsAllocated();
+      if ( numberBitsAllocated == 0 || numberBitsAllocated == 12 )
+      {
+         numberBitsAllocated = 16;
+      }
+
+      size_t decSize = HeaderInternal->GetXSize()
+                    * HeaderInternal->GetYSize() 
+                    * HeaderInternal->GetZSize()
+                    * ( numberBitsAllocated / 8 )
+                    * HeaderInternal->GetSamplesPerPixel();
+      size_t rgbSize = decSize;
+      if( HeaderInternal->HasLUT() )
+         rgbSize = decSize * 3;
+
       switch(WriteMode)
       {
-         case WMODE_NATIVE :
-            break;
-         case WMODE_DECOMPRESSED :
-            if(GetImageDataRawSize()!=ImageDataSize)
+         case WMODE_RAW :
+            if( decSize!=PixelWriteConverter->GetUserDataSize() )
             {
-               std::cerr<<"RAW : "<<GetImageDataRawSize()<<" / "<<ImageDataSize<<std::endl;
+               dbg.Verbose(0, "File::CheckWriteIntegrity: Data size is incorrect (Raw)");
+               //std::cerr << "File::CheckWriteIntegrity: Data size is incorrect (Raw)\n"
+               //          << decSize << " / " << PixelWriteConverter->GetUserDataSize() << "\n";
                return false;
             }
             break;
          case WMODE_RGB :
-            if(GetImageDataSize()!=ImageDataSize)
+            if( rgbSize!=PixelWriteConverter->GetUserDataSize() )
             {
-               std::cerr<<"RGB : "<<GetImageDataSize()<<" / "<<ImageDataSize<<std::endl;
+               dbg.Verbose(0, "File::CheckWriteIntegrity: Data size is incorrect (RGB)");
+               //std::cerr << "File::CheckWriteIntegrity: Data size is incorrect (RGB)\n"
+               //          << decSize << " / " << PixelWriteConverter->GetUserDataSize() << "\n";
                return false;
             }
             break;
@@ -492,20 +636,7 @@ bool File::CheckWriteIntegrity()
    return true;
 }
 
-void File::SetWriteToNative()
-{
-   if(Pixel_Data)
-   {
-      BinEntry* pixel = CopyBinEntry(GetHeader()->GetGrPixel(),GetHeader()->GetNumPixel());
-      pixel->SetValue(GDCM_BINLOADED);
-      pixel->SetBinArea(Pixel_Data,false);
-      pixel->SetLength(ImageDataSize);
-
-      Archive->Push(pixel);
-   }
-}
-
-void File::SetWriteToDecompressed()
+void File::SetWriteToRaw()
 {
    if(HeaderInternal->GetNumberOfScalarComponents()==3 && !HeaderInternal->HasLUT())
    {
@@ -513,30 +644,23 @@ void File::SetWriteToDecompressed()
    } 
    else
    {
-      ValEntryphotInt = CopyValEntry(0x0028,0x0004);
+      ValEntry *photInt = CopyValEntry(0x0028,0x0004);
       if(HeaderInternal->HasLUT())
       {
          photInt->SetValue("PALETTE COLOR ");
-         photInt->SetLength(14);
       }
       else
       {
          photInt->SetValue("MONOCHROME1 ");
-         photInt->SetLength(12);
       }
 
-      BinEntry* pixel = CopyBinEntry(GetHeader()->GetGrPixel(),GetHeader()->GetNumPixel());
+      PixelWriteConverter->SetReadData(PixelReadConverter->GetRaw(),
+                                       PixelReadConverter->GetRawSize());
+
+      BinEntry *pixel = CopyBinEntry(GetHeader()->GetGrPixel(),GetHeader()->GetNumPixel());
       pixel->SetValue(GDCM_BINLOADED);
-      if(Pixel_Data)
-      {
-         pixel->SetBinArea(Pixel_Data,false);
-         pixel->SetLength(ImageDataSize);
-      }
-      else
-      {
-         pixel->SetBinArea(PixelConverter->GetDecompressed(),false);
-         pixel->SetLength(PixelConverter->GetDecompressedSize());
-      }
+      pixel->SetBinArea(PixelWriteConverter->GetData(),false);
+      pixel->SetLength(PixelWriteConverter->GetDataSize());
 
       Archive->Push(photInt);
       Archive->Push(pixel);
@@ -547,38 +671,33 @@ void File::SetWriteToRGB()
 {
    if(HeaderInternal->GetNumberOfScalarComponents()==3)
    {
-      PixelConverter->BuildRGBImage();
+      PixelReadConverter->BuildRGBImage();
       
-      ValEntryspp = CopyValEntry(0x0028,0x0002);
+      ValEntry *spp = CopyValEntry(0x0028,0x0002);
       spp->SetValue("3 ");
-      spp->SetLength(2);
 
-      ValEntryplanConfig = CopyValEntry(0x0028,0x0006);
+      ValEntry *planConfig = CopyValEntry(0x0028,0x0006);
       planConfig->SetValue("0 ");
-      planConfig->SetLength(2);
 
-      ValEntryphotInt = CopyValEntry(0x0028,0x0004);
+      ValEntry *photInt = CopyValEntry(0x0028,0x0004);
       photInt->SetValue("RGB ");
-      photInt->SetLength(4);
 
-      BinEntry* pixel = CopyBinEntry(GetHeader()->GetGrPixel(),GetHeader()->GetNumPixel());
-      pixel->SetValue(GDCM_BINLOADED);
-      if(Pixel_Data)
+      if(PixelReadConverter->GetRGB())
       {
-         pixel->SetBinArea(Pixel_Data,false);
-         pixel->SetLength(ImageDataSize);
+         PixelWriteConverter->SetReadData(PixelReadConverter->GetRGB(),
+                                          PixelReadConverter->GetRGBSize());
       }
-      else if(PixelConverter->GetRGB())
+      else // Raw data
       {
-         pixel->SetBinArea(PixelConverter->GetRGB(),false);
-         pixel->SetLength(PixelConverter->GetRGBSize());
-      }
-      else // Decompressed data
-      {
-         pixel->SetBinArea(PixelConverter->GetDecompressed(),false);
-         pixel->SetLength(PixelConverter->GetDecompressedSize());
+         PixelWriteConverter->SetReadData(PixelReadConverter->GetRaw(),
+                                          PixelReadConverter->GetRawSize());
       }
 
+      BinEntry *pixel = CopyBinEntry(GetHeader()->GetGrPixel(),GetHeader()->GetNumPixel());
+      pixel->SetValue(GDCM_BINLOADED);
+      pixel->SetBinArea(PixelWriteConverter->GetData(),false);
+      pixel->SetLength(PixelWriteConverter->GetDataSize());
+
       Archive->Push(spp);
       Archive->Push(planConfig);
       Archive->Push(photInt);
@@ -597,17 +716,14 @@ void File::SetWriteToRGB()
       // samples per pixels = 1 (in the read file)
       if(HeaderInternal->GetBitsAllocated()==24) 
       {
-         ValEntrybitsAlloc = CopyValEntry(0x0028,0x0100);
+         ValEntry *bitsAlloc = CopyValEntry(0x0028,0x0100);
          bitsAlloc->SetValue("8 ");
-         bitsAlloc->SetLength(2);
 
-         ValEntrybitsStored = CopyValEntry(0x0028,0x0101);
+         ValEntry *bitsStored = CopyValEntry(0x0028,0x0101);
          bitsStored->SetValue("8 ");
-         bitsStored->SetLength(2);
 
-         ValEntryhighBit = CopyValEntry(0x0028,0x0102);
+         ValEntry *highBit = CopyValEntry(0x0028,0x0102);
          highBit->SetValue("7 ");
-         highBit->SetLength(2);
 
          Archive->Push(bitsAlloc);
          Archive->Push(bitsStored);
@@ -616,7 +732,7 @@ void File::SetWriteToRGB()
    }
    else
    {
-      SetWriteToDecompressed();
+      SetWriteToRaw();
    }
 }
 
@@ -641,10 +757,44 @@ void File::RestoreWrite()
    Archive->Restore(0x0028,0x1203);
 }
 
+void File::SetWriteFileTypeToACR()
+{
+   Archive->Push(0x0002,0x0010);
+}
+
+void File::SetWriteFileTypeToExplicitVR()
+{
+   std::string ts = Util::DicomString( 
+      Document::GetTransferSyntaxValue(ExplicitVRLittleEndian).c_str() );
+
+   ValEntry *tss = CopyValEntry(0x0002,0x0010);
+   tss->SetValue(ts);
+
+   Archive->Push(tss);
+}
+
+void File::SetWriteFileTypeToImplicitVR()
+{
+   std::string ts = Util::DicomString(
+      Document::GetTransferSyntaxValue(ImplicitVRLittleEndian).c_str() );
+
+   ValEntry *tss = CopyValEntry(0x0002,0x0010);
+   tss->SetValue(ts);
+
+   Archive->Push(tss);
+}
+
+void File::RestoreWriteFileType()
+{
+   Archive->Restore(0x0002,0x0010);
+}
+
 void File::SetWriteToLibido()
 {
-   ValEntry *oldRow = dynamic_cast<ValEntry *>(HeaderInternal->GetDocEntryByNumber(0x0028, 0x0010));
-   ValEntry *oldCol = dynamic_cast<ValEntry *>(HeaderInternal->GetDocEntryByNumber(0x0028, 0x0011));
+   ValEntry *oldRow = dynamic_cast<ValEntry *>
+                (HeaderInternal->GetDocEntryByNumber(0x0028, 0x0010));
+   ValEntry *oldCol = dynamic_cast<ValEntry *>
+                (HeaderInternal->GetDocEntryByNumber(0x0028, 0x0011));
    
    if( oldRow && oldCol )
    {
@@ -662,18 +812,38 @@ void File::SetWriteToLibido()
       Archive->Push(newRow);
       Archive->Push(newCol);
    }
+
+   ValEntry *libidoCode = CopyValEntry(0x0008,0x0010);
+   libidoCode->SetValue("ACRNEMA_LIBIDO_1.1");
+   Archive->Push(libidoCode);
+}
+
+void File::SetWriteToNoLibido()
+{
+   ValEntry *recCode = dynamic_cast<ValEntry *>
+                (HeaderInternal->GetDocEntryByNumber(0x0008,0x0010));
+   if( recCode )
+   {
+      if( recCode->GetValue() == "ACRNEMA_LIBIDO_1.1" )
+      {
+         ValEntry *libidoCode = CopyValEntry(0x0008,0x0010);
+         libidoCode->SetValue("");
+         Archive->Push(libidoCode);
+      }
+   }
 }
 
-void File::RestoreWriteFromLibido()
+void File::RestoreWriteOfLibido()
 {
    Archive->Restore(0x0028,0x0010);
    Archive->Restore(0x0028,0x0011);
+   Archive->Restore(0x0008,0x0010);
 }
 
-ValEntryFile::CopyValEntry(uint16_t group,uint16_t element)
+ValEntry *File::CopyValEntry(uint16_t group,uint16_t element)
 {
-   DocEntryoldE = HeaderInternal->GetDocEntryByNumber(group, element);
-   ValEntrynewE;
+   DocEntry *oldE = HeaderInternal->GetDocEntryByNumber(group, element);
+   ValEntry *newE;
 
    if(oldE)
    {
@@ -688,10 +858,10 @@ ValEntry* File::CopyValEntry(uint16_t group,uint16_t element)
    return(newE);
 }
 
-BinEntryFile::CopyBinEntry(uint16_t group,uint16_t element)
+BinEntry *File::CopyBinEntry(uint16_t group,uint16_t element)
 {
-   DocEntryoldE = HeaderInternal->GetDocEntryByNumber(group, element);
-   BinEntrynewE;
+   DocEntry *oldE = HeaderInternal->GetDocEntryByNumber(group, element);
+   BinEntry *newE;
 
    if(oldE)
    {
@@ -703,30 +873,55 @@ BinEntry* File::CopyBinEntry(uint16_t group,uint16_t element)
       newE = GetHeader()->NewBinEntryByNumber(group,element);
    }
 
-
    return(newE);
 }
 
 //-----------------------------------------------------------------------------
-// Private
+// Protected
 /**
- * \brief Set the pixel datas in the good entry of the Header
+ * \brief Factorization for various forms of constructors.
  */
-void File::SetPixelData(uint8_t* data)
+void File::Initialise()
 {
-   GetHeader()->SetEntryByNumber( GDCM_BINLOADED,
-      GetHeader()->GetGrPixel(), GetHeader()->GetNumPixel());
+   WriteMode = WMODE_RAW;
+   WriteType = ExplicitVR;
+
+   PixelReadConverter = new PixelReadConvert;
+   PixelWriteConverter = new PixelWriteConvert;
+   Archive = new DocEntryArchive( HeaderInternal );
+
+   if ( HeaderInternal->IsReadable() )
+   {
+      PixelReadConverter->GrabInformationsFromHeader( HeaderInternal );
+   }
+}
 
-   // Will be 7fe0, 0010 in standard case
-   DocEntry* currentEntry = GetHeader()->GetDocEntryByNumber(GetHeader()->GetGrPixel(), GetHeader()->GetNumPixel());
-   if ( currentEntry )
+uint8_t *File::GetRaw()
+{
+   uint8_t *raw = PixelReadConverter->GetRaw();
+   if ( ! raw )
    {
-      if ( BinEntry* binEntry = dynamic_cast<BinEntry *>(currentEntry) )
-         // Flag is to false because datas are kept in the gdcmPixelConvert
-         binEntry->SetBinArea( data, false );
+      // The Raw image migth not be loaded yet:
+      std::ifstream *fp = HeaderInternal->OpenFile();
+      PixelReadConverter->ReadAndDecompressPixelData( fp );
+      if(fp) 
+         HeaderInternal->CloseFile();
+
+      raw = PixelReadConverter->GetRaw();
+      if ( ! raw )
+      {
+         dbg.Verbose(0, "File::GetRaw: read/decompress of "
+                        "pixel data apparently went wrong.");
+         return 0;
+      }
    }
+
+   return raw;
 }
 
+//-----------------------------------------------------------------------------
+// Private
+
 //-----------------------------------------------------------------------------
 } // end namespace gdcm