]> Creatis software - creaBruker.git/commitdiff
Use shorter names for 'Study description'
authorjean-pierre roux <jean-pierre.roux@creatis.insa-lyon.fr>
Mon, 14 Jun 2010 08:48:45 +0000 (08:48 +0000)
committerjean-pierre roux <jean-pierre.roux@creatis.insa-lyon.fr>
Mon, 14 Jun 2010 08:48:45 +0000 (08:48 +0000)
lib/src1/bruker2dicom.cxx

index bf947e7bee044ac51b33d73586e2d226edfe51bb..b1e939267edf00570b4358100149a86862cd1f40 100644 (file)
@@ -221,7 +221,10 @@ void Bruker2Dicom::DealWithSingleStudyDirectory (GDCM_NAME_SPACE::DirListType &f
       subject_study_name = "defStudyName";         // Why not? 
    }
     // subject_name is already in 'Patient Name' 
-    strStudyDescr = /*subject_name + "." + */ subject_study_name + "." + subject_entry + "." + subject_position + "." + subject_date;
+    if (0) // We don't realy need so much information for 'Study'
+       strStudyDescr = /*subject_name + "." + */ subject_study_name + "." + subject_entry + "." + subject_position + "." + subject_date;
+    else
+       strStudyDescr = /*subject_name + "." + */ subject_study_name /*+ "." + subject_entry + "." + subject_position + "." + subject_date*/;
 
     // creation directory : 'nom de la Study'
 
@@ -1196,7 +1199,7 @@ void Bruker2Dicom::dealWithCarto(GDCM_NAME_SPACE::DirListType &fileNames, int NX
 int Bruker2Dicom::CheckUserDirectory(std::string &userDirName)
 {
    GDCM_NAME_SPACE::DirList dirList(userDirName, false, true); // DON'T get recursively the list of files
-     
+
    GDCM_NAME_SPACE::DirListType fileNames;
    fileNames = dirList.GetFilenames();
 
@@ -1219,7 +1222,7 @@ int Bruker2Dicom::CheckUserDirectory(std::string &userDirName)
        acqpFound= true; // user passes a 'serie', not a 'study' 
     else
        acqpFound = false; // user passed a 'non study' directory; Hope it's a 'set of studies' directory!
-   
+
    int type;    
    if  (subjectFound )  type = 1; // user passed a 'study
    else if  (acqpFound) type = 2; // user passed a 'serie' 
@@ -1281,7 +1284,7 @@ void Bruker2Dicom::cleanString(std::string &s)
       l--;
       s = s.substr(0, l);
    }   
-   
+
    if (s[0] == '<')      
       s= s.substr(1,l-2);
    std::string repChar("_");   
@@ -1364,7 +1367,7 @@ std::vector<BrukerImage> Bruker2Dicom::CreateImageSet ( )
          br_acqp.SetLoopStructure();
          std::vector<int> tempVect                      = br_acqp.GetLoopStructure() ;
          std::map<std::string, BrukerFieldData> map     = br_acqp.GetBrukerHeaderMap();
-        
+
          bool result                                    = br_acqp.ObjectVaryingProperties.init(map,tempVect);
 
          if (result == false)
@@ -1387,7 +1390,6 @@ std::vector<BrukerImage> Bruker2Dicom::CreateImageSet ( )
 
  // Just for checking
  /*
-
          std::vector<std::vector <double> > imageOrientation;
          std::vector <double> imagePosition; 
          for(int i=0;i<brukerImageList.size();i++)
@@ -1438,7 +1440,7 @@ void Bruker2Dicom::MakeDicomImage(unsigned char *tabPixels,
 
    GDCM_NAME_SPACE::File *file;
    file = GDCM_NAME_SPACE::File::New();       
-      
+
   // Set the image size
    str.str(""); 
    str << X;
@@ -1468,7 +1470,7 @@ void Bruker2Dicom::MakeDicomImage(unsigned char *tabPixels,
 
   // Set the pixel representation // 0/1 , 0=unsigned
    file->InsertEntryString("1",0x0028,0x0103, "US"); // Pixel Representation
-   
+
   // Set the samples per pixel // 1:Grey level, 3:RGB
    file->InsertEntryString("1",0x0028,0x0002, "US"); // Samples per Pixel
 
@@ -1476,25 +1478,25 @@ void Bruker2Dicom::MakeDicomImage(unsigned char *tabPixels,
    str.str("");
    str << spacingX << "\\" << spacingY;
    file->InsertEntryString(str.str(),0x0028,0x0030, "DS"); // Pixel Spacing     
-   
+
  //   0018 0050 DS 1 Slice Thickness 
    str.str("");    
    str << sliceDistance;
    file->InsertEntryString(str.str(),0x0018,0x0050, "DS"); 
-   
+
 //    0020 0011 IS 1 Series Number
    str.str("");    
    str << serieNumber;
    file->InsertEntryString(str.str(),0x0020,0x0011, "IS");
-      
+
 //    0020|0013 [IS]  [Instance Number] 
    instanceNumber++;
    str.str("");    
    str << instanceNumber;
    file->InsertEntryString(str.str(),0x0020,0x0013, "IS");
-     
+
   // 1.2.840.10008.5.1.4.1.1.4.1 : Enhanced MR Image Storage
- //  file->InsertEntryString("1.2.840.10008.5.1.4.1.1.4.1" , 0x0002, 0x0002, "UI");  // [Media Storage SOP Class UID]
 //  file->InsertEntryString("1.2.840.10008.5.1.4.1.1.4.1" , 0x0002, 0x0002, "UI");  // [Media Storage SOP Class UID]
   // file->InsertEntryString("1.2.840.10008.5.1.4.1.1.4.1" , 0x0008, 0x0016, "UI");  // [SOP Class UID]
 
 // OK : MR is NOT multiframe, but I want just a quick an dirty solution
@@ -1509,7 +1511,7 @@ void Bruker2Dicom::MakeDicomImage(unsigned char *tabPixels,
 
    file->InsertEntryString(studyUID, 0x0020, 0x000d, "UI");
    file->InsertEntryString(serieUID, 0x0020, 0x000e, "UI");
-   
+
 //  0008 0020 DA 1 Study Date
 //  0008 0030 TM 1 Study Time
 
@@ -1542,7 +1544,7 @@ std::cout << "charImageOrientation  " <<
                               imageSet[imgNum].getRotationMatrixRPS2XYZ()[1][0],
                               imageSet[imgNum].getRotationMatrixRPS2XYZ()[1][1],
                               imageSet[imgNum].getRotationMatrixRPS2XYZ()[1][2] ) ;
-    
+
    file->InsertEntryString(charImageOrientation,0x0020,0x0037, "DS");
 
 // 0020 0032 DS 3 Image Position (Patient) 
@@ -1552,13 +1554,12 @@ std::cout << "charImageOrientation  " <<
                              imageSet[imgNum].getTranslationVectorRPS2XYZ()[0], 
                              imageSet[imgNum].getTranslationVectorRPS2XYZ()[1],
                              imageSet[imgNum].getTranslationVectorRPS2XYZ()[2]);
-  
+
    file->InsertEntryString(charImagePosition,0x0020,0x0032, "DS");  //0020 0032 DS 3 Image Position (Patient) 
-   
-   
+
    // Private creator
    /*
-   
+
 0029|0010 [LO]   [Private Creator] [CREATIS HEADER ]
 0029|0011 [LO]   [Private Creator] [CREATIS for BRUKER to DICOM]
 0029|1031 [LO]   [] [4.0.7571167 ]
@@ -1592,46 +1593,46 @@ std::cout << "charImageOrientation  " <<
 > 138 1600
 >
    */
-   
+
    file->InsertEntryString("CREATIS HEADER", 0x0029,0x0010, "LO");
    file->InsertEntryString("CREATIS FOR BRUKER-TO-DICOM", 0x0029,0x0011, "LO");
-   
+
    BrukerFieldData brf_method            = br_method.GetFieldData("Method");
    const std::vector<std::string> method = brf_method.GetStringValue();
 
   // std::string method = br_method.GetFieldData("Method").GetStringValue();    
    file->InsertEntryString(method[0], 0x0029,0x0100, "LO");
-   
+
    std::cout << " method[" << method[0]  << "]" << std::endl;
    if (method[0] == "DtiEpi" || method[0] == "DtiStandard")
    {
      // const std::vector<int> NDiffDir       = br_method.GetFieldData("PVM_DwNDiffDir").GetIntValue();
       BrukerFieldData brf_NDiffDir    = br_method.GetFieldData("PVM_DwNDiffDir");
       const std::vector<int> NDiffDir = brf_NDiffDir.GetIntValue();
-      
+
       //const std::vector<int> NDiffExpEach   = br_method.GetFieldData("PVM_DwNDiffExpEach").GetIntValue();
       BrukerFieldData brf_NDiffExpEach      = br_method.GetFieldData("PVM_DwNDiffExpEach");
       const std::vector<int> NDiffExpEach   =brf_NDiffExpEach.GetIntValue();
-      
+
 //      const std::vector<int> AoImages     = br_method.GetFieldData("PVM_DwAoImages").GetIntValue();
       BrukerFieldData brf_AoImages        = br_method.GetFieldData("PVM_DwAoImages");
       const std::vector<int> AoImages     = brf_AoImages.GetIntValue();
-      
+
       BrukerFieldData brf_DwBvalEach       = br_method.GetFieldData("PVM_DwBvalEach");
       const std::vector<double> DwBvalEach = brf_DwBvalEach.GetDoubleValue();
-      
+
  std::cout << "nb directions : " <<  NDiffDir[0] << " nb valeurs de b par direction : " << NDiffExpEach[0] 
            << " nb images for b=0 : " << AoImages[0] << std::endl;
-          
+
        str.str(""); 
        str << NDiffDir[0];
-       file->InsertEntryString(str.str(),0x0029,0x0101,"US"); // directions
+       file->InsertEntryString(str.str(),0x0029,0x0101,"US"); // nb de directions
        str.str(""); 
        str << NDiffExpEach[0];
-       file->InsertEntryString(str.str(),0x0029,0x0102,"US"); // valeurs de b par direction       
+       file->InsertEntryString(str.str(),0x0029,0x0102,"US"); // nb valeurs de b par direction       
        str.str(""); 
        str << AoImages[0];
-       file->InsertEntryString(str.str(),0x0029,0x0103,"US"); // image pour b=0
+       file->InsertEntryString(str.str(),0x0029,0x0103,"US"); // nb images pour b=0
 
        str.str("");
        for (unsigned int i=0; i<NDiffExpEach[0]; i++) {
@@ -1660,7 +1661,7 @@ std::cout << "charImageOrientation  " <<
    GDCM_NAME_SPACE::FileHelper *fileH;
    fileH = GDCM_NAME_SPACE::FileHelper::New(file);
    fileH->SetContentType(contentType);   
-    
+
    // cast is just to avoid warnings (*no* conversion is performed)
    //fileH->SetImageData((uint8_t *)img,int(maxX*maxY)*sizeof(uint16_t)); // troubles when maxX, mayY are *actually* float!