import time;
def elapsed(var):
t = time.clock();
r = eval(var)
return (r,time.clock()-t)
epsilon = float(sys.argv[1]);
print ("Loading first set of GeoLife Dataset");
A=trajcomp.geolife(1,"/home/martin/datasets/geolife/")
handles = range(A,trajcomp.size());
print "Found " + str(trajcomp.size())
# list of lists is suboptimal, but does not depend on numpy
for i in handles:
for j in handles:
(v_real,t_real) = elapsed("trajcomp.frechet(%d,%d,-1)" % (i,j))
A = trajcomp.douglas_peucker(i,epsilon)
B = trajcomp.douglas_peucker(j,epsilon)
# trajcomp.summary(i);
# trajcomp.summary(j);
# trajcomp.summary(A);
# trajcomp.summary(B);
(v_sim,t_sim) = elapsed("trajcomp.frechet(%d,%d,-1)" % (A,B))
print str(epsilon) + " " + str(i) +" "+ str(j)+" " + str (t_real) + " " + str(t_sim) + " " +str(v_real)+" " + str(v_sim) +" " + str(v_real / v_sim)
import trajcomp;
from operator import itemgetter
import matplotlib.pyplot as plt
from matplotlib.patches import Rectangle
A=trajcomp.load("../../data/prague1.dat");
print "Reference is " + str(A)
# list of lists is suboptimal, but does not depend on numpy
for i in range(1,150):
epsilon = i/10.0;
epsilon = epsilon * epsilon;
print epsilon;
B = trajcomp.douglas_peucker(A,epsilon);
trajcomp.summary(B);
C =trajcomp.get(B)
X = map(itemgetter(0), C) ;
Y = map(itemgetter(1),C);
plt.clf();
fig, ax1 = plt.subplots()
ax1.plot(X, Y, '-o')
ax1.set_xlabel('X')
ax1.set_ylabel('Y')
ax2 = ax1.twinx()
plt.ylim(ymax=(15*15),ymin=0);
space=abs(max(X)-min(X))/20;
ax2.bar(max(X)+2*space,epsilon,width=space);
