风玫瑰图

library(openair)
mydata = read.csv(file='C:\\Users\\yutian\\Desktop\\wind.csv')#读取csv文件
#mydata$date = as.POSIXct(strptime(mydata$time,format='%Y/%m/%d %H:%M'))#告诉时间格式
#windRose(mydata)
#polarPlot(mydata,pollutant = "M_flux_total",x = "ws",wd = "wd")
#png(file='C:\\con.png',unit='cm',width = 10,height = 8,res=300)

#p=polarPlot(mydata,pollutant = "mBC",x = "ws",wd = "wd",key.header ='',
#          key.footer=expression(ug/m^3),main='BC mass conc.',key.position='right',
#          auto.text=FALSE,force.positive = FALSE)

p=polarPlot(mydata,pollutant = "M_flux_total",x = "ws",wd = "wd",key.header ='',
            key.footer=expression(ng/m^2/s),main='BC mass flux',key.position='right',
            auto.text=FALSE,force.positive = FALSE)
#dev.off()

#statistic = 'nwr'
#expression(ng/m^2/s)