]> Creatis software - gdcm.git/blobdiff - Example/RawToDicom.cxx
Normalization
[gdcm.git] / Example / RawToDicom.cxx
index e9fa0f66c386c78f0034dacb5244c5f6d52e9222..da5f00b2f964465a4dcb95c2051a412f59033cef 100755 (executable)
@@ -3,8 +3,8 @@
   Program:   gdcm
   Module:    $RCSfile: RawToDicom.cxx,v $
   Language:  C++
-  Date:      $Date: 2007/03/23 15:01:47 $
-  Version:   $Revision: 1.9 $
+  Date:      $Date: 2007/09/20 11:15:54 $
+  Version:   $Revision: 1.14 $
                                                                                 
   Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
   l'Image). All rights reserved. See Doc/License.txt or
@@ -61,8 +61,6 @@ void ConvertSwapZone(int pixelSize, void *Raw, size_t RawSize)
    }
 }
 
-
-
 int main(int argc, char *argv[])
 {
    START_USAGE(usage)
@@ -72,19 +70,23 @@ int main(int argc, char *argv[])
    "                   fileout=outputFileName                                 ",
    "                   rows=nb of Rows                                        ",
    "                   lines=nb of Lines,                                     ",
+   "                   [frames = nb of Frames] //defaulted to 1               ",
    "                   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]                         ",
+   "                   [samples = {1|3}}       //(1:Gray,3:RGB) defaulted to 1",
+   "                   [monochrome1]                                          ",
+   "                   [studyid = ] [patientname = Patient's name]            ",
    "                   [debug]                                                ",
    "                                                                          ",
-   "      debug      : developper wants to run the program in 'debug mode'    ",
+   "  monochrome1 = user wants MONOCHROME1 photom. interp. (0=white)          ", 
+   "  studyUID   : *aware* user wants to add the serie                        ",
+   "                                             to an already existing study ",
+   "  debug      : developper wants to run the program in 'debug mode'        ",
    FINISH_USAGE
    
 
-   // Initialize Arguments Manager   
-   gdcm::ArgMgr *am= new gdcm::ArgMgr(argc, argv);
+   // ------------ Initialize Arguments Manager ----------------  
+   GDCM_NAME_SPACE::ArgMgr *am= new GDCM_NAME_SPACE::ArgMgr(argc, argv);
   
    if (argc == 1 || am->ArgMgrDefined("usage") )
    {
@@ -107,9 +109,23 @@ int main(int argc, char *argv[])
    int l = am->ArgMgrDefined("l");
       
    char *pixelType = am->ArgMgrWantString("pixeltype", usage);
-   
+
+   bool monochrome1 = ( 0 != am->ArgMgrDefined("monochrome1") );
+      
    if (am->ArgMgrDefined("debug"))
-      gdcm::Debug::DebugOn();
+      GDCM_NAME_SPACE::Debug::DebugOn();
+   bool userDefinedStudy = ( 0 != am->ArgMgrDefined("studyUID") );
+
+   const char *studyUID;
+   if (userDefinedStudy)
+      studyUID  = am->ArgMgrGetString("studyUID");  
+
+   // not described *on purpose* in the Usage !
+   bool userDefinedSerie = ( 0 != am->ArgMgrDefined("serieUID") );       
+
+   const char *serieUID;
+   if(userDefinedSerie)
+      serieUID = am->ArgMgrGetString("serieUID");
 
    /* if unused Param we give up */
    if ( am->ArgMgrPrintUnusedLabels() )
@@ -123,6 +139,7 @@ int main(int argc, char *argv[])
 
    // ----------- End Arguments Manager ---------
    
+ /// \TODO Deal with all the images of a directory
   
  // Read the Raw file  
    std::ifstream *Fp = new std::ifstream(inputFileName, std::ios::in | std::ios::binary);
@@ -134,7 +151,7 @@ int main(int argc, char *argv[])
       return 0;
    }  
 
-   bool bigEndian = gdcm::Util::IsCurrentProcessorBigEndian();
+   bool bigEndian = GDCM_NAME_SPACE::Util::IsCurrentProcessorBigEndian();
 
    std::string strPixelType(pixelType);
    int pixelSign;
@@ -175,7 +192,21 @@ int main(int argc, char *argv[])
       std::cout << "Wrong 'pixeltype' (" << strPixelType << ")" << std::endl;
       return 1;
    }
+   std::string strStudyUID;
+   std::string strSerieUID;
+
+   if (userDefinedStudy)
+      strSerieUID =  studyUID;
+   else
+      strStudyUID =  GDCM_NAME_SPACE::Util::CreateUniqueUID();
    
+   if (userDefinedStudy)
+     strSerieUID =  serieUID;
+   else
+      strStudyUID =  GDCM_NAME_SPACE::Util::CreateUniqueUID();  
+      
+        
    int dataSize =  nX*nY*nZ*pixelSize*samplesPerPixel;
    uint8_t *pixels = new uint8_t[dataSize];
    
@@ -184,38 +215,30 @@ int main(int argc, char *argv[])
    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();
-     
-   
-   // If you want to use this program as a template to create
-   // a 'Single Study UID - Single Serie UID' file set
-   // keep the following lines out of the loop
+   GDCM_NAME_SPACE::File *fileToBuild = fileH->GetFile();
 
    // 'Study Instance UID'
    // The user is allowed to create his own Study, 
    //          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(); 
-   fileToBuild->InsertEntryString(studyUID, 0x0020,0x000d,"UI");
+
+   fileToBuild->InsertEntryString(strStudyUID,0x0020,0x000d,"UI");  //  Study UID   
 
    // 'Serie Instance UID'
    // The user is allowed to create his own Series, 
    // 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();    
-   fileToBuild->InsertEntryString(serieUID, 0x0020,0x000e,"UI");   
-   // end of 'keep out of loop lines  
+   
+   fileToBuild->InsertEntryString(strSerieUID,0x0020,0x000e,"UI");  //  Serie UID
    
    std::ostringstream str;
 
@@ -252,8 +275,9 @@ int main(int argc, char *argv[])
    str.str("");
    str << samplesPerPixel;
    
-// If you deal with a Serie of images, it up to you to tell gdcm,
-// for each image, what are the values of :
+// If you deal with a Serie of images, as slices of a volume,
+// it up to you to tell gdcm, for each image, what are the values of :
+// 
 // 0020 0032 DS 3 Image Position (Patient)
 // 0020 0037 DS 6 Image Orientation (Patient)
 
@@ -261,8 +285,11 @@ int main(int argc, char *argv[])
 
    if (strlen(patientName) != 0)
       fileToBuild->InsertEntryString(patientName,0x0010,0x0010, "PN"); // Patient's Name
-   
-   
+    
+   //  0=white  
+   if(monochrome1)
+      fileH->SetPhotometricInterpretationToMonochrome1();
+     
 // Set the image Pixel Data
    fileH->SetImageData(pixels,dataSize);