Skip to content
Prev 324813 / 398502 Next

Map Antarctica

You have mixed your x with your y, try this:

 map("worldHires","Antarctica",ylim=c(-90,-60),xlim=c(-180,180),col="gray90",fill=TRUE)

I tend not to use maps/mapdata if I can help it (there are problems as
you can see, though there is support for projections with mapproject,
see ?map), if you want to use an alternative try this.

library(maptools)
data(wrld_simpl)
## select out just the continent
antarctica <-  wrld_simpl[wrld_simpl$NAME == "Antarctica", ]
plot(antarctica)
degAxis(1)
degAxis(2)

For a better display of this part of the world

library(rgdal)
## define a sensible projection
pr <- "+proj=laea +lat_0=-90 +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"
antarctica.laea <- spTransform(antarctica, CRS(pr))


plot(antarctica.laea)

Removing the dud coordinates at -90 is a bit more difficult, but can
be done (the first plot by map() has a similar problem).

There is a mailing list R-Sig-Geo which is more topical for questions
like this if you want to explore further.

Cheers, Mike.
On Wed, Jun 5, 2013 at 11:00 PM, ejb <ejb274 at psu.edu> wrote: