Tôi muốn chiếu bản đồ này trong dự đoán robinson:
library(ggmap)
world <- map_data("world")
ggplot() + geom_path(data = world,
aes(long, lat, group = group))
Và tôi muốn thay đổi phép chiếu thành "Robinson" (theo lời khuyên từ câu trả lời cho câu hỏi trước đây của tôi: Bản đồ khu vực khí hậu toàn cầu từ Wikipedia sử dụng phép chiếu nào?
Tôi đã có một thời gian khó khăn để tìm một triển khai mặc định của phép chiếu này, tôi đã tìm ra cách sau để sử dụng proj4
thư viện:
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")
Tôi đã thử một số cách tiếp cận, bao gồm:
# 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)
nhưng không ai trong số này làm việc Đây có phải là một lỗi đánh máy, hoặc tôi thiếu một cái gì đó cơ bản về phương pháp này?