]> Creatis software - gdcm.git/blobdiff - src/gdcmDocument.cxx
Add the NO_SHADOWSEQ load mode
[gdcm.git] / src / gdcmDocument.cxx
index 0147abbe90c9d4b0434b6cfa21f4f9350043bbba..ce90658f9349684f408f35519382dfe928e24ae3 100644 (file)
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: gdcmDocument.cxx,v $
   Language:  C++
-  Date:      $Date: 2005/03/22 11:29:37 $
-  Version:   $Revision: 1.233 $
+  Date:      $Date: 2005/06/29 15:58:33 $
+  Version:   $Revision: 1.252 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -35,8 +35,6 @@
 namespace gdcm 
 {
 //-----------------------------------------------------------------------------
-// Refer to Document::CheckSwap()
-//const unsigned int Document::HEADER_LENGTH_TO_READ = 256;
 
 // Refer to Document::SetMaxSizeLoadEntry()
 const unsigned int Document::MAX_SIZE_LOAD_ELEMENT_VALUE = 0xfff; // 4096
@@ -47,8 +45,7 @@ const unsigned int Document::MAX_SIZE_LOAD_ELEMENT_VALUE = 0xfff; // 4096
 
 /**
  * \brief This default constructor neither loads nor parses the file. 
- *        You should then invoke \ref Document::SetFileName 
- *        and \ref Document::Load.
+ *        You should then invoke \ref Document::Load.
  *         
  */
 Document::Document() 
@@ -60,15 +57,17 @@ Document::Document()
    Initialize();
    SwapCode = 1234;
    Filetype = ExplicitVR;
+   // Load will set it to true if sucessfull
    Group0002Parsed = false;
-   LoadMode = 0x00000000; // Load everything
+   IsDocumentAlreadyLoaded = false;
+   LoadMode = 0x00000000; // default : load everything, later
 }
 
 /**
  * \brief   Constructor (not to break the API) 
- * @param   filename 'Document' (File or DicomDir) to be opened for parsing
+ * @param   fileName 'Document' (File or DicomDir) to be open for parsing
  */
-Document::Document( std::string const &filename )
+Document::Document( std::string const &fileName )
          :ElementSet(-1) 
 {
    Fp = 0;
@@ -80,7 +79,10 @@ Document::Document( std::string const &filename )
    Group0002Parsed = false;
    LoadMode = 0x00000000; // Load everything
 
-   Load(filename); 
+   // Load will set it to true if sucessfull
+   IsDocumentAlreadyLoaded = false;
+
+   Load(fileName);
 }
 /**
  * \brief   Canonical destructor.
@@ -96,18 +98,34 @@ Document::~Document ()
 
 /**
  * \brief   Loader  
- * @param   filename 'Document' (File or DicomDir) to be opened for parsing
+ * @param   fileName 'Document' (File or DicomDir) to be open for parsing
+ * @return false if file cannot be open or no swap info was found,
+ *         or no tag was found.
  */
-void Document::Load( std::string const &filename ) 
+bool Document::Load( std::string const &fileName ) 
 {
-   Filename = filename;
-
    // We should clean out anything that already exists.
+   // Check IsDocumentAlreadyLoaded to be sure.
+   if ( IsDocumentAlreadyLoaded )
+   {
+      gdcmWarningMacro( "A file was already parsed inside this "
+                        << "gdcm::Document (previous name was: "
+                        << Filename.c_str() << ". New name is :"
+                        << fileName );
+     // todo : clean out the 'Document'
+     // Should we call ClearEntry() on the parent object ?!?
+   }
+
+   Filename = fileName;
 
    Fp = 0;
    if ( !OpenFile() )
    {
-      return;
+      // warning already performed in OpenFile()
+      //gdcmWarningMacro( "Unable to open as an ACR/DICOM file: "
+      //                 << Filename.c_str() );
+      Filetype = Unknown;
+      return false;
    }
 
    Group0002Parsed = false;
@@ -119,22 +137,44 @@ void Document::Load( std::string const &filename )
 
    Fp->seekg(0, std::ios::beg);
 
-   CheckSwap();
+   // CheckSwap returns a boolean 
+   // (false if no swap info of any kind was found)
+   if (! CheckSwap() )
+   {
+      gdcmWarningMacro( "Neither a DICOM V3 nor an ACR-NEMA file: " 
+                   << Filename.c_str());
+      CloseFile(); 
+      return false;      
+    }
+
    long beg = Fp->tellg();      // just after DICOM preamble (if any)
 
    lgt -= beg;                  // remaining length to parse    
 
    ParseDES( this, beg, lgt, false); // Loading is done during parsing
 
+   if ( IsEmpty() )
+   { 
+      gdcmWarningMacro( "No tag in internal hash table for: "
+                        << Filename.c_str());
+      CloseFile(); 
+      return false;
+   }
+   IsDocumentAlreadyLoaded = true;
+
    Fp->seekg( 0, std::ios::beg);
    
    // Load 'non string' values
       
    std::string PhotometricInterpretation = GetEntryValue(0x0028,0x0004);   
-   if( PhotometricInterpretation == "PALETTE COLOR " )
+   if ( PhotometricInterpretation == "PALETTE COLOR " )
    {
-      LoadEntryBinArea(0x0028,0x1200);  // gray LUT   
-      /// FIXME FIXME FIXME
+   // FIXME
+   // Probabely this line should be outside the 'if'
+   // Try to find an image sample holding a 'gray LUT'
+      LoadEntryBinArea(0x0028,0x1200);  // gray LUT
+   
+      /// FIXME
       /// The tags refered by the three following lines used to be CORRECTLY
       /// defined as having an US Value Representation in the public
       /// dictionary. BUT the semantics implied by the three following
@@ -159,9 +199,25 @@ void Document::Load( std::string const &filename )
       LoadEntryBinArea(0x0028,0x1222);
       // Segmented Blue  Palette Color LUT Data
       LoadEntryBinArea(0x0028,0x1223);
-   } 
+   }
    //FIXME later : how to use it?
-   LoadEntryBinArea(0x0028,0x3006);  //LUT Data (CTX dependent) 
+   SeqEntry *modLutSeq = GetSeqEntry(0x0028,0x3000);
+   if ( modLutSeq !=0 )
+   {
+      SQItem *sqi= modLutSeq->GetFirstSQItem();
+      if ( sqi != 0 )
+      {
+         BinEntry *b = sqi->GetBinEntry(0x0028,0x3006);
+         if ( b != 0 )
+         {
+            if ( b->GetLength() != 0 )
+            {
+               LoadEntryBinArea(b);    //LUT Data (CTX dependent)
+            }   
+        }
+     }      
+   }
 
    CloseFile(); 
   
@@ -188,6 +244,8 @@ void Document::Load( std::string const &filename )
          SetValEntry(rows   , 0x0028, 0x0011);
    }
    // --- End of ACR-LibIDO kludge --- 
+
+   return true;
 }
 
 /**
@@ -227,16 +285,15 @@ bool Document::SetShaDict(DictKey const &dictName)
 }
 
 /**
- * \brief  This predicate, based on hopefully reasonable heuristics,
- *         decides whether or not the current Document was properly parsed
- *         and contains the mandatory information for being considered as
- *         a well formed and usable Dicom/Acr File.
- * @return true when Document is the one of a reasonable Dicom/Acr file,
- *         false otherwise. 
+ * \brief  This predicate tells us whether or not the current Document 
+ *         was properly parsed and contains at least *one* Dicom Element
+ *         (and nothing more, sorry).
+ * @return false when we're 150 % sure it's NOT a Dicom/Acr file,
+ *         true otherwise. 
  */
 bool Document::IsReadable()
 {
-   if( Filetype == Unknown)
+   if ( Filetype == Unknown )
    {
       gdcmWarningMacro( "Wrong filetype");
       return false;
@@ -258,7 +315,7 @@ bool Document::IsReadable()
 bool Document::IsDicomV3()
 {
    // Checking if Transfer Syntax exists is enough
-   // Anyway, it's to late check if the 'Preamble' was found ...
+   // Anyway, it's too late check if the 'Preamble' was found ...
    // And ... would it be a rich idea to check ?
    // (some 'no Preamble' DICOM images exist !)
    return GetDocEntry(0x0002, 0x0010) != NULL;
@@ -419,60 +476,67 @@ std::ifstream *Document::OpenFile()
       return 0;
    }
 
-   if(Fp)
+   if ( Fp )
    {
       gdcmWarningMacro( "File already open: " << Filename.c_str());
       CloseFile();
    }
 
    Fp = new std::ifstream(Filename.c_str(), std::ios::in | std::ios::binary);
-   if( ! *Fp )
+   if ( ! *Fp )
    {
-      gdcmDebugMacro( "Cannot open file: " << Filename.c_str());
+   // Don't user gdcmErrorMacro :
+   // a spurious message will appear when you use, for instance 
+   // gdcm::FileHelper *fh = new gdcm::FileHelper( outputFileName );
+   // to create outputFileName.
+      gdcmWarningMacro( "Cannot open file: " << Filename.c_str());
       delete Fp;
       Fp = 0;
       return 0;
+      //exit(1); // No function is allowed to leave the application instead
+                 // of warning the caller
    }
  
    uint16_t zero = 0;
    Fp->read((char*)&zero, (size_t)2);
-   if( Fp->eof() )
+   if ( Fp->eof() )
    {
       CloseFile();
       return 0;
    }
  
-   //ACR -- or DICOM with no Preamble; may start with a Shadow Group --
-   if( 
+   //-- ACR or DICOM with no Preamble; may start with a Shadow Group --
+   if 
        zero == 0x0001 || zero == 0x0100 || zero == 0x0002 || zero == 0x0200 ||
        zero == 0x0003 || zero == 0x0300 || zero == 0x0004 || zero == 0x0400 ||
        zero == 0x0005 || zero == 0x0500 || zero == 0x0006 || zero == 0x0600 ||
        zero == 0x0007 || zero == 0x0700 || zero == 0x0008 || zero == 0x0800 )
    {
-      std::string msg 
-         = Util::Format("ACR/DICOM with no preamble: (%04x)\n", zero);
+      std::string msg = Util::Format(
+        "ACR/DICOM starting at the begining of the file:(%04x)\n", zero);
       gdcmWarningMacro( msg.c_str() );
       return Fp;
    }
  
-   //DICOM
+   //-- DICOM --
    Fp->seekg(126L, std::ios::cur);
    char dicm[4] = {' ',' ',' ',' '};
    Fp->read(dicm,  (size_t)4);
-   if( Fp->eof() )
+   if ( Fp->eof() )
    {
       CloseFile();
       return 0;
    }
-   if( memcmp(dicm, "DICM", 4) == 0 )
+   if ( memcmp(dicm, "DICM", 4) == 0 )
    {
       HasDCMPreamble = true;
       return Fp;
    }
+
+   // -- Neither ACR/No Preamble Dicom nor DICOMV3 file
    CloseFile();
-   gdcmWarningMacro( "Not DICOM/ACR (missing preamble)" << Filename.c_str());
+   gdcmWarningMacro( "Neither ACR/No Preamble Dicom nor DICOMV3 file: "
+                      << Filename.c_str()); 
    return 0;
 }
 
@@ -482,7 +546,7 @@ std::ifstream *Document::OpenFile()
  */
 bool Document::CloseFile()
 {
-   if( Fp )
+   if ( Fp )
    {
       Fp->close();
       delete Fp;
@@ -546,7 +610,7 @@ void Document::LoadEntryBinArea(uint16_t group, uint16_t elem)
       return;
 
    BinEntry *binElement = dynamic_cast<BinEntry *>(docElement);
-   if( !binElement )
+   if ( !binElement )
       return;
 
    LoadEntryBinArea(binElement);
@@ -559,11 +623,11 @@ void Document::LoadEntryBinArea(uint16_t group, uint16_t elem)
  */
 void Document::LoadEntryBinArea(BinEntry *elem) 
 {
-   if(elem->GetBinArea())
+   if (elem->GetBinArea() )
       return;
 
    bool openFile = !Fp;
-   if(openFile)
+   if ( openFile )
       OpenFile();
 
    size_t o =(size_t)elem->GetOffset();
@@ -571,14 +635,14 @@ void Document::LoadEntryBinArea(BinEntry *elem)
 
    size_t l = elem->GetLength();
    uint8_t *a = new uint8_t[l];
-   if( !a )
+   if ( !a )
    {
       gdcmWarningMacro( "Cannot allocate BinEntry content");
       return;
    }
 
    Fp->read((char*)a, l);
-   if( Fp->fail() || Fp->eof())
+   if ( Fp->fail() || Fp->eof() )
    {
       delete[] a;
       return;
@@ -586,7 +650,7 @@ void Document::LoadEntryBinArea(BinEntry *elem)
 
    elem->SetBinArea(a);
 
-   if(openFile)
+   if ( openFile )
       CloseFile();
 }
 
@@ -598,7 +662,7 @@ void Document::LoadEntryBinArea(BinEntry *elem)
  */
 void Document::LoadDocEntrySafe(DocEntry *entry)
 {
-   if(Fp)
+   if ( Fp )
    {
       long PositionOnEntry = Fp->tellg();
       LoadDocEntry(entry);
@@ -618,11 +682,11 @@ bool Document::operator<(Document &document)
    // Patient Name
    std::string s1 = GetEntryValue(0x0010,0x0010);
    std::string s2 = document.GetEntryValue(0x0010,0x0010);
-   if(s1 < s2)
+   if (s1 < s2)
    {
       return true;
    }
-   else if( s1 > s2 )
+   else if ( s1 > s2 )
    {
       return false;
    }
@@ -648,7 +712,7 @@ bool Document::operator<(Document &document)
          {
             return true;
          }
-         else if( s1 > s2 )
+         else if ( s1 > s2 )
          {
             return false;
          }
@@ -661,7 +725,7 @@ bool Document::operator<(Document &document)
             {
                return true;
             }
-            else if( s1 > s2 )
+            else if ( s1 > s2 )
             {
                return false;
             }
@@ -687,7 +751,7 @@ uint16_t Document::ReadInt16()
    {
       throw FormatError( "Document::ReadInt16()", " file error." );
    }
-   if( Fp->eof() )
+   if ( Fp->eof() )
    {
       throw FormatError( "Document::ReadInt16()", "EOF." );
    }
@@ -709,7 +773,7 @@ uint32_t Document::ReadInt32()
    {
       throw FormatError( "Document::ReadInt32()", " file error." );
    }
-   if( Fp->eof() )
+   if ( Fp->eof() )
    {
       throw FormatError( "Document::ReadInt32()", "EOF." );
    }
@@ -746,19 +810,19 @@ int Document::ComputeGroup0002Length( FileType filetype )
    {
       gr = entry->GetGroup();
 
-      if( gr == 0x0002 )
+      if ( gr == 0x0002 )
       {
          found0002 = true;
 
-         if( entry->GetElement() != 0x0000 )
+         if ( entry->GetElement() != 0x0000 )
          {
             vr = entry->GetVR();
  
-            if( filetype == ExplicitVR )
+            if ( filetype == ExplicitVR )
             {
-               if ( (vr == "OB") || (vr == "OW") || (vr == "SQ") ) 
+               if ( (vr == "OB") || (vr == "OW") || (vr == "SQ") || (vr == "UT") 
                {
-                  // explicit VR AND OB, OW, SQ : 4 more bytes
+                  // explicit VR AND OB, OW, SQ, UT : 4 more bytes
                   groupLength +=  4;
                }
             }
@@ -854,7 +918,7 @@ void Document::ParseDES(DocEntrySet *set, long offset,
             }
 
             LoadDocEntry( newBinEntry );
-            if( !set->AddEntry( newBinEntry ) )
+            if ( !set->AddEntry( newBinEntry ) )
             {
               //Expect big troubles if here
               //delete newBinEntry;
@@ -896,7 +960,7 @@ void Document::ParseDES(DocEntrySet *set, long offset,
                      if ( strLgrGroup != GDCM_UNFOUND)
                      {
                         lgrGroup = atoi(strLgrGroup.c_str());
-                        Fp->seekg(lgrGroup , std::ios::cur);
+                        Fp->seekg(lgrGroup, std::ios::cur);
                         used = false;
                         continue;
                      }
@@ -904,7 +968,7 @@ void Document::ParseDES(DocEntrySet *set, long offset,
                }
              }
 
-            if( !set->AddEntry( newValEntry ) )
+            if ( !set->AddEntry( newValEntry ) )
             {
               // If here expect big troubles
               // delete newValEntry; //otherwise mem leak
@@ -913,13 +977,13 @@ void Document::ParseDES(DocEntrySet *set, long offset,
 
             if (delimitor)
             {
-               if(!used)
+               if ( !used )
                   delete newDocEntry;
                break;
             }
-            if ( !delim_mode && ((long)(Fp->tellg())-offset) >= l_max)
+            if ( !delim_mode && ((long)(Fp->tellg())-offset) >= l_max )
             {
-               if(!used)
+               if ( !used )
                   delete newDocEntry;
                break;
             }
@@ -930,7 +994,8 @@ void Document::ParseDES(DocEntrySet *set, long offset,
       }
       else
       {
-         // VR = "SQ"
+         /////////////////////// SeqEntry :  VR = "SQ"
+
          unsigned long l = newDocEntry->GetReadLength();          
          if ( l != 0 ) // don't mess the delim_mode for zero-length sequence
          {
@@ -943,9 +1008,21 @@ void Document::ParseDES(DocEntrySet *set, long offset,
               delim_mode = false;
             }
          }
-
-         if ( (LoadMode & NO_SEQ) && ! delim_mode ) // User asked to skip SQ
+        if ( (LoadMode & NO_SHADOWSEQ) && ! delim_mode )
+        { 
+           // User asked to skip SQ only if they belong to Shadow group
+           if ( newDocEntry->GetGroup()%2 != 0 )
+           {
+              Fp->seekg( l, std::ios::cur);
+              used = false;
+              continue;  
+           } 
+        } 
+         if ( (LoadMode & NO_SEQ) && ! delim_mode ) 
          {
+           // User asked to skip *any* SeQuence
             Fp->seekg( l, std::ios::cur);
             used = false;
             continue;
@@ -979,20 +1056,20 @@ void Document::ParseDES(DocEntrySet *set, long offset,
                      newDocEntry->GetOffset(),
                      l, delim_mode);
          }
-         if( !set->AddEntry( newSeqEntry ) )
+         if ( !set->AddEntry( newSeqEntry ) )
          {
             used = false;
          }
 
          if ( !delim_mode && ((long)(Fp->tellg())-offset) >= l_max)
          {
-            if( !used )
+            if ( !used )
                delete newDocEntry;
             break;
          }
       }
 
-      if( !used )
+      if ( !used )
          delete newDocEntry;
    }
 }
@@ -1018,7 +1095,7 @@ void Document::ParseSQ( SeqEntry *seqEntry,
          // FIXME Should warn user
          break;
       }
-      if( delim_mode )
+      if ( delim_mode )
       {
          if ( newDocEntry->IsSequenceDelimitor() )
          {
@@ -1085,7 +1162,7 @@ void Document::LoadDocEntry(DocEntry *entry)
    //          (fffe e000) tells us an Element is beginning
    //          (fffe e00d) tells us an Element just ended
    //          (fffe e0dd) tells us the current SeQuence just ended
-   if( group == 0xfffe )
+   if ( group == 0xfffe )
    {
       // NO more value field for SQ !
       return;
@@ -1107,17 +1184,16 @@ void Document::LoadDocEntry(DocEntry *entry)
    {
       if (BinEntry *binEntryPtr = dynamic_cast< BinEntry* >(entry) )
       {  
-         //s << "gdcm::NotLoaded (BinEntry)";
          s << GDCM_NOTLOADED;
-         s << " Address:" << (long)entry->GetOffset();
-         s << " Length:"  << entry->GetLength();
+         s << " Ad.:" << (long)entry->GetOffset();
+         s << " x(" << std::hex << entry->GetOffset() << ")";
+         s << std::dec;
+         s << " Lgt:"  << entry->GetLength();
          s << " x(" << std::hex << entry->GetLength() << ")";
          binEntryPtr->SetValue(s.str());
       }
-      // Be carefull : a BinEntry IS_A ValEntry ... 
       else if (ValEntry *valEntryPtr = dynamic_cast< ValEntry* >(entry) )
       {
-        // s << "gdcm::NotLoaded. (ValEntry)";
          s << GDCM_NOTLOADED;  
          s << " Address:" << (long)entry->GetOffset();
          s << " Length:"  << entry->GetLength();
@@ -1151,7 +1227,7 @@ void Document::LoadDocEntry(DocEntry *entry)
       uint32_t NewInt;
       int nbInt;
       // When short integer(s) are expected, read and convert the following 
-      // n *two characters properly i.e. consider them as short integers as
+      // (n * 2) characters properly i.e. consider them as short integers as
       // opposed to strings.
       // Elements with Value Multiplicity > 1
       // contain a set of integers (not a single one)       
@@ -1200,11 +1276,11 @@ void Document::LoadDocEntry(DocEntry *entry)
    str[length] = '\0'; //this is only useful when length is odd
    // Special DicomString call to properly handle \0 and even length
    std::string newValue;
-   if( length % 2 )
+   if ( length % 2 )
    {
       newValue = Util::DicomString(str, length+1);
       gdcmWarningMacro("Warning: bad length: " << length <<
-                       ",For string :" <<  newValue.c_str()); 
+                       " For string :" <<  newValue.c_str()); 
       // Since we change the length of string update it length
       //entry->SetReadLength(length+1);
    }
@@ -1218,12 +1294,17 @@ void Document::LoadDocEntry(DocEntry *entry)
    {
       if ( Fp->fail() || Fp->eof())
       {
-         gdcmWarningMacro("Unread element value");
+         if ( Fp->fail() )
+            gdcmWarningMacro("--> fail");
+
+         gdcmWarningMacro("Unread element value " << valEntry->GetKey() 
+                          << " lgt : " << valEntry->GetReadLength() 
+                          << " at " << std::hex << valEntry->GetOffset());
          valEntry->SetValue(GDCM_UNREAD);
          return;
       }
 
-      if( vr == "UI" )
+      if ( vr == "UI" )
       {
          // Because of correspondance with the VR dic
          valEntry->SetValue(newValue);
@@ -1235,7 +1316,9 @@ void Document::LoadDocEntry(DocEntry *entry)
    }
    else
    {
-      gdcmErrorMacro( "Should have a ValEntry, here !");
+      gdcmWarningMacro("Should have a ValEntry, here ! " << valEntry->GetKey() 
+                          << " lgt : " << valEntry->GetReadLength() 
+                          << " at " << std::hex << valEntry->GetOffset());
    }
 }
 
@@ -1251,7 +1334,7 @@ void Document::FindDocEntryLength( DocEntry *entry )
    
    if ( Filetype == ExplicitVR && !entry->IsImplicitVR() ) 
    {
-      if ( vr == "OB" || vr == "OW" || vr == "SQ" || vr == "UN" ) 
+      if ( vr == "OB" || vr == "OW" || vr == "SQ" || vr == "UT" || vr == "UN" ) 
       {
          // The following reserved two bytes (see PS 3.5-2003, section
          // "7.1.2 Data element structure with explicit vr", p 27) must be
@@ -1401,7 +1484,7 @@ std::string Document::FindDocEntryVR()
    Fp->read (vr, (size_t)2);
    vr[2] = 0;
 
-   if( !CheckDocEntryVR(vr) )
+   if ( !CheckDocEntryVR(vr) )
    {
       Fp->seekg(positionOnEntry, std::ios::beg);
       return GDCM_UNKNOWN;
@@ -1449,14 +1532,14 @@ std::string Document::GetDocEntryValue(DocEntry *entry)
       // Elements with Value Multiplicity > 1
       // contain a set of short integers (not a single one) 
    
-      if( vr == "US" || vr == "SS" )
+      if ( vr == "US" || vr == "SS" )
       {
          uint16_t newInt16;
 
          nbInt = length / 2;
          for (int i=0; i < nbInt; i++) 
          {
-            if( i != 0 )
+            if ( i != 0 )
             {
                s << '\\';
             }
@@ -1472,14 +1555,14 @@ std::string Document::GetDocEntryValue(DocEntry *entry)
       // as usual for standard multivaluated filels
       // Elements with Value Multiplicity > 1
       // contain a set of integers (not a single one) 
-      else if( vr == "UL" || vr == "SL" )
+      else if ( vr == "UL" || vr == "SL" )
       {
          uint32_t newInt32;
 
          nbInt = length / 4;
          for (int i=0; i < nbInt; i++) 
          {
-            if( i != 0)
+            if ( i != 0)
             {
                s << '\\';
             }
@@ -1606,7 +1689,7 @@ void Document::FixDocEntryFoundLength(DocEntry *entry,
    // Allthough not recent many such GE corrupted images are still present
    // on Creatis hard disks. Hence this fix shall remain when such images
    // are no longer in use (we are talking a few years, here)...
-   // Note: XMedCom probably uses such a trick since it is able to read
+   // Note: XMedCon probably uses such a trick since it is able to read
    //       those pesky GE images ...
    if ( foundLength == 13)
    {
@@ -1636,12 +1719,12 @@ void Document::FixDocEntryFoundLength(DocEntry *entry,
    //////// We encountered a 'delimiter' element i.e. a tag of the form 
    // "fffe|xxxx" which is just a marker. Delimiters length should not be
    // taken into account.
-   else if( gr == 0xfffe )
+   else if ( gr == 0xfffe )
    {    
      // According to the norm, fffe|0000 shouldn't exist. BUT the Philips
      // image gdcmData/gdcm-MR-PHILIPS-16-Multi-Seq.dcm happens to
      // causes extra troubles...
-     if( entry->GetElement() != 0x0000 )
+     if ( entry->GetElement() != 0x0000 )
      {
         foundLength = 0;
      }
@@ -1704,12 +1787,7 @@ bool Document::IsDocEntryAnInteger(DocEntry *entry)
  *         true  when we hope ours assuptions are OK
  */
 bool Document::CheckSwap()
-{
-   // The only guaranted way of finding the swap code is to find a
-   // group tag since we know it's length has to be of four bytes i.e.
-   // 0x00000004. Finding the swap code in then straigthforward. Trouble
-   // occurs when we can't find such group...
-   
+{   
    uint32_t  s32;
    uint16_t  s16;
        
@@ -1724,7 +1802,7 @@ bool Document::CheckSwap()
    Fp->read(deb, 256);
    
    char *entCur = deb + 128;
-   if( memcmp(entCur, "DICM", (size_t)4) == 0 )
+   if ( memcmp(entCur, "DICM", (size_t)4) == 0 )
    {
       gdcmWarningMacro( "Looks like DICOM Version3 (preamble + DCM)" );
       
@@ -1748,11 +1826,11 @@ bool Document::CheckSwap()
       // even if elem 0002,0010 (Transfer Syntax) tells us the file is
       // *Implicit* VR  (see former 'gdcmData/icone.dcm')
       
-      if( memcmp(entCur, "UL", (size_t)2) == 0 ||
-          memcmp(entCur, "OB", (size_t)2) == 0 ||
-          memcmp(entCur, "UI", (size_t)2) == 0 ||
-          memcmp(entCur, "CS", (size_t)2) == 0 )  // CS, to remove later
-                                                  // when Write DCM *adds*
+      if ( memcmp(entCur, "UL", (size_t)2) == 0 ||
+           memcmp(entCur, "OB", (size_t)2) == 0 ||
+           memcmp(entCur, "UI", (size_t)2) == 0 ||
+           memcmp(entCur, "CS", (size_t)2) == 0 )  // CS, to remove later
+                                                   // when Write DCM *adds*
       // FIXME
       // Use Document::dicom_vr to test all the possibilities
       // instead of just checking for UL, OB and UI !? group 0000 
@@ -1780,19 +1858,56 @@ bool Document::CheckSwap()
       
       // Position the file position indicator at first tag 
       // (i.e. after the file preamble and the "DICM" string).
-      Fp->seekg(0, std::ios::beg);
+
+      Fp->seekg(0, std::ios::beg); // FIXME : Is it usefull?
+
       Fp->seekg ( 132L, std::ios::beg);
       return true;
-   } // End of DicomV3
+   } // ------------------------------- 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).
-   gdcmWarningMacro( "Not a DICOM Version3 file");
+
+   gdcmWarningMacro( "Not a Kosher DICOM Version3 file (no preamble)");
+
    Fp->seekg(0, std::ios::beg);
 
+   // Let's check 'No Preamble Dicom File' :
+   // Should start with group 0x0002
+   // and be Explicit Value Representation
+
+   s16 = *((uint16_t *)(deb));
+   SwapCode = 0;     
+   switch ( s16 )
+   {
+      case 0x0002 :
+         SwapCode = 1234;
+         entCur = deb + 4;
+         break;
+      case 0x0200 :
+         SwapCode = 4321;
+         entCur = deb + 6;
+    } 
+
+   if ( SwapCode != 0 )
+   {
+      if ( memcmp(entCur, "UL", (size_t)2) == 0 ||
+           memcmp(entCur, "OB", (size_t)2) == 0 ||
+           memcmp(entCur, "UI", (size_t)2) == 0 ||
+           memcmp(entCur, "SH", (size_t)2) == 0 ||
+           memcmp(entCur, "AE", (size_t)2) == 0 ||
+           memcmp(entCur, "OB", (size_t)2) == 0 )
+         {
+            Filetype = ExplicitVR;
+            gdcmWarningMacro( "Group 0002 : Explicit Value Representation");
+            return true;
+          }
+    }
+// ------------------------------- End of 'No Preamble' DicomV3 -------------
+
    // Our next best chance would be to be considering a 'clean' ACR/NEMA file.
-   // By clean we mean that the length of the first tag is written down.
+   // By clean we mean that the length of the first group is written down.
    // If this is the case and since the length of the first group HAS to be
    // four (bytes), then determining the proper swap code is straightforward.
 
@@ -1801,7 +1916,6 @@ bool Document::CheckSwap()
    // representation of a 32 bits integer. Hence the following dirty
    // trick :
    s32 = *((uint32_t *)(entCur));
-
    switch( s32 )
    {
       case 0x00040000 :
@@ -1875,7 +1989,8 @@ bool Document::CheckSwap()
  */
 void Document::SwitchByteSwapCode() 
 {
-   gdcmWarningMacro( "Switching Byte Swap code from "<< SwapCode);
+   gdcmWarningMacro( "Switching Byte Swap code from "<< SwapCode
+                     << " at :" <<std::hex << Fp->tellg() );
    if ( SwapCode == 1234 ) 
    {
       SwapCode = 4321;
@@ -1946,37 +2061,38 @@ DocEntry *Document::ReadNextDocEntry()
    std::string vr = FindDocEntryVR();
    std::string realVR = vr;
 
-   if( vr == GDCM_UNKNOWN)
+   if ( vr == GDCM_UNKNOWN )
    {
       if ( elem == 0x0000 ) // Group Length
          realVR = "UL";     // must be UL
       else
       {
          DictEntry *dictEntry = GetDictEntry(group,elem);
-         if( dictEntry )
+         if ( dictEntry )
             realVR = dictEntry->GetVR();
       }
    }
 
    DocEntry *newEntry;
-   if( Global::GetVR()->IsVROfSequence(realVR) )
+   if ( Global::GetVR()->IsVROfSequence(realVR) )
       newEntry = NewSeqEntry(group, elem);
-   else if( Global::GetVR()->IsVROfStringRepresentable(realVR) )
+   else if ( Global::GetVR()->IsVROfStringRepresentable(realVR) )
       newEntry = NewValEntry(group, elem,vr);
    else
       newEntry = NewBinEntry(group, elem,vr);
 
-   if( vr == GDCM_UNKNOWN )
+   if ( vr == GDCM_UNKNOWN )
    {
-      if( Filetype == ExplicitVR )
+      if ( Filetype == ExplicitVR )
       {
          // We thought this was explicit VR, but we end up with an
          // implicit VR tag. Let's backtrack.
          if ( newEntry->GetGroup() != 0xfffe )
          { 
             std::string msg;
-            msg = Util::Format("Entry (%04x,%04x) should be Explicit VR\n", 
-                          newEntry->GetGroup(), newEntry->GetElement());
+            int offset = Fp->tellg();
+            msg = Util::Format("Entry (%04x,%04x) at %x should be Explicit VR\n", 
+                          newEntry->GetGroup(), newEntry->GetElement(), offset );
             gdcmWarningMacro( msg.c_str() );
           }
       }
@@ -2001,7 +2117,7 @@ DocEntry *Document::ReadNextDocEntry()
 
 /**
  * \brief   Handle broken private tag from Philips NTSCAN
- *          where the endianess is being switch to BigEndian for no
+ *          where the endianess is being switched to BigEndian for no
  *          apparent reason
  * @return  no return
  */