X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmElValSet.cxx;h=d716e6af4afa7782926199f484256a50146216dd;hb=885580b981648efd2b861bb92873c43d09757f27;hp=9101d26f97ab2541b25bdf4a44110a7620d27318;hpb=e8be02129ece7a8c7dd90307a8c83e1ef5d157b6;p=gdcm.git diff --git a/src/gdcmElValSet.cxx b/src/gdcmElValSet.cxx index 9101d26f..d716e6af 100644 --- a/src/gdcmElValSet.cxx +++ b/src/gdcmElValSet.cxx @@ -1,8 +1,14 @@ -#include "gdcmlib.h" +// gdcmElValSet.cxx + +#include "gdcm.h" #include "gdcmUtil.h" +TagElValueHT & ElValSet::GetTagHt(void) { + return tagHt; +} + void ElValSet::Add(ElValue * newElValue) { - tagHt[newElValue->GetKey()] = newElValue; + tagHt [newElValue->GetKey()] = newElValue; NameHt[newElValue->GetName()] = newElValue; } @@ -15,7 +21,7 @@ void ElValSet::Print(ostream & os) { os << "[" << tag->second->GetName() << "]"; os << "[" << tag->second->GetVR() << "]" << endl; } -} +} void ElValSet::PrintByName(ostream & os) { for (TagElValueNameHT::iterator tag = NameHt.begin(); @@ -23,27 +29,240 @@ void ElValSet::PrintByName(ostream & os) { ++tag){ os << tag->first << ": "; os << "[" << tag->second->GetValue() << "]"; - os << "[" << tag->second->GetKey() << "]"; + os << "[" << tag->second->GetKey() << "]"; os << "[" << tag->second->GetVR() << "]" << endl; } } -string ElValSet::GetElValue(guint32 group, guint32 element) { +ElValue* ElValSet::GetElementByNumber(guint32 group, guint32 element) { TagKey key = gdcmDictEntry::TranslateToKey(group, element); if ( ! tagHt.count(key)) - return "UNFOUND"; + return (ElValue*)0; if (tagHt.count(key) > 1) - dbg.Verbose(0, "ElValSet::GetElValue", + dbg.Verbose(0, "ElValSet::GetElementByNumber", + "multiple entries for this key (FIXME) !"); + return tagHt.find(key)->second; +} + +ElValue* ElValSet::GetElementByName(string TagName) { + if ( ! NameHt.count(TagName)) + return (ElValue*)0; + if (NameHt.count(TagName) > 1) + dbg.Verbose(0, "ElValSet::GetElement", + "multipe entries for this key (FIXME) !"); + return NameHt.find(TagName)->second; +} + +string ElValSet::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", "multiple entries for this key (FIXME) !"); return tagHt.find(key)->second->GetValue(); } -string ElValSet::GetElValue(string TagName) { +string ElValSet::GetElValueByName(string TagName) { if ( ! NameHt.count(TagName)) - return "UNFOUND"; + return "gdcm::Unfound"; if (NameHt.count(TagName) > 1) dbg.Verbose(0, "ElValSet::GetElValue", "multipe entries for this key (FIXME) !"); return NameHt.find(TagName)->second->GetValue(); } +int ElValSet::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", + "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); +} + +int ElValSet::SetElValueByName(string content, string TagName) { + if ( ! NameHt.count(TagName)) + return 0; + if (NameHt.count(TagName) > 1) { + dbg.Verbose(0, "ElValSet::SetElValueByName", + "multipe entries for this key (FIXME) !"); + return 0; + } + NameHt.find(TagName)->second->SetValue(content); + NameHt.find(TagName)->second->SetLength(strlen(content.c_str())); + return(1); +} + + +int ElValSet::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, "ElValSet::SetElValueLengthByNumber", + "multiple entries for this key (FIXME) !"); + return (0); + } + // m à j LgrElem + tagHt[key]->SetLength(l); + return(1); +} + + +int ElValSet::SetElValueLengthByName(guint32 l, string TagName) { + if ( ! NameHt.count(TagName)) + return 0; + if (NameHt.count(TagName) > 1) { + dbg.Verbose(0, "ElValSet::SetElValueByName", + "multipe entries for this key (FIXME) !"); + return 0; + } + NameHt.find(TagName)->second->SetLength(l); + return(1); +} + + +int ElValSet::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_int32; + guint16 val_int16; + 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; + 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 + + 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 + + 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){ + + // Question : + // peut-on se passer des affectations? + // - passer l'adresse du resultat d'une fonction (???) + // - acceder au champ sans passer par un accesseur ? + + 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 + + if (vr == "US" || vr == "SS") { + val_int16 = atoi(val); + ptr = &val_int16; + 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; + } + + // 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); +}