testDicomReader.cpp

Go to the documentation of this file.

#include <iostream>
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/io/readers/DicomReader.h"

#include "ConfigTest.h"


using namespace std;
using namespace DGtal;

// Functions for testing class DicomReader.

bool testDicomReader()
{
  //Default image selector = STLVector
  typedef ImageContainerBySTLVector<Z3i::Domain,  unsigned char > Image3D;

  std::string filename = testPath + "samples/dicomSample/1629.dcm";
  Image3D image = DicomReader< Image3D >::importDicom( filename );

  trace.info() << image <<endl;

  unsigned int nbVal=0, nbPos = 0;
  for ( Image3D::ConstIterator it=image.begin(), itend=image.end() ; it != itend ; ++it )
  {
    nbVal++;
          if ( (*it) > 0 ) nbPos++;
  }

  trace.info() << "Number of points with (val>0) = " << nbVal << endl;

  return nbVal==2130048 && nbPos==296030;
}


bool testIOException()
{
  //Default image selector = STLVector
  typedef ImageContainerBySTLVector<Z3i::Domain,  int > Image3D;

  std::string filename = testPath + "samples/null.dcm";
  try {
          Image3D image = DicomReader< Image3D >::importDicom( filename );
  }
  catch(exception& e) {
          trace.info() << "Exception catched. Message : " << e.what()<<endl;
    return true;
  }

  return false;
}

// Standard services - public :

int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class DicomReader" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
        trace.info() << " " << argv[ i ];
  trace.info() << endl;

  bool res = testDicomReader() && testIOException();
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();

  return !res;
}
//                                                                           //