testDistanceTransformationND.cpp File Reference

#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/images/ImageSelector.h"
#include "DGtal/geometry/volumes/distance/ExactPredicateLpSeparableMetric.h"
#include "DGtal/geometry/volumes/distance/DistanceTransformation.h"
#include "DGtal/io/colormaps/HueShadeColorMap.h"
#include "DGtal/io/colormaps/GrayscaleColorMap.h"
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/images/SimpleThresholdForegroundPredicate.h"

Go to the source code of this file.

Functions
bool	testDistanceTransformND ()
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

David Coeurjolly (david.nosp@m..coe.nosp@m.urjol.nosp@m.ly@l.nosp@m.iris..nosp@m.cnrs.nosp@m..fr ) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), CNRS, France

Date

2010/12/24

Functions for testing class DistanceTransformND.

This file is part of the DGtal library.

Definition in file testDistanceTransformationND.cpp.

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 122 of file testDistanceTransformationND.cpp.

{
  trace.beginBlock ( "Testing class DistanceTransformND" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
 
  bool res = testDistanceTransformND(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), DGtal::Trace::info(), testDistanceTransformND(), and DGtal::trace.

testDistanceTransformND()

bool testDistanceTransformND

(

)

Example of a test. To be completed.

Definition at line 56 of file testDistanceTransformationND.cpp.

{
  unsigned int nbok = 0;
  unsigned int nb = 0;
 
  trace.beginBlock ( "Testing dT dim=5 ..." );
 
  typedef SpaceND<5> TSpace;
  typedef TSpace::Point Point;
  typedef HyperRectDomain<TSpace> Domain;
  TSpace::Integer t[5] = {0,0,0,0,0};
  Point a ( t );
  TSpace::Integer t2[5] = {15,15,15,15,15};
  Point b ( t2 );
  TSpace::Integer t3[5] = {3,3,3,3,3};
  Point c ( t3 );
  Point d;
 
  typedef ImageSelector<Domain, unsigned int>::Type Image;
  Domain domain(a,b);
  Image image(domain);
 
  //We create an object image with a signle background point (set to 0)
  for (Image::Iterator it=image.begin(),itend=image.end(); it!=itend; ++it)
    *it = 128;
  image.setValue( c , 0 );
 
  typedef SimpleThresholdForegroundPredicate<Image> Predicate;
  Predicate aPredicate(image,0);
 
  typedef ExactPredicateLpSeparableMetric<TSpace, 2> L2Metric;
  L2Metric l2;
  DistanceTransformation<TSpace,Predicate,L2Metric> dt(&domain,&aPredicate, &l2 );
 
  //We check the result
  bool res=true;
  for(Domain::ConstIterator itDom = domain.begin(), itDomend = domain.end();
      itDom != itDomend; ++itDom)
    {
      //distance from the point to the seed
      d = (*itDom) - c;
      L2Metric::RawValue norm2=0;
      for(Point::Iterator itd=d.begin(), itdend=d.end(); itd!=itdend; ++itd)
        norm2+= (*itd)*(*itd);
 
      if ( dt.metric()->rawDistance( (*itDom), dt.getVoronoiSite(*itDom) ) != norm2)
  {
    trace.error()<<"Error at "<<(*itDom)
                 << ": expected="<<norm2<<" and computed="
                 <<dt.metric()->rawDistance( (*itDom), dt.getVoronoiSite(*itDom) )
                 <<endl;
  res=false;
  }
    }
  nbok += res ? 1 : 0;
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "true == true" << std::endl;
  trace.endBlock();
 
  return nbok == nb;
}

References DGtal::HyperRectDomain< TSpace >::begin(), DGtal::Trace::beginBlock(), domain, DGtal::HyperRectDomain< TSpace >::end(), DGtal::Trace::endBlock(), DGtal::Trace::error(), DGtal::Trace::info(), and DGtal::trace.

Referenced by main().