Openstreetmap maps in Python

I figured out a simple way to present data on a map in Python.
Openstreetmap server often experiences heavy loads (everyone wants access to these tiles, including me). A little more code to handle denials:
Combined with matplotlib's hexbin one can create nice density plots.

By the way: hexbin has an amazing parameter C (lacking in every other density plot function I've encountered). According to documentation:

If C is specified, it specifies values at the coordinate (x[i],y[i]). These values are accumulated for each hexagonal bin and then reduced according to reduce_C_function, which defaults to numpy’s mean function (np.mean).

 My data is counts (of Geiger-Mueller counter) per minute versus geographic position. I want to plot an average counts per minute value in each bin, not total, as other systems I've tried would have me do. Some would even limit me to number of data entries.