[R] Create a map

fub2011 a.dobler at fu-berlin.de
Tue Oct 11 12:47:25 CEST 2011


hi,

I've just figured out how to plot (more) up-to-date poilitical borders in R
in an easy way:

You can get the outline file Borders_MWDB3 from the NASA panoply site
http://www.giss.nasa.gov/tools/panoply/overlays/

Then read it into R, get the indices with a jump over the -180/180 E line,
remove the first point of the part "on the other side" and finally plot the
boarders as line:

bord <- read.table("Borders_MWDB3.cno",sep=",",na.strings=9999,fill=T)
around <- which(abs(diff(bord[,1]))>180)+1
bord[around,] <- NA
plot(bord,type="l",xlab="degrees east",ylab="degrees north")

hope this works.

--
View this message in context: http://r.789695.n4.nabble.com/Create-a-map-tp3689877p3893581.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list