如何在 R 中的地图和栅格中实现等地球投影？

Question

我想知道是否有人知道如何在 R 中为地图和栅格实现 Equal Earth Projection？

我找到了这个教程：

https://rud.is/b/2018/09/28/quick-hit-using-the-new-equal-earth-projection-in-r/

但我无法复制它。似乎“+proj=eqearth”不可用，但现在看来应该可用。

library(rgdal)
library(maps)
library(maptools)


RobinsonProj="+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
EqualEarthProj= "+proj=eqearth +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
LatLonProj="+proj=longlat +datum=WGS84"

worldmap <- maps::map("world", plot=F,interior = F,wrap=c(-180,180))   #create map of the world
worldmapLines <- map2SpatialLines(worldmap, proj4string=CRS(LatLonProj))

#reproject
RobinsonMap <- spTransform(worldmapLines, CRS(RobinsonProj))
LatLonMap <- spTransform(worldmapLines, CRS(LatLonProj))
EqualEarthMap <- spTransform(worldmapLines, CRS(EqualEarthProj))

Answer 1

那篇博文的开头是“我注意到 proj4 已更新到 5.2”

但我假设你在 windows 上，当你这样做时

library(rgdal)

你得到

# Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]

这就解释了。还不是 5.2 版。