6

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

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

编辑:这是我现在使用的。基本绘图代码取自 Matplotlib 示例之一

# 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()

由此产生的情节是

在此处输入图像描述

4

1 回答 1

4

我设法把他的径向轴倒置了。我不得不重新映射半径,以匹配新轴:

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。您必须对正在绘制的所有变量使用映射函数。

于 2012-09-12T16:02:33.153 回答