等值线世界地图

我读了很多帖子和文章，但不断出现错误。我正在尝试制作一个等值线？使用我从全球恐怖主义数据库获得的数据绘制的世界地图。我想根据击杀次数或仅根据该国家/地区的攻击次数对国家进行着色。我现在不在乎。由于拥有数据的国家太多，因此绘制任何图表来显示这些数据是不合理的。

非常感谢您的帮助，如果我没有正确询问这个问题，我真诚地道歉，我正在学习这个网站的规则。

我的代码（到目前为止..）

library(maps)
library(ggplot2)
map("world")
world<- map_data("world")
gtd<- data.frame(gtd)
names(gtd)<- tolower(names(gtd))
gtd$country_txt<- tolower(rownames(gtd))
demo<- merge(world, gts, sort=FALSE, by="country_txt")

在 gtd 数据框中，国家列的名称是“country_txt”，所以我想我会使用它，但在 fix.by(by.x, x) 中出现错误：“by”必须指定唯一有效的列

如果那行得通，我会像我在一些网站上看到的那样进行绘图。 老实说，我已经在这方面工作了很长时间，并且阅读了很多代码/其他类似问题/网站/R 手册等。当涉及到 R 时，我会很乐意寻求一些帮助，我会接受我的无能。

像这样的东西吗？这是一个使用的解决方案rgdal and ggplot。我很久以前就放弃了使用 R 基础来做这类事情。

library(rgdal)        # for readOGR(...)
library(RColorBrewer) # for brewer.pal(...)
library(ggplot2)
setwd(" < directory with all files >")

gtd        <- read.csv("globalterrorismdb_1213dist.csv")
gtd.recent <- gtd[gtd$iyear>2009,]
gtd.recent <- aggregate(nkill~country_txt,gtd.recent,sum)
world      <- readOGR(dsn=".",
                      layer="world_country_admin_boundary_shapefile_with_fips_codes")

countries <- world@data
countries <- cbind(id=rownames(countries),countries)
countries <- merge(countries,gtd.recent, 
                   by.x="CNTRY_NAME", by.y="country_txt", all.x=T)
map.df <- fortify(world)
map.df <- merge(map.df,countries, by="id")
ggplot(map.df, aes(x=long,y=lat,group=group)) +
  geom_polygon(aes(fill=nkill))+
  geom_path(colour="grey50")+
  scale_fill_gradientn(name="Deaths",
                       colours=rev(brewer.pal(9,"Spectral")),
                       na.value="white")+
  coord_fixed()+labs(x="",y="")

全球恐怖主义数据库有多个版本。我使用了可用的完整数据集here http://www.start.umd.edu/gtd/contact/，然后对 > 2009 年进行子集化。因此，此地图显示了 2010-01-01 至 2013-01-01 期间按国家/地区划分的恐怖主义造成的总死亡人数（此来源提供的最后数据）。这些文件可作为 MS Excel 下载，我将其转换为 csv 以导入到 R 中。

The 世界地图 http://geocommons.com/overlays/33578.zip可以作为 shapefile 从地理共享网站 http://geocommons.com/.

制作分区统计图的棘手部分是将数据与正确的多边形（国家/地区）相关联。这通常是一个四步过程：

1. 在 shapefile 属性表中查找一个映射（无双关语）到数据中相应字段的字段。在本例中，shapefile 中的字段“CNTRY_NAME”似乎映射到 gtd 数据库中的字段“country_txt”。
2. 创建多边形 ID（存储在属性表的行名称中）和 CNTRY_NAME 字段之间的关联。
3. 使用 CNTRY_NAME 和country_txt 将结果与您的数据合并。
4. 将其结果与使用创建的数据框合并fortify(map)- 这将多边形与死亡联系起来（nkill).