%matplotlib inline
import matplotlib.pyplot as plt
import pysal as ps
import numpy as np
from pysal.contrib.viz import mapping as maps
import cartopy.io.img_tiles as cimgt
import cartopy.crs as ccrs

shp = ps.open('boroughs.shp')
bg = cimgt.OSM()

extent = [shp.bbox[0], shp.bbox[2], shp.bbox[1], shp.bbox[3]]
boros = maps.map_poly_shp(shp)
boros.set_transform(ccrs.Geodetic())
boros.set_facecolor('none')

f = plt.figure(figsize=(10, 8))

ax = plt.axes(projection=bg.crs)
ax.set_extent(extent)
ax.add_image(bg, 11)

ax.add_collection(boros)

f.suptitle('London boroughs, all together')

plt.show()

f = plt.figure(figsize=(9, 9))

for i, boro in enumerate(shp):
    extent = [boro.bbox[0], boro.bbox[2], boro.bbox[1], boro.bbox[3]]
    boroP = maps.map_poly_shp([boro], setup=False)
    boroP.set_transform(ccrs.PlateCarree())
    boroP.set_facecolor('k')
    
    ax = plt.subplot(6, 6, i+1, projection=ccrs.OSGB())
    ax.set_extent(extent)

    ax.add_collection(boroP)
    ax.outline_patch.set_visible(False)

    f.add_axes(ax)

f.suptitle('London boroughs, one by one')
plt.show()

from IPython.core.display import HTML
def css_styling():
    styles = open("styles.css", "r").read()
    return HTML(styles)
css_styling()