在matplotlib中将垂直线添加到Date格式化的时间序列

Question

我正在尝试将红色垂直线添加到时间序列图中，其中x轴格式为％Y-％m-％d。我想添加该行的日期是2013-05-14。只需在＆＃34; plt.show（）＆＃34;之前添加该行：

plt.axvline(x=2013-05-14)

或：

plt.axvline(x='2013-05-14')

返回错误：

RuntimeError: RRuleLocator estimated to generate 23972 ticks from 0044-05-12 23:59:59.999990+00:00 to 2013-06-07 00:00:00.000010+00:00: exceeds Locator.MAXTICKS * 2 (2000) 

这是函数，它运行良好：

 def time_series(self):
    fig = plt.figure(figsize=(20, 20), frameon = False)
    ax1 = fig.add_subplot(3, 1, 1)


    d_dates, d_flux_n2o, d_sem_n2o = np.loadtxt('%stime_series/dynamic.csv' %self.dir['r_directory'],  delimiter = ',', unpack=True, converters={0: mdates.strpdate2num('%Y-%m-%d')})

    ax1.set_xlabel('Date')
    ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d'))
    ax1.xaxis.set_major_locator(mdates.MonthLocator())
    ax1.xaxis.set_minor_locator(mdates.DayLocator())
    plt.gcf().autofmt_xdate()

    ax1.errorbar(d_dates, d_flux_n2o, yerr=d_sem_n2o, fmt="y-", linewidth=1.5, label = 'Biodynamic')
    ax1.legend(loc = 0)

    plt.show()

Answer 1

您必须为axvline方法提供数值，而不是字符串。您可以通过定义将日期的sting表示转换为datenums的转换器来实现此目的。您的np.loadtxt方法调用中已经有一个这样的转换器。如果将其定义为函数，则可以在加载数据时使用它，也可以在单个日期字符串中使用它。

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


def time_series(self):
    fig = plt.figure(figsize=(20, 20), frameon = False)
    ax1 = fig.add_subplot(3, 1, 1)

    # Converter to convert date strings to datetime objects
    conv = np.vectorize(mdates.strpdate2num('%Y-%m-%d'))


    d_dates, d_flux_n2o, d_sem_n2o = np.loadtxt('%stime_series/dynamic.csv' %self.dir['r_directory'],  delimiter = ',', unpack=True, converters={0: conv})

    ax1.errorbar(d_dates, d_flux_n2o, yerr=d_sem_n2o, fmt="y-", linewidth=1.5, label = 'Biodynamic')
    ax1.legend(loc = 0)
    ax1.axvline(conv('2013-05-14'), color='r', zorder=0)

    ax1.set_xlabel('Date')
    ax1.xaxis.set_major_formatter(mdates.DateFormatter('%Y-%m-%d'))
    ax1.xaxis.set_major_locator(mdates.MonthLocator())
    ax1.xaxis.set_minor_locator(mdates.DayLocator())
    plt.gcf().autofmt_xdate()

    plt.show()