2013-07-18 1 views
3

Я пытаюсь добавить красную вертикальную линию к графику временного ряда, где ось x отформатирована как% Y-% m-% d. Дата, когда я хочу добавить строку, - 2013-05-14. Просто добавив строку перед "plt.show()":Добавление вертикальной строки в формат даты в формате matplotlib

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

ответ

6

Вы должны отправить axvline метод числовое значение, а не строка. Вы можете достичь этого, указав конвертер, который преобразует представления с указанием дат в 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() 

enter image description here