knitr::opts_chunk$set( comment = "#>", collapse = TRUE, warning = FALSE, message = FALSE )
Get sea ice data at ftp://sidads.colorado.edu/DATASETS/NOAA/G02135/shapefiles.
library('rnoaa') library('plyr')
urls <- sapply(seq(1979, 1990, 1), function(x) seaiceeurls(yr = x, mo = 'Feb', pole = 'S')) out <- lapply(urls, seaice) head(out[[1]])
urls <- seaiceeurls(mo = 'Apr', pole = 'N', yr = 1990) out <- seaice(urls) library('ggplot2') ggplot(out, aes(long, lat, group = group)) + geom_polygon(fill = "steelblue") + theme_ice()
urls <- seaiceeurls(mo = 'Apr', pole = 'N')[1:10] out <- lapply(urls, seaice) names(out) <- seq(1979, 1988, 1) df <- ldply(out) library('ggplot2') ggplot(df, aes(long, lat, group = group)) + geom_polygon(fill = "steelblue") + theme_ice() + facet_wrap(~.id)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.