导入数据:
library(openair)
setwd("文件存储地址如:C:/Program Files/RStudio/R")
Sys.setlocale("LC_TIME","English")
mydata=read.csv("文件名.csv",header=TRUE)
mydata[mydata==-99]<-NA
mydata$date<-as.POSIXct(strptime(mydata$date,format="%Y/%m/%d",tz="GMT"))
head(mydata)
tail(mydata)
做污染日历图
calendarPlot(mydata,pollutant="level",year=2019,month=c(4),breaks=c(-1,0,1,2,3,4,5,6),labels=c("无效","优","良","轻度污染","中度污染","重度污染","严重污染"),cols=c("gray","limegreen","yellow","orange","red","darkmagenta","darkred"),main="某地区某时间空气质量等级月历",w.shift=2)
给你一个例子,我把数据源文件中AQI无效的都改成了-1,做的是4月份的污染日历图,效果如下: