testL1LengthEstimator.cpp

#include <iostream>
#include "DGtal/base/Common.h"

#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/kernel/sets/DigitalSetSelector.h"
#include "DGtal/topology/KhalimskySpaceND.h"
#include "DGtal/topology/SurfelAdjacency.h"
#include "DGtal/topology/SurfelNeighborhood.h"


#include "DGtal/geometry/curves/GridCurve.h"
#include "DGtal/geometry/2d/L1LengthEstimator.h"
#include "DGtal/geometry/2d/TwoStepLocalLengthEstimator.h"
#include "DGtal/geometry/2d/MLPLengthEstimator.h"
#include "DGtal/geometry/2d/FPLengthEstimator.h"
#include "DGtal/geometry/2d/DSSLengthEstimator.h"

#include "ConfigTest.h"



using namespace std;
using namespace DGtal;

// Functions for testing Length Estimator classes.

bool testL1LengthEstimator(std::string &filename)
{

  trace.info() << "Reading GridCurve " << endl;
  
  ifstream instream; // input stream
  instream.open (filename.c_str(), ifstream::in);

  GridCurve<KhalimskySpaceND<2> > c; //grid curve
  c.initFromVectorStream(instream);

  GridCurve<KhalimskySpaceND<2> >::ArrowsRange ra = c.getArrowsRange(); //range
  L1LengthEstimator<  GridCurve<KhalimskySpaceND<2> >::ArrowsRange::ConstIterator > l1length;
  TwoStepLocalLengthEstimator<  GridCurve<KhalimskySpaceND<2> >::ArrowsRange::ConstIterator > locallength(1.0,sqrt(2.0));
    
  l1length.init(1, ra.begin(), ra.end(), c.isClosed());
  trace.info() << "L1 length (h=1) = "<< l1length.eval()<<std::endl;

  l1length.init(10, ra.begin(), ra.end(), c.isClosed());
  trace.info() << "L1 length (h=10) = "<< l1length.eval()<<std::endl;

  TwoStepLocalLengthEstimator<  GridCurve<KhalimskySpaceND<2> >::ArrowsRange::ConstIterator > locallength(1.0,sqrt(2.0));
 
  localength.init(1, ra.begin(), ra.end(), c.isClosed());
  trace.info() << "Local length (h=1) = "<< localength.eval()<<std::endl;

    
  GridCurve<KhalimskySpaceND<2> >::PointsRange rp = c.getPointsRange(); //range
  MLPLengthEstimator<  GridCurve<KhalimskySpaceND<2> >::PointsRange::ConstIterator > MLPlength;
    
  MLPlength.init(1, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "MLP Length (h=1) = "<< MLPlength.eval()<<std::endl;

  MLPlength.init(10, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "MLP Length (h=10) = "<< MLPlength.eval()<<std::endl;

  FPLengthEstimator<  GridCurve<KhalimskySpaceND<2> >::PointsRange::ConstIterator > FPlength;
    
  FPlength.init(1, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "FP Length (h=1) = "<< FPlength.eval()<<std::endl;

  FPlength.init(10, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "FP Length (h=10) = "<< FPlength.eval()<<std::endl;

  DSSLengthEstimator<  GridCurve<KhalimskySpaceND<2> >::PointsRange::ConstIterator > DSSlength;
    
  DSSlength.init(1, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "DSS Length (h=1) = "<< DSSlength.eval()<<std::endl;

  DSSlength.init(10, rp.begin(), rp.end(), c.isClosed());
  trace.info() << "DSS Length (h=10) = "<< DSSlength.eval()<<std::endl;

  return true;
}



// Standard services - public :

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

  std::string sinus2D4 = testPath + "samples/sinus2D4.dat";



  bool res = testL1LengthEstimator(sinus2D4); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //