Тихоокеанская проекция Робинсона с ggplot в R

Question

Я пытаюсь создать тихо-ориентированную проекцию Робинсона в R с помощью ggplot...

преобразование шейп-файла в Robinson CRS с использованием spTransform и установка долготы на 150 градусов не работает, и fortify ошибки.

как это:

world_robin <- spTransform(world2, CRS("+proj=robin +lon_0=150 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))
plot(world_robin, col = "khaki", bg = "azure2")

И, кажется, нет способа сделать это позже в процессе установки проекции для ggmap: т.е.

coord_map("mollweide", orientation=c(90, 0, 150)) 

это кажется лучшим, но это не совсем там.

Answer 1

Вам понадобится GDAL / rgdal для этого и 180 против 150, но вы должны быть в состоянии экстраполировать. ГИС ТАК FTW!

library(ggplot2)
library(ggthemes)
library(sp)
library(rgdal)

# assumes you are in the ne_110m... directory
# split the world and stitch it back together again

system("ogr2ogr world_part1.shp ne_110m_admin_0_countries.shp -clipsrc -180 -90 0 90")
system("ogr2ogr world_part2.shp ne_110m_admin_0_countries.shp -clipsrc 0 -90 180 90")
system('ogr2ogr world_part1_shifted.shp world_part1.shp -dialect sqlite -sql "SELECT ShiftCoords(geometry,360,0), admin FROM world_part1"')
system("ogr2ogr world_0_360_raw.shp world_part2.shp")
system("ogr2ogr -update -append world_0_360_raw.shp world_part1_shifted.shp -nln world_0_360_raw")

world <- readOGR("ne_110m_admin_0_countries/world_0_360_raw.shp", "world_0_360_raw")
world_robin <- spTransform(world, CRS("+proj=robin +lon_0=180 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))

world_dat <- fortify(world_robin)

gg <- ggplot()
gg <- gg + geom_map(data=world_dat, map=world_dat,
                    aes(x=long, y=lat, map_id=id),
                    fill="khaki", color="black", size=0.25)
gg <- gg + coord_equal()
gg <- gg + theme_map()
gg <- gg + theme(plot.background=element_rect(fill="azure2"))
gg

Ты можешь использовать gdalUtils для этого также, но он просто вызывает системные двоичные файлы, и вы должны указать полигоны отсечения с хорошо известными строками (так что system вызов на несколько символов короче для строки).

К вашему сведению: это означает, что вы должны сдвинуть и спроецировать все точки / фигуры, прежде чем строить с помощью ggplot2.