1
我正在使用okmesonet
包獲取有關降雨量的數據。我試過使用這個軟件包中的avgokmts
來計算每天的降雨量,但我得到的是非感性的值。avgokmts從ok mesonet數據返回不正確的最大雨量
獲取降雨數據諾曼,OK(以5分鐘間隔毫米累計降雨超過一天)
library(okmesonet)
rainDat <- okmts(begintime="2016-06-21 00:00:00", endtime="2016-07-04 00:00:00",
station="NRMN", variables="RAIN", localtime=TRUE)
計算每它返回這些值一天
avgokmts(rainDat, by="day", metric="max")
最大雨
STID STNM DAY MONTH YEAR RAIN Time Date
1 NRMN 121 21 06 2016 0.00 23:55:00 2016-06-22
2 NRMN 121 22 06 2016 0.25 23:55:00 2016-06-23
3 NRMN 121 23 06 2016 59.70 23:55:00 2016-06-24
4 NRMN 121 24 06 2016 0.00 23:55:00 2016-06-25
5 NRMN 121 25 06 2016 0.00 23:55:00 2016-06-26
6 NRMN 121 26 06 2016 0.00 23:55:00 2016-06-27
7 NRMN 121 27 06 2016 0.00 23:55:00 2016-06-28
8 NRMN 121 28 06 2016 0.00 23:55:00 2016-06-29
9 NRMN 121 29 06 2016 0.00 23:55:00 2016-06-30
10 NRMN 121 30 06 2016 28.19 23:55:00 2016-07-01
11 NRMN 121 01 07 2016 0.00 23:55:00 2016-07-02
12 NRMN 121 02 07 2016 0.51 23:55:00 2016-07-03
13 NRMN 121 03 07 2016 0.00 23:55:00 2016-07-04
14 NRMN 121 04 07 2016 0.00 00:00:00 2016-07-04
但是,這些降雨量值非常明顯地與降雨量不匹配,如下圖所示(高峯期出現在6月27日和7月3日)。
plot(rainDat$TIME, rainDat$RAIN, xlab="Date", ylab="Cumulative Daily Rain (mm)")
爲什麼不avgokmts
在這種情況下工作?我是如何調用函數的錯誤?有沒有其他方法可以使用這個數據集來計算每日降雨量?