带有Python的2d极坐标直方图

如何解决带有Python的2d极坐标直方图

我正在尝试制作二维极坐标直方图。我的数据是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是面元边缘。但我不确定如何进行实际的归一化：

任何建议都会很有帮助

解决方法

在maptlotlib中使用极坐标图之前，必须转换辐射角。