转载

【R】用 ggplot2 绘制漂亮的分级统计地图

最近我一直尝试利用R绘制地图，我从网上找到了上百种不同的实现方法，然而其中却没有适用于我的数据的方法。最终，我从以下几个博客【1】中找到了灵感。我在整合这些资源的基础上，通过不断的试验和修正得到了一个较好地解决方案。这个方案就是本篇博文的主要内容。

本篇博文中展示了如何利用 ggplot2 来绘制分级统计地图，同时还介绍了如何更改图例、颜色等参数指标，以及如何导出图像文件。

数据预处理

绘制分级统计地图需要一些软件包，你最好确认你的电脑中已经安装并加载了它们。我们利用 maptools 库中的 readShapeSpatial() 函数来读取形状数据。形状数据是储存了经纬度等地图信息的空间向量数据。好消息是，你可以免费下载这些数据！比如，你可以在这个网站中【2】下载全球各个国家的形状数据。

我的项目是对印度的数据进行绘图分析，所以我下载了印度的形状数据并将它们导入到 R 中。

【R】用 ggplot2 绘制漂亮的分级统计地图

我们所导入的形状数据是一个空间类型的数据，我们可以检查它的变量名来观测它包含哪些内容。比如我们可以输出各个邦的名称：

【R】用 ggplot2 绘制漂亮的分级统计地图

现在我们已经得到目标数据。接下来我将会编造一些数据，你也可以导入 csv 格式的数据或者从你估计的模型中提取数据。最重要的一点是，新导入的数据中 id 号码必须和形状数据中的 id 一致，因为我们等会要合并这两个数据集。

【R】用 ggplot2 绘制漂亮的分级统计地图

[译者注]此处代码有误，原数据的ID编码为1287-1321，所以“id=1:num.states”应改为“id=1286+(1:num.states)”。

现在我们需要整合形状数据和目标数据集。首先，我们利用 ggplot 中的 fortify() 函数将形状数据转换成数据框格式的数据。该函数可以根据特定的识别变量将数据转化成数据框格式。

【R】用 ggplot2 绘制漂亮的分级统计地图

[译者注]此处显示的 ID 数据与原数据集的 ID 数据不一致。

我们可以看到这是一个常见的数据框格式的数据，其中列向量包含了经纬度数据。

接下来我们根据 ID 数据合并两个数据集，并保证空间数据集中没有缺失值。更重要的是，我们需要根据变量 order 对数据集进行排序。

【R】用 ggplot2 绘制漂亮的分级统计地图

基础绘图

我们准备利用 ggplot()、geom_polygon() 和 coord_map() 函数来绘图。

【R】用 ggplot2 绘制漂亮的分级统计地图

高阶绘图

上一部分只是绘图的开始，我们还可以通过调整其他参数使图形变得更漂亮。首先，我们可以利用 ggmap() 库中的 theme_nothing() 来删除背景和网格线，同时还需要设置 legend 的参数绘制图例。当然，我们还可以更改图例和整个图像的标题。

接下来，我们将配色改为“YlGn”（从黄色到绿色过渡）。你可以通过 display.brewer.all() 函数来观察具体配色情况。

我们也可以利用 scales 库中的 pretty_breaks() 函数来获取最佳分割点。如果你不喜欢 pretty_breaks() 的方法，你还可以通过 scale_fill_gradient() 来实现：scale_fill_gradient(name=”var”,limits=c(0,100),low=”white”,high=”red”)

【R】用 ggplot2 绘制漂亮的分级统计地图