+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" );
+ if (! testFILE )
+ {
+ uint8_t* testedImageData = tested->GetImageData(); // Kludge
+ (void)testedImageData;
+
+ tested->SetWriteModeToRGB();
+ tested->WriteDcmExplVR( referenceFileName );
+ }
+ else
+ {
+ fclose( testFILE );
+ }
+
+ ////// 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;
+ }
+
+ ////// Step 3b:
+ std::cout << "3b...";
+ int testedDataSize = tested->GetImageDataSize();
+ uint8_t* testedImageData = tested->GetImageData();
+
+ int referenceDataSize = reference->GetImageDataSize();
+ uint8_t* referenceImageData = reference->GetImageData();
+
+ if (testedDataSize != referenceDataSize)
+ {
+ std::cout << " Failed" << std::endl
+ << " pixel areas lengths differ: "
+ << testedDataSize << " # " << referenceDataSize
+ << std::endl;
+ delete tested;
+ delete reference;
+ return 1;
+ }
+
+ if (int res = memcmp(testedImageData, referenceImageData,
+ testedDataSize) != 0 )
+ {
+ (void)res;
+ std::cout << " Failed" << std::endl
+ << " pixel 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;
+}
+