我正在尝试做一个二维极坐标直方图。我的数据是 2 个 np.arrays：
方位角 = [方位角 -180 到 +180 度]
天顶 = [天顶角 0 到 90 度]
这是我到目前为止的代码：
# define binning
rbins = np.linspace(0, 90,10)
abins = np.linspace(-180, 180, 30)
#calculate histogram
hist, _, _ = np.histogram2d(azim, zenith, bins=(abins, rbins))
A, R = np.meshgrid(abins, rbins)
# plot
fig, ax = plt.subplots(subplot_kw=dict(projection="polar"))
pc = ax.pcolormesh(A, R, hist.T, cmap="magma_r")
fig.colorbar(pc)
ax.grid(True)
plt.show()
我得到的情节：
看起来好像它的一部分丢失了，我不知道为什么。我也想将天顶角归一化为立体角，而不是归一化为一个。归一化的公式是这样的，其中 θ1，θ2 是 bin 边缘。但我不确定如何进行实际的标准化：
任何建议都会很有帮助