From 02899c034acb29e15219663f18bfee580fdfb266 Mon Sep 17 00:00:00 2001
From: malaterre <malaterre>
Date: Fri, 10 Dec 2004 15:50:04 +0000
Subject: [PATCH] ENH: Adding an example that only work on linux, not on Win32.
 It shows how to create a DICOM image from scratch, we are only allowed to
 copy the string based entry from a previous DICOM

---
 Example/CMakeLists.txt      |  1 +
 Example/TestFromScratch.cxx | 94 +++++++++++++++++++++++++++++++++++++
 2 files changed, 95 insertions(+)
 create mode 100644 Example/TestFromScratch.cxx

diff --git a/Example/CMakeLists.txt b/Example/CMakeLists.txt
index 328f2f03..3eb21a2f 100644
--- a/Example/CMakeLists.txt
+++ b/Example/CMakeLists.txt
@@ -15,6 +15,7 @@ SET(EXAMPLE_SOURCES
   PrintDicomDir
   TestCopyDicom
   TestChangeHeader
+  TestFromScratch
   TestReadWriteReadCompare
   TestDcm2Acr
   TestWrite
diff --git a/Example/TestFromScratch.cxx b/Example/TestFromScratch.cxx
new file mode 100644
index 00000000..27fd2627
--- /dev/null
+++ b/Example/TestFromScratch.cxx
@@ -0,0 +1,94 @@
+/*=========================================================================
+                                                                                
+  Program:   gdcm
+  Module:    $RCSfile: TestFromScratch.cxx,v $
+  Language:  C++
+  Date:      $Date: 2004/12/10 15:50:04 $
+  Version:   $Revision: 1.1 $
+                                                                                
+  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 "gdcmDictEntry.h"
+#include "gdcmDocEntry.h"
+#include "gdcmValEntry.h"
+#include "gdcmDebug.h"
+
+// The whole purpose of this example is to create a DICOM from scratch
+// well almost the only thing we allow use is to copy the string entry
+// From the previous image read.
+
+int main(int argc, char *argv[])
+{
+   if( argc < 2 )
+   {
+      std::cerr << "Usage: " << argv[0] << " Image.dcm" << std::endl;
+      return 1;
+   }
+
+   // Doesn't seems to do anything:
+   dbg.SetDebug(-1);
+   // Doesn't link properly:
+   //gdcm::Debug::GetReference().SetDebug(1);
+
+   std::string filename = argv[1];
+   //gdcm::File *f1 = new gdcm::File( "/home/malaterre/Creatis/gdcmData/012345.002.050.dcm" );
+   gdcm::File *f1 = new gdcm::File( filename );
+   gdcm::Header *h1 = f1->GetHeader();
+
+   int dataSize = f1->GetImageDataSize();
+   std::cout << "DataSize:      " << dataSize << std::endl;
+   // Since we know the image is 16bits:
+   uint8_t* imageData = f1->GetImageData();
+ 
+   // Hopefully default to something
+   gdcm::Header *h2 = new gdcm::Header();
+
+   const gdcm::TagDocEntryHT &nameHt = h1->GetTagHT();
+   //const gdcm::TagDocEntryHT &nameHt = h2->GetTagHT(); seems to be empty ??
+   for (gdcm::TagDocEntryHT::const_iterator tag = nameHt.begin(); tag != nameHt.end(); ++tag)
+   {
+      //Copy string only:
+      if (tag->second->GetVR().find("SQ") == 0)
+      {
+         // skip DICOM SeQuence, otherwise following cast will crash
+         continue;
+      }
+      const std::string &name = tag->second->GetName();
+      const std::string &value = ((gdcm::ValEntry*)(tag->second))->GetValue();
+      if( name != "unkn"
+       && value.find( "gdcm::NotLoaded" ) != 0
+       && value.find( "gdcm::Binary" )    != 0
+       && value.find( "gdcm::Loaded" )    != 0 )
+      {
+        //std::cout << name << "," << value << std::endl;
+        gdcm::DictEntry *dictEntry = h2->GetPubDict()->GetDictEntryByName(name);
+        h2->ReplaceOrCreateByNumber( value, dictEntry->GetGroup(), dictEntry->GetElement());
+      }
+   }
+
+   h2->ReplaceOrCreateByNumber( imageData, dataSize, 0x7fe0, 0x0010, "PXL" );
+   h2->Print( std::cout );
+
+   gdcm::File *f2 = new gdcm::File( h2 );
+ 
+   // SetImageData does not work because you need to have called
+   // somwhow ReplaceOrCreateByNumber before with the proper group/element
+   f2->SetImageData(imageData, dataSize);
+
+   f2->WriteDcmExplVR( "output.dcm" );
+
+   delete f1;
+   delete f2;
+   delete h2;
+
+   return 0;
+}
-- 
2.45.1