如何重新投影地图，如“wrld_simpl”？

Question

要绘制全局地图，我只需使用：

library(maptools)
data(wrld_simpl)
plot(wrld_simpl)

似乎绘制的地图的投影是WGS84。如何将投影更改为另一个投影。

感谢任何提示

Answer 1

使用包spTransform中的rgdal：

library(rgdal)
library(maptools)
data(wrld_simpl)
wrld_transf <- spTransform(wrld_simpl, CRS(" +proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +a=6371228 +b=6371228 +units=m +no_defs "))
plot(wrld_transf)