0

我有timeseries50 年的每日降雨量。我正在对冬季(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]

我有兴趣获得满足上述条件的一年中所有 10 月每日降雨量的总和> 1mm。11 月份也是如此。

4

0 回答 0