X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Testing%2FTestWriteSimple.cxx;h=5e981eee707f219dd9ec66348f38f53f765bcfd4;hb=d9f8b4b69e332d38e9c883ba60f40dc858883de9;hp=273865090d82abc3f4aae7e42515dc5e083fbe3e;hpb=ca1c8a4541a05aaec4fbe7bb3b23e18cbe2d31bd;p=gdcm.git diff --git a/Testing/TestWriteSimple.cxx b/Testing/TestWriteSimple.cxx index 27386509..5e981eee 100644 --- a/Testing/TestWriteSimple.cxx +++ b/Testing/TestWriteSimple.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: TestWriteSimple.cxx,v $ Language: C++ - Date: $Date: 2005/07/27 03:02:12 $ - Version: $Revision: 1.33 $ + Date: $Date: 2007/09/28 14:20:22 $ + Version: $Revision: 1.53 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -25,6 +25,8 @@ #include "gdcmFile.h" #include "gdcmFileHelper.h" #include "gdcmDebug.h" +#include "gdcmGlobal.h" +#include "gdcmDictSet.h" #include #include @@ -35,8 +37,8 @@ typedef struct int sizeY; // Size Y of the image int sizeZ; // Size Z of the image int components; // Number of components for a pixel - int componentSize; // Component size (in bits : 8, 16) - int componentUse ; // Component size (in bits) + int componentSize; // Component size (in bits : 8, 16) // Bits Allocated + int componentUse ; // Component size (in bits) // Bits Stored int sign; // Sign of components char writeMode; // Write mode // - 'a' : ACR @@ -45,84 +47,110 @@ typedef struct } Image; Image Images [] = { - {256, 256, 1, 1, 8, 8, 0, 'a'}, - {256, 256, 1, 1, 8, 8, 0, 'e'}, - {256, 256, 1, 1, 8, 8, 0, 'i'}, - - {512, 256, 1, 1, 8, 8, 0, 'a'}, - {512, 256, 1, 1, 8, 8, 0, 'e'}, - {512, 256, 1, 1, 8, 8, 0, 'i'}, - - {256, 512, 1, 1, 8, 8, 0, 'a'}, - {256, 512, 1, 1, 8, 8, 0, 'e'}, - {256, 512, 1, 1, 8, 8, 0, 'i'}, - - {256, 512, 1, 1, 16, 16, 0, 'a'}, - {256, 512, 1, 1, 16, 16, 0, 'e'}, - {256, 512, 1, 1, 16, 16, 0, 'i'}, - {256, 512, 1, 1, 16, 16, 0, 'a'}, - {256, 512, 1, 1, 16, 16, 0, 'e'}, - {256, 512, 1, 1, 16, 16, 0, 'i'}, - - {512, 256, 10, 1, 8, 8, 0, 'a'}, - {512, 256, 10, 1, 8, 8, 0, 'e'}, - {512, 256, 10, 1, 8, 8, 0, 'i'}, - {512, 256, 10, 3, 8, 8, 0, 'a'}, - {512, 256, 10, 3, 8, 8, 0, 'e'}, - {512, 256, 10, 3, 8, 8, 0, 'i'}, - - {256, 256, 1, 1, 8, 8, 1, 'a'}, - {256, 256, 1, 1, 8, 8, 1, 'e'}, - {256, 256, 1, 1, 8, 8, 1, 'i'}, - - {512, 256, 1, 1, 8, 8, 1, 'a'}, - {512, 256, 1, 1, 8, 8, 1, 'e'}, - {512, 256, 1, 1, 8, 8, 1, 'i'}, - - {256, 512, 1, 1, 8, 8, 1, 'a'}, - {256, 512, 1, 1, 8, 8, 1, 'e'}, - {256, 512, 1, 1, 8, 8, 1, 'i'}, - - {256, 512, 1, 1, 16, 16, 1, 'a'}, - {256, 512, 1, 1, 16, 16, 1, 'e'}, - {256, 512, 1, 1, 16, 16, 1, 'i'}, - {256, 512, 1, 1, 16, 16, 1, 'a'}, - {256, 512, 1, 1, 16, 16, 1, 'e'}, - {256, 512, 1, 1, 16, 16, 1, 'i'}, - - {512, 256, 10, 1, 8, 8, 1, 'a'}, - {512, 256, 10, 1, 8, 8, 1, 'e'}, - {512, 256, 10, 1, 8, 8, 1, 'i'}, - {512, 256, 10, 3, 8, 8, 1, 'a'}, - {512, 256, 10, 3, 8, 8, 1, 'e'}, - {512, 256, 10, 3, 8, 8, 1, 'i'}, +// these ones are use to check further oddities. + + {63, 127, 1, 1, 16, 8, 0 ,'e'}, + {63, 127, 1, 1, 16, 9, 0 ,'e'}, + {63, 127, 1, 1, 16, 10, 0 ,'e'}, + {63, 127, 1, 1, 16, 11, 0 ,'e'}, + {63, 127, 1, 1, 16, 12, 0 ,'e'}, + {63, 127, 1, 1, 16, 13, 0 ,'e'}, + {63, 127, 1, 1, 16, 14, 0 ,'e'}, + {63, 127, 1, 1, 16, 15, 0 ,'e'}, + {63, 127, 1, 1, 16, 16, 0 ,'e'}, + + {63, 127, 1, 1, 32, 32, 0 ,'e'}, // Pixel VR should be OL? + + {128, 128, 1, 1, 8, 8, 0, 'e'}, + {256, 128, 1, 1, 8, 8, 0, 'a'}, + {128, 128, 1, 1, 8, 8, 0, 'i'}, + + {128, 128, 1, 1, 8, 8, 0, 'a'}, + {256, 128, 1, 1, 8, 8, 0, 'i'}, + + + {256, 128, 1, 1, 8, 8, 0, 'e'}, + {128, 128, 1, 1, 16, 16, 0, 'e'}, + {128, 256, 1, 1, 16, 16, 0, 'e'}, + + {128, 256, 1, 1, 8, 8, 0, 'e'}, + {128, 256, 1, 1, 8, 8, 0, 'i'}, + {128, 256, 1, 1, 8, 8, 0, 'a'}, + + {128, 256, 1, 1, 16, 16, 0, 'i'}, + {128, 256, 1, 1, 16, 16, 0, 'i'}, + {128, 256, 1, 1, 16, 16, 0, 'a'}, + {128, 256, 1, 1, 16, 16, 0, 'a'}, + + {256, 128, 10, 1, 8, 8, 0, 'e'}, + {256, 128, 10, 3, 8, 8, 0, 'e'}, + {256, 128, 10, 3, 8, 8, 0, 'i'}, + {256, 128, 10, 1, 8, 8, 0, 'i'}, + {256, 128, 10, 1, 8, 8, 0, 'a'}, + {256, 128, 10, 3, 8, 8, 0, 'a'}, + + {128, 128, 1, 1, 8, 8, 1, 'e'}, + {128, 128, 1, 1, 8, 8, 1, 'i'}, + {128, 128, 1, 1, 8, 8, 1, 'a'}, + + {256, 128, 1, 1, 8, 8, 1, 'e'}, + {256, 128, 1, 1, 8, 8, 1, 'i'}, + {256, 128, 1, 1, 8, 8, 1, 'a'}, + + {128, 256, 1, 1, 8, 8, 1, 'a'}, + {128, 256, 1, 1, 8, 8, 1, 'e'}, + {128, 256, 1, 1, 8, 8, 1, 'i'}, + + {128, 256, 1, 1, 16, 16, 1, 'e'}, + {128, 256, 1, 1, 16, 16, 1, 'e'}, + {128, 256, 1, 1, 16, 16, 1, 'i'}, + {128, 256, 1, 1, 16, 16, 1, 'i'}, + {128, 256, 1, 1, 16, 16, 1, 'a'}, + {128, 256, 1, 1, 16, 16, 1, 'a'}, + + {256, 128, 10, 1, 8, 8, 1, 'e'}, + {256, 128, 10, 1, 8, 8, 1, 'i'}, + {256, 128, 10, 1, 8, 8, 1, 'a'}, + + {256, 128, 10, 3, 8, 8, 1, 'e'}, + {256, 128, 10, 3, 8, 8, 1, 'i'}, + {256, 128, 10, 3, 8, 8, 1, 'a'}, {0, 0, 1, 1, 8, 8, 0, 'i'} // to find the end }; + +const unsigned int MAX_NUMBER_OF_DIFFERENCE = 10; + int WriteSimple(Image &img) { + + std::cout << "======================= WriteSimple =========(begin of processing current image)" << std::endl; std::ostringstream fileName; fileName.str(""); fileName << "TestWriteSimple"; -// Step 1 : Create the header of the image +// Step 1 : Create an empty FileHelper + std::cout << " 1..."; - gdcm::File *fileToBuild = new gdcm::File(); + GDCM_NAME_SPACE::FileHelper *fileH = GDCM_NAME_SPACE::FileHelper::New(); + + // Get the (empty) image header. + GDCM_NAME_SPACE::File *fileToBuild = fileH->GetFile(); std::ostringstream str; // Set the image size str.str(""); str << img.sizeX; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0011); // Columns + fileToBuild->InsertEntryString(str.str(),0x0028,0x0011,"US"); // Columns str.str(""); str << img.sizeY; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0010); // Rows + fileToBuild->InsertEntryString(str.str(),0x0028,0x0010,"US"); // Rows if(img.sizeZ>1) { str.str(""); str << img.sizeZ; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0008); // Number of Frames + fileToBuild->InsertEntryString(str.str(),0x0028,0x0008, "IS"); // Number of Frames } fileName << "-" << img.sizeX << "-" << img.sizeY << "-" << img.sizeZ; @@ -130,44 +158,62 @@ int WriteSimple(Image &img) // Set the pixel type str.str(""); str << img.componentSize; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0100); // Bits Allocated + fileToBuild->InsertEntryString(str.str(),0x0028,0x0100,"US"); // Bits Allocated str.str(""); str << img.componentUse; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0101); // Bits Stored + fileToBuild->InsertEntryString(str.str(),0x0028,0x0101,"US"); // Bits Stored str.str(""); - str << ( img.componentSize - 1 ); - fileToBuild->InsertValEntry(str.str(),0x0028,0x0102); // High Bit + str << ( img.componentUse - 1 ); + fileToBuild->InsertEntryString(str.str(),0x0028,0x0102,"US"); // High Bit // Set the pixel representation str.str(""); str << img.sign; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0103); // Pixel Representation + fileToBuild->InsertEntryString(str.str(),0x0028,0x0103,"US"); // Pixel Representation fileName << "-" << img.componentSize; if(img.sign == 0) fileName << "U"; else fileName << "S"; - + + fileName << "-" << img.componentSize << "-" << img.componentUse; + + fileToBuild->InsertEntryString("0",0x0008,0x0000,"UL"); // Should be removed + // except for ACR switch (img.writeMode) { case 'a' : - fileName << ".ACR"; break; + fileName << ".ACR"; + break; case 'e' : - fileName << ".EXPL"; break; + fileName << ".EXPL"; + break; case 'i' : - fileName << ".IMPL"; break; -} + fileName << ".IMPL"; + break; + } + if(img.componentSize == 32) + { + // Create a Private DataElement; VR =OL + fileToBuild->InsertEntryString("gdcm test OL",0x0009,0x0010,"LO"); + uint32_t binArea[4]; + binArea[0] = 0x01234567; + binArea[1] = 0x89abcdef; + binArea[2] = 0x2468ace0; + binArea[2] = 0xfdb97531; + fileToBuild->InsertEntryBinArea((uint8_t *)binArea,16,0x0009,0x0100,"OL"); + } + + + std::cout << "[" << fileName.str() << "]..."; // Set the samples per pixel str.str(""); str << img.components; - fileToBuild->InsertValEntry(str.str(),0x0028,0x0002); // Samples per Pixel - - - + fileToBuild->InsertEntryString(str.str(),0x0028,0x0002,"US"); // Samples per Pixel // Step 2 : Create the output image std::cout << "2..."; @@ -200,15 +246,14 @@ int WriteSimple(Image &img) } } -// Step 3 : Create the file of the image +// Step 3 : Set the image Pixel Data std::cout << "3..."; - gdcm::FileHelper *fileH = new gdcm::FileHelper(fileToBuild); fileH->SetImageData(imageData,size); // Step 4 : Set the writting mode and write the image - std::cout << "4..."; fileH->SetWriteModeToRaw(); + std::cout << "4'..."; switch (img.writeMode) { case 'a' : // Write an ACR file @@ -224,45 +269,65 @@ int WriteSimple(Image &img) break; default : - std::cout << "Failed\n" + std::cout << "Failed for [" << fileName.str() << "]\n" << " Write mode '"<Delete(); delete[] imageData; return 1; } + std::cout << std::endl; + //fileToBuild->Print(); + if( !fileH->Write(fileName.str()) ) { - std::cout << "Failed\n" - << " File in unwrittable\n"; + std::cout << "Failed for [" << fileName.str() << "]\n" + << " File is unwrittable" << std::endl; + + fileH->Delete(); - delete fileH; - delete fileToBuild; delete[] imageData; return 1; } // Step 5 : Read the written image - std::cout << "5..."; - gdcm::FileHelper *reread = new gdcm::FileHelper( fileName.str() ); - //gdcm::FileHelper *reread = new gdcm::FileHelper( ); + std::cout << "5..." << std::endl; + // old form. + //GDCM_NAME_SPACE::FileHelper *reread = new GDCM_NAME_SPACE::FileHelper( fileName.str() ); + // Better use : + GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New( ); + f->SetLoadMode(GDCM_NAME_SPACE::LD_ALL); + f->SetFileName( fileName.str() ); + //reread->SetFileName( fileName.str() ); - //reread->Load(); + //reread->SetLoadMode(GDCM_NAME_SPACE::LD_ALL); // Load everything + // Possible values are + // GDCM_NAME_SPACE::LD_ALL, + // GDCM_NAME_SPACE::LD_NOSEQ, + // GDCM_NAME_SPACE::LD_NOSHADOW, + // GDCM_NAME_SPACE::LD_NOSEQ|GDCM_NAME_SPACE::LD_NOSHADOW, + // GDCM_NAME_SPACE::LD_NOSHADOWSEQ + + f->Load(); + // reread->Load(); + GDCM_NAME_SPACE::FileHelper *reread = GDCM_NAME_SPACE::FileHelper::New( f ); + + //reread->Print(); + if( !reread->GetFile()->IsReadable() ) { std::cerr << "Failed" << std::endl - << "Could not read written image : " << fileName << std::endl; - delete fileToBuild; - delete fileH; - delete reread; + << "Could not read written image : " << fileName.str() << std::endl; + fileToBuild->Delete(); + fileH->Delete(); + reread->Delete(); delete[] imageData; return 1; } // Step 6 : Compare to the written image - std::cout << "6..."; + std::cout << "6..." << std::endl; size_t dataSizeWritten = reread->GetImageDataSize(); uint8_t *imageDataWritten = reread->GetImageData(); @@ -273,20 +338,20 @@ int WriteSimple(Image &img) << " File type differ: " << fileH->GetWriteType() << " # " << reread->GetFile()->GetFileType() << std::endl; - delete fileToBuild; - delete fileH; - delete reread; + fileToBuild->Delete(); + fileH->Delete(); + reread->Delete(); delete[] imageData; return 1; } - + std::cout << "6.1..." << std::endl; // Test the image size if (fileToBuild->GetXSize() != reread->GetFile()->GetXSize() || fileToBuild->GetYSize() != reread->GetFile()->GetYSize() || fileToBuild->GetZSize() != reread->GetFile()->GetZSize()) { - std::cout << "Failed" << std::endl + std::cout << "Failed for [" << fileName.str() << "]" << std::endl << " X Size differs: " << "X: " << fileToBuild->GetXSize() << " # " << reread->GetFile()->GetXSize() << " | " @@ -294,48 +359,61 @@ int WriteSimple(Image &img) << reread->GetFile()->GetYSize() << " | " << "Z: " << fileToBuild->GetZSize() << " # " << reread->GetFile()->GetZSize() << std::endl; - delete fileToBuild; - delete fileH; - delete reread; + fileToBuild->Delete(); + fileH->Delete(); + reread->Delete(); delete[] imageData; return 1; } - + std::cout << "6.2..." << std::endl; // Test the data size if (size != dataSizeWritten) { std::cout << "Failed" << std::endl << " Pixel areas lengths differ: " << size << " # " << dataSizeWritten << std::endl; - delete fileToBuild; - delete fileH; - delete reread; + fileToBuild->Delete(); + fileH->Delete(); + reread->Delete(); delete[] imageData; return 1; } - - // Test the data's content +std::cout << "6.3..." << std::endl; + // Test the data content if ( memcmp(imageData, imageDataWritten, size) !=0 ) { - std::cout << "Failed" << std::endl + std::cout << fileName.str() << " Failed " << std::endl << " Pixel differ (as expanded in memory)." << std::endl; - delete fileToBuild; - delete fileH; - delete reread; + std::cout << " list of the first " << MAX_NUMBER_OF_DIFFERENCE + << " pixels differing (pos : test - ref) :" + << std::endl; + unsigned int i; + unsigned int j; + for(i=0, j=0;iDelete(); + fileH->Delete(); + reread->Delete(); delete[] imageData; return 1; } - - std::cout << "OK" << std::endl; - - delete fileToBuild; - delete fileH; - delete reread; +std::cout << "======================= OK" << std::endl; + fileH->Delete(); + reread->Delete(); delete[] imageData; - return 0; } @@ -349,15 +427,16 @@ int TestWriteSimple(int argc, char *argv[]) return 1; } - gdcm::Debug::DebugOn(); - + // GDCM_NAME_SPACE::Debug::DebugOn(); + int ret=0; int i=0; while( Images[i].sizeX>0 && Images[i].sizeY>0 ) { - std::cout << "Test n :" << i; + std::cout << "Test n :" << i <