import json from geopy.distance import vincenty import sys data = [] with open(sys.argv[1]) as f: for line in f: tmp = json.loads("{" + line.split('{',1)[-1]) if "lat" in tmp: data.append(tmp) a=(0,0) t=0 for d in data: if a != (0,0): dist = vincenty(a, (d['lat'], d['lon'])).km if dist > 1: #print(dist) delta=int(d['tst'])-int(t) if float(dist) == 0: continue if delta == 0: continue speed=(float(dist)/delta)*3600 if speed > 500: #print(speed) print("sed -i.bak '/" + str(d['tst']) + "/d' " + sys.argv[1]) a = (d['lat'], d['lon']) t = int(d['tst'])