2
繪製全球地圖,我簡單地使用:如何重新投影諸如「wrld_simpl」之類的地圖?
library(maptools)
data(wrld_simpl)
plot(wrld_simpl)
看來,繪製地圖的投影WGS84
。我怎樣才能改變投影到另一個投影。
感謝您的任何暗示
繪製全球地圖,我簡單地使用:如何重新投影諸如「wrld_simpl」之類的地圖?
library(maptools)
data(wrld_simpl)
plot(wrld_simpl)
看來,繪製地圖的投影WGS84
。我怎樣才能改變投影到另一個投影。
感謝您的任何暗示
從包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)
由於wrld_simpl的x和y座標是在-90至90和-180到180的範圍,而它不是在wrld_transf或任何其他投影不是'+ proj = latlon'的情況。要知道新地圖類型'bbox(wrld_transf)' – plannapus
的座標範圍,否則您必須重新投影您的柵格! – plannapus
如果你不熟悉重新投影背後的數學,我強烈建議不要重投影,並堅持經典的緯度/經度投影。 – plannapus