featureson.h

/*Copyright 2014 Francisco Alvaro

 This file is part of SESHAT.

    SESHAT is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    SESHAT 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 SESHAT.  If not, see <http://www.gnu.org/licenses/>.

This file is a modification of the online features original software
covered by the following copyright and permission notice:

*/
/*
    Copyright (C) 2006,2007 Moisés Pastor <mpastorg@dsic.upv.es>

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU 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/>.
*/
#ifndef FEATURES_H
#define FEATURES_H

#include <math.h>
#include <iostream>
#include <iomanip>
#include <vector>
#include <values.h>
#include "online.h"

#define MAXNUMHATS 200
#define OFFSET_INS 20

using namespace std;

class frame {
  public:
    double x,y,dx,dy,ax,ay,k;

    void print(ostream & fd);
    int get_fr_dim();

    double getFea(int i){
      switch (i){
      case 0: return x;
      case 1: return y;
      case 2: return dx;
      case 3: return dy;
      case 4: return ax;
      case 5: return ay;
      case 6: return k;
      default:
	fprintf(stderr, "Error: getFea(%d)\n", i);
	exit(-1);
      }
    }
};

class sentenceF {
  public:
    string transcrip;
    int n_frames;
    frame * frames;
    
    sentenceF();
    ~sentenceF();

    bool data_plot(ostream & fd);
    bool print(ostream & fd);

    void calculate_features(sentence &s);

  private:
    vector<PointR> normalizaAspect(vector<Point> & puntos);
    void calculate_derivatives(vector<PointR> & points, bool norm=true);
    void calculate_kurvature();
};


#endif