X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=src%2FgdcmElValSet.cxx;h=cfdae3cfe9b9fb5b4a29fc59fa79ba7e0955cd46;hb=d5e17ffd70201ed72ebb59a2d0c09c8eac177c52;hp=dc2093e28b01c671d058a73985a97eac8587f9eb;hpb=551b84c1dc82294c4fc7356da0b5b1b3d74f7645;p=gdcm.git

diff --git a/src/gdcmElValSet.cxx b/src/gdcmElValSet.cxx
index dc2093e2..cfdae3cf 100644
--- a/src/gdcmElValSet.cxx
+++ b/src/gdcmElValSet.cxx
@@ -1,7 +1,12 @@
 // gdcmElValSet.cxx
 
-#include "gdcm.h"
 #include "gdcmUtil.h"
+#include "gdcmElValSet.h"
+
+
+#include <vector>
+static void Tokenize (const string& str, vector<string>& tokens, const string& delimiters = " ");
+
 
 TagElValueHT & ElValSet::GetTagHt(void) {
 	return tagHt;
@@ -142,8 +147,13 @@ int ElValSet::Write(FILE * _fp) {
 	guint32 lgr;
 	const char * val;
 	string vr;
-	guint32 val_int32;
-	guint16 val_int16;
+	guint32 val_uint32;
+	gint32  val_int32;
+	guint16 val_uint16;
+	gint16  val_int16;
+	
+	vector<string> tokens;
+	
 	void *ptr;
 	char str_lgrCalcGroupe[10];
 	
@@ -166,6 +176,8 @@ int ElValSet::Write(FILE * _fp) {
 	//
 	// cf : code IdDcmWriteFile dans libido/src/dcmwrite.c
 		
+	if (0)  // Risque de pb dans le calcul des lgr de chaque groupe. On le saute pour le moment!
+	
 	for (TagElValueHT::iterator tag = tagHt.begin();
 		  tag != tagHt.end();
 		  ++tag){
@@ -189,7 +201,7 @@ int ElValSet::Write(FILE * _fp) {
 			}	
 			
 			if (deja) {
-				sprintf(str_lgrCalcGroupe,"%d",lgrCalcGroupe);
+				//sprintf(str_lgrCalcGroupe,"%d",lgrCalcGroupe);
 				elemZPrec->SetValue(str_lgrCalcGroupe);
 				lgrCalcGroupe = 0;
 			}
@@ -208,7 +220,6 @@ int ElValSet::Write(FILE * _fp) {
 			lgrCalcGroupe += 2 + 2 + 4 + elem->GetLength();  // Gr + Num + Lgr + LgrElem 
 			
 			printf ("lgrCalcGroupe apres %d\n",lgrCalcGroupe);
-
 		}		
 	}
 	
@@ -221,20 +232,15 @@ int ElValSet::Write(FILE * _fp) {
 		
 	// restent à tester les echecs en écriture (apres chaque fwrite)
 	
-	for (TagElValueHT::iterator tag = tagHt.begin();
-		  tag != tagHt.end();
-		  ++tag){
-
-		// Question :
-		// peut-on se passer des affectations?
-		// - passer l'adresse du resultat d'une fonction (???)
-		// - acceder au champ sans passer par un accesseur ?
+	for (TagElValueHT::iterator tag2 = tagHt.begin();
+		  tag2 != tagHt.end();
+		  ++tag2){
 		
-		gr =  tag->second->GetGroup();
-		el =  tag->second->GetElement();
-		lgr = tag->second->GetLength();
-		val = tag->second->GetValue().c_str();
-		vr =  tag->second->GetVR();
+		gr =  tag2->second->GetGroup();
+		el =  tag2->second->GetElement();
+		lgr = tag2->second->GetLength();
+		val = tag2->second->GetValue().c_str();
+		vr =  tag2->second->GetVR();
 			
 		fwrite ( &gr,(size_t)2 ,(size_t)1 ,_fp); 	//group
 		fwrite ( &el,(size_t)2 ,(size_t)1 ,_fp); 	//element
@@ -243,19 +249,47 @@ int ElValSet::Write(FILE * _fp) {
 		
 		// si on n'est pas en IMPLICIT VR voir pb (lgr  + VR)
 		
-		fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp); 			//lgr
+		fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp); 		//lgr
+		
+		tokens.erase(tokens.begin(),tokens.end());
+		Tokenize (tag2->second->GetValue(), tokens, "\\");
+		
+		//printf ("%04x %04x [%s] : [%s]\n",gr, el, vr.c_str(), val);
+		//if (tokens.size() > 1) { printf ("size : %d\n",tokens.size());}
+		
 		
 		if (vr == "US" || vr == "SS") {
+			/*
 			val_int16 = atoi(val);
 			ptr = &val_int16;
 			fwrite ( ptr,(size_t)2 ,(size_t)1 ,_fp);	
 			continue;
+			*/
+			
+			for (unsigned int i=0; i<tokens.size();i++) {
+				val_uint16 = atoi(tokens[i].c_str());		
+				ptr = &val_uint16;
+				fwrite ( ptr,(size_t)2 ,(size_t)1 ,_fp);
+			}
+			continue;
+			
 		}
 		if (vr == "UL" || vr == "SL") {
+			/*
 			val_int32 = atoi(val);
 			ptr = &val_int32;
 			fwrite ( ptr,(size_t)4 ,(size_t)1 ,_fp);	
 			continue;
+			*/
+			
+			
+			for (unsigned int i=0; i<tokens.size();i++) {
+				val_uint32 = atoi(tokens[i].c_str());		
+				ptr = &val_uint32;
+				fwrite ( ptr,(size_t)4 ,(size_t)1 ,_fp);
+			}
+			continue;				
+			
 		}	
 		
 		// Les pixels ne sont pas chargés dans l'element !
@@ -266,3 +300,173 @@ int ElValSet::Write(FILE * _fp) {
 		
 	return(1);
 }
+
+
+
+
+int ElValSet::WriteAcr(FILE * _fp) {
+
+
+// ATTENTION : fonction non terminée (commitée a titre de precaution)
+// ATTENTION : fusioner le code avec celui de lValSet::Write
+
+
+	guint16 gr, el;
+	guint32 lgr;
+	const char * val;
+	string vr;
+	guint32 val_uint32;
+	gint32  val_int32;
+	guint16 val_uint16;
+	gint16  val_int16;
+	
+	vector<string> tokens;
+	
+	void *ptr;
+	char str_lgrCalcGroupe[10];
+	
+	//string implicitVRTransfertSyntax = "1.2.840.10008.1.2"; // supprime par rapport à Write
+	
+	// Utilisées pour le calcul Group Length
+	int deja = 0;
+	guint32 lgrCalcGroupe=0;
+	ElValue *elem, *elemZ, *elemZPrec;
+	guint16 grCourant = 0;
+	
+	// Question :
+	// Comment pourrait-on tester si on est TrueDicom ou non ,
+	// (FileType est un champ de gdcmHeader ...)
+	//
+
+	// On parcourt la table pour recalculer la longueur des 'elements 0x0000'
+	// au cas ou un tag ai été ajouté par rapport à ce qui a été lu
+	// dans l'image native
+	//
+	// cf : code IdDcmWriteFile dans libido/src/dcmwrite.c
+		
+	if (1)  // Risque de pb dans le calcul des lgr de chaque groupe. On le saute pour le moment!
+	
+	for (TagElValueHT::iterator tag = tagHt.begin();
+		  tag != tagHt.end();
+		  ++tag){
+
+		elem = tag->second;
+		//printf("gr %04x el %04x lgr %d\n",elem->GetGroup(), elem->GetElement(), elem->GetLength());
+	
+		if ( (elem->GetGroup() != grCourant) &&	  
+			 (elem->GetGroup() != 0xfffe)	) { 	// On arrive sur un nv Groupe
+			 
+		//printf("Nouv Groupegr %04x el %04x \n",elem->GetGroup(), elem->GetElement());
+
+			elemZ = elem; 
+			 
+			if(elemZ->GetElement() != 0x0000) { 	// pas d'element 'Lgr groupe'
+				// On crée
+			 	gdcmDictEntry * tagZ = new gdcmDictEntry(grCourant, 0x0000, "UL");
+				elemZ = new ElValue(tagZ); // on le cree
+				elemZ->SetLength(4);
+				Add(elemZ);	 				 // On l'accroche à sa place	
+			}	
+			
+			if (deja) {
+				//sprintf(str_lgrCalcGroupe,"%d",lgrCalcGroupe);
+				elemZPrec->SetValue(str_lgrCalcGroupe);
+				lgrCalcGroupe = 0;
+			}
+			deja = 1;
+			
+			lgrCalcGroupe =  12; //2 + 2 + 4 + 4; // Gr + Num + Lgr + LgrGroupe
+			//printf ("lgrCalcGroupe %d\n",lgrCalcGroupe);
+			
+			elemZPrec = elemZ;
+	 		grCourant = elem->GetGroup();
+							
+		} else {		// On n'EST PAS sur un nv Groupe
+		
+			//printf ("lgrCalcGroupe avant : %d LgrElem %d\n",lgrCalcGroupe,elem->GetLength());
+
+			lgrCalcGroupe += 2 + 2 + 4 + elem->GetLength();  // Gr + Num + Lgr + LgrElem 
+			
+			//printf ("lgrCalcGroupe apres %d\n",lgrCalcGroupe);
+		}		
+	}
+	
+	// Si on fait de l'implicit VR little Endian 
+	// (pour moins se fairche sur processeur INTEL)
+	// penser a forcer le TRANSFERT SYNTAX UID
+	
+	// supprime par rapport à Write		
+	//SetElValueByNumber(implicitVRTransfertSyntax, 0x0002, 0x0010);	
+	//SetElValueLengthByNumber(18, 0x0002, 0x0010);  // Le 0 de fin de chaine doit etre stocké, dans ce cas	
+		
+	// restent à tester les echecs en écriture (apres chaque fwrite)
+	
+	for (TagElValueHT::iterator tag2 = tagHt.begin();
+		  tag2 != tagHt.end();
+		  ++tag2){
+
+		gr =  tag2->second->GetGroup();
+						// saut des groupes speciaux DICOM V3
+		if (gr < 0x0008) continue; 	// ajouté par rapport à Write
+						// saut des groupes impairs
+		if (gr %2) 	 continue; 	// ajouté par rapport à Write
+		
+		el =  tag2->second->GetElement();
+		lgr = tag2->second->GetLength();
+		val = tag2->second->GetValue().c_str();
+		vr =  tag2->second->GetVR();
+			
+		fwrite ( &gr,(size_t)2 ,(size_t)1 ,_fp); 	//group
+		fwrite ( &el,(size_t)2 ,(size_t)1 ,_fp); 	//element
+				
+		// si on n'est pas en IMPLICIT VR voir pb (lgr  + VR)
+		
+		fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp); 		//lgr
+		
+		tokens.erase(tokens.begin(),tokens.end());
+		Tokenize (tag2->second->GetValue(), tokens, "\\");
+		
+		if (vr == "US" || vr == "SS") {
+
+			for (unsigned int i=0; i<tokens.size();i++) {
+				val_uint16 = atoi(tokens[i].c_str());		
+				ptr = &val_uint16;
+				fwrite ( ptr,(size_t)2 ,(size_t)1 ,_fp);
+			}
+			continue;
+			
+		}
+		if (vr == "UL" || vr == "SL") {
+			for (unsigned int i=0; i<tokens.size();i++) {
+				val_uint32 = atoi(tokens[i].c_str());		
+				ptr = &val_uint32;
+				fwrite ( ptr,(size_t)4 ,(size_t)1 ,_fp);
+			}
+			continue;				
+			
+		}	
+		
+		// Les pixels ne sont pas chargés dans l'element !
+		if ((gr == 0x7fe0) && (el == 0x0010) ) break;
+
+		fwrite ( val,(size_t)lgr ,(size_t)1 ,_fp); //valeur Elem
+	}
+		
+	return(1);
+}
+
+
+
+
+// mettre ça dans une bibliothèque d'utilitaires ?
+// ca peut servir
+
+static void Tokenize (const string& str, vector<string>& tokens, const string& delimiters = " ") {
+	string::size_type lastPos = str.find_first_not_of(delimiters,0);
+	string::size_type pos     = str.find_first_of(delimiters,lastPos);
+	while (string::npos != pos || string::npos != lastPos) {
+		tokens.push_back(str.substr(lastPos, pos - lastPos));
+		lastPos = str.find_first_not_of(delimiters, pos);
+		pos = str.find_first_of(delimiters, lastPos);
+	}
+}