我们在csv文件中包含以下数据,我们希望从中创建极坐标图。我们将使用这个R包 - openair来创建极坐标图。
timestamp humandate NOppb
1 1412877113 09/10/2014 13:51 19
2 1412876508 09/10/2014 13:41
3 1412876508 09/10/2014 13:41
4 1412877118 09/10/2014 13:51 17
....
但是,我们缺少使用polarPlot()
,
# Load package.
library("openair")
polarPlot(dat, pollutant = "NOppb", na.rm = TRUE)
结果:
Can't find the variable(s) wd ws
Error in checkPrep(mydata, vars, type, remove.calm = FALSE) :
它需要 wd 和 ws 的列,以获取我们没有的风向和速度。
我被告知我们可以从wunderground's api中提取这些缺失的数据,但问题是:
如何从wunderground的api中提取数据以匹配上面数据的每一行?
天气数据按小时计算和记录,但我们的数据不会按小时记录,因为您可以在上面看到。那怎么会匹配呢?
任何想法我该怎么办?
答案 0 :(得分:2)
openair
套餐可以方便地访问英国空气质量监测站数据,包括伦敦的几个站点。这些数据将自动包括风速和风向(ws和wd)。此功能由openair的importAURN
和importKCL
函数提供。
使用这些功能之一从您站点附近的监控站下载您感兴趣的一段时间内的每小时数据集,并按日期(时间戳)将其与您的数据合并。在openair中的时间戳(日期列)是一个POSIXct
日期,通常是整整一个小时。您需要使用POSIXct
将时间戳或humandate转换为as.POSIXct
,并将结果列命名为date
。然后round
您的日期到最近的整个小时,然后将date
与AURN数据集合并。
然后,您可以根据风数据制作极坐标图,甚至将污染物测量结果与城市自己的监测站数据进行比较。
我对伦敦的具体电台一无所知,但在开放手册中的importAURN
和importKCL
功能下阅读此内容,或在装载露天广场后在R中提供帮助。请参阅CRAN上的openair,或github上的最新更新(https://github.com/davidcarslaw/openair)。 ps:露天作者是伦敦空气质量专家。