arcgis中如何导出奥维可识别的图_基于ArcGIS和R语言的世界疫情分析

“随着新冠肺炎疫情发展，目前国内疫情得到了初步控制，世界范围的疫情越来越严重，在2月28日，世卫组织将新冠肺炎疫情风险从“高”，调整为“非常高”，新冠肺炎已经开始在世界范围内扩散，防控新冠肺炎，成为了全人类的共同任务 ”

本文使用ArcGIS和R语言，对世界新冠肺炎疫情进行分析。

01

ArcGIS制作世界疫情分布图

—

先看一下成图，以新增确诊病例为例，制作世界各国新增确诊病例分布图。

然后看一下整理的数据，数据分成两部分，疫情属性数据和世界矢量数据，疫情表格中，有“国家”字段，和矢量数据中“CH_NAME”字段对应，可以通过挂接的方式，将两组数据连接起来。

数据挂接完成后，就可以对数据进行符号化了，在这里采用了分级设色的方法，对每个国家的新增病例进行符号化。

然后对国家名称和新增病例数进行标注，具体方法不再赘述，请参阅：

如何使用ArcGIS制作疫情分布图

如何制作地图？以ArcGIS勘测定界成果制作说明为例

这样，在ArcGIS中的制图工作就完成了。

02

R语言进行世界疫情分析

—

重点在这一部分。

首先加载所需的一些程序包

读取EXCEL文件用的readxl包

时间序列用的lubridate包

绘图神器ggplot2包

格网grid包

数据处理dplyr包

library(readxl)library(lubridate)library(ggplot2)library(grid)library(dplyr)

读取数据，由于数据有csv格式的，也有xls格式的，分成两段读取，一部分读csv，一部分读xls，在这里我将csv文件和xls文件分别放在了worldDay1和worldDay2文件夹下。为了生成时间序列，我最后使用ymd函数生成了一个R语言可以识别的时间序列。

a = list.files("worldDay1")                                             #list.files命令将input文件夹下所有文件名输入adir = paste("./worldDay1/",a,sep="")                                      #用paste命令构建路径变量dirn = length(dir)