title: "bomage" author: "Michael Sumner" date: "December 21, 2015" output: html_document
```{r setup, include=FALSE} knitr::opts_chunk$set(echo = TRUE)
## Get images from the BOM
```R
r <- raster("dev/www.bom.gov.au/charts_data/IDY20001/current/windarrow/10m/2015-12-21/IDY20001.windarrow-10m.066.png")
a <- regions("windAntarctica")
projection(r) <- a$proj
extent(r) <- a$extent
plot(r)
library(maptools)
data(wrld_simpl)
library(rgdal)
w <- spTransform(wrld_simpl, projection(r))
plot(w, add = TRUE)
ex <- domath(pts, rawxy)
gdalvrt(r, a_ullr = c(xmin(ex), ymax(ex), xmax(ex), ymin(ex)),
a_srs = "\"+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0\"")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.