StackOverflow2013

Note that there are some explanatory texts on larger screens.

plurals

PO
text
Body
copied!<p>Here it follows two functions: <code>hist2d_bubble</code> and <code>hist3d_bubble</code>; that may fit for your purpose:</p> <p><img src="https://i.stack.imgur.com/48K32.png" alt="enter image description here"></p> <pre><code>import numpy as np import matplotlib.pyplot as pyplot from mpl_toolkits.mplot3d import Axes3D def hist2d_bubble(x_data, y_data, bins=10): ax = np.histogram2d(x_data, y_data, bins=bins) xs = ax[1] ys = ax[2] points = [] for (i, j), v in np.ndenumerate(ax[0]): points.append((xs[i], ys[j], v)) points = np.array(points) fig = pyplot.figure() sub = pyplot.scatter(points[:, 0],points[:, 1], color='black', marker='o', s=128*points[:, 2]) sub.axes.set_xticks(xs) sub.axes.set_yticks(ys) pyplot.ion() pyplot.grid() pyplot.show() return points, sub def hist3d_bubble(x_data, y_data, z_data, bins=10): ax1 = np.histogram2d(x_data, y_data, bins=bins) ax2 = np.histogram2d(x_data, z_data, bins=bins) ax3 = np.histogram2d(z_data, y_data, bins=bins) xs, ys, zs = ax1[1], ax1[2], ax3[1] smart = np.zeros((bins, bins, bins),dtype=int) for (i1, j1), v1 in np.ndenumerate(ax1[0]): if v1 == 0: continue for k2, v2 in enumerate(ax2[0][i1]): v3 = ax3[0][k2][j1] if v1 == 0 or v2 == 0 or v3 == 0: continue num = min(v1, v2, v3) smart[i1, j1, k2] += num v1 -= num v2 -= num v3 -= num points = [] for (i, j, k), v in np.ndenumerate(smart): points.append((xs[i], ys[j], zs[k], v)) points = np.array(points) fig = pyplot.figure() sub = fig.add_subplot(111, projection='3d') sub.scatter(points[:, 0], points[:, 1], points[:, 2], color='black', marker='o', s=128*points[:, 3]) sub.axes.set_xticks(xs) sub.axes.set_yticks(ys) sub.axes.set_zticks(zs) pyplot.ion() pyplot.grid() pyplot.show() return points, sub </code></pre> <p>The two figures above were created using:</p> <pre><code>temperature = [4, 3, 1, 4, 6, 7, 8, 3, 1] radius = [0, 2, 3, 4, 0, 1, 2, 10, 7] density = [1, 10, 2, 24, 7, 10, 21, 102, 203] import matplotlib matplotlib.rcParams.update({'font.size':14}) points, sub = hist2d_bubble(radius, density, bins=4) sub.axes.set_xlabel('radius') sub.axes.set_ylabel('density') points, sub = hist3d_bubble(temperature, density, radius, bins=4) sub.axes.set_xlabel('temperature') sub.axes.set_ylabel('density') sub.axes.set_zlabel('radius') </code></pre> <p>Related:</p> <p><a href="https://stackoverflow.com/questions/1721273/howto-bin-series-of-float-values-into-histogram-in-python">Howto bin series of float values into histogram in Python?</a></p> <p><a href="https://stackoverflow.com/questions/8437788/how-to-correctly-generate-a-3d-histogram-using-numpy-or-matplotlib-built-in-func">How to correctly generate a 3d histogram using numpy or matplotlib built in functions in python?</a></p> <p><a href="https://stackoverflow.com/questions/6832168/2d-histogram-with-python">2D histogram with Python</a></p>

Querying!

Guidance

An individual column

Larger individual text columns get their own page to allow for proper reading.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload