imageSetDT.cpp

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <algorithm>
#include <boost/algorithm/minmax_element.hpp>

#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
 
#include "DGtal/base/BasicFunctors.h"
#include "DGtal/kernel/BasicPointPredicates.h"
#include "DGtal/kernel/sets/DigitalSetInserter.h"
 
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/images/ImageHelper.h"
#include "DGtal/geometry/volumes/distance/DistanceTransformation.h"
#include "DGtal/images/IntervalForegroundPredicate.h"
 
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/io/readers/PGMReader.h"
#include "DGtal/io/colormaps/HueShadeColorMap.h"
#include "DGtal/io/colormaps/GrayscaleColorMap.h"
 
#include "ConfigExamples.h"
 

using namespace DGtal;
 
int main()
{
  typedef DGtal::ImageContainerBySTLVector< Z2i::Domain, unsigned char> Image;
  typedef DGtal::GrayscaleColorMap<unsigned char> Gray;
 

  std::string filename =  examplesPath + "samples/contourS.pgm";
  Image image = DGtal::PGMReader<Image>::importPGM(filename);
  DGtal::trace.info() << "Imported image: "<<image<<std::endl;
 

  DGtal::Board2D aBoard;
  aBoard << image.domain();
  aBoard.saveSVG("imageDomainTuto.svg");
  aBoard.clear();
  Display2DFactory::drawImage<Gray>(aBoard, image, (unsigned char)0, (unsigned char)255);
  aBoard.saveEPS("imageDomainTuto2.eps");
 
 
  Z2i::DigitalSet mySet(image.domain());
  DigitalSetInserter<Z2i::DigitalSet> inserter(mySet);
  setFromImage(image, inserter, 1, 135);
  aBoard.clear();
  aBoard << mySet.domain()
         << mySet;
  aBoard.saveEPS("imageDomainTuto2bis.eps");
 

  typedef functors::IntervalForegroundPredicate<Image> Binarizer;
  Binarizer b(image,1, 135);
  typedef DGtal::DistanceTransformation<Z2i::Space, Binarizer, Z2i::L2Metric> DTL2;
  DTL2 dt(&image.domain(),&b, &Z2i::l2Metric );
 

  DTL2::Value maxDT = (*boost::first_max_element(dt.constRange().begin(),
                                                 dt.constRange().end(), std::less<DTL2::Value>() ));
  typedef DGtal::HueShadeColorMap<DTL2::Value,2> HueTwice;
 
  aBoard.clear();
  Display2DFactory::drawImage<HueTwice>(aBoard, dt, (DTL2::Value)0,
                                        (DTL2::Value)maxDT);
  aBoard.saveEPS("imageDomainTuto3.eps");
 
  return 0;
 
}