Anonymous avatar Anonymous committed 0523d88

added velocity.py

Comments (0)

Files changed (2)

analysis/dat2kml.py

 		continue
 	alt = cols[3]
 	lat, long = cols[2].split(",")
+	if not float(lat) and not float(long):
+		continue
 	latlong+= "%s,%s,%s\n" % (long, lat, alt)
 
 print kml_template % (latlong,)

analysis/velocity.py

+#!/usr/bin/env python
+
+import sys
+
+if len(sys.argv) != 2:
+	print "usage: %s <datfile>"
+	sys.exit(2)
+
+latlong = ""
+fp = open(sys.argv[1],"r")
+count = 0
+for line in fp:
+	cols = line.split(" ")
+	if cols[1] != 'gps:':
+		continue
+
+	lat, long = cols[2].split(",")
+	if not float(lat) and not float(long):
+		count = 0
+		continue
+		
+	alt = float(cols[3])
+	ts = float(cols[4]) / 1000.0
+	if not count:
+		start_alt = alt
+		start_ts = ts
+	count+= 1
+
+	if count == 60:
+		vel = (alt - start_alt) / (ts - start_ts)
+		print vel
+		count = 0
+
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.