![[MacOS]](martinkoechy.de/Made_w_MacOS.gif)
A map of the world in cylindrical equal area projection with true scale at 10°N, the North Pole tilted to 50°N, 170°W and the map rotated 140° counter-clockwise.
The map is produced using R.
library("maps")
library("mapproj")
library("mapdata")
par(mar=c(1,0.1, 1, 0.1))
all=map(projection="cylequalarea", parameters=10, orientation=c(50, -170, -140),
wrap=TRUE, interior=FALSE, plot=FALSE, namesonly=TRUE)
# 50, -170, -140 # land left, Oceania right
# 50, -160, 100 # land right, Oceania left
#countries
par(fg="black", xpd=FALSE)
map(projection="", wrap=TRUE, bg="skyblue1", fill=TRUE, col="grey90",
interior=FALSE, fg="grey90")
map(region="Saudi",projection="", fill=TRUE, col="grey90", add=TRUE)
#surface
par(fg="grey90", bg="skyblue1")
map(region=all[which(all!="Lesotho")], projection="", exact=TRUE, wrap=TRUE,
fg="black", fill=TRUE, interior=FALSE, col="grey90", add=FALSE)
map(region="Saudi",projection="", fill=TRUE, col="grey90", interior=FALSE,
add=TRUE)
par(fg="black")
map(region=c("Caspian", "Aral", "Lake", "USSR:Lake Baikal"),projection="",
wrap=TRUE, fill=TRUE, col="skyblue1", interior=FALSE, boundary=TRUE, add=TRUE,
lwd=0.5)
map(region=all[which(all!="Lesotho")], projection="", exact=TRUE, wrap=TRUE,
fill=FALSE, boundary=TRUE, interior=FALSE, col="black", add=TRUE, lwd=0.5)
NAm=as.data.frame(map(regions=c("USA","Canada"), exact=TRUE, plot=FALSE)[1:2])
GreatLakes=na.omit(NAm[NAm$x>-92 & NAm$x< -75 & NAm$y>40 & NAm$y<49,])
polygon(mapproject(GreatLakes[!is.na(NAm$x),]), col="skyblue1", lwd=0.5)
my.rivers<-map('rivers', plot=FALSE, add=TRUE)
lines(mapproject(my.rivers), col="blue", lwd=0.5)
for (lon in seq(-160,180, by=20)){
m=mapproject(rep(lon,81),seq(-80,80,by=2))
m.min=which.min(m$x)
ma=data.frame(
x=m$x[c(1:(m.min-1),NA,m.min:length(m$x))],
y=m$y[c(1:(m.min-1),NA,m.min:length(m$x))]
)
lines(ma, col="yellow")
}
for (lat in c(-80, -66, -23, 0, 23, 66, 80)){
m=mapproject(seq(-175,180, by=5),rep(lat,72))
m.min=which.min(m$x)
m.max=which.max(m$x)
ma=data.frame(
x=m$x[c(m.min:length(m$x), 1:(m.min-1))],
y=m$y[c(m.min:length(m$x), 1:(m.min-1))]
)
lines(ma, col="yellow")
}
rect(-3.05,-1,3.05,1)
Page updated 2012-06-18.