[Dcmlib] TestFromScratch
Mathieu Malaterre
mathieu.malaterre at kitware.com
Fri Dec 10 16:40:42 CET 2004
Correction, ca ne marche que sous linux. Je l'ajoute dans Exemple
seulement. Des que c'est fixer ca passera dans Tests.
Mathieu
Mathieu Malaterre wrote:
> Woohoo !
>
> J'ai trouver, il faut faire:
>
> h2->ReplaceOrCreateByNumber(imageData,dataSize,0x7fe0, 0x0010,"PXL");
>
> puis
>
> f2->SetImageData(imageData, dataSize);
>
> Et ca marche !
>
> Je l'ajoute dans les exemples ou pas ?
>
> Mathieu
> Ps: ca n'a plus bcp de sens si on a vrai writer, mais c'est pour gdcm
> 0.6 dans ITK. Donc la question c'est est-ce qu'on veut conserver ce
> comportement on bien les user ne sont pas senser faire aussi bas niveau ?
>
> Mathieu Malaterre wrote:
>
>> ben y'a vraiment plus rien qui marche ou c moi qui a besoin de cafe ?
>>
>>
>
> ------------------------------------------------------------------------
>
> #include "gdcmHeader.h"
> #include "gdcmFile.h"
> #include "gdcmDictEntry.h"
> #include "gdcmDocEntry.h"
> #include "gdcmValEntry.h"
> #include "gdcmDebug.h"
>
> // The whole purpose of this example is to create a DICOM from scratch
> // well almost the only thing we allow use is to copy the string entry
> // From the previous image read.
>
> int TestFromScratch(int argc, char *argv[])
> {
> if( argc < 2 )
> {
> std::cerr << "Usage: " << argv[0] << " Image.dcm" << std::endl;
> return 1;
> }
>
> // Doesn't seems to do anything:
> dbg.SetDebug(-1);
> // Doesn't link properly:
> //gdcm::Debug::GetReference().SetDebug(1);
>
> std::string filename = argv[1];
> //gdcm::File *f1 = new gdcm::File( "/home/malaterre/Creatis/gdcmData/012345.002.050.dcm" );
> gdcm::File *f1 = new gdcm::File( filename );
> gdcm::Header *h1 = f1->GetHeader();
>
> int dataSize = f1->GetImageDataSize();
> std::cout << "DataSize: " << dataSize << std::endl;
> // Since we know the image is 16bits:
> uint8_t* imageData = f1->GetImageData();
>
> // Hopefully default to something
> gdcm::Header *h2 = new gdcm::Header();
>
> const gdcm::TagDocEntryHT &nameHt = h1->GetTagHT();
> //const gdcm::TagDocEntryHT &nameHt = h2->GetTagHT(); seems to be empty ??
> for (gdcm::TagDocEntryHT::const_iterator tag = nameHt.begin(); tag != nameHt.end(); ++tag)
> {
> //Copy string only:
> if (tag->second->GetVR().find("SQ") == 0)
> {
> // skip DICOM SeQuence, otherwise following cast will crash
> continue;
> }
> std::string temp = ((gdcm::ValEntry*)(tag->second))->GetValue();
> const std::string &name = tag->second->GetName();
> const std::string &value = ((gdcm::ValEntry*)(tag->second))->GetValue();
> if( name != "unkn"
> && value.find( "gdcm::NotLoaded" ) != 0
> && value.find( "gdcm::Binary" ) != 0
> && value.find( "gdcm::Loaded" ) != 0 )
> {
> //std::cout << name << "," << temp << std::endl;
> gdcm::DictEntry *dictEntry = h2->GetPubDict()->GetDictEntryByName(name);
> h2->ReplaceOrCreateByNumber( temp, dictEntry->GetGroup(), dictEntry->GetElement());
> }
> }
>
> h2->ReplaceOrCreateByNumber( imageData, dataSize, 0x7fe0, 0x0010, "PXL" );
> h2->Print( std::cout );
>
> gdcm::File *f2 = new gdcm::File( h2 );
>
> // SetImageData does not work because you need to have called
> // somwhow ReplaceOrCreateByNumber before with the proper group/element
> f2->SetImageData(imageData, dataSize);
>
> f2->WriteDcmExplVR( "output.dcm" );
>
> delete f1;
> delete f2;
> delete h2;
>
> return 0;
> }
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Dcmlib mailing list
> Dcmlib at creatis.insa-lyon.fr
> http://www.creatis.insa-lyon.fr/mailman/listinfo/dcmlib
More information about the Dcmlib
mailing list