ixp-voronoi-cells.py

#!/usr/bin/python
### Samer Lahoud http://samer.lahoud.fr

import json
import os
import matplotlib.pyplot as plt
from scipy.spatial import Voronoi, voronoi_plot_2d
import numpy as np
import mplleaflet

# Load up the geojson data
filename = ('buildings.geojson')
with open(filename) as f:
    gj = json.load(f)

# Grab the coordinates (longitude, latitude) from the features, which we
# know are Points
xy = np.array([feat['geometry']['coordinates'] for feat in gj['features']])

vor_polygon = Voronoi(xy)
voronoi_plot_2d(vor_polygon, show_vertices = False)

mapfile = 'ixp-voronoi-cells.html'
# Create the map. Save the file to html
mplleaflet.show(path=mapfile)