我试图使用新的choroplethr包根据一些指标(在工作中经历过事故的人口的百分比)制作欧洲国家的等值线图。 这是一个可重现的例子,用于演示我需要解决的两个问题。
a<-c(4.1,2.5,0.4,6.4,1.4,1.8,3.8,1.3,2.3,8.4,5.2,1.9,0.8,1.5,2.1,1.2,3.8,1.4,3.1,0.8,4.0,1.3,4.8,2.6,2.8,2.3,3.1,2.5)
target<-c("austria","belgium","bulgaria","switzerland","cyprus","czech republic","denmark",
"estonia","spain","finland","france","greece","croatia","hungary","ireland","italy",
"lithuania","luxembourg","latvia","norway","poland","portugal","romania","sweden",
"slovenia","slovakia","turkey","united kingdom")
datas<-data.frame(region=target,value=a)
datas$region<-as.character(datas$region)
install.packages("choroplethr")
install.packages("choroplethrMaps")
library(choroplethr);library(choroplethrMaps)
data(country.map)
data(country.regions)
country_choropleth(datas,legend="%",num_colors=1,zoom=target)
1)我遇到的第一个深刻问题是属于某些国家(例如法国)的这些额外区域,并且不要让地图完全放大欧洲。
我给变量参数提供了数据框中存在的国家的整个向量。
如果我试图country_choropleth(datas,legend="%",num_colors=1,zoom="switzerland")
然后它只会放大瑞士
问题来自法国等国家
country_choropleth(datas,legend="%",num_colors=1,zoom="france")
因此,我需要一种方法来解除属于欧洲国家但不属于欧洲地区的区域。有什么想法吗?
2)我遇到的第二个问题是在country.regions$region
内部存在与包一起使用的可用区域,缺少马耳他。它在地图上只有一个印迹,但我仍然无法跳过它。
有没有办法添加一个国家?
答案 0 :(得分:2)
你应该提出这两个问题。您可以使用其他ggplot
层/美学来增加等值线pkg的输出。要“放大”并使用接近正确的投影,您可以执行以下操作:
library(ggplot2)
gg <- country_choropleth(datas,legend="%",num_colors=1,zoom=target)
gg <- gg + xlim(-31.266001, 39.869301)
gg <- gg + ylim(27.636311, 81.008797)
gg <- gg + coord_map("lambert", lat0=27.636311, lat1=81.008797)
gg
(并且您可以使用边界框坐标来进行更紧凑的操作)\