6

我正在尝试从使用熊猫创建的时间序列图中将绘图的 xlimits 作为 python 日期时间对象。使用ax.get_xlim()将轴限制返回为 a numpy.float64,我不知道如何将数字转换为可用的日期时间。

import pandas 
from matplotlib import dates
import matplotlib.pyplot as plt
from datetime import datetime
from numpy.random import randn

ts = pandas.Series(randn(10000), index=pandas.date_range('1/1/2000',
    periods=10000, freq='H')) 
ts.plot()
ax = plt.gca()

ax.set_xlim(datetime(2000,1,1))
d1, d2 = ax.get_xlim()
print "%s(%s) to %s(%s)" % (d1, type(d1), d2, type(d2))

print "Using matplotlib: %s" % dates.num2date(d1)
print "Using datetime: %s" % datetime.fromtimestamp(d1)

返回:

262968.0 (<type 'numpy.float64'>) to 272967.0 (<type 'numpy.float64'>)
Using matplotlib: 0720-12-25 00:00:00+00:00
Using datetime: 1970-01-03 19:02:48

根据pandas timeseries docs,pandas 使用 numpy.datetime64 dtype。我正在使用熊猫版本'0.9.0'。

我使用get_xlim()的是直接访问 pandas 系列,因为当用户在绘图区域中移动时,我使用 xlim_changed回调来做其他事情。

破解以获取可用值

对于上面的示例,限制从 Epoch 开始以小时为单位返回。所以我可以转换为自 Epoch 以来的秒数time.gmtime(),并使用它来获得可用的地方,但这仍然感觉不对。

In [66]: d1, d2 = ax.get_xlim()

In [67]: time.gmtime(d1*60*60)
Out[67]: time.struct_time(tm_year=2000, tm_mon=1, tm_mday=1, tm_hour=0, tm_min=0, tm_sec=0, tm_wday=5, tm_yday=1, tm_isdst=0)  
4

2 回答 2

6

matplotlib.dates的当前行为:

datetime 对象被转换为浮点数,表示自 0001-01-01 UTC 以来的天数加上 1。例如,0001-01-01, 06:00 是 1.25,而不是 0.25。辅助函数 date2num()、num2date() 和 drange() 用于促进轻松转换为日期时间和数字范围。

pandas.tseries.converter.PandasAutoDateFormatter() 似乎建立在此之上,所以:

x = pandas.date_range(start='01/01/2000', end='01/02/2000')
plt.plot(x, x)
matplotlib.dates.num2date(plt.gca().get_xlim()[0])

给出:

datetime.datetime(2000, 1, 1, 0, 0, tzinfo=<matplotlib.dates._UTC object at 0x7ff73a60f290>)
于 2015-04-13T20:20:51.573 回答
5
# First convert to pandas Period
period = pandas.tseries.period.Period(ordinal=int(d1), freq=ax.freq)
# Then convert to pandas timestamp
ts = period.to_timestamp()
# Then convert to date object
dt = ts.to_datetime()
于 2012-11-16T20:28:22.530 回答