在Matplotlib极坐标图上设置径向轴 [英] Set radial axis on Matplotlib polar plots

问题描述

我正在极坐标图上绘制方位角-仰角曲线，其中仰角是径向分量.默认情况下，Matplotlib 绘制从中心的 0 到周边的 90 的径向值.我想反转它，所以 90 度是在中心.我尝试通过调用ax.set_ylim(90,0)来设置限制，但这导致抛出LinAlgError异常.ax是通过调用add_axes获得的axes对象.

I'm plotting an azimuth-elevation curve on a polar plot where the elevation is the radial component. By default, Matplotlib plots the radial value from 0 in the center to 90 on the perimeter. I want to reverse that so 90 degrees is at the center. I tried setting the limits with a call to ax.set_ylim(90,0) but this results in a LinAlgError exception being thrown. ax is the axes object obtained from a call to add_axes.

可以这样做吗?如果可以，我该怎么办?

Can this be done and, if so, what must I do?

这是我现在使用的.基本绘图代码来自Matplotlib示例之一

Here is what I'm using now. The basic plotting code was taken from one of the Matplotlib examples

# radar green, solid grid lines
rc('grid', color='#316931', linewidth=1, linestyle='-')
rc('xtick', labelsize=10)
rc('ytick', labelsize=10)

# force square figure and square axes looks better for polar, IMO
width, height = matplotlib.rcParams['figure.figsize']
size = min(width, height)
# make a square figure
fig = figure(figsize=(size, size))
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], projection='polar', axisbg='#d5de9c')

# Adjust radius so it goes 90 at the center to 0 at the perimeter (doesn't work)
#ax.set_ylim(90, 0)

# Rotate plot so 0 degrees is due north, 180 is due south

ax.set_theta_zero_location("N")

obs.date = datetime.datetime.utcnow()
az,el = azel_calc(obs, ephem.Sun())
ax.plot(az, el, color='#ee8d18', lw=3)
obs.date = datetime.datetime.utcnow()
az,el = azel_calc(obs, ephem.Moon())
ax.plot(az, el, color='#bf7033', lw=3)

ax.set_rmax(90.)
grid(True)

ax.set_title("Solar Az-El Plot", fontsize=10)
show()

由此产生的情节是

推荐答案

我设法将他的径向轴倒置了.为了匹配新轴，我不得不重新映射半径:

I managed to put he radial axis inverted. I had to remap the radius, in order to match the new axis:

fig = figure()
ax = fig.add_subplot(1, 1, 1, polar=True)

def mapr(r):
   """Remap the radial axis."""
   return 90 - r

r = np.arange(0, 90, 0.01)
theta = 2 * np.pi * r / 90

ax.plot(theta, mapr(r))
ax.set_yticks(range(0, 90, 10))                   # Define the yticks
ax.set_yticklabels(map(str, range(90, 0, -10)))   # Change the labels

请注意，这只是一个小技巧，轴的中心仍为0，周长为90.您必须对要绘制的所有变量使用映射函数.

Note that is just a hack, the axis is still with the 0 in the center and 90 in the perimeter. You will have to use the mapping function for all the variables that you are plotting.

这篇关于在Matplotlib极坐标图上设置径向轴的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！