Есть ли способ создания двумерной версии сюжета scipy.stats.skewnorm?

Заимствовать код непосредственно из этих двух сайтов:

Документы Skewnorm

Hist3D демо

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import skewnorm

fig = plt.figure(figsize=(15,15))
ax = fig.add_subplot(111, projection='3d')

#generate data from skewnorm dist
x,y = skewnorm.rvs(.2, size=1000),skewnorm.rvs(.2, size=1000)
hist, xedges, yedges = np.histogram2d(x, y, bins=4, range=[[0, 4], [0, 4]])

xpos, ypos = np.meshgrid(xedges[:-1] + 0.25, yedges[:-1] + 0.25)
xpos = xpos.flatten('F')
ypos = ypos.flatten('F')
zpos = np.zeros_like(xpos)

dx = 0.5 * np.ones_like(zpos)
dy = dx.copy()
dz = hist.flatten()

ax.bar3d(xpos, ypos, zpos, dx, dy, dz, color='b', zsort='average')

plt.show()

[ 113]