matplotlibの年月をティックマークとして使用してx軸を設定する方法は？

Question

私は次元[年、月]を持つnetcdfファイルから2D変数XVARを持っています。私はフラット化されたXVAR（長さnyear * nmonthの1D配列）をプロットし、x軸を次のように設定したい：主要なティックの年とマイナーティックの月。難しいのは、毎月のステップで1dアレイを作成する方法を知らないことです。私が使用できるmonthdeltaメソッドはありません（ただし、毎月の日数が異なるため理由が分かります）。

デルタ=ですか？デルタ= relativedelta.relativedelta（月= 1）を試しましたが、「オブジェクトに属性「total_seconds」がありません」というエラーが表示されました。これは完全に理解できません。

import numpy as np 
from netCDF4 import Dataset 
import matplotlib.pyplot as plt 
import matplotlib.dates as mdates 
from mpl_toolkits.basemap import Basemap 
from datetime import date, timedelta 

ncfile = Dataset('filepath',mode='r') 
XVAR4d = ncfile.variables['XVAR'][:] 
XVAR2d = np.nanmean(XVAR4d,axis=(2,3)).flatten() 
yrs = ncfile.variables['YEAR'] 

stt = date(np.min(yrs),1,1) 
end = date(np.max(yrs)+1,1,1) 
delta = ? 
dates = mdates.drange(stt,end,delta) 

years = mdates.YearLocator() # every year 
months = mdates.MonthLocator() # every month 
yearsFmt = mdates.DateFormatter('%Y') 

fig = plt.figure() 
ax1 = fig.add_subplot(211) 
ax1.xaxis.set_major_locator(years) 
ax1.xaxis.set_major_formatter(yearsFmt) 
ax1.xaxis.set_minor_locator(months) 
ax1.set_xlim(stt,end) 
ax1.plot(dates,xvar2d,c='r') 

Answer 1

私は2番目のコメントのアイデアが気に入ったので、実際の提案された回答に変えています。

代わりのdrangeを使用して、日付を自分で作成します。drange`は `datetime.timedelta`を取る`よう

totalMonths = 12*(np.max(yrs) - np.min(yrs)+1) 
dates = mdates.date2num([date(np.min(yrs)+(i//12),i%12+1,1) for i in range(totalMonths)])