R语言与地图（四）

Hello,

       大家好！欢迎收听 每周五分钟 与您分享 一点儿数据分析 的那些事。

作为用R语言绘制地图的最后一期，本周，我们将讲解如何利用 REmap 包以交互的形式来绘制迁徙图和热力图。

首先，简单介绍一下REmap包的三个函数：

(1). remapB(): 通过调用百度地图API绘制一个可以通过鼠标进行伸缩的地图；

(2). remapC(): 绘制迁徙图，路线图等；

(3). remapH(): 绘制中心辐射的热力图，应用于气象、人口密度、海拔测绘等领域。

本文，我们将通过模拟的数据用remapC() 和remapH()来说明某假期华师大闵行校区的小伙伴游玩的区域分布。

1. 模拟产生数据：

set.seed(20161126)

pop <- runif(19,0,1)

stu<-data.frame( district = mapNames('shanghai'), values = round(pop*100) )

stu<-stu[order(stu[,2],decreasing =T),]

head(stu)

如图所示，我们得知小伙伴最想去的前五个地区是黄浦区，浦东新区，崇明县，杨浦区和普陀区。接下来，我们用迁徙图和热力图加以展示。

2.绘制迁徙图：

       library(devtools)

        install_github('lchiffon/REmap') 

        library(REmap)

        markLine_data <- data.frame(origin=rep("华东师范大学闵行校区",5), destination=stu[1:5,1] )

        markPoint_data <- as.character(unlist(markLine_data[,2]))

        remapC(stu, maptype = '上海',

          title="模拟的迁徙图", subtitle="迁徙人数前5的目标省份",

          theme=get_theme( theme="Dark", titleColor = "#FFFFFF"),

          maxdata = 100,mindata=0,

          color=c('#FF0000','#C0C0C0'),   #颜色渐变:红－>灰

          markLineData=markLine_data,   

          markLineTheme=markLineControl( symbolSize = 6, effect =T, color="white", lineWidth=4, lineType="dashed" ),

          markPointData=markPoint_data,

          markPointTheme=markPointControl( symbol = 'heart', symbolSize=12, effect=T, effectType = 'bounce', color="white")

        )

3. 热力图：

library(baidumap)

ll<-matrix(0,19,2)

for(i in 1:19)  ll[i,]<- get_city_coord(as(stu$district[i],'character')) # 比 getCoordinate 精度高

data<-cbind(stu,ll)[,c(3,4,2)]

theme <- get_theme(theme ="none", titleColor ="red", backgroundColor ="white", borderColor ="blue", labelShow = T)

remapH(data, title = '模拟热力图',

maptype = '上海',

theme = theme,

blurSize = 60,

color = "darkred",

minAlpha= 1,

opacity = 1

)

4.其他工具：

        展现和分析时空数据的流行工具还包括：echarts、ArcGIS 和SaTScan等软件。

        附：color及所对应的十六进制编码图(百度里找的)：

          Okay, that's all. 谢谢您的收听，咱们下期再见，再见，再见～

如有兴趣，可以关注一下微信公众号：