目前我有以下脚本,用于生成方位角/半径数据的极坐标图。"R1"只是一个由ArcGIS表中的[方位角,倾角]值列表派生而来的值列表。
import matplotlib.pyplot as plt
import numpy as np
for(a,r) in R1:
angles.append(a)
radius.append(90-r)
theta = np.radians(angles)
r = radius
ax = plt.subplot(111,polar=True)
ax.plot(theta, r, color='black', ls='-', linewidth=1)
ax.fill(theta,r,'0.75') ## should I use ax.fill_betweenx() ?
ax.set_theta_zero_location('N')
ax.set_theta_direction(-1)
ax.set_rmax(90)
ax.set_rmin(0)
ax.set_yticks(range(0,90,10))
yLabel=['90','','','60','','','30','','','']
ax.set_yticklabels(yLabel)
ax.grid(True)
plt.show()
目前,这将创建以下绘图:
如何“反转”填充,使灰色部分变为白色,白色部分变为灰色?
我尝试了ax.fill_betweenx(theta,90,r,color='0.75'),但没有成功。我已经花了一些时间与此作斗争,但都无济于事。
非常感谢任何帮助或建议!
如果有任何方法可以使此更清晰,请在评论中让我知道。
theta = np.linspace(0, 2*np.pi, 100)
和r = 2 + np.sin(theta)
来演示。 - tacaswellset_yticklabels
是危险的,应该避免使用。 - tacaswellset_yticklabels
要好一些。为了看到危险,请更改r限制或刻度位置。标签将完全与数据脱钩。 - tacaswell