在 R 中绘制 1990 年前的世界地图
Plot pre-1990 world map in R
我正在处理 1991 年之前的一些全球数据,因此在苏联、南斯拉夫和捷克斯洛伐克分裂之前。我想使用 rworldmap 或地图绘制数据,但该包似乎只有现代世界地图可以轻松访问。所有 1991 年之前的国家/地区都显示为空白,并且边界将其 post-1991 对应国家分开。
此代码生成历史地图:
if (requireNamespace("mapdata", quietly=TRUE) && packageVersion("mapdata") >= "2.3")
{map("mapdata::worldLores", fill = TRUE, col = 1:10)}
编辑:此外,根据下面的有用评论,历史地图 shapefile 很容易从以下位置获得:
library(cshapes)
cshp.data<-cshp(as.Date("1990-01-01"))
plot(cshp.data)
但我不知道是否可以将其与 rworldmap 函数结合使用...或者我是否必须弄清楚如何使用 maps 包,这似乎工作方式不同。 (或者也许有一个 ggplot 解决方案?)
我目前使用的rworldmap代码(获取现代地图)是:
#make example data including Soviet Union
country <- as.vector(c("Afghanistan","Australia","Iceland","Soviet Union",
"Zimbabwe"))
value <- as.vector(c(5,10,100,10,50))
df<-data.frame(country,value)
#make map
map1 <- joinCountryData2Map(df, joinCode = "NAME", nameJoinColumn =
"country")
mapCountryData( map1, addLegend=F, catMethod="fixedWidth",
nameColumnToPlot="value" )
#...Soviet Union is blank
啊哈,有一个使用 mapdata 包中的旧地图的 ggplot 解决方案:
library(ggplot2)
library(dplyr)
library(mapdata)
df<-data.frame(country=c("Afghanistan","Australia","Iceland","USSR","Zimbabwe"),
value=c(5,10,100,10,50),stringsAsFactors=FALSE)
WorldData <- map_data('worldLores') #use the old map
WorldData <- fortify(WorldData)
mapped <- ggplot() +
geom_map(data=WorldData, map=WorldData,
aes(x=long, y=lat, group=group, map_id=region),
fill="white", colour="#7f7f7f", size=0.5) +
geom_map(data=df, map=WorldData,
aes(fill=value, map_id=country),
colour="#7f7f7f", size=0.5)
mapped
(映射代码借自 post,加油@hrbrmstr)
我正在处理 1991 年之前的一些全球数据,因此在苏联、南斯拉夫和捷克斯洛伐克分裂之前。我想使用 rworldmap 或地图绘制数据,但该包似乎只有现代世界地图可以轻松访问。所有 1991 年之前的国家/地区都显示为空白,并且边界将其 post-1991 对应国家分开。
此代码生成历史地图:
if (requireNamespace("mapdata", quietly=TRUE) && packageVersion("mapdata") >= "2.3")
{map("mapdata::worldLores", fill = TRUE, col = 1:10)}
编辑:此外,根据下面的有用评论,历史地图 shapefile 很容易从以下位置获得:
library(cshapes)
cshp.data<-cshp(as.Date("1990-01-01"))
plot(cshp.data)
但我不知道是否可以将其与 rworldmap 函数结合使用...或者我是否必须弄清楚如何使用 maps 包,这似乎工作方式不同。 (或者也许有一个 ggplot 解决方案?)
我目前使用的rworldmap代码(获取现代地图)是:
#make example data including Soviet Union
country <- as.vector(c("Afghanistan","Australia","Iceland","Soviet Union",
"Zimbabwe"))
value <- as.vector(c(5,10,100,10,50))
df<-data.frame(country,value)
#make map
map1 <- joinCountryData2Map(df, joinCode = "NAME", nameJoinColumn =
"country")
mapCountryData( map1, addLegend=F, catMethod="fixedWidth",
nameColumnToPlot="value" )
#...Soviet Union is blank
啊哈,有一个使用 mapdata 包中的旧地图的 ggplot 解决方案:
library(ggplot2)
library(dplyr)
library(mapdata)
df<-data.frame(country=c("Afghanistan","Australia","Iceland","USSR","Zimbabwe"),
value=c(5,10,100,10,50),stringsAsFactors=FALSE)
WorldData <- map_data('worldLores') #use the old map
WorldData <- fortify(WorldData)
mapped <- ggplot() +
geom_map(data=WorldData, map=WorldData,
aes(x=long, y=lat, group=group, map_id=region),
fill="white", colour="#7f7f7f", size=0.5) +
geom_map(data=df, map=WorldData,
aes(fill=value, map_id=country),
colour="#7f7f7f", size=0.5)
mapped
(映射代码借自