r绘制中国地图 并展示流行病学数据_光环大数据培训机构.pdfVIP

光环大数据--大数据培训人工智能培训 R 绘制中国地图 并展示流行病学数据_光环大数据培训机构 流行病学的数据讲究“三间分布”，即人群分布、时间分布和空间分布。其 中的“空间分布”最好是在地图上展示，才比较清楚。R 软件集统计分析与高级 绘图于大成，是最适合做这项工作了。关于地图的绘制过程，谢益辉、邱怡轩和 陈丽云等人都早有文章讲述，开R 地图中文教程之先河。由于目前指导毕业论文 用到，因此研究了一下。本来因为网上教程很多，曾打消了写些文字的计划，但 怡轩版主鼓励说“教程者众，整合者鲜”，所以才战胜拖延症，提起拙笔综述整 合一下，并对DIY 统计GIS 地图提出了一点自己的想法。 1 地图GIS 数据的来源与R 绘制软件包 中国地图 GIS 数据的官方数据可以在国家基础地理信息中心的网站 （）里面可以免费下载。官方公开的数据包括：地图 数据，及居住地、交通、河流等辅助数据。今年6 月开始，官方正组织开始制作 新版数据。老数据暂时无法下载，读者要自行百度搜索，本文以旧版数据为例。 旧版地图数据中部分地名和地市区划已经过时，使用时需注意。 地图数据有 4 个压缩文件：bou1_4m.zip、bou2_4m.zip、bou3_4m.zip 和 bou4_4m.zip。bou 代表边界的意思，数字 1~4 代表国家、省、市、县的4 级行 政划分；4m 代表比例是 400 万分之一，这个比例的图形是公开的。每个文件解 压缩后含有两类文件：以字母p 结尾的表示多边形数据，用来绘制区域；以字母 l 结尾的文件是线形数据，用来绘制边界。但是老版数据中，市级数据中缺少绘 制区域的多边形数据，让市级分布图的绘制稍麻烦一些，新版中也许会有改进。 用R 绘制地图比较简单。比如画一下全国范围的区域，可以用如下代码： library(maptools)mydat = 光环大数据 光环大数据--大数据培训人工智能培训 readShapePoly(maps/bou1/bou1_4p.shp)plot(mydat) 但是，可以看出这样绘制的地图的形状有些扁平。这是因为，在绘图的过程 中，默认把经度和纬度作为普通数据，均匀平等对待，绘制在笛卡尔坐标系上造 成的。其实，地球的球面图形如何映射到平面图上，在地理学上是有一系列不同 的专业算法的。地图不应该画在普通的笛卡尔坐标系上，而是要画在地理学专业 的坐标系上。在这一点上，R 的ggplot2 包提供了专门的coord_map()函数。所 以推荐R 的ggplot2 包来绘制地图。 library(ggplot2)mymap = ggplot(data = fortify(mydat)) + fill = NA) + theme_grey()print(mymap + coord_map()) 这次中国地图的形状与百度地图一样了。 ggplot2包的coord_map 函数默认的映射类型是mercator。如果有其他需要， 可以使用其他的映射类型来绘制地图，如： mymap + coord_map(projection = azequidistant) coord_map 函数的映射类型及其含义可以通过下列代码查询帮助，一般我们 用默认的就可以。 library(mapproj)?mapproject 光环大数据 光环大数据--大数据培训人工智能培训 2 GIS 地图的数据结构及省市地图的绘制 GIS 地图有很多种存储格式，其中 shapefile 格式（.shp）可以通过R 的 maptools 包打开。其他格式可以去R 官网查询相应的软件包。 地图数据基本可以分为点、线、面三种数据，在maptools 包内分别有对应 的函数来读取（readShapePoints、readShapeLines 和 readShap