我想用ggplot::borders()
展示太平洋。我的问题是,我无法弄清楚如何向东和西方的美洲和亚洲展示澳大利亚和亚洲。
ggplot() + borders() + coord_cartesian(xlim = c(-290, -70))
在结果图中,亚洲和澳大利亚应该展示,但他们不会。
我想按照以下方式做点什么:
ggplot() + borders() + borders(x = long - 360)
获取这样的地图:
有人可以帮我解决这个问题吗?有关maps
和ggmap
个套餐的答案,但我找不到使用ggplot2::borders()
的答案。
任何答案的一个重要方面是lat-lon值是正确的。挑战在于:对于传统的经度坐标,太平洋的视图将从增加的正值(E半球)开始,然后在反射(W半球)处切换到减小的负值。如何在R中优先borders()
?
答案 0 :(得分:1)
类似的东西:
library(sp)
library(maps)
library(maptools)
library(ggplot2)
library(ggthemes)
world <- map("world", fill=TRUE, col="transparent", plot=FALSE)
worldSpP <- map2SpatialPolygons(world, world$names, CRS("+proj=longlat +ellps=WGS84"))
worldSpP <- worldSpP[-grep("Antarctica", row.names(worldSpP)),]
worldSpP <- worldSpP[-grep("Ghana", row.names(worldSpP)),]
worldSpP <- worldSpP[-grep("UK:Great Britain", row.names(worldSpP)),]
worldSpPnr <- nowrapRecenter(worldSpP)
world_map <- fortify(worldSpPnr)
gg <- ggplot()
gg <- gg + geom_map(data=world_map, map=world_map,
aes(x=long, y=lat, map_id=id),
color="black", fill="white", size=0.25)
gg <- gg + coord_map()
gg <- gg + theme_map()
gg
答案 1 :(得分:1)
我尝试缩放它以更贴近您发布的内容(下面的代码)。感谢this answer
kohske# install.packages("mapdata", dependencies = TRUE)
# install.packages("ggplot2", dependencies = TRUE)
library(ggplot2)
library(mapdata)
mp1 <- fortify(map(fill=TRUE, plot=FALSE))
mp2 <- mp1
mp2$long <- mp2$long + 360
mp2$group <- mp2$group + max(mp2$group) + 1
mp <- rbind(mp1, mp2)
ggplot(aes(x = long, y = lat, group = group), data = mp) +
geom_path() +
scale_x_continuous(limits = c(110, 300)) +
scale_y_continuous(limits = c(-50, 70))