Я хочу спроецировать эту карту в проекции Робинсона:
library(ggmap)
world <- map_data("world")
ggplot() + geom_path(data = world,
aes(long, lat, group = group))
И я хотел бы изменить проекцию на «Робинсон» (следуя совету из ответа на мой предыдущий вопрос: какую проекцию использует карта региона глобального климата из Википедии?
Мне было трудно найти реализацию этой проекции по умолчанию, я разработал следующее для использования proj4
библиотеки:
library(proj4)
robinson <- project(cbind(world$long, world$lat),
proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
Я пробовал несколько подходов, в том числе:
# using ggmap::get.map()
get_map("world", projection = mapprojection(robinson))
# using ggplot2::coord_map
coord_map(projection = robinson)
# and sp::coordinates:
library(sp)
coordinates(world) <- ~ lat + long
gridded(world) <- TRUE # returns error
proj4string(world) <- CRS(robinson)
но ничего из этого не работает. Это опечатка, или я упустил что-то фундаментальное об этом методе?