空汙視覺化

說明書

http://www.openair-project.org/PDF/OpenAir_Manual.pdf

安裝套件 OpenAir

install.packages("openair") //安裝OpenAir套件

library(openair) //開啟套件

OpenAir內有一內建資料mydata，是空汙資料，可供練習

Airdata <- as.data.frame(mydata)

View(Airdata)

建立日曆熱圖

calendarPlot(Airdata, pollutant = "pm25" , year = 2004) //PM2.5在2004年的日曆熱圖

data(mydata) ## make sure openair 'mydata' loaded fresh

calendarPlot(mydata, pollutant = "pm10", annotate = "value", lim =50,

cols = "Purples", col.lim = c("black", "orange"),

layout = c(4, 3) , year = 2004)

calendarPlot(mydata, pollutant = "o3", year = 2003, annotate = "ws")

timePlot(Airdata , pollutant = c("pm25" , "o3" , "ws"))

timeVariation(Airdata , pollutant = c("pm10" , "pm25" , "o3"))

polarFreq(mydata, pollutant = "pm25", type = "year", statistic = "weighted.mean", min.bin = 2)