创建极坐标图但缺少数据

Create polar plots but missing data

我们在 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) : 

它需要 wdws 的风向和风速列,我们没有。

有人告诉我我们可以从 wunderground's api 中提取这些丢失的数据,但问题是:

  1. 如何从 wunderground 的 api 中提取数据来匹配我们上面的每一行数据?

  2. 天气数据似乎是每小时测量和记录一次,但我们的数据并非每小时记录一次,如您在上面所见。那么这将如何匹配?

有什么办法吗?

openair 软件包可以轻松访问英国空气质量监测站数据,包括伦敦的多个监测站。这些数据将自动包括风速和风向(ws 和 wd)。此功能由 openair 的 importAURNimportKCL 函数提供。

使用这些功能之一从您站点附近的监测站下载您感兴趣的时间段内的每小时数据集,并按日期(时间戳)将其与您的数据合并。 Timestamps(日期栏)在露天是一个POSIXct日期,通常为整整一个小时。您需要使用 as.POSIXct 将您的时间戳或人类日期转换为 POSIXct,并将结果列命名为 date。然后 round 你的日期精确到整小时,然后 date 与 AURN 数据集合并。

然后您可以根据风数据制作极坐标图,甚至可以将污染物测量值与城市自己的监测站数据进行比较。

我对伦敦的特定车站一无所知,但在 openair 手册的 importAURNimportKCL 功能下阅读了相关信息,或者在加载 openair 后在 R 中提供帮助。在 CRAN 上查看 openair,或在 github (https://github.com/davidcarslaw/openair) 上查看最新更新。 ps: 作者是伦敦空气质量专家。