0
 library("openair")
 library("plyr")

我有以下行代码,我需要用罗盘值(即 N=1000、NE-6859、W=9585 等)总结 WindRose 图的结果。该功能基于包openair。我不知道如何从图中提取每个风向桨的频率值。

windRose(viento,"WS_ms_Avg","WindDir",type="mes", paddle=FALSE,layout=c(3,4))

4

1 回答 1

1

这是你要找的吗?

tmp_1 <- windRose(mydata)

head(tmp_1$data)

##                           default Interval1 Interval2 Interval3 Interval4
## 1 01 January 1998 to 23 June 2005 0.9918888  4.560834  6.322132  6.838161
## 2 01 January 1998 to 23 June 2005 0.4665894  2.679027  4.344535  5.113944
## 3 01 January 1998 to 23 June 2005 0.4264195  2.479722  4.339900  4.956354
## 4 01 January 1998 to 23 June 2005 0.5561993  2.399382  3.589031  4.018540
## 5 01 January 1998 to 23 June 2005 0.6056392  2.717652  4.105060  4.863654
## 6 01 January 1998 to 23 June 2005 0.7122441  3.797605  7.341831 10.246427
##    wd calm panel.fun mean.wd freqs
## 1  30  0.1      4.49  -118.3  4426
## 2  60  0.1      4.49  -118.3  3310
## 3  90  0.1      4.49  -118.3  3208
## 4 120  0.1      4.49  -118.3  2601
## 5 150  0.1      4.49  -118.3  3148
## 6 180  0.1      4.49  -118.3  6632


tmp_2 <- windRose(mydata, type="year")

head(tmp_2$data)

##   year Interval1 Interval2 Interval3 Interval4  wd calm panel.fun mean.wd
## 1 1998 1.5962554  6.204993  7.057129  7.093135  30  0.2      4.38  -110.4
## 2 1998 0.5160826  3.300528  4.176668  4.200672  60  0.2      4.38  -110.4
## 3 1998 0.2760442  1.500240  2.472396  2.652424  90  0.2      4.38  -110.4
## 4 1998 0.5160826  2.244359  3.192511  3.480557 120  0.2      4.38  -110.4
## 5 1998 0.7561210  2.568411  3.552568  3.972636 150  0.2      4.38  -110.4
## 6 1998 0.9481517  3.984638  7.333173  9.841575 180  0.2      4.38  -110.4
##   freqs
## 1   591
## 2   350
## 3   221
## 4   290
## 5   331
## 6   820
于 2015-08-26T03:47:23.073 回答