+int InternalTest(std::string const & filename,
+ std::string const & referenceFileName )
+{
+ std::cout << " Testing: " << filename << std::endl;
+
+ ////// Step 1:
+ std::cout << " 1...";
+ gdcm::File* tested = new gdcm::File( filename );
+ if( !tested->GetHeader()->IsReadable() )
+ {
+ std::cout << " Failed" << std::endl
+ << " Image not gdcm compatible:"
+ << filename << std::endl;
+ delete tested;
+ return 1;
+ }
+
+ ////// Step 2:
+ ////// Check for existence of reference baseline dicom file:
+ std::cout << "2...";
+
+ //FILE* testFILE = fopen( referenceFileName.c_str(), "r" );
+ std::ifstream testFILE( referenceFileName.c_str() );
+ if (! testFILE )
+ {
+ uint8_t* testedImageData = tested->GetImageData(); // Kludge
+ (void)testedImageData;
+
+ tested->SetWriteModeToRGB();
+ tested->WriteDcmExplVR( referenceFileName );
+ }
+ else
+ {
+ //fclose( testFILE );
+ }
+ testFILE.close();
+
+ ////// Step 3a:
+ ////// When reference file is not gdcm readable test is failed:
+ std::cout << "3a...";
+
+ gdcm::File* reference = new gdcm::File( referenceFileName );
+ if( !reference->GetHeader()->IsReadable() )
+ {
+ std::cout << " Failed" << std::endl
+ << " reference image "
+ << referenceFileName
+ << " is not gdcm compatible." << std::endl;
+ delete tested;
+ delete reference;
+ return 1;
+ }
+
+ std::string PixelType = reference->GetHeader()->GetPixelType();
+
+ ////// Step 3b:
+ std::cout << "3b...";
+ int testedDataSize = tested->GetImageDataSize();
+ uint8_t* testedImageData = tested->GetImageData();
+
+ int referenceDataSize = reference->GetImageDataSize();
+ uint8_t* referenceImageData = reference->GetImageData();
+
+ // Test the image size
+ if (tested->GetHeader()->GetXSize() != reference->GetHeader()->GetXSize() ||
+ tested->GetHeader()->GetYSize() != reference->GetHeader()->GetYSize() ||
+ tested->GetHeader()->GetZSize() != reference->GetHeader()->GetZSize())
+ {
+ std::cout << "Failed" << std::endl
+ << " Size differs: "
+ << "X: " << tested->GetHeader()->GetXSize() << " # "
+ << reference->GetHeader()->GetXSize() << " | "
+ << "Y: " << tested->GetHeader()->GetYSize() << " # "
+ << reference->GetHeader()->GetYSize() << " | "
+ << "Z: " << tested->GetHeader()->GetZSize() << " # "
+ << reference->GetHeader()->GetZSize() << std::endl;
+ delete reference;
+ delete tested;
+ return 1;
+ }
+
+ // Test the data size
+ if (testedDataSize != referenceDataSize)
+ {
+ std::cout << " Failed" << std::endl
+ << " pixel ("
+ << PixelType
+ <<") areas lengths differ: "
+ << testedDataSize << " # " << referenceDataSize
+ << std::endl;
+ delete tested;
+ delete reference;
+ return 1;
+ }
+
+ // Test the data content
+ if (int res = memcmp(testedImageData, referenceImageData,
+ testedDataSize) != 0 )
+ {
+ (void)res;
+ std::cout << " Failed" << std::endl
+ << " pixel ("
+ << PixelType
+ << ") differ (as expanded in memory)."
+ << std::endl;
+ delete tested;
+ delete reference;
+ return 1;
+ }
+
+ //////////////// Clean up:
+ delete tested;
+ delete reference;
+
+ std::cout << "OK." << std::endl;
+
+ return 0;
+}
+