X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FTestCopyDicom.cxx;h=24400816bab017ec246e08936efd8cd927a762b4;hb=b8ca604b2c1c17790f19ac61b63dfb5b550ac6c0;hp=f7139f88615de71ca356d631ef773183736a5f2a;hpb=5dd999e6394b31ea3b040e8cd71563dc1951c23c;p=gdcm.git

diff --git a/Example/TestCopyDicom.cxx b/Example/TestCopyDicom.cxx
index f7139f88..24400816 100644
--- a/Example/TestCopyDicom.cxx
+++ b/Example/TestCopyDicom.cxx
@@ -1,10 +1,30 @@
+/*=========================================================================
+                                                                                
+  Program:   gdcm
+  Module:    $RCSfile: TestCopyDicom.cxx,v $
+  Language:  C++
+  Date:      $Date: 2004/11/17 10:20:06 $
+  Version:   $Revision: 1.13 $
+                                                                                
+  Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
+  l'Image). All rights reserved. See Doc/License.txt or
+  http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
+                                                                                
+     This software is distributed WITHOUT ANY WARRANTY; without even
+     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+     PURPOSE.  See the above copyright notices for more information.
+                                                                                
+=========================================================================*/
 #include "gdcmHeader.h"
 #include "gdcmFile.h"
 #include "gdcmDocument.h"
 #include "gdcmValEntry.h"
+#include "gdcmBinEntry.h"
 
 #ifndef _WIN32
-#include <unistd.h>
+#include <unistd.h> //for access, unlink
+#else
+#include <io.h> //for _access
 #endif
 
 // return true if the file exists
@@ -46,74 +66,87 @@ int main(int argc, char* argv[])
       return 1;
    }
 
-   if( FileExists( argv[2] ) )
-   {
-      std::cerr << "Don't try to cheat, I am removing the file anyway" << std::endl;
-      if( !RemoveFile( argv[2] ) )
+// don't modify identation in order to let this source xdiffable with ../Test
+
+      std::string filename = argv[1];
+      std::string output = argv[2];
+
+      if( FileExists( output.c_str() ) )
       {
-         std::cerr << "Ouch, the file exist, but I cannot remove it" << std::endl;
-         return 1;
+         std::cerr << "Don't try to cheat, I am removing the file anyway" << std::endl;
+         if( !RemoveFile( output.c_str() ) )
+         {
+            std::cerr << "Ouch, the file exist, but I cannot remove it" << std::endl;
+            return 1;
+         }
       }
-   }
-   gdcmFile *original = new gdcmFile( argv[1] );
+      gdcm::File *original = new gdcm::File( filename );
    
-   std::cout << "--- Original ----------------------" << std::endl;
-   //original->GetHeader()->Print();
+      std::cout << "--- Original ----------------------" << std::endl;
+      //original->GetHeader()->Print();
    
-   gdcmFile *copy = new gdcmFile( argv[2] );
+      gdcm::File *copy = new gdcm::File( output );
+
+      const gdcm::TagDocEntryHT & Ht = original->GetHeader()->GetTagHT();
 
-   //First of all copy the header field by field
+      size_t dataSize = original->GetImageDataSize();
+      uint8_t* imageData = original->GetImageData();
   
-   // Warning :Accessor gdcmElementSet::GetEntry() should not exist 
-   //It was commented out by Mathieu, that was a *good* idea
-   // (the user does NOT have to know the way we implemented the Header !)
-    
-   TagDocEntryHT & Ht = original->GetHeader()->GetEntry();   
-   
-   for (TagDocEntryHT::iterator tag = Ht.begin(); tag != Ht.end(); ++tag)
-   {
-      if (tag->second->GetVR() == "SQ") //to skip pb of SQ recursive exploration
-         continue;
-
-      std::string value = ((gdcmValEntry*)(tag->second))->GetValue();
-      //According to JPR I should also skip those:
-//      if (tag->second->GetVR() == "unkn") //to skip pb of SQ recursive exploration
-//         continue;
-
-//      if( value.find( "gdcm::NotLoaded" ) == 0 )
-//         continue;
-
-      if( value.find( "gdcm::Loaded" ) == 0 )
-         continue;
-         
-      //std::cerr << "Reading: " << tag->second->GetVR() << std::endl;
-      //tag->second->Print(); std::cout << std::endl;
-    
-      std::cerr << "Reading: " << value  << std::endl;
-
-      // Well ... Should have dynamic cast here 
-      copy->GetHeader()->ReplaceOrCreateByNumber( 
-                                 value,
-                                 tag->second->GetGroup(), 
-                                 tag->second->GetElement() );
-   
-     // todo : Setting Offset to 0 to avoid further missprint 
-   }
+      //First of all copy the header field by field
+  
+      // Warning :Accessor gdcmElementSet::GetEntry() should not exist 
+      // It was commented out by Mathieu, that was a *good* idea
+      // (the user does NOT have to know the way we implemented the Header !)
+      // Waiting for a 'clean' solution, I keep the method ...JPRx
+
+      gdcm::DocEntry* d;
+
+      for (gdcm::TagDocEntryHT::const_iterator tag = Ht.begin(); tag != Ht.end(); ++tag)
+      {
+         d = tag->second;
+         d->Print(); std::cout << std::endl;
+         if ( gdcm::BinEntry* b = dynamic_cast<gdcm::BinEntry*>(d) )
+         {              
+            copy->GetHeader()->ReplaceOrCreateByNumber( 
+                                 b->GetBinArea(),
+                                 b->GetLength(),
+                                 b->GetGroup(), 
+                                 b->GetElement(),
+                                 b->GetVR() );
+         }
+         else if ( gdcm::ValEntry* v = dynamic_cast<gdcm::ValEntry*>(d) )
+         {   
+            copy->GetHeader()->ReplaceOrCreateByNumber( 
+                                 v->GetValue(),
+                                 v->GetGroup(), 
+                                 v->GetElement(),
+                                 v->GetVR() ); 
+         }
+         else
+         {
+          // We skip pb of SQ recursive exploration
+          //std::cout << "Skipped Sequence " 
+          //          << "------------- " << d->GetVR() << " "<< std::hex
+          //          << d->GetGroup() << " " << d->GetElement()
+          //  << std::endl;    
+         }
+      }
+
+
+
 
-   size_t dataSize = original->GetImageDataSize();
-   void *imageData = original->GetImageData();
 
-   //copy->GetImageData();
-   copy->SetImageData(imageData, dataSize);
+      //copy->GetImageData();
+      //copy->SetImageData(imageData, dataSize);
 
-   std::cout << "--- Copy ----------------------" << std::endl;
-   std::cout <<std::endl << "DO NOT care about Offset"  <<std::endl<<std::endl;; 
-   //copy->GetHeader()->Print();
-   std::cout << "--- ---- ----------------------" << std::endl;
+      std::cout << "--- Copy ----------------------" << std::endl;
+      std::cout <<std::endl << "DO NOT care about Offset"  <<std::endl<<std::endl;; 
+      copy->GetHeader()->Print();
+      std::cout << "--- ---- ----------------------" << std::endl;
    
-   copy->WriteDcmExplVR( argv[2] );
+      copy->WriteDcmExplVR( output );
 
-   return 0;
+      return 0;
 }