在R中绘制带有重新调整的Mollweide投影的栅格

时间:2016-10-31 17:29:52

标签: r raster map-projections

我想在Mollweide投影中生成一个世界地图,但以经度150为中心。我想在上面绘制一个世界的光栅和矢量地图。

library(maptools)
library(raster)
library(rgdal)

# generate some dummy data
data(wrld_simpl)
ras <- raster(ext=extent(wrld_simpl), res=c(2,2))
values(ras) <- runif(ncell(ras))
ras <- mask(ras, wrld_simpl, inverse=TRUE)

# Here is the map unprojected, without recentering
plot(ras)
plot(wrld_simpl, add=TRUE, col='black')

# now I transform to Mollweide 
mollproj <- '+proj=moll +lon_0=150 +ellps=WGS84'

# project raster
ext <- projectExtent(ras, mollproj)
rasMoll <- projectRaster(ras, to=ext)

# project vector map
wrldMoll <- spTransform(wrld_simpl, CRS(mollproj))

# plot
plot(rasMoll)
plot(wrldMoll, add=TRUE)

enter image description here

这里有几个问题。地图不是完整的,矢量地图有水平线,并且有超出世界范围的光栅的浮动部分。

有关如何使其发挥作用的任何建议吗?

谢谢!

0 个答案:

没有答案