[R-sig-Geo] How to overlay world map using non-standard projection
Robert Hijmans
r.hijmans at gmail.com
Thu Jan 20 19:22:31 CET 2011
>Problem is that I've been unable to overlay a map of the world. I've tried
to use the map() function, and >I've generated the map I want (right
projection and centered correctly) with the call:
>library(maps)
>world <-map(projection = "cylequalarea", parameter=30, interior=FALSE,
wrap=TRUE, ylim=c(->60,60),xlim=c(-180,180), orientation=c(90,-205,0))
>But when I try to overlay this map on the plot of my transformed
SpatialGrid I find that the scale of both >maps is completely different.
Their coordinates don't match:
Perhaps you need to provide more parameters to map() (ellipsoid, datum?), as
it appears that it now uses something equivalent to an ellipsoid with radii
of about ~1 meter
Robert
--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/How-to-overlay-world-map-using-non-standard-projection-tp5941681p5944824.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
More information about the R-sig-Geo
mailing list