2014-01-20 27 views

回答

3

從包rgdal使用spTransform

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) 
+1

由於wrld_simpl的x和y座標是在-90至90和-180到180的範圍,而它不是在wrld_transf或任何其他投影不是'+ proj = latlon'的情況。要知道新地圖類型'bbox(wrld_transf)' – plannapus

+0

的座標範圍,否則您必須重新投影您的柵格! – plannapus

+0

如果你不熟悉重新投影背後的數學,我強烈建議不要重投影,並堅持經典的緯度/經度投影。 – plannapus

相關問題