Doug Nychka’s fields package in R is really excellent. Among its nice functions is a command world() that draws an outline of the continents. I wrote him today inquiring whether there was an option in which land could be colored. Unfortunately there wasn’t. However by the end of the day, there was. Doug sent me a function world.color() which yielded the following pretty map.
It will be included in the next release of fields, but in the meantime CA readers can download the function here
library(fields) #load the package
source(“http://data.climateaudit.org/scripts/world.color.R”) #download function
You can handle this like the world() function where you can use xlim and ylim delimiters as in a plot function; you can use the add=T option if you want to overlay on a plot of points to make a location map. Thanks, Doug.