Avoid confusion due to missing values

[gdcm.git] / Example / RawToDicom.cxx

diff --git a/Example/RawToDicom.cxx b/Example/RawToDicom.cxx
index d6c40f731fa9fc51d680543c3432e1026d34bba9..62710fc7e0f3670476fee4e4b02249f2a30906e2 100755 (executable)
--- a/Example/RawToDicom.cxx
+++ b/Example/RawToDicom.cxx
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: RawToDicom.cxx,v $
   Language:  C++
-  Date:      $Date: 2006/07/17 13:27:57 $
-  Version:   $Revision: 1.6 $
+  Date:      $Date: 2007/05/23 14:18:04 $
+  Version:   $Revision: 1.10 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -13,7 +13,7 @@
      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.
-                                                                                
+                                      
 =========================================================================*/
 
 /**
@@ -30,6 +30,39 @@
 #include <iostream>
 #include <sstream>
 
+void ConvertSwapZone(int pixelSize, void *Raw, size_t RawSize);
+
+void ConvertSwapZone(int pixelSize, void *Raw, size_t RawSize)
+{
+   unsigned int i;    
+   if ( pixelSize == 2 )
+   {
+      uint16_t *im16 = (uint16_t*)Raw;
+      for( i = 0; i < RawSize / 2; i++ )
+      {
+         im16[i]= (im16[i] >> 8) | (im16[i] << 8 );
+      }     
+   }
+   else if ( pixelSize == 4 )
+   {
+      uint32_t s32;
+      uint16_t high;
+      uint16_t low;
+      uint32_t *im32 = (uint32_t*)Raw;
+
+      for( i = 0; i < RawSize / 4; i++ )
+      {
+         low     = im32[i] & 0x0000ffff; // 3412
+         high    = im32[i] >> 16;
+         s32     = low;
+         im32[i] = ( s32 << 16 ) | high;
+      }
+      
+   }
+}
+
+
+
 int main(int argc, char *argv[])
 {
    START_USAGE(usage)
@@ -40,6 +73,7 @@ int main(int argc, char *argv[])
    "                   rows=nb of Rows                                        ",
    "                   lines=nb of Lines,                                     ",
    "                   pixeltype={8U|8S|16U|16S|32U|32S}                      ",
+   "                   [{b|l}] b:BigEndian,l:LittleEndian default : l         ",
    "                   [frames = nb of Frames] //defaulted to 1               ",
    "                   [samples = {1|3}}       //defaulted to 1(1:Gray,3:RGB) ",
    "                   [patientname = Patient's name]                         ",
@@ -50,7 +84,7 @@ int main(int argc, char *argv[])
    
 
    // Initialize Arguments Manager   
-   gdcm::ArgMgr *am= new gdcm::ArgMgr(argc, argv);
+   GDCM_NAME_SPACE::ArgMgr *am= new GDCM_NAME_SPACE::ArgMgr(argc, argv);
   
    if (argc == 1 || am->ArgMgrDefined("usage") )
    {
@@ -69,11 +103,13 @@ int main(int argc, char *argv[])
    int nZ = am->ArgMgrGetInt("frames", 1);
    int samplesPerPixel = am->ArgMgrGetInt("samples", 1);
    
-   
+   int b = am->ArgMgrDefined("b");
+   int l = am->ArgMgrDefined("l");
+      
    char *pixelType = am->ArgMgrWantString("pixeltype", usage);
    
    if (am->ArgMgrDefined("debug"))
-      gdcm::Debug::DebugOn();
+      GDCM_NAME_SPACE::Debug::DebugOn();
 
    /* if unused Param we give up */
    if ( am->ArgMgrPrintUnusedLabels() )
@@ -98,6 +134,8 @@ int main(int argc, char *argv[])
       return 0;
    }  
 
+   bool bigEndian = GDCM_NAME_SPACE::Util::IsCurrentProcessorBigEndian();
+
    std::string strPixelType(pixelType);
    int pixelSign;
    int pixelSize;
@@ -105,32 +143,32 @@ int main(int argc, char *argv[])
    if (strPixelType == "8S")
    {
       pixelSize = 1;
-      pixelSign = 0;
+      pixelSign = 1;
    }
    else if (strPixelType == "8U")
    {
       pixelSize = 1;
-      pixelSign = 1;
+      pixelSign = 0;
    }
    else if (strPixelType == "16S")
    {
       pixelSize = 2;
-      pixelSign = 0; 
+      pixelSign = 1; 
    }   
    else if (strPixelType == "16U")
    {
       pixelSize = 2;
-      pixelSign = 1;
+      pixelSign = 0;
    }      
    else if (strPixelType == "32S")
    {
       pixelSize = 4;
-      pixelSign = 0;
+      pixelSign = 1;
    }   
    else if (strPixelType == "32U")
    {
       pixelSize = 4;
-      pixelSign = 1;
+      pixelSign = 0;
    }
    else
    {
@@ -142,14 +180,19 @@ int main(int argc, char *argv[])
    uint8_t *pixels = new uint8_t[dataSize];
    
    Fp->read((char*)pixels, (size_t)dataSize);
-  
-
+     
+   if ( pixelSize !=1 && ( (l && bigEndian) || (b && ! bigEndian) ) )
+   {  
+      ConvertSwapZone(pixelSize, pixels, dataSize);   
+   }
+   
+   
 // Create an empty FileHelper
 
-   gdcm::FileHelper *fileH = gdcm::FileHelper::New();
+   GDCM_NAME_SPACE::FileHelper *fileH = GDCM_NAME_SPACE::FileHelper::New();
  
  // Get the (empty) image header.  
-   gdcm::File *fileToBuild = fileH->GetFile();
+   GDCM_NAME_SPACE::File *fileToBuild = fileH->GetFile();
      
    
    // If you want to use this program as a template to create
@@ -161,7 +204,7 @@ int main(int argc, char *argv[])
    //          keeping the same 'Study Instance UID' for various images
    // The user may add images to a 'Manufacturer Study',
    //          adding new Series to an already existing Study
-   std::string studyUID =  gdcm::Util::CreateUniqueUID(); 
+   std::string studyUID =  GDCM_NAME_SPACE::Util::CreateUniqueUID(); 
    fileToBuild->InsertEntryString(studyUID, 0x0020,0x000d,"UI");
 
    // 'Serie Instance UID'
@@ -169,7 +212,7 @@ int main(int argc, char *argv[])
    // keeping the same 'Serie Instance UID' for various images
    // The user shouldn't add any image to a 'Manufacturer Serie'
    // but there is no way no to prevent him for doing that
-   std::string serieUID =  gdcm::Util::CreateUniqueUID();    
+   std::string serieUID =  GDCM_NAME_SPACE::Util::CreateUniqueUID();    
    fileToBuild->InsertEntryString(serieUID, 0x0020,0x000e,"UI");   
  
    // end of 'keep out of loop lines