X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Testing%2FTestCopyDicom.cxx;h=5b35c48bd57c22bf2cc8dc365f98fa28d81a4b01;hb=7eb058ef98b0003b379bf0db33e5c4c673b7f1e6;hp=e67a9926ce380aebb25487a81fa6963cd9a313e3;hpb=152194ed48de20a5189b610bd8d3d26fea95c2ce;p=gdcm.git diff --git a/Testing/TestCopyDicom.cxx b/Testing/TestCopyDicom.cxx index e67a9926..5b35c48b 100644 --- a/Testing/TestCopyDicom.cxx +++ b/Testing/TestCopyDicom.cxx @@ -2,6 +2,7 @@ #include "gdcmFile.h" #include "gdcmDocument.h" #include "gdcmValEntry.h" +#include "gdcmBinEntry.h" //Generated file: #include "gdcmDataImages.h" @@ -40,7 +41,7 @@ bool RemoveFile(const char* source) // Here we load a gdcmFile and then try to create from scratch a copy of it, // copying field by field the dicom image -int TestCopyDicom(int , char* []) +int TestCopyDicom(int argc, char* argv[]) { int i =0; int retVal = 0; //by default this is an error @@ -55,7 +56,7 @@ int TestCopyDicom(int , char* []) if( FileExists( output.c_str() ) ) { - std::cerr << "Don't try to cheat, I am removing the file anyway" << std::endl; + // std::cerr << "Don't try to cheat, I am removing the file anyway" << std::endl; if( !RemoveFile( output.c_str() ) ) { std::cerr << "Ouch, the file exist, but I cannot remove it" << std::endl; @@ -66,50 +67,68 @@ int TestCopyDicom(int , char* []) gdcmFile *original = new gdcmFile( filename ); gdcmFile *copy = new gdcmFile( output ); + TagDocEntryHT & Ht = original->GetHeader()->GetEntry(); + //First of all copy the header field by field - TagNameHT & nameHt = original->GetHeader()->GetPubDict()->GetEntriesByName(); - for (TagNameHT::iterator tag = nameHt.begin(); tag != nameHt.end(); ++tag) + + // Warning :Accessor gdcmElementSet::GetEntry() should not exist + // It was commented out by Mathieu, that was a *good* idea + // (the user does NOT have to know the way we implemented the Header !) + // Waiting for a 'clean' solution, I keep the method ...JPRx + + gdcmDocEntry* d; + + for (TagDocEntryHT::iterator tag = Ht.begin(); tag != Ht.end(); ++tag) { - if (tag->second->GetVR() == "SQ") //to skip pb of SQ recursive exploration - continue; - - //According to JPR I should also skip those: -// if (tag->second->GetVR() == "unkn") //to skip pb of SQ recursive exploration -// continue; - - //no clue what to do with this one - //std::cerr << "Reading: " << tag->second->GetVR() << std::endl; - std::string value = ((gdcmValEntry*)(tag->second))->GetValue(); -// if( value.find( "gdcm::NotLoaded" ) == 0 ) -// continue; - - //no clue what to do with this one - if( value.find( "gdcm::Loaded" ) == 0 ) - continue; - - copy->GetHeader()->ReplaceOrCreateByNumber( - value, - tag->second->GetGroup(), - tag->second->GetElement() ); + d = tag->second; + if ( gdcmBinEntry* b = dynamic_cast(d) ) + { + copy->GetHeader()->ReplaceOrCreateByNumber( + b->GetVoidArea(), + b->GetLength(), + b->GetGroup(), + b->GetElement(), + b->GetVR() ); + } + else if ( gdcmValEntry* v = dynamic_cast(d) ) + { + copy->GetHeader()->ReplaceOrCreateByNumber( + v->GetValue(), + v->GetGroup(), + v->GetElement(), + v->GetVR() ); + } + else + { + // We skip pb of SQ recursive exploration + //std::cout << "Skipped Sequence " + // << "------------- " << d->GetVR() << " "<< std::hex + // << d->GetGroup() << " " << d->GetElement() + // << std::endl; + } } size_t dataSize = original->GetImageDataSize(); void *imageData = original->GetImageData(); copy->SetImageData(imageData, dataSize); - //original->GetHeader()->SetImageDataSize(dataSize); + original->GetHeader()->SetImageDataSize(dataSize); + + copy->WriteDcmExplVR( output ); + + delete original; + delete copy; - //copy->GetHeader()->PrintEntry(); + copy = new gdcmFile( output ); - //copy->WriteDcmExplVR( output ); - //Is the file written still gdcm parsable ? - gdcmFile check( output ); - retVal += !check.GetHeader()->IsReadable(); + if ( !copy->GetHeader()->IsReadable() ) + { + retVal +=1; + std::cout << output << " Failed" << std::endl; + } + i++; } - return retVal; } - -