Recorded location of three CapMetro vehicles on 01/20/15. Click on an individual marker to see the time at that position.

11106@2015-01-20T14:57:00-06:00<---->2015-01-20T15:31:00-06:00.geojson

5003@2015-01-20T14:42:00-06:00<---->2015-01-20T14:57:00-06:00.geojson

5053@2015-01-20T14:42:00-06:00<---->2015-01-20T15:04:00-06:00.geojson

example.py

#!/usr/bin/env python
import rethinkdb as r
import geojson as gj
import arrow

import sys


def make_feature(vehicle, id):
    point = gj.Point(tuple(vehicle['location']['coordinates']))
    time = arrow.get(vehicle['timestamp']).to('America/Chicago')
    f = {
        'geometry': point,
        'properties': {'time': time.isoformat()},
        'id': id
    }
    return gj.Feature(**f)


def make_featcoll(vehicle_id, limit):
    docs = get_vehicle_data(vehicle_id, limit)

    features = []
    for ind, doc in enumerate(docs):
        features.append(make_feature(doc, ind))

    feat_coll = gj.FeatureCollection(features)

    start = arrow.get(docs[0]['timestamp']).to('America/Chicago').isoformat()
    end = arrow.get(docs[-1]['timestamp']).to('America/Chicago').isoformat()
    fname = '{vid}@{start}<-->{end}.geojson'.format(vid=vehicle_id, start=start, end=end)
    with open(fname, 'w') as f:
        gj.dump(feat_coll, f)


def get_vehicle_data(vehicle_id, limit):
    return list(r.table('vehicle_position') \
                .filter(r.row['vehicle_id'] == vehicle_id) \
                .order_by('timestamp') \
                .limit(int(limit)).run())


if __name__ == '__main__':
    if len(sys.argv) < 3:
        print 'Usage: explore.py <vehicle_id> <limit>'
        sys.exit(-1)

    r.connect('localhost', 28015).repl()
    make_featcoll(sys.argv[1], sys.argv[2])