X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Testing%2FTestAllReadCompareDicom.cxx;h=7a7ab5c6f4f79ba630bb947d523798702d42d610;hb=ead1ddd3fcdebedc60b1f24c59a1892a7a2da2df;hp=6f6748410935582f53b1a8ab508d3e4d9b443b07;hpb=0ced3dcc91cb31e52f9756b8bd1f62a3690f883c;p=gdcm.git diff --git a/Testing/TestAllReadCompareDicom.cxx b/Testing/TestAllReadCompareDicom.cxx index 6f674841..7a7ab5c6 100644 --- a/Testing/TestAllReadCompareDicom.cxx +++ b/Testing/TestAllReadCompareDicom.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: TestAllReadCompareDicom.cxx,v $ Language: C++ - Date: $Date: 2005/03/30 15:30:33 $ - Version: $Revision: 1.31 $ + Date: $Date: 2005/08/30 15:13:07 $ + Version: $Revision: 1.47 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -18,9 +18,11 @@ #include "gdcmDirList.h" #include "gdcmFile.h" #include "gdcmFileHelper.h" +#include "gdcmGlobal.h" +#include "gdcmTS.h" +#include "gdcmDebug.h" #include -#include //Generated file: #include "gdcmDataImages.h" @@ -43,16 +45,16 @@ public: ~TestFile(void); bool IsReadable(void) {return readable;} - unsigned int GetXSize(void) {return sizeX;}; - void SetXSize(unsigned int size) {sizeX = size;}; - unsigned int GetYSize(void) {return sizeY;}; - void SetYSize(unsigned int size) {sizeY = size;}; - unsigned int GetZSize(void) {return sizeZ;}; - void SetZSize(unsigned int size) {sizeZ = size;}; - unsigned int GetScalarSize(void) {return scalarSize;}; - void SetScalarSize(unsigned int size) {scalarSize = size;}; - unsigned int GetNumberOfComponents(void) {return components;}; - void SetNumberOfComponents(unsigned int size) {components = size;}; + int GetXSize(void) {return sizeX;}; + void SetXSize(int size) {sizeX = size;}; + int GetYSize(void) {return sizeY;}; + void SetYSize(int size) {sizeY = size;}; + int GetZSize(void) {return sizeZ;}; + void SetZSize(int size) {sizeZ = size;}; + int GetScalarSize(void) {return scalarSize;}; + void SetScalarSize(int size) {scalarSize = size;}; + int GetNumberOfComponents(void) {return components;}; + void SetNumberOfComponents(int size) {components = size;}; unsigned long GetDataSize(void) {return GetLineSize()*sizeY*sizeZ;} uint8_t *GetData(void) {return data;} @@ -75,30 +77,46 @@ private: bool WriteFileHeader(std::ofstream *fp); bool WriteFileData(std::ofstream *fp); - uint8_t ReadInt8 (std::ifstream *fp); - uint16_t ReadInt16(std::ifstream *fp); - uint32_t ReadInt32(std::ifstream *fp); - void WriteInt8 (std::ofstream *fp,uint8_t data); - void WriteInt16(std::ofstream *fp,uint16_t data); - void WriteInt32(std::ofstream *fp,uint32_t data); + uint8_t ReadInt8 (std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) + throw( std::ios::failure ); +#else + ; +#endif + uint16_t ReadInt16(std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) + throw( std::ios::failure ); +#else + ; +#endif + uint32_t ReadInt32(std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) + throw( std::ios::failure ); +#else + ; +#endif + void WriteInt8 (std::ofstream *fp,uint8_t value); + void WriteInt16(std::ofstream *fp,uint16_t value); + void WriteInt32(std::ofstream *fp,uint32_t value); std::string fileName; bool readable; - unsigned int sizeX; - unsigned int sizeY; - unsigned int sizeZ; - unsigned int scalarSize; - unsigned int components; + int sizeX; + int sizeY; + int sizeZ; + uint16_t scalarSize; + uint16_t components; uint8_t *data; int swapCode; static const unsigned int HEADER_SIZE; }; +const unsigned int MAX_NUMBER_OF_DIFFERENCE = 10; const unsigned int TestFile::HEADER_SIZE = 20; -TestFile::TestFile(void) +TestFile::TestFile() { fileName = ""; readable=false; @@ -233,7 +251,7 @@ bool TestFile::ReadFileHeader(std::ifstream *fp) sizeY = ReadInt32(fp); // Size Y sizeZ = ReadInt32(fp); // Size Z scalarSize = ReadInt16(fp)/8; // bits per scalar - components = ReadInt16(fp); // Number of components + components = ReadInt16(fp); // Number of components return(true); } @@ -279,7 +297,7 @@ bool TestFile::WriteFileHeader(std::ofstream *fp) WriteInt32(fp,sizeY); // Size Y WriteInt32(fp,sizeZ); // Size Z WriteInt16(fp,scalarSize*8); // bits per scalar - WriteInt16(fp,components); // number of components + WriteInt16(fp,components); // number of components return(true); } @@ -292,26 +310,34 @@ bool TestFile::WriteFileData(std::ofstream *fp) } uint8_t TestFile::ReadInt8 (std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) throw( std::ios::failure ) +#endif { uint8_t g; fp->read ((char*)&g, (size_t)1); +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) if ( fp->fail() ) throw std::ios::failure( "TestFile::ReadInt8() - file error." ); if( fp->eof() ) throw std::ios::failure( "TestFile::ReadInt8() - EOF." ); +#endif return g; } uint16_t TestFile::ReadInt16(std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) throw( std::ios::failure ) +#endif { uint16_t g; fp->read ((char*)&g, (size_t)2); +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) if ( fp->fail() ) throw std::ios::failure( "TestFile::ReadInt16() - file error." ); if( fp->eof() ) throw std::ios::failure( "TestFile::ReadInt16() - EOF." ); +#endif #if defined(GDCM_WORDS_BIGENDIAN) g = ( g << 8 | g >> 8 ); @@ -320,14 +346,18 @@ uint16_t TestFile::ReadInt16(std::ifstream *fp) } uint32_t TestFile::ReadInt32(std::ifstream *fp) +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) throw( std::ios::failure ) +#endif { uint32_t g; fp->read ((char*)&g, (size_t)4); +#if !(__GNUC__==2 && __GNUC_MINOR__<=96) if ( fp->fail() ) throw std::ios::failure( "TestFile::ReadInt32() - file error." ); if( fp->eof() ) throw std::ios::failure( "TestFile::ReadInt32() - EOF." ); +#endif #if defined(GDCM_WORDS_BIGENDIAN) g = ( (g<<24) | ((g<<8) & 0x00ff0000) | @@ -336,26 +366,26 @@ uint32_t TestFile::ReadInt32(std::ifstream *fp) return g; } -void TestFile::WriteInt8 (std::ofstream *fp,uint8_t data) +void TestFile::WriteInt8 (std::ofstream *fp,uint8_t value) { - fp->write((char*)&data, (size_t)1); + fp->write((char*)&value, (size_t)1); } -void TestFile::WriteInt16(std::ofstream *fp,uint16_t data) +void TestFile::WriteInt16(std::ofstream *fp,uint16_t value) { #if defined(GDCM_WORDS_BIGENDIAN) - data = ( data << 8 | data >> 8 ); + value = ( value << 8 | value >> 8 ); #endif - fp->write((char*)&data, (size_t)2); + fp->write((char*)&value, (size_t)2); } -void TestFile::WriteInt32(std::ofstream *fp,uint32_t data) +void TestFile::WriteInt32(std::ofstream *fp,uint32_t value) { #if defined(GDCM_WORDS_BIGENDIAN) - data = ( (data<<24) | ((data<<8) & 0x00ff0000) | - ( (data>>8) & 0x0000ff00) | (data>>24) ); + value = ( (value<<24) | ((value<<8) & 0x00ff0000) | + ( (value>>8) & 0x0000ff00) | (value>>24) ); #endif - fp->write((char*)&data, (size_t)4); + fp->write((char*)&value, (size_t)4); } int InternalTest(std::string const &filename, @@ -366,26 +396,34 @@ int InternalTest(std::string const &filename, ////// Step 1: std::cout << "1..."; - gdcm::FileHelper *tested = new gdcm::FileHelper( filename ); - if( !tested->GetFile()->IsReadable() ) + + // new style + gdcm::File *f = new gdcm::File(); + f->SetLoadMode ( gdcm::LD_ALL ); // Load everything + f->SetFileName( filename ); + f->Load(); + + if( !f->IsReadable() ) { std::cout << " Failed" << std::endl << " Image not gdcm compatible:" << filename << std::endl; - delete tested; + delete f; return 1; } - + gdcm::FileHelper *tested = new gdcm::FileHelper( f ); + ////// Step 2: ////// Check for existence of reference baseline dicom file: std::cout << "2..."; TestFile *reference = new TestFile(); - reference->Load(referenceFileName); - if(!reference->IsReadable()) + std::ifstream refFile(referenceFileName.c_str(), + std::ios::binary|std::ios::in); + if(!refFile) { std::cout << " Failed" << std::endl - << " Image not BMP compatible:" + << " Image not found:" << referenceFileName << std::endl; reference->SetXSize(tested->GetFile()->GetXSize()); reference->SetYSize(tested->GetFile()->GetYSize()); @@ -395,6 +433,8 @@ int InternalTest(std::string const &filename, reference->SetData(tested->GetImageData()); reference->Write(referenceFileName); } + else + refFile.close(); reference->Load(referenceFileName); if(!reference->IsReadable()) @@ -404,6 +444,7 @@ int InternalTest(std::string const &filename, << filename << std::endl; delete reference; delete tested; + delete f; return 1; } @@ -431,6 +472,7 @@ int InternalTest(std::string const &filename, << reference->GetZSize() << std::endl; delete reference; delete tested; + delete f; return 1; } @@ -446,7 +488,7 @@ int InternalTest(std::string const &filename, << reference->GetNumberOfComponents() << std::endl; delete reference; delete tested; - return 1; + delete f; } // Test the data size @@ -465,6 +507,7 @@ int InternalTest(std::string const &filename, << std::endl; delete tested; delete reference; + delete f; return 1; } @@ -473,19 +516,44 @@ int InternalTest(std::string const &filename, testedDataSize) != 0 ) { (void)res; + std::string ts = tested->GetFile()->GetTransferSyntax(); + std::cout << " Failed" << std::endl << " pixel (" << PixelType << ") differ (as expanded in memory)." + << std::endl + << " compression : " + << gdcm::Global::GetTS()->GetValue(ts) << std::endl; + + std::cout << " list of the first " << MAX_NUMBER_OF_DIFFERENCE + << " pixels differing (pos : test - ref) :" << std::endl; + int i; + unsigned int j; + for(i=0, j=0;i= 3 ) { // The test is specified a specific filename, use it instead of looping // over all images @@ -502,7 +573,7 @@ int TestAllReadCompareDicom(int argc, char *argv[]) const std::string reference = argv[2]; return InternalTest( input, reference ); } - else if ( argc > 3 || argc == 2 ) + else if ( argc > 4 || argc == 2 ) { std::cerr << " Usage: " << argv[0] << " (no arguments needed)." << std::endl; @@ -522,23 +593,18 @@ int TestAllReadCompareDicom(int argc, char *argv[]) std::cout << " step 1: parse the image (as gdcmFile) and call" << " IsReadable(). " << std::endl; - std::cout << " step 2: find in GDCM_DATA_ROOT/BaselineDicom/filename.dcm" - << std::endl - << " (with format DICOM V3, explicit Value" - << "Representation)" - << std::endl; - std::cout << " step 3a: when image NOT found on step 2, write " + std::cout << " step 2: find in GDCM_DATA_ROOT/BaselineDicom/filename.tst" << std::endl - << " GDCM_DATA_ROOT/BaselineDicom/filename.dcm" + << " special internal file format containing the" << std::endl - << " (with format DICOM V3, explicit Value" - << "Representation)" + << " caracteristic of the image and the pixel datas " + << "(uncompressed). This file is written if it's not found." << std::endl; - std::cout << " step 3b: when image found on step 2, and when IsReadable()" + std::cout << " step 3: compare the DICOM image with the reference image" << std::endl - << " compare it (in memory with memcmp) with the" + << " (.tst file). The test is made on the caracteristics" << std::endl - << " image we are testing (the one of step 1). " + << " of the image and the pixel datas" << std::endl << std::endl; int i = 0; @@ -548,7 +614,7 @@ int TestAllReadCompareDicom(int argc, char *argv[]) ////// Check for existence of reference baseline directory std::string baseLineDir = GDCM_DATA_ROOT; - baseLineDir += "/BaselineDicom/"; + baseLineDir += "/BaselineDicom"; if( !gdcm::DirList::IsDirectory(baseLineDir) ) { @@ -565,6 +631,7 @@ int TestAllReadCompareDicom(int argc, char *argv[]) filename += "/"; filename += gdcmDataImages[i]; + baseLineDir += '/'; std::string referenceFileName = baseLineDir + gdcmDataImages[i++]; std::string::size_type slash_pos = referenceFileName.rfind( "." ); if( slash_pos != std::string::npos )