Chcę wyświetlić tę mapę w projekcji Robinsona:
library(ggmap)
world <- map_data("world")
ggplot() + geom_path(data = world,
aes(long, lat, group = group))
I chciałbym zmienić projekcję na „Robinson” (zgodnie z radą z odpowiedzi na moje poprzednie pytanie: Jakiej projekcji używa mapa globalnego regionu klimatycznego z Wikipedii?
Trudno mi było znaleźć domyślną implementację tej projekcji, pracowałem nad użyciem proj4
biblioteki:
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")
Próbowałem wielu podejść, w tym:
# 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)
ale żadna z nich nie działa. Czy to literówka, czy brakuje mi czegoś fundamentalnego w tej metodzie?