From 60ccac3fe785b813cca84c15814e8adc92e46116 Mon Sep 17 00:00:00 2001 From: jpr Date: Thu, 30 Aug 2007 11:29:03 +0000 Subject: [PATCH] Add a JPEG200-writter test Right now, it fails on 12 images (amongst 150) --- Testing/TestReadWriteJPEG2000ReadCompare.cxx | 286 +++++++++++++++++++ 1 file changed, 286 insertions(+) create mode 100755 Testing/TestReadWriteJPEG2000ReadCompare.cxx diff --git a/Testing/TestReadWriteJPEG2000ReadCompare.cxx b/Testing/TestReadWriteJPEG2000ReadCompare.cxx new file mode 100755 index 00000000..9a097e0d --- /dev/null +++ b/Testing/TestReadWriteJPEG2000ReadCompare.cxx @@ -0,0 +1,286 @@ +/*========================================================================= + + Program: gdcm + Module: $RCSfile: TestReadWriteJPEG2000ReadCompare.cxx,v $ + Language: C++ + Date: $Date: 2007/08/30 11:29:03 $ + Version: $Revision: 1.1 $ + + Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de + l'Image). All rights reserved. See Doc/License.txt or + http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details. + + This software is distributed WITHOUT ANY WARRANTY; without even + the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR + PURPOSE. See the above copyright notices for more information. + +=========================================================================*/ +#include "gdcmFile.h" +#include "gdcmFileHelper.h" +#include "gdcmDebug.h" + +#include "gdcmGlobal.h" +#include "gdcmTS.h" + +//Generated file: +#include "gdcmDataImages.h" + + +const unsigned int MAX_NUMBER_OF_DIFFERENCE = 10; + +int nb_of_success2000___; +int nb_of_failure2000___; +int nb_of_diffPM12000___; +static int CompareInternalJPEG2000(std::string const &filename, std::string const &output) +{ + std::cout << "----------------------------------------------------------------------" << std::endl + << " Testing: " << filename << std::endl; + + //////////////// Step 1 (see above description): + + GDCM_NAME_SPACE::File *file = GDCM_NAME_SPACE::File::New( ); + file->SetFileName( filename ); + file->Load (); + if( !file->IsReadable() ) + { + std::cout << "Failed" << std::endl + << "Test::TestReadWriteJPEG2000ReadCompare: Image not gdcm compatible:" + << filename << std::endl; + file->Delete(); + nb_of_failure2000___++; + return 1; + } + std::cout << " step 1..."; + + //////////////// Step 2: + GDCM_NAME_SPACE::FileHelper *filehelper = GDCM_NAME_SPACE::FileHelper::New( file ); + int dataSize = filehelper->GetImageDataRawSize(); + uint8_t *imageData = filehelper->GetImageDataRaw(); //EXTREMELY IMPORTANT + // Sure, it is : It's up to the user to decide if he wants to + // GetImageData or if he wants to GetImageDataRaw + // (even if we do it by setting a flag, *he* will have to decide) + + //filehelper->SetImageData(imageData, dataSize); + + filehelper->SetContentType(GDCM_NAME_SPACE::UNMODIFIED_PIXELS_IMAGE); // lossless compression : pixels reain unimpared + filehelper->SetWriteModeToRaw(); + filehelper->SetWriteTypeToJPEG2000( ); + + filehelper->SetUserData(imageData,dataSize); // This one ensures the compression + filehelper->Write( output ); + + std::cout << "2..."; + + //////////////// Step 3: + GDCM_NAME_SPACE::File *fileout = GDCM_NAME_SPACE::File::New(); + fileout->SetFileName( output ); + fileout->Load(); + + if( !fileout->IsReadable() ) + { + std::cout << "Failed" << std::endl + << "Test::TestReadWriteJPEG2000ReadCompare: Could not parse the newly " + << "written image:" << filename << std::endl; + file->Delete(); + filehelper->Delete(); + fileout->Delete(); + nb_of_failure2000___++; + return 1; + } + + GDCM_NAME_SPACE::FileHelper *reread = GDCM_NAME_SPACE::FileHelper::New( fileout ); + + std::cout << "3..."; + // For the next step: + int dataSizeWritten = reread->GetImageDataRawSize(); + uint8_t *imageDataWritten = reread->GetImageDataRaw(); + + //////////////// Step 4: + // Test the image size + if (file->GetXSize() != reread->GetFile()->GetXSize() || + file->GetYSize() != reread->GetFile()->GetYSize() || + file->GetZSize() != reread->GetFile()->GetZSize()) + { + std::cout << "Failed" << std::endl + << " Size differs: " + << "X: " << file->GetXSize() << " # " + << reread->GetFile()->GetXSize() << " | " + << "Y: " << file->GetYSize() << " # " + << reread->GetFile()->GetYSize() << " | " + << "Z: " << file->GetZSize() << " # " + << reread->GetFile()->GetZSize() << std::endl; + file->Delete(); + filehelper->Delete(); + fileout->Delete(); + reread->Delete(); + nb_of_failure2000___++; + return 1; + } + + // Test the data size + // beware of odd length Pixel Element! + int dataSizeFixed = dataSize + dataSize%2; + int dataSizeWrittenFixed = dataSizeWritten + dataSizeWritten%2; + + if (dataSizeFixed != dataSizeWrittenFixed) + { + std::cout << "Failed" << std::endl + << " Pixel areas lengths differ: " + << dataSize << " # " << dataSizeWritten << std::endl; + file->Delete(); + filehelper->Delete(); + fileout->Delete(); + reread->Delete(); + nb_of_failure2000___++; + return 1; + } + + // Test the data content + + if (memcmp(imageData, imageDataWritten, dataSizeFixed) !=0) + { + std::string PixelType = filehelper->GetFile()->GetPixelType(); + std::string ts = filehelper->GetFile()->GetTransferSyntax(); + + unsigned int j =0; + for(int i1=0; i1 2) + { + j=1; + break; + } + } + + if (j!=0) + { + std::cout << " Failed" << std::endl + << " pixel (" + << PixelType + << ") differ (as expanded in memory)." + << std::endl + << " compression : " + << GDCM_NAME_SPACE::Global::GetTS()->GetValue(ts) << std::endl; + + std::cout << " list of the first " << MAX_NUMBER_OF_DIFFERENCE + << " pixels differing (pos : original - written) :" + << std::endl; + + for(int i=0, j=0; i 2) + { + std::cout << std::hex << "(" << i << " : " + << std::hex << (int)(imageData[i]) << " - " + << std::hex << (int)(imageDataWritten[i]) << ") " + << std::dec; + ++j; + } + } + std::cout << std::endl; + + file->Delete(); + filehelper->Delete(); + fileout->Delete(); + reread->Delete(); + nb_of_failure2000___++; + return 1; + } + else + { + std::cout << " some pixels" << std::endl + << " (" + << PixelType + << ") differ +/-1 (as expanded in memory)." + << std::endl + << " compression : " + << GDCM_NAME_SPACE::Global::GetTS()->GetValue(ts) << std::endl; + + std::cout << " list of the first " << MAX_NUMBER_OF_DIFFERENCE + << " pixels differing (pos : original - written) :" + << std::endl; + + for(int i=0, j=0; iDelete(); + filehelper->Delete(); + fileout->Delete(); + reread->Delete(); + + return 0; +} + +// ------------------------------------------------------------------------------------------- + +int TestReadWriteJPEG2000ReadCompare(int argc, char *argv[]) +{ + int result = 0; + nb_of_success2000___ = 0; + nb_of_failure2000___ = 0; + nb_of_diffPM12000___ = 0; + + if (argc == 4) + GDCM_NAME_SPACE::Debug::DebugOn(); + + if (argc >= 3) + { + const std::string input = argv[1]; + const std::string output = argv[2]; + result += CompareInternalJPEG2000(input, output); + } + else if( argc > 4 || argc == 2 ) + { + std::cout << "Please read the manual" << std::endl; + } + else + { + std::cout<< "Test::TestReadWriteJPEG2000ReadCompare: description " << std::endl; + std::cout << " For all images in gdcmData (and not blacklisted in " + "Test/CMakeLists.txt)" << std::endl; + std::cout << " apply the following multistep test: " << std::endl; + std::cout << " step 1: parse the image (as gdcmFile) and call" + << " IsReadable(). " << std::endl; + std::cout << " step 2: write the corresponding image in JPEG2000 DICOM V3 " + << "with explicit Value Representation " << std::endl + << " in temporary file " + << "TestReadWriteJPEG2000ReadCompare.dcm." << std::endl; + std::cout << " step 3: read the image written on step2 and call " + << " IsReadable(). " << std::endl; + std::cout << " step 4: compare (in memory with memcmp) that the two " + << "images " << std::endl + << " match (as expanded by gdcm)." << std::endl; + + int i = 0; + while( gdcmDataImages[i] != 0 ) + { + std::string filename = GDCM_DATA_ROOT; + filename += "/"; + filename += gdcmDataImages[i++]; + result += CompareInternalJPEG2000(filename, "TestReadWriteJPEG2000ReadCompare.dcm"); + } + } + std::cout << "==================================" << std::endl; + std::cout << "nb of success " << nb_of_success2000___ << std::endl; + std::cout << "nb of failure " << nb_of_failure2000___ << std::endl; + std::cout << "nb of diff+/-1 " << nb_of_diffPM12000___ << std::endl; + return result; +} -- 2.45.1