数据地图系列8|R语言版数据地图(下)

今天要跟大家分享的是数据地图系列的第八篇——R语言版数据地图(下),分省(市级)热力地图。

步骤与昨天分享的中国热力地图步骤基本一致,只是需要调用的数据文件和需要自定义的指标文件略有本不同。

R语言系统环境配置:

R version 3.2.3

RStudio Version 0.99.484

在正式开始之前,必须确保你的R语言环境中已经安装以下包:

ggplot2

plyr

maptools

如果还未安装,需要先安装:install.packages("ggplot2","plyr","maptools")

下面是正式作图步骤:

1、加载所需要的R包:

library(ggplot2)

library(plyr)

library(maptools)

2、读入作图数据:

CHN_adm2 <- readShapePoly("c:/rstudy/CHN_adm/CHN_adm2.shp") 读入数据文件

3、整合与转化数据

CHN_adm2_1 <- fortify(CHN_adm2) #转化为数据框

x <- CHN_adm2@data #读取行政信息

xs <- data.frame(x,id=seq(0:344)-1) #总共345行

china_map_data <- join(CHN_adm2_1, xs, type = "full") #合并形状数据与行政

4、抽取安徽省的整合数据

anhui <-subset(china_map_data,NAME_1=="Anhui") #取出安徽省的子集

5、读入指标数据

mydata=read.csv("c:/rstudy/anhui.csv") #打开自定义的指标文件。

以上指标数据需要你自己定义,使用excel文件创建,指标数据输入完成之后,另存为.CSV格式(存至c:/rstudy/下,这样便于引用,当然这些指标数据文件在昨天分享的数据文档里已经为你创建好了,你只需要修改数据就可以了)

anhui_data <- join(anhui, mydata, type="full") #合并地图数据与业务数据

6、运行代码:

ggplot(anhui_data, aes(x = long, y = lat, group = group,fill=zhibiao)) +

geom_polygon(colour="grey40") +

scale_fill_gradient(low="white",high="steelblue") +

coord_map("polyconic") +

theme(

panel.grid = element_blank(),

panel.background = element_blank(),

axis.text = element_blank(),

axis.ticks = element_blank(),

axis.title = element_blank()

) #【无标签】

以上代码与前一篇所讲到的代码几乎一样,只是具体数据不一样。

midpos <- function(x) mean(range(x,na.rm=TRUE))

centres <- ddply(anhui_data,.(city),colwise(midpos,.(long,lat)))

ggplot(anhui_data,aes(long,lat)) +

geom_polygon(aes(group=group,fill=zhibiao),colour="black") +

scale_fill_gradient(low="white",high="steelblue") +

geom_text(aes(label=city),data=centres) +

theme(

panel.grid = element_blank(),

panel.background = element_blank(),

axis.text = element_blank(),

axis.ticks = element_blank(),

axis.title = element_blank()

)#【带标签】

以上代码可以输出结果带有市级行政名称的带标签数据地图来。

当然,其实我们只需要修改以上所有代码的引用文件名称,可以跑出来任意一个省份的市级热力地图。

(前提是,你已经定义好每一个省份的指标数据)

事实上,我已经搜集整理了全国27个省级行政单位的热力地图代码(TXT)和数据文件(主要是指标文件需要自己定义,而且所有省份的指标文件已经在昨天分享的数据文件夹里了,地图数据文件都是引用的同一个数据源)。

如果你能够获取这些代码文件,你只需要用excel打开对应省份.CSV文件,替换,继续存为.CSV就可以直接跑我给你的代码了,直接出图(最好与我代码里给的路径一致)

原文发布于微信公众号 - 数据小魔方(datamofang)

原文发表时间:2016-06-20

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏SAP最佳业务实践

想学FM系列(13)-SAP FM模块:预算结构(4)-预算结构配置-层次结构

4)激活多级预算结构 ? 指定财务管理范围激活多级预算结构。如果激活后,在某一年度想要不激活多层级预算结构,需加一条配置,从该这一年度不激活,“激活”参数不选即...

4417
来自专栏Crossin的编程教室

【Pygame 第1课】 hello pygame

我们已经把python的基本内容讲得差不多了,所以从今天起,尝试一下新的方面:pygame -- 用python来写游戏。 pygame是一个python的游戏...

35510
来自专栏美团技术团队

人工智能在线特征系统中的数据存取技术

一、在线特征系统 主流互联网产品中,不论是经典的计算广告、搜索、推荐,还是垂直领域的路径规划、司机派单、物料智能设计,建立在人工智能技术之上的策略系统已经深入到...

8985
来自专栏UML

什么是实体关系图(ERD)?

实体 - 关系(ER)图(也称为ERD或ER模型)是Peter最初在1976年提出的经典且流行的概念数据模型。它是系统内不同实体的视觉表示以及它们如何相互关联。...

3645
来自专栏数据小魔方

数据地图系列7|R语言版(上)

今天跟大家分享的是数据地图系列的第七篇——使用R语言制作热力数据地图! 也许很多小伙伴儿对于R语言还很陌生,感觉很神秘。 确实,R语言的数据地图需要使用很长的代...

3894
来自专栏玉树芝兰

如何用VOSviewer分析CNKI关键词共现?

用VOSviewer尝试CNKI中文文献关键词共现(keyword co-occurence)分析时,你可能会踩到一个大坑。本文帮助你绕开这个坑,或是从坑里爬出...

1252
来自专栏FreeBuf

双尾蝎(APT-C-23):伸向巴以两国的毒针

? 摘要 2016 年 5 月起至今,双尾蝎组织(APT-C-23)对巴勒斯坦教育机构、军事机构等重要领域展开了有组织、有计划、有针对性的长时间不间断攻击。 ...

28310
来自专栏CSDN技术头条

介绍两款大数据清洗工具——DataWrangler、Google Refine

在进行数据分析和可视化之前,经常需要先“清洗”数据。这意味着什么?可能有些词条列表里是“New York City”,而其他人写成“New York,NY”。然...

6529
来自专栏吉浦迅科技

Titan V做计算真的这么不靠谱么?

今年3月份有一篇文章在高性能计算这个领域算是引起了一个不小的波动英伟达的 Titan V GPU 计算故障:2+2=4,呃=4.1,不,=4.3.....

6972
来自专栏逍遥剑客的游戏开发

一个困扰我一个多星期的Nebula3的BUG

1703

扫码关注云+社区

领取腾讯云代金券