// $Id: gdcmElValSet.cxx,v 1.20 2003/03/13 17:00:10 jpr Exp $ #include "gdcmUtil.h" #include "gdcmElValSet.h" TagElValueHT & gdcmElValSet::GetTagHt(void) { return tagHt; } void gdcmElValSet::Add(gdcmElValue * newElValue) { tagHt [newElValue->GetKey()] = newElValue; NameHt[newElValue->GetName()] = newElValue; } void gdcmElValSet::Print(ostream & os) { for (TagElValueHT::iterator tag = tagHt.begin(); tag != tagHt.end(); ++tag){ os << tag->first << ": "; os << "[" << tag->second->GetValue() << "]"; os << "[" << tag->second->GetName() << "]"; os << "[" << tag->second->GetVR() << "]" << endl; } } 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->GetVR() << "]" << endl; } } gdcmElValue* gdcmElValSet::GetElementByNumber(guint32 group, guint32 element) { TagKey key = gdcmDictEntry::TranslateToKey(group, element); if ( ! tagHt.count(key)) return (gdcmElValue*)0; if (tagHt.count(key) > 1) dbg.Verbose(0, "gdcmElValSet::GetElementByNumber", "multiple entries for this key (FIXME) !"); return tagHt.find(key)->second; } gdcmElValue* gdcmElValSet::GetElementByName(string TagName) { if ( ! NameHt.count(TagName)) return (gdcmElValue*)0; if (NameHt.count(TagName) > 1) dbg.Verbose(0, "gdcmElValSet::GetElement", "multipe entries for this key (FIXME) !"); return NameHt.find(TagName)->second; } 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, "gdcmElValSet::GetElValueByNumber", "multiple entries for this key (FIXME) !"); return tagHt.find(key)->second->GetValue(); } 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, "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); } int gdcmElValSet::SetElValueByName(string content, string TagName) { if ( ! NameHt.count(TagName)) return 0; if (NameHt.count(TagName) > 1) { dbg.Verbose(0, "gdcmElValSet::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 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 gdcmElValSet::SetElValueLengthByName(guint32 l, string TagName) { if ( ! NameHt.count(TagName)) return 0; if (NameHt.count(TagName) > 1) { dbg.Verbose(0, "gdcmElValSet::SetElValueByName", "multipe entries for this key (FIXME) !"); return 0; } 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 (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 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