testMeasure.cpp

#include <cstdio>
#include <cmath>
#include <iostream>
#include <fstream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/math/MeasureOfStraightLines.h"


using namespace DGtal;
using namespace std;



void testUnitSquare()
{
    vector<double> a;
    vector<double> b;

    a.push_back(0);
    b.push_back(0);
    a.push_back(1);
    b.push_back(0);
    a.push_back(1);
    b.push_back(1);
    a.push_back(0);
    b.push_back(1);

    MeasureOfStraightLines measure;


    trace.info() << "Measure of the Straight of Lines of the unit square = " << measure.computeMeasure(a,b)<< std::endl;
    trace.emphase() <<"Expected value = 0.707107 (sqrt(2)/2)"<<endl;
}


void testUnitSquareCentroid()
{
    vector<double> a;
    vector<double> b;

    a.push_back(0);
    b.push_back(0);
    a.push_back(1);
    b.push_back(0);
    a.push_back(1);
    b.push_back(1);
    a.push_back(0);
    b.push_back(1);

    MeasureOfStraightLines measure;

    trace.info() << "Centroid measure of the unit square = (" << measure.computeCentroidA(a,b)
    << ","<<measure.computeCentroidB(a,b)<<")"<<std::endl;
    trace.emphase() <<"Expected value = (0.4142,0.5)"<<endl;
}


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