X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmElValSet.cxx;h=871c1aebddae0da10667a66de74c1155b9bed030;hb=58eb3ab71df4f87544d370c045386e7ced745da7;hp=70f21fc15ffc47f8ae25538ca06f97651e3aee99;hpb=d0e0404925f8a3dc6f771ed598c75f9dea95912d;p=gdcm.git diff --git a/src/gdcmElValSet.cxx b/src/gdcmElValSet.cxx index 70f21fc1..871c1aeb 100644 --- a/src/gdcmElValSet.cxx +++ b/src/gdcmElValSet.cxx @@ -1,18 +1,19 @@ -// gdcmElValSet.cxx +// $Id: gdcmElValSet.cxx,v 1.21 2003/03/14 14:26:01 jpr Exp $ -#include "gdcm.h" #include "gdcmUtil.h" +#include "gdcmElValSet.h" -TagElValueHT & ElValSet::GetTagHt(void) { + +TagElValueHT & gdcmElValSet::GetTagHt(void) { return tagHt; } -void ElValSet::Add(ElValue * newElValue) { - tagHt[newElValue->GetKey()] = newElValue; +void gdcmElValSet::Add(gdcmElValue * newElValue) { + tagHt [newElValue->GetKey()] = newElValue; NameHt[newElValue->GetName()] = newElValue; } -void ElValSet::Print(ostream & os) { +void gdcmElValSet::Print(ostream & os) { for (TagElValueHT::iterator tag = tagHt.begin(); tag != tagHt.end(); ++tag){ @@ -21,111 +22,434 @@ void ElValSet::Print(ostream & os) { os << "[" << tag->second->GetName() << "]"; os << "[" << tag->second->GetVR() << "]" << endl; } -} - - -int ElValSet::Write(FILE * _fp) { - - guint16 gr, el; - guint32 lgr; - const char * val; - // restent à tester les echecs en écriture - 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 ? - gr = tag->second->GetGroup(); - el = tag->second->GetElement(); - lgr = tag->second->GetLength(); - val = tag->second->GetValue().c_str(); - - fwrite ( &gr,(size_t)2 ,(size_t)1 ,_fp); //group - fwrite ( &el,(size_t)2 ,(size_t)1 ,_fp); //element - //fwrite ( tag->second->GetVR(),(size_t)2 ,(size_t)1 ,_fp); //VR - - // voir pb lgr + VR - fwrite ( &lgr,(size_t)4 ,(size_t)1 ,_fp); //lgr - // voir pb des int16 et int32 : les identifier, les convertir, modifier la longueur - fwrite ( val,(size_t)lgr ,(size_t)1 ,_fp); //valeur Elem - } - return(1); -} -void ElValSet::PrintByName(ostream & os) { +void gdcmElValSet::PrintByName(ostream & os) { for (TagElValueNameHT::iterator tag = NameHt.begin(); tag != NameHt.end(); ++tag){ os << tag->first << ": "; os << "[" << tag->second->GetValue() << "]"; - os << "[" << tag->second->GetKey() << "]"; + os << "[" << tag->second->GetKey() << "]"; os << "[" << tag->second->GetVR() << "]" << endl; } } -ElValue* ElValSet::GetElementByNumber(guint32 group, guint32 element) { +gdcmElValue* gdcmElValSet::GetElementByNumber(guint32 group, guint32 element) { TagKey key = gdcmDictEntry::TranslateToKey(group, element); if ( ! tagHt.count(key)) - return (ElValue*)0; + return (gdcmElValue*)0; if (tagHt.count(key) > 1) - dbg.Verbose(0, "ElValSet::GetElementByNumber", + dbg.Verbose(0, "gdcmElValSet::GetElementByNumber", "multiple entries for this key (FIXME) !"); return tagHt.find(key)->second; } -ElValue* ElValSet::GetElementByName(string TagName) { +gdcmElValue* gdcmElValSet::GetElementByName(string TagName) { if ( ! NameHt.count(TagName)) - return (ElValue*)0; + return (gdcmElValue*)0; if (NameHt.count(TagName) > 1) - dbg.Verbose(0, "ElValSet::GetElement", + dbg.Verbose(0, "gdcmElValSet::GetElement", "multipe entries for this key (FIXME) !"); return NameHt.find(TagName)->second; } -string ElValSet::GetElValueByNumber(guint32 group, guint32 element) { +string gdcmElValSet::GetElValueByNumber(guint32 group, guint32 element) { TagKey key = gdcmDictEntry::TranslateToKey(group, element); if ( ! tagHt.count(key)) return "gdcm::Unfound"; if (tagHt.count(key) > 1) - dbg.Verbose(0, "ElValSet::GetElValueByNumber", + dbg.Verbose(0, "gdcmElValSet::GetElValueByNumber", "multiple entries for this key (FIXME) !"); return tagHt.find(key)->second->GetValue(); } -int ElValSet::SetElValueByNumber(string content, guint32 group, guint32 element) { +string gdcmElValSet::GetElValueByName(string TagName) { + if ( ! NameHt.count(TagName)) + return "gdcm::Unfound"; + if (NameHt.count(TagName) > 1) + dbg.Verbose(0, "gdcmElValSet::GetElValue", + "multipe entries for this key (FIXME) !"); + return NameHt.find(TagName)->second->GetValue(); +} + +int gdcmElValSet::SetElValueByNumber(string content, + guint32 group, guint32 element) { TagKey key = gdcmDictEntry::TranslateToKey(group, element); if ( ! tagHt.count(key)) return 0; if (tagHt.count(key) > 1) { - dbg.Verbose(0, "ElValSet::SetElValueByNumber", + dbg.Verbose(0, "gdcmElValSet::SetElValueByNumber", "multiple entries for this key (FIXME) !"); return (0); - } - + } tagHt[key]->SetValue(content); + + // Question : m à j LgrElem ? + tagHt[key]->SetLength(strlen(content.c_str())); + + // FIXME should we really update the element length ? + tagHt[key]->SetLength(content.length()); + return(1); } -string ElValSet::GetElValueByName(string TagName) { +int gdcmElValSet::SetElValueByName(string content, string TagName) { if ( ! NameHt.count(TagName)) - return "gdcm::Unfound"; - if (NameHt.count(TagName) > 1) - dbg.Verbose(0, "ElValSet::GetElValue", + return 0; + if (NameHt.count(TagName) > 1) { + dbg.Verbose(0, "gdcmElValSet::SetElValueByName", "multipe entries for this key (FIXME) !"); - return NameHt.find(TagName)->second->GetValue(); + return 0; + } + NameHt.find(TagName)->second->SetValue(content); + NameHt.find(TagName)->second->SetLength(strlen(content.c_str())); + return(1); +} + + +int gdcmElValSet::SetElValueLengthByNumber(guint32 l, + guint32 group, guint32 element) { + TagKey key = gdcmDictEntry::TranslateToKey(group, element); + if ( ! tagHt.count(key)) + return 0; + if (tagHt.count(key) > 1) { + dbg.Verbose(0, "gdcmElValSet::SetElValueLengthByNumber", + "multiple entries for this key (FIXME) !"); + return (0); + } + // FIXME JPR: comments in English please ! + // m à j LgrElem + tagHt[key]->SetLength(l); + return(1); } -int ElValSet::SetElValueByName(string content, string TagName) { + +int gdcmElValSet::SetElValueLengthByName(guint32 l, string TagName) { if ( ! NameHt.count(TagName)) return 0; if (NameHt.count(TagName) > 1) { - dbg.Verbose(0, "ElValSet::SetElValue", + dbg.Verbose(0, "gdcmElValSet::SetElValueByName", "multipe entries for this key (FIXME) !"); return 0; } - NameHt.find(TagName)->second->SetValue(content); + NameHt.find(TagName)->second->SetLength(l); + return(1); +} + + +int gdcmElValSet::Write(FILE * _fp) { + +// ATTENTION : fonction non terminée (commitée a titre de precaution) + + guint16 gr, el; + guint32 lgr; + const char * val; + string vr; + guint32 val_uint32; + gint32 val_int32; + guint16 val_uint16; + gint16 val_int16; + + vector tokens; + + void *ptr; + char str_lgrCalcGroupe[10]; + + string implicitVRTransfertSyntax = "1.2.840.10008.1.2"; + + // Utilisées pour le calcul Group Length + int deja = 0; + guint32 lgrCalcGroupe=0; + gdcmElValue *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 (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){ + + 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 gdcmElValue(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 + + 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(); + 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 + + //fwrite ( vr,(size_t)2 ,(size_t)1 ,_fp); //VR + + // 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, "\\"); + + //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; + + 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; + gdcmElValue *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 (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){ + + 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 gdcmElValue(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