# testing equal-area projections test in R
library(mapdata)
x11()
map("world")
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "default (rectangular?)")
x11()
map("world", projection="rectangular", parameter=0)
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "rectangular, parameter=0, wrap=FALSE")
x11()
map("world", projection="rectangular", parameter=0, wrap=TRUE)
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "rectangular, parameter=0, wrap=TRUE")
x11()
map("world", projection="mercator")
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Mercator, wrap=FALSE")
x11()
map("world", projection="mercator", wrap=TRUE)
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Mercator, wrap=TRUE")
x11()
map("world", projection="mercator", wrap=c(0,360))
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Mercator, wrap=c(0,360)")
x11()
map("world", projection="bonne", parameter=0, wrap=c(0,360))
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Bonne, parameter=0, wrap=c(0,360)")
x11()
map("world", projection="bonne", parameter=45, wrap=c(0,360))
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Bonne, parameter=45, wrap=c(0,360)")
x11()
map("world", projection="mollweide", wrap=c(0,360))
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Mollweide, wrap=c(0,360)")
x11()
map("world", projection="azequalarea")
gpa <- par()$usr
text((gpa[1]+gpa[2])/2, gpa[3], pos=3, "Azequalarea")
