Como converter histograma 3D em mapa de calor

A solução mais fácil funciona via seaborn (uma extensão do matplotlib especializada em visualização estatística):

import matplotlib.pyplot as plt
import seaborn as sns
import numpy as np

# Fixing random state for reproducibility
np.random.seed(19680801)

x, y = np.random.rand(2, 100) * 4
sns.histplot(x=x, y=y, bins=4, binrange=[[0, 4], [0, 4]])

Com matplotlib puro, uma possibilidade com plt.imshow():

import matplotlib.pyplot as plt
import numpy as np

# Fixing random state for reproducibility
np.random.seed(19680801)

x, y = np.random.rand(2, 100) * 4
hist, xedges, yedges = np.histogram2d(x, y, bins=4, range=[[0, 4], [0, 4]])

plt.imshow(hist.T, origin='lower', extent=[xedges[0], xedges[-1], yedges[0], yedges[-1]], cmap='Blues')
plt.xticks(xedges)
plt.yticks(yedges)
plt.grid()
plt.show()