Saya ingin memproyeksikan peta ini dalam proyeksi robinson:
library(ggmap)
world <- map_data("world")
ggplot() + geom_path(data = world,
aes(long, lat, group = group))
Dan saya ingin mengubah proyeksi menjadi "Robinson" (mengikuti saran dari jawaban atas pertanyaan saya sebelumnya: Proyeksi apa yang digunakan peta wilayah iklim global dari Wikipedia?
Saya kesulitan menemukan implementasi default dari proyeksi ini, saya mengerjakan yang berikut untuk menggunakan proj4
perpustakaan:
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")
Saya telah mencoba sejumlah pendekatan, termasuk:
# 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)
tapi tidak ada yang berhasil. Apakah ini salah ketik, atau saya melewatkan sesuatu yang mendasar tentang metode ini?