mtp.cc

   1
   2 ///////////////////////////////////////////////////////////////////////////
   3 // This program is free software: you can redistribute it and/or modify  //
   4 // it under the terms of the version 3 of the GNU General Public License //
   5 // as published by the Free Software Foundation.                         //
   6 //                                                                       //
   7 // This program is distributed in the hope that it will be useful, but   //
   8 // WITHOUT ANY WARRANTY; without even the implied warranty of            //
   9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU      //
  10 // General Public License for more details.                              //
  11 //                                                                       //
  12 // You should have received a copy of the GNU General Public License     //
  13 // along with this program. If not, see <http://www.gnu.org/licenses/>.  //
  14 //                                                                       //
  15 // Written by and Copyright (C) Francois Fleuret                         //
  16 // Contact <francois.fleuret@idiap.ch> for comments & bug reports        //
  17 ///////////////////////////////////////////////////////////////////////////
  18
  19 // Multi-Tracked Path
  20
  21 // Takes the graph description file as input and produces a dot file.
  22
  23 // EXAMPLE: ./mtp ./graph2.txt  | dot -T pdf -o- | xpdf -
  24
  25 #include <iostream>
  26 #include <fstream>
  27 #include <stdlib.h>
  28
  29 using namespace std;
  30
  31 #include "tracker.h"
  32
  33 //////////////////////////////////////////////////////////////////////
  34
  35 scalar_t detection_score(int true_label, scalar_t flip_noise) {
  36   if((true_label > 0) == (drand48() < flip_noise)) {
  37     return   1.0 + 0.2 * (drand48() - 0.5);
  38   } else {
  39     return - 1.0 + 0.2 * (drand48() - 0.5);
  40   }
  41 }
  42
  43 int main(int argc, char **argv) {
  44   int nb_locations = 5;
  45   int nb_time_steps = 20;
  46   int motion_amplitude = 1;
  47
  48   Tracker *tracker = new Tracker(nb_time_steps, nb_locations);
  49
  50   for(int l = 0; l < nb_locations; l++) {
  51     for(int k = 0; k < nb_locations; k++) {
  52       tracker->set_allowed_motion(l, k, abs(l - k) <= motion_amplitude);
  53     }
  54     tracker->set_as_exit(0, 1);
  55     tracker->set_as_entrance(0, 1);
  56   }
  57
  58   tracker->build_graph();
  59
  60   for(int r = 0; r < 10; r++) {
  61     cout << "* ROUND " << r << endl;
  62
  63     // We generate synthetic detection scores at location
  64     // nb_locations/2, with 10% false detection (FP or FN)
  65
  66     for(int t = 0; t < nb_time_steps; t++) {
  67       for(int l = 0; l < nb_locations; l++) {
  68         tracker->set_detection_score(t, l, detection_score(-1, 0.9));
  69       }
  70       tracker->set_detection_score(t, nb_locations/2, detection_score(1, 0.9));
  71     }
  72
  73     tracker->track();
  74
  75     for(int t = 0; t < tracker->nb_trajectories(); t++) {
  76       cout << "TRAJECTORY "
  77            << t
  78            << " [starting " << tracker->trajectory_entrance_time(t) << "]";
  79       for(int u = 0; u < tracker->trajectory_duration(t); u++) {
  80         cout << " " << tracker->trajectory_location(t, u);
  81       }
  82       cout << endl;
  83     }
  84   }
  85
  86   delete tracker;
  87
  88   exit(EXIT_SUCCESS);
  89 }