testDicomReader.cpp File Reference

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

Include dependency graph for testDicomReader.cpp:

Go to the source code of this file.

Functions
template<typename Image3D >
bool	testDicomReader ()
static std::vector< std::string >	getFirstDicomSerieFileNames (const std::string &path)
template<typename Image3D >
bool	testDicomReaderFromDirectory ()
template<typename PixelType >
bool	testSpatialInformation ()
bool	testIOException ()
int	main (int argc, char **argv)

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author

Adrien Krähenbühl (adrie.nosp@m.n.kr.nosp@m.ahenb.nosp@m.uhl@.nosp@m.loria.nosp@m..fr ) LORIA (CNRS, UMR 7503), Université de Lorraine, France

Date

2013/10/14

Functions for testing class DicomReader.

This file is part of the DGtal library.

Definition in file testDicomReader.cpp.

Function Documentation

getFirstDicomSerieFileNames()

static std::vector<std::string> getFirstDicomSerieFileNames ( const std::string &  path )

static

Definition at line 75 of file testDicomReader.cpp.

{
  typedef itk::GDCMSeriesFileNames NamesGeneratorType;
  NamesGeneratorType::Pointer nameGenerator = NamesGeneratorType::New();
  nameGenerator->SetUseSeriesDetails( true );
  nameGenerator->SetDirectory( path );
  
  typedef itk::GDCMSeriesFileNames::SeriesUIDContainerType SeriesIdContainer;
  const SeriesIdContainer & seriesUID = nameGenerator->GetSeriesUIDs();
  
  if (! seriesUID.empty() )
  {
    return nameGenerator->GetFileNames( *(seriesUID.begin()) );
  }
  return std::vector<std::string>();
}

Referenced by testDicomReaderFromDirectory(), and testSpatialInformation().

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 164 of file testDicomReader.cpp.

{
  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<ImageContainerBySTLVector<Z3i::Domain, unsigned char> >()
    && testDicomReader<ImageContainerBySTLVector<Z3i::Domain, uint16_t> >()
    && testDicomReader<ImageContainerByITKImage<Z3i::Domain, unsigned char> >()
    && testDicomReader<ImageContainerByITKImage<Z3i::Domain, uint16_t> >()
    && testDicomReaderFromDirectory<ImageContainerBySTLVector<Z3i::Domain, unsigned char> >()
    && testDicomReaderFromDirectory<ImageContainerByITKImage<Z3i::Domain, uint16_t> >()
    && testSpatialInformation<unsigned char>()
    && testSpatialInformation<uint16_t>()
    && testIOException();
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
 
  return !res;
}

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), DGtal::Trace::info(), testDicomReader(), testDicomReaderFromDirectory(), testIOException(), and DGtal::trace.

testDicomReader()

template<typename Image3D >

bool testDicomReader

(

)

Example of a test. To be completed.

Definition at line 53 of file testDicomReader.cpp.

{
  const std::string filename = testPath + "samples/dicomSample/1629.dcm";
  Image3D image = DicomReader< Image3D >::importDicom( filename );
 
  trace.info() << image <<endl;
 
  unsigned int nbVal=0, nbPos = 0;
  typename Image3D::ConstRange r = image.constRange();
  for ( typename Image3D::ConstRange::ConstIterator it=r.begin(), itend=r.end() ; it != itend ; ++it )
  {
    nbVal++;
    if ( (*it) > 0 ) nbPos++;
  }
 
  trace.info() << "Number of points with (val>0) = " << nbVal << endl;
 
  return nbVal==2130048 && nbPos==296030;
}

References DGtal::Trace::info(), and DGtal::trace.

Referenced by main().

testDicomReaderFromDirectory()

template<typename Image3D >

bool testDicomReaderFromDirectory

(

)

Definition at line 93 of file testDicomReader.cpp.

{
  const std::string path = testPath + "samples/dicomSample";
  const std::vector<std::string> fileNames = getFirstDicomSerieFileNames( path );
  Image3D image = DicomReader< Image3D >::importDicomSeries( fileNames );
 
  trace.info() << image <<endl;
 
  unsigned int nbVal=0, nbPos = 0;
  typename Image3D::ConstRange r = image.constRange();
  for ( typename Image3D::ConstRange::ConstIterator it=r.begin(), itend=r.end() ; it != itend ; ++it )
  {
    nbVal++;
    if ( (*it) > 0 ) nbPos++;
  }
 
  trace.info() << "Number of points with (val>0) = " << nbVal << endl;
 
  return nbVal==2130048 && nbPos==296030;
}

References getFirstDicomSerieFileNames(), DGtal::Trace::info(), and DGtal::trace.

Referenced by main().

testIOException()

bool testIOException

(

)

Definition at line 142 of file testDicomReader.cpp.

{
  //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;
}

References DGtal::Trace::info(), and DGtal::trace.

Referenced by main().

testSpatialInformation()

template<typename PixelType >

bool testSpatialInformation

(

)

Definition at line 116 of file testDicomReader.cpp.

{
  const std::string path = testPath + "samples/dicomSample";
 
  const std::vector<std::string> fileNames = getFirstDicomSerieFileNames( path );
 
  typedef ImageContainerByITKImage<Z3i::Domain,  PixelType> Image3D;
  Image3D image = DicomReader< Image3D >::importDicomSeries( fileNames );
  typename Image3D::ITKImagePointer dgtal_itk = image.getITKImagePointer();
 
  typedef itk::Image<PixelType, 3> ItkImage;
  typedef itk::ImageSeriesReader<ItkImage> ItkReader;
  typename ItkReader::Pointer reader = ItkReader::New();
  reader->SetFileNames( fileNames );
 
  reader->Update();
  
  typename ItkImage::Pointer itk = reader->GetOutput();
 
  const bool ok1 = ( dgtal_itk->GetSpacing() == itk->GetSpacing() );
  const bool ok2 = ( dgtal_itk->GetOrigin() == itk->GetOrigin() );
  const bool ok3 = ( dgtal_itk->GetDirection() == itk->GetDirection() );
  return ok1 && ok2 && ok3;
}

References getFirstDicomSerieFileNames().