# NOT RUN {
Month <- SO2$month + (SO2$year - 1990) * 12
Year <- SO2$year + (SO2$month - 0.5) / 12
Location <- cbind(SO2$longitude, SO2$latitude)
back <- I
if (require(maps)) {
mapxy <- map('world', plot = FALSE,
xlim = range(SO2$longitude), ylim = range(SO2$latitude))
back <- function() map(mapxy, add = TRUE)
}
rp.plot4d(Location, Year, SO2$logSO2, col.palette = rev(heat.colors(12)),
background.plot = back)
# }
Run the code above in your browser using DataLab