testMeaningfulScaleAnalysis.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtalCatch.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/math/Profile.h"
#include "DGtal/math/MeaningfulScaleAnalysis.h"
 
 
using namespace std;
using namespace DGtal;
 
// Functions for testing class MeaningfulScaleAnalysis.
 
 
 
 
struct LogFct{
  double operator()(const double &a) const {
    return log(a);
  }
};
  
 
TEST_CASE( "Testing MeaningfulScaleAnalysis" )
{
  
  Profile<LogFct> sp (Profile<LogFct>::MEAN);
  sp.init(6);
  
  SECTION("Testing noise level detect of Profile")
    {
      sp.addValue(0,22);
      sp.addValue(1,15);
      sp.addValue(2,8);
      sp.addValue(3,17);
      sp.addValue(4,7);
      sp.addValue(5,2);      
 
      MeaningfulScaleAnalysis<Profile<LogFct>> msa(sp);
      std::vector< std::pair<unsigned int, unsigned int> > interval;      
      msa.computeMeaningfulScales(interval, 1);
      msa.getSlopeFromMeaningfulScales(0,0,2);
      msa.lowerBoundedNoiseLevel(0,10,2,2,2);
      unsigned int n = msa.noiseLevel();
      REQUIRE( interval[0].first == 1 );
      REQUIRE( interval[0].second == 3 );
      REQUIRE( interval[1].first == 4 );
      REQUIRE( n == 1 );
    }
 
 
 
}