#include "gdcmImageReader.h"
#include "gdcmImage.h"
#include "gdcmWriter.h"
#include "gdcmAttribute.h"
#include "gdcmSystem.h"
#include <iostream>
#include <fstream>
int main(int argc, char *argv[])
{
if( argc < 3 )
{
std::cerr << argv[0] << " input1.dcm input2.dcm" << std::endl;
return 1;
}
const char *filename1 = argv[1];
const char *filename2 = argv[2];
gdcm::ImageReader reader1;
reader1.SetFileName( filename1 );
if( !reader1.Read() )
{
std::cerr << "Could not read: " << filename1 << std::endl;
return 1;
}
gdcm::ImageReader reader2;
reader2.SetFileName( filename2 );
if( !reader2.Read() )
{
std::cerr << "Could not read: " << filename2 << std::endl;
return 1;
}
std::cout << "Both files can be read and looks like DICOM" << std::endl;
size_t s1 = gdcm::System::FileSize(filename1);
size_t s2 = gdcm::System::FileSize(filename2);
if( s1 != s2 )
{
std::cout << "Size mismatch: " << s1 << " != " << s2 << std::endl;
return 1;
}
else
{
std::cout << "Size match: " << s1 << " = " << s2 << std::endl;
}
std::ifstream is1( filename1 );
char *buffer1 = new char[s1];
is1.read(buffer1, s1);
std::ifstream is2( filename2 );
char *buffer2 = new char[s2];
is2.read(buffer2, s2);
assert( s1 == s2 );
if( memcmp(buffer1, buffer2, s1 ) == 0 )
{
std::cout << "memcmp succeed ! File are bit identical" << std::endl;
}
else
{
std::cout << "memcmp failed!" << std::endl;
}
std::multiset<char> set1( buffer1, buffer1 + s1 );
std::multiset<char> set2( buffer2, buffer2 + s2 );
if( set1 == set2 )
{
std::cout << "set1 == set2. Byte histogram seems valid" << std::endl;
}
else
{
std::cout << "set1 != set2" << std::endl;
}
delete[] buffer1;
delete[] buffer2;
return 0;
}