From a716f958bd2a6105ac06997fa9509ebaad6fc285 Mon Sep 17 00:00:00 2001
From: malaterre <malaterre>
Date: Tue, 8 Nov 2005 16:32:20 +0000
Subject: [PATCH] ENH: Saving demo code

---
 Example/WriteDicomAsMPEG.cxx | 108 +++++++++++++++++++++++++++++++++++
 1 file changed, 108 insertions(+)
 create mode 100644 Example/WriteDicomAsMPEG.cxx

diff --git a/Example/WriteDicomAsMPEG.cxx b/Example/WriteDicomAsMPEG.cxx
new file mode 100644
index 00000000..6ee8a2fb
--- /dev/null
+++ b/Example/WriteDicomAsMPEG.cxx
@@ -0,0 +1,108 @@
+/*=========================================================================
+                                                                                
+  Program:   gdcm
+  Module:    $RCSfile: WriteDicomAsMPEG.cxx,v $
+  Language:  C++
+  Date:      $Date: 2005/11/08 16:32:20 $
+  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 "gdcmFile.h"
+#include "gdcmFileHelper.h"
+#include "gdcmUtil.h"
+ 
+
+// Open a dicom file and compress it as JPEG stream
+int main(int argc, char *argv[])
+{
+  const int xsize = 352; //352x240 DirectClass 248kb 0.000u 0:01
+  const int ysize = 240;
+  const int zsize = 10;
+  const int samplesPerPixel = 3;
+  (void)argc;
+  (void)argv;
+
+   // Step 1 : Create the header of the image
+
+   gdcm::File *fileToBuild = gdcm::File::New();
+   std::ostringstream str;
+
+   // Set the image size
+   str.str("");
+   str << xsize;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0011); // Columns
+   str.str("");
+   str << ysize;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0010); // Rows
+
+   if(zsize>1)
+   {
+      str.str("");
+      str << zsize;
+      fileToBuild->InsertEntryString(str.str(),0x0028,0x0008); // Number of Frames
+   }
+
+   // Set the pixel type
+   str.str("");
+   str << 8; //img.componentSize;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0100); // Bits Allocated
+
+   str.str("");
+   str << 8; //img.componentUse;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0101); // Bits Stored
+
+   str.str("");
+   str << 7; //( img.componentSize - 1 );
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0102); // High Bit
+
+   // Set the pixel representation
+   str.str("");
+   str << 0; //img.sign;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0103); // Pixel Representation
+
+   // Set the samples per pixel
+   str.str("");
+   str << samplesPerPixel; //img.components;
+   fileToBuild->InsertEntryString(str.str(),0x0028,0x0002); // Samples per Pixel
+
+// Step 2 : Create the output image
+//   std::cout << "2...";
+//   if( img.componentSize%8 > 0 )
+//   {
+//      img.componentSize += 8-img.componentSize%8;
+//   }
+   size_t size = xsize * ysize * zsize
+               * samplesPerPixel /* * img.componentSize / 8*/;
+
+   uint8_t *imageData = new uint8_t[size];
+   gdcm::FileHelper *fileH = gdcm::FileHelper::New(fileToBuild);
+   //fileH->SetImageData(imageData,size);
+   std::ifstream mpeg("/tmp/ts.mpg");
+   mpeg.seekg(0, std::ios::end);
+   size = mpeg.tellg();
+   std::cerr << "Size MPEG:" << size << std::endl;
+   mpeg.seekg(0, std::ios::beg);
+   mpeg.read((char*)imageData, size);
+   fileH->SetImageData(imageData, size);
+   fileH->SetWriteTypeToJPEG(  );
+   std::string outfilename = "/tmp/ts.dcm";
+   if( !fileH->Write(outfilename) )
+     {
+     std::cerr << "Badddd" << std::endl;
+     }
+
+   fileToBuild->Delete();
+   fileH->Delete();
+
+   return 0;
+}
+
-- 
2.48.1