Как реализовать равные проекции Земли на карте и растровых в R?

Question

Интересно, кто-нибудь знает, как реализовать Равную проекцию Земли для карты и растра в R?

Я нашел этот урок:

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.