私はtimeseries
50年間毎日降雨があります。冬季 (10 月から 11 月) の季節傾向分析を行っています。
import scikits.timeseries as ts
import datetime as dt
start_date = dt.datetime(1954,1,1)
end_date = dt.datetime(2003,12,31)
dateRange = ts.date_array(start_date,end_date,freq='D')
dly_series = ts.time_series(dates=dateRange,freq='D',data=rain)
# Autumn season OctNov
oct_series = dly_series[dly_series.month == 10]
nov_series = dly_series[dly_series.month == 11]
# October November months daily series
on_dly_series = ts.concatenate((oct_series,nov_series))
loc = np.where(on_dly_series.data >= 1.0)[0]
corr_dates = on_dly_series.dates[loc]
上の条件 を満たす 1 年の 10 月のすべての日降水量の合計を取得したいと考えてい> 1mm
ます。11月も同様。