import subprocess
subprocess.call('pip install mapbox'.split())

nyc_shoreline_geojson_url = 'https://data.cityofnewyork.us/api/geospatial/2qj2-cctx?method=export&format=GeoJSON'

import requests
response = requests.get(nyc_shoreline_geojson_url)
response_json = response.json()

from shapely.geometry import GeometryCollection, shape
features = response_json['features']
shoreline_geometry = GeometryCollection([shape(f['geometry']) for f in features])

shoreline_geometry.simplify(0.01)

# Get coastline
# Get senior centers
# Get subway stations