“随着新冠肺炎疫情发展,目前国内疫情得到了初步控制,世界范围的疫情越来越严重,在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)