X-Git-Url: https://fleuret.org/cgi-bin/gitweb/gitweb.cgi?a=blobdiff_plain;f=mtp.cc;h=2f9e7a98f593227e1d4d84be31d012a2865f235a;hb=8f7e0196dd3c89cb4e68732d6812a13d43bdfbc7;hp=979915b57b9b384207648847ac25f9cd8925b34c;hpb=535b8ede7d6646ea17135da10c8a16416524f6e9;p=mtp.git diff --git a/mtp.cc b/mtp.cc index 979915b..2f9e7a9 100644 --- a/mtp.cc +++ b/mtp.cc @@ -25,62 +25,37 @@ using namespace std; #include "tracker.h" -////////////////////////////////////////////////////////////////////// +int main(int argc, char **argv) { -scalar_t detection_score(int true_label, scalar_t flip_noise) { - if((true_label > 0) == (drand48() < flip_noise)) { - return 1.0 + 0.2 * (drand48() - 0.5); - } else { - return - 1.0 + 0.2 * (drand48() - 0.5); + if(argc < 2) { + cerr << argv[0] << " " << endl; + exit(EXIT_FAILURE); } -} -int main(int argc, char **argv) { - int nb_locations = 6; - int nb_time_steps = 5; - int motion_amplitude = 1; - - Tracker *tracker = new Tracker(nb_time_steps, nb_locations); - - for(int l = 0; l < nb_locations; l++) { - for(int k = 0; k < nb_locations; k++) { - tracker->allowed_motion[l][k] = abs(l - k) <= motion_amplitude; - } - tracker->entrances[0] = 1; - tracker->exits[nb_locations - 1] = 1; - } + Tracker tracker; - tracker->build_graph(); + ifstream in_tracker(argv[1]); - // We generate synthetic detection scores at location - // nb_locations/2, with 5% false detection (FP or FN) + if(in_tracker.good()) { - for(int t = 0; t < nb_time_steps; t++) { - for(int l = 0; l < nb_locations; l++) { - tracker->detection_score[t][l] = detection_score(-1, 0.95); - } - tracker->detection_score[t][nb_locations/2] = detection_score(1, 0.95); - } + tracker.read(&in_tracker); + tracker.build_graph(); + tracker.track(); - tracker->track(); + ofstream out_traj("result.trj"); + tracker.write_trajectories(&out_traj); + cout << "Wrote result.trj" << endl; - for(int t = 0; t < tracker->nb_trajectories(); t++) { - cout << "TRAJECTORY " - << t - << " [starting " << tracker->trajectory_entrance_time(t) << "]"; - for(int u = 0; u < tracker->trajectory_duration(t); u++) { - cout << " " << tracker->trajectory_location(t, u); - } - cout << endl; - } + ofstream out_dot("graph.dot"); + tracker.print_graph_dot(&out_dot); + cout << "Wrote graph.dot" << endl; - { - ofstream dot("graph.dot"); - tracker->print_graph_dot(&dot); - cout << "Wrote graph.dot." << endl; - } + } else { - delete tracker; + cerr << "Can not open " << argv[1] << endl; + exit(EXIT_FAILURE); + + } exit(EXIT_SUCCESS); }