def fromTimeDataUpdateBaseMap(self, form, end):
# 删除掉以前的数据
a = ["../temp/matched_trips_directory_1", "../temp/newTripOut"]
for folder in a:
utils.clean_and_mkdir(folder)
trips = Trip_get(form, end, self.trace).load_trip_from_db()
self.match_trip_from_db(self.init_graphdb_filename, trips,
self.matched_trips_directory_1)
# 将中断的点提取出来 作为 下一步的输入
r.getNewTrip(self.matched_trips_directory_1, self.newTripOut, 65)
if (not os.listdir(self.newTripOut)):
print("\n 没有更新")
os._exit(0)
trips = TripLoader.load_all_trips(self.newTripOut)
self.createMap(trips, self.cell_size, self.gaussian_blur,
self.init_graphdb_filename, 3, self.graphdb_filename, 1)
cv.SaveImage("raw_data.png", lines)
print "done."
#print "Intensity map acquired."
sys.stdout.write("Smoothing... ")
sys.stdout.flush()
# # create the mask and compute the contour
cv.Smooth(themap, themap, cv.CV_GAUSSIAN, gaussian_blur, gaussian_blur)
cv.SaveImage("kde.png", themap)
print "done."
print "\nKDE generation complete."
if __name__ == '__main__':
opts,args = getopt.getopt(sys.argv[1:],"c:b:p:h")
for o,a in opts:
if o == "-c":
cell_size=int(a)
elif o == "-b":
gaussian_blur = int(a)
elif o == "-p":
trips_path = str(a)
elif o == "-h":
print "Usage: kde.py [-c <cell_size>] [-b <gaussian_blur_size>] [-p <trips_path>] [-h]\n"
sys.exit()
k = KDE()
k.create_kde_with_trips(TripLoader.load_all_trips(trips_path))
