Program: gdcm
Module: $RCSfile: TestAllReadCompareDicom.cxx,v $
Language: C++
- Date: $Date: 2004/11/24 10:23:46 $
- Version: $Revision: 1.17 $
+ Date: $Date: 2005/01/11 15:58:36 $
+ Version: $Revision: 1.23 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
=========================================================================*/
#include "gdcmHeader.h"
#include "gdcmFile.h"
-#include <stdio.h>
+
+#include <iostream>
+#include <fstream>
//Generated file:
#include "gdcmDataImages.h"
{
std::cout << " Testing: " << filename << std::endl;
+ ////// Step 1:
+ std::cout << " 1...";
gdcm::File* tested = new gdcm::File( filename );
if( !tested->GetHeader()->IsReadable() )
{
- std::cout << " Image not gdcm compatible:"
+ 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" );
+ //FILE* testFILE = fopen( referenceFileName.c_str(), "r" );
+ std::ifstream testFILE( referenceFileName.c_str() );
if (! testFILE )
{
- ////// Step 3a:
uint8_t* testedImageData = tested->GetImageData(); // Kludge
(void)testedImageData;
+ tested->SetWriteModeToRGB();
tested->WriteDcmExplVR( referenceFileName );
- std::cerr << " Creating reference baseline file :" << std::endl
- << " " << referenceFileName
- << std::endl;
- delete tested;
- //delete (char*)testedImageData;
- return 0;
}
else
{
- fclose( testFILE );
+ //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: reference image " << std::endl
- << " " << referenceFileName <<std::endl
- << " is not gdcm compatible." << std::endl;
+ std::cout << " Failed" << std::endl
+ << " reference image "
+ << referenceFileName
+ << " is not gdcm compatible." << std::endl;
delete tested;
delete reference;
return 1;
}
- ////// Step 3b:
+ 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 areas lengths differ: "
+ << " pixel ("
+ << PixelType
+ <<") areas lengths differ: "
<< testedDataSize << " # " << referenceDataSize
<< std::endl;
delete tested;
return 1;
}
+ // Test the data content
if (int res = memcmp(testedImageData, referenceImageData,
testedDataSize) != 0 )
{
(void)res;
std::cout << " Failed" << std::endl
- << " pixel differ (as expanded in memory)."
+ << " pixel ("
+ << PixelType
+ << ") differ (as expanded in memory)."
<< std::endl;
delete tested;
delete reference;
delete tested;
delete reference;
- std::cout << " Passed..." << std::endl;
+ std::cout << "OK." << std::endl;
return 0;
}