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);
import trajcomp;
from operator import itemgetter
import matplotlib.pyplot as plt
from matplotlib.patches import Rectangle
print ("Loading first set of GeoLife Dataset");
A=trajcomp.geolife(15,"/home/martin/datasets/geolife/")
handles = range(A,trajcomp.size());
print "Found " + str(trajcomp.size())
xlimit = [39.689602,40.122410]
ylimit = [116.105789,116.670021]
# list of lists is suboptimal, but does not depend on numpy
for i in handles:
A=trajcomp.get(i);
trajcomp.summary(i);
X = map(itemgetter(0), A) ;
Y = map(itemgetter(1),A);
#plt.clf()
plt.plot(X, Y, '-')
# Grossraum
# plt.xlim([39.75,40.05]);
# plt.ylim([116.116,116.667]);
plt.xlim(xlimit);
plt.ylim(ylimit)
#plt.show()
plt.savefig(("frame-%04d.png" % i), bbox_inches='tight')
