使用 matplotlib 绘制年度数据的极坐标图

Question

我正在尝试在 matplotlib 中将一整年的数据绘制为极坐标图，但无法找到任何此类示例。我设法根据this thread 转换了 pandas 的日期，但我无法将头绕在（字面上）y 轴或 theta 周围。

这是我走了多远：

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
import pandas as pd

times = pd.date_range("01/01/2016", "12/31/2016")
rand_nums = np.random.rand(len(times),1)
df = pd.DataFrame(index=times, data=rand_nums, columns=['A'])

ax = plt.subplot(projection='polar')
ax.set_theta_direction(-1)
ax.set_theta_zero_location("N")
ax.plot(mdates.date2num(df.index.to_pydatetime()), df['A'])
plt.show()

这给了我这个情节：

缩小日期范围以了解发生了什么
times = pd.date_range("01/01/2016", "01/05/2016")我得到了这个情节：

我收集到该系列的开始在 90 到 135 之间，但我怎样才能“重新映射”它，以便我的年份日期范围在北原点开始和结束？

Answer 1

极坐标图的角度范围以辐射线为单位，即[0, 2π]。 因此，需要将日期范围标准化为完整的圆圈。

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
import pandas as pd

times = pd.date_range("01/01/2016", "12/31/2016")
rand_nums = np.random.rand(len(times),1)
df = pd.DataFrame(index=times, data=rand_nums, columns=['A'])

ax = plt.subplot(projection='polar')
ax.set_theta_direction(-1)
ax.set_theta_zero_location("N")

t = mdates.date2num(df.index.to_pydatetime())
y = df['A']
tnorm = (t-t.min())/(t.max()-t.min())*2.*np.pi
ax.fill_between(tnorm,y ,0, alpha=0.4)
ax.plot(tnorm,y , linewidth=0.8)
plt.show()