2

http://www.r-bloggers.com/mapping-the-worlds-biggest-airlines/にあるチュートリアルに従おうとしています。チュートリアルに従ってプロットすることはできますが、ロビンソン図法にも変換したいと思います。readOGR コマンドを実行した後、ne_10m_populated_places.shp でエラーが発生する可能性があることに気付きました。誰かが私が空間形状で何をしなければならないかを理解するのを手伝ってくれませんか?

library(ggplot2)
library(maps)
library(rgeos)
library(maptools)
gpclibPermit()
library(rgdal)

worldmapin<- readShapePoly("ne_10m_admin_0_countries.shp")
worldmap_ogr <- readOGR(".", "ne_10m_admin_0_countries")
worldRobinson <- spTransform(worldmap_ogr, CRS("+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))

urbanareasin<- readShapePoly("ne_10m_urban_areas.shp")
urbanareas_ogr <- readOGR(".", "ne_10m_urban_areas")
urbanareasRobinson <- spTransform(urbanareas_ogr, CRS("+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))
simp<- gSimplify(urbanareasRobinson, 10, topologyPreserve=T)
urbanareas<-fortify(simp)

placesin<- readShapeSpatial("ne_10m_populated_places")
placesin_ogr <- readOGR(".", "ne_10m_populated_places")
placesinRobinson <- spTransform(placesin_ogr, CRS("+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"))

base<- ggplot(worldmapin)
wrld<-c(geom_polygon(aes(long,lat,group=group), size = 0.1, colour= "#090D2A", fill="#FFFFFF", alpha=1, data=worldmap))urb<- c(geom_polygon(aes(long,lat,group=group), size = 0.1, colour= "#FCFFF1", fill="#FCFFF1", alpha=1, data=urbanareas))
places<-geom_text(data=placesinRobinson, aes(LONGITUDE, LATITUDE, label = NAME), size=1)

base+wrld+urb+places
4

0 に答える 0