こんにちは、netCDFファイルから抽出したいくつかのデータから次のプロットを作成しました。 Plot created。夜間配列を定数で除算して曜日を作る
# import netCDF
import matplotlib.pyplot as plt
import numpy as np
import netCDF4
from pylab import rcParams
f = netCDF4.Dataset('data', 'r')
v = f.variables['store_Bio'][0:27, 0:23, 0]
# Reshape array
r = np.reshape(v, 621)
# plot asthetics
plt.plot(r, linestyle='-', linewidth=5.0, c='c')
plt.xlabel('Date')
plt.ylabel('Average Biogenic Emissions (kg/hour/per cell area)')
plt.title('Average Biogenic Emissions March 2013')
plt.ylim(0, 350)
plt.xlim(0, 621)
# Annotate plot
plt.annotate('Daytime Maximum', xy=(72, 140), xytext=(24, 200),
arrowprops=dict(facecolor='black', shrink=0.05,
width=0.1, headwidth=10))
# change plot size
rcParams['figure.figsize'] = 25, 5
# change ticks
plt.xticks(range(0, 621, 12), [str(i) for i in range(0, 621, 12)])
plt.yticks(range(0, 350, 50), [str(i) for i in range(0, 350, 50)])
plt.show()
は、しかし、私はx軸は日数の代わりに、時間を表示し、他のすべてが同じままであるように、私は24で日付だけを分割する可能性がどのように思っていた:私は、次のスクリプトを使用していました。どんな助けでも大歓迎です!
日付は 'r'の最初の列にあると思いますので、' date = r [: 、0]/24.'、 'plot(date、r [:、1])'となります。 –
コメントありがとう!私はこれを試して、 '配列のインデックスが多すぎます'というエラーを取得します。 –