Any hints on how to visualize the density difference between two hexbin plot, each with x-y (2D) data?

Each set has roughly 300,000 points.

The range in x and y values for each data set are roughly similar but with slightly different density:

range x,x1: -1.222 to 3.656

range y,y1: 13.191, 18.150

So the steps:

(1) Produce the two hexbin maps:

fig1=hexbin(c_b204_jmk,c_b204_k,C = None, gridsize = 100, bins = None, xscale = ‘linear’, yscale = ‘linear’, cmap=None, norm=None, vmin=None, vmax=None, alpha=None, linewidths=None, edgecolors=‘none’,reduce_C_function = np.mean, mincnt=None, marginals=False)

fig2=hexbin(c_b211_jmk,c_b211_k,C = None, gridsize = 100, bins = None, xscale = ‘linear’, yscale = ‘linear’, cmap=None, norm=None, vmin=None, vmax=None, alpha=None, linewidths=None, edgecolors=‘none’,reduce_C_function = np.mean, mincnt=None, marginals=False)

(2) Determine the difference in the hexbin counts, with

diff_hex=fig2.get_array()-fig1.get_array()

(3) BUT when i try to plot the diff hexbin map

fig3=hexbin(c_b211_jmk,c_b211_k,C = diff_hex, gridsize = 100, bins = None, xscale = ‘linear’, yscale = ‘linear’, cmap=None, norm=None, vmin=None, vmax=None, alpha=None, linewidths=None, edgecolors=‘none’,reduce_C_function = np.mean, mincnt=None, marginals=False)

I obviously get a:

IndexError: index out of bounds

Because:

len(c_b211_jmk) = len(c_b211_k) != len(diff_hex),

Since diff_hex are the binned counts.

(4) Any simple way to get around this, to plot the hexbin difference counts on top of the the hexbin (c_b211_jmk,c_b211_k) distribution?

thanks in advance & with best regards,

- Sebastian