geojsoncompat.py

import json
from math import radians, cos, sin, asin, sqrt
import string


def distBetween(lon1, lat1, lon2, lat2):
  lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])
  # haversine formula
  dlon = lon2 - lon1
  dlat = lat2 - lat1
  a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
  c = 2 * asin(sqrt(a))
  r = 3956  # Radius of earth in kilometers. Use 6371 for kilometers.
  return(c * r)

def filterByRadiusAndTime(dict, lat, long, radius, past_day):
  python_obj = dict.get('result').get('records')
  geoJsonDicts = {"type": "FeatureCollection"}
  arrayOfObjects = []

  for object in python_obj:
    if (object["lat"] == None or object["long"] == None or object["month"] == None or object["year"] == None or object["street"] == None or object["occurred_on_date"] == None):
      continue

    lat1 = float(object["lat"])
    long1 = float(object["long"])
    month = object["month"]
    year = object["year"]
    date = object["occurred_on_date"]

    if (distBetween(long1, lat1, long, lat) <= radius) and (date > past_day):
      arrayOfObjects.append({
          "geometry": {
              "type": "Point",
              "coordinates": [long1, lat1]
          },
          "type": "Feature",
          "properties": {
              "Offense description": object["offense_description"][0] + object["offense_description"][1:].lower(),
              "Date": month + "/" + date[8:10] + "/" + year,
              "Time": date[11:],
              "District": object["district"],
              "Street": object["street"].title()
          }
      })
      geoJsonDicts["features"] = arrayOfObjects
  return (geoJsonDicts)