如何用openair包按月份展示风玫瑰图的变化？

先处理数据 #

首先得把日期列转成R能识别的格式，不然openair没法正确按月份分组。用lubridate包的mdy()函数就行，因为你的日期是月/日/年的格式：

library(openair)
library(lubridate)

# 转换日期格式（假设你的数据框名为nitrogen）
nitrogen$Data <- mdy(nitrogen$Data)

按月份绘制风速风向玫瑰图 #

openair的windRose()自带分组功能，只要加个type="month"参数，就能直接生成12个月的分面风玫瑰图，不用手动拆分数据：

windRose(nitrogen,
         type = "month",  # 指定按月份分面
         key = list(header="Wind Rose Acri", footer="wind speed",
                    plot.style = c("ticks", "border"),
                    fit = "all", height = 1,
                    space = "top"),
         layout = c(3, 4))  # 设置为3行4列，刚好放下12个月的图

按月份绘制NO₂污染玫瑰图 #

同样的逻辑，给pollutionRose()加上type="month"参数即可：

pollutionRose(nitrogen, 
              pollutant = "no2",
              type = "month",
              layout = c(3, 4))

额外提示 #

• 如果想区分不同年份的同月（比如对比2015年1月和2016年1月的差异），可以把type设为c("year", "month")，会按年+月的组合分组展示
• layout参数可按需调整行列数，比如不想图太挤可以设为c(2,6)

内容的提问来源于stack exchange，提问作者Elena Iakimova