%matplotlib inline
import matplotlib.pyplot as plt
plt.rcParams['figure.figsize'] = (10.0, 6.0)

# Downloaded from http://www.volcano.si.edu/list_volcano_excel.cfm
# and converted to CSV.
from matplotlib.pylab import csv2rec
gvp = csv2rec('/downloads/GVP_Volcano_List.csv', skiprows=1)
print gvp.size, gvp.longitude, gvp.latitude

import cartopy.crs as ccrs
ax = plt.axes(projection=ccrs.PlateCarree())
plt.plot(gvp.longitude, gvp.latitude, '^r', transform=ccrs.Geodetic())
ax.stock_img()
plt.show()

def draw_volcanoes(projection):
    ax = plt.axes(projection=projection)
    plt.plot(gvp.longitude, gvp.latitude, '^r', transform=ccrs.Geodetic())
    ax.stock_img()
    return ax

draw_volcanoes(ccrs.PlateCarree(central_longitude=180))
plt.show()

draw_volcanoes(ccrs.InterruptedGoodeHomolosine())
plt.show()