## import AQ data and add meta data to data frame
aq <- importAURN(site = c("kc1", "my1"), year = 2009)
meta <- importMeta(source = "aurn")
aq <- merge(aq, meta, by = "site")
## aggregate first before adding meta data (useful for many sites/years)
aq <- importAURN(site = c("kc1", "my1"), year = 2009)
meta <- importMeta(source = "aurn")
## calculate annual means
annual <- timeAverage(aq, avg.time = "year")
annual <- merge(annual, meta, by = "site")
## make a GoogleMapsPlot
GoogleMapsPlot(annual, pollutant = "no2")
Run the code above in your browser using DataLab