绘图杂记【5】R语言绘制热力地图（重庆热力地图）

最新推荐文章于 2022-12-01 10:44:41 发布

卖山楂啦prss

最新推荐文章于 2022-12-01 10:44:41 发布

阅读量7.5k

点赞数 10

分类专栏：数据可视化

本文链接：https://blog.csdn.net/qq_42374697/article/details/110476955

版权

数据可视化专栏收录该内容

28 篇文章 25 订阅

订阅专栏

注：文章中出现的地图不可以进行文章、期刊等出版，仅限学习使用。

安装和加载包

install.packages(mapdata)
install.packages(maptools)
install.packages(ggplot2)
install.packages(plyr)
install.packages(rgdal)
install.packages(mapproj)


library(ggplot2)
library(maps)
library(mapdata)
library(maptools)
library(plyr)
library(rgdal)
library(mapproj)

快速绘图

library(maps)
map()  # 世界地图


map("state") # 美国各个州


map("county") # 各个州更细分

读取 shp 文件

绘制热力图

china_map <- rgdal::readOGR("bou2_4p.shp")
plot(china_map)

在这里插入图片描述
图像过于扁平，需要进行调整

ggplot(
	china_map,
	aes(x=long,y=lat,group=group)) +
geom_polygon(
	fill="white",
	colour="grey") +
coord_map("polyconic")

在这里插入图片描述
或者

ggplot(
	china_map, 
	aes(x = long, y = lat,group = group)) +
geom_polygon(
	fill = "white", 
	colour = "black") +
coord_map("polyconic") +
theme(
   panel.grid = element_blank(),
   panel.background = element_blank(),
   axis.text = element_blank(),
   axis.ticks =element_blank(),
   axis.title = element_blank(),
   legend.position = c(0.2, 0.3)
)

在这里插入图片描述
准备加入数据

x <- china_map@data #读取行政信息
xs <- data.frame(x,id=seq(0:924)-1) #含岛屿共925个形状
china_map1 <- fortify(china_map) #转化为数据框
china_map_data <- join(china_map1, xs, type = "full")      #合并两个数据框

library(readxl)
a <- read_excel("C:/Users/Administrator/Desktop/a.xlsx")
a

数据
在这里插入图片描述

下面读取业务指标数据，并与地图数据合并：

china_data <- join(china_map_data, mydata, type="full") #合并两个数据框

绘制地图

ggplot(
	china_data, 
	aes(x = long, y = lat, group = group, fill = number)) +
geom_polygon(colour="grey40") +
scale_fill_gradient(low="#99FF99",high="#FF6666") +  #指定渐变填充色，可使用RGB
coord_map("polyconic") #指定投影方式为polyconic，获得常见视角中国地图

在这里插入图片描述
改进

去掉图中的背景色、坐标轴、经纬线等

在这里插入图片描述

ggplot(
	china_data, 
	aes(x = long, y = lat, group = group, fill = number)) +
geom_polygon(colour="grey40") +
scale_fill_gradient(low="#99FF99",high="#FF6666") +  #指定渐变填充色，可使用RGB
coord_map("polyconic")+ #指定投影方式为polyconic，获得常见视角中国地图
#清除不需要的元素
theme(
   panel.grid=element_blank(),
   panel.background=element_blank(),
   axis.text=element_blank(),
   axis.ticks=element_blank(),
   axis.title=element_blank(),
   legend.position = c(0.9,0.6)
)

接下来处理标签问题：

创建各城市中心地理坐标：

在这里插入图片描述
chinaprovincecity

province_city <- read.csv("C:/Users/Administrator/Desktop/chinaprovincecity.csv")  #读取省会城市坐标

ggplot(china_data,aes(long,lat))+
geom_polygon(aes(group=group,fill=number),colour="grey60")+
scale_fill_gradient(low="#99FF99",high="#FF6666") +
coord_map("polyconic") +
geom_text(aes(x = jd,y = wd,label = province), data =province_city,colour="black",size=3)+
theme(
         panel.grid = element_blank(),
        panel.background = element_blank(),
         axis.text = element_blank(),
         axis.ticks = element_blank(),
         axis.title = element_blank()
)

按照同样的方法，绘制重庆地图

地图数据

cq <- rgdal::readOGR("重庆.shp")
plot(cq)

在这里插入图片描述

# 有可能乱码，看一下那些字段乱码了
head(cq@data)
# 逐一修改
Encoding(cq@data$name) <- "UTF-8" ##转换格式为UTF-8

# 再检查一次
head(cq@data)

data1 <- cq@data         
data2 <- data.frame(id=row.names(data1),data1)

cq_ <- fortify(cq) #转化为数据框
cq_map_data <- join(cq_, data2, type = "full") 
head(cq_map_data)

准备填充数据

library(readxl)
library(xlsx)
# 把城市名称保存为excel表，这样就不用了自己写城市名，或者如果业务数据本身有城市名，也可以跳过这步
write.xlsx(cq@data$name,'b.xls',row.names = F)

如果需要这一步，就把数据按城市名对应的添加

例如

在这里插入图片描述

如果不需要这一步，可以用连接的方式，合并数据

这里我使用 group 字段进行合并

#读取业务数据
mydata1 <- read_excel("C:/Users/Administrator/Desktop/b - 副本.xlsx")
mydata1

数据为虚构
在这里插入图片描述

#读取业务指标数据，并与地图数据合并
cq_data <- join(cq_map_data, mydata1, type="full") #合并两个数据框

在这里插入图片描述
开始绘图

# 绘图
windowsFonts(myFont = windowsFont("微软雅黑")) 
ggplot(
	cq_data, 
	aes(x = long, y = lat, group = group, fill = number)) +
geom_polygon(colour="grey40") +ggtitle("2020年重庆xxx数据分布")+
scale_fill_gradient(low="#99FF99",high="#FF6666") +
theme(plot.title=element_text(family="myFont",size=19,hjust = 0.5), 
)

# 主题设置https://blog.csdn.net/u014801157/article/details/24372531

在这里插入图片描述
关于颜色

参考：R语言的颜色应用

颜色数据集RColorBrewer

library(RColorBrewer)
brewer.pal.info

在这里插入图片描述

使用：

brewer.pal(n, name)，n为要生成多少个颜色，name对应上面图片中颜色条的名字

比如：brewer.pal(11,'Spectral')

在这里插入图片描述

自定义渐变色板 colorRampPalette()

用法：
1.输入调色板的主要颜色，返回函数
2.在函数中输入想要获取的颜色数量，返回颜色数值

library(scales) #调用scales库中的show_col函数
mypalette<- colorRampPalette(c("#FF66CC" , "#FF6600"))
mycolors<- mypalette(12)#12种颜色，返回颜色数值
show_col(mycolors)

在这里插入图片描述

mypalette<- colorRampPalette(c("#FF66CC" , "#FF6600",'#66CCFF'))
mycolors<- mypalette(16)#16种颜色，返回颜色数值
show_col(mycolors)

在这里插入图片描述
在这里，可以使用scale_fill_gradientn来设置梯度颜色

在这里插入图片描述

清除一些元素

# 绘图
windowsFonts(myFont = windowsFont("微软雅黑")) 
ggplot(
	cq_data, 
	aes(x = long, y = lat, group = group, fill = number)) +
geom_polygon(colour="grey40") +ggtitle("2020年重庆xxx数据分布")+
scale_fill_gradientn(colours = colorRampPalette(brewer.pal(5,'GnBu'))(24))+
theme(plot.title=element_text(family="myFont",size=19,hjust = 0.5), 
,   panel.grid=element_blank(),
   panel.background=element_blank(),
   axis.text=element_blank(),
   axis.ticks=element_blank(),
   axis.title=element_blank())

在这里插入图片描述

接下来处理标签问题：

创建各城市中心地理坐标：

midpos <- function(data1) mean(range(data1,na.rm=TRUE))
centres <- ddply(cq_data,.(name),colwise(midpos,.(long,lat)))
# 绘图
windowsFonts(myFont = windowsFont("微软雅黑")) 
ggplot(
	cq_data, 
	aes(x = long, y = lat)) +
geom_polygon(aes(group=group,fill=number),colour="white")+
scale_fill_gradientn(colours = colorRampPalette(brewer.pal(5,'GnBu'))(24))+
coord_map("polyconic") +
geom_text(aes(label=name),size =3,family="myFont",fontface="plain",data=centres)+ 
ggtitle("2020年重庆xxx数据分布")+
theme(plot.title=element_text(family="myFont",size=19,hjust = 0.5), 
,   panel.grid=element_blank(),
   panel.background=element_blank(),
   axis.text=element_blank(),
   axis.ticks=element_blank(),
   axis.title=element_blank())

在这里插入图片描述
不过这个加了标签后，有点丑…

添加指北针

用tmap程序包绘制带指北针、图例和经纬网的采样点地图

这里我用就简单的方式添加

library(ggsn)     ## 提供指北针
library(legendMap) ## 提供比例尺

midpos <- function(data1) mean(range(data1,na.rm=TRUE))
centres <- ddply(cq_data,.(name),colwise(midpos,.(long,lat)))
# 绘图
windowsFonts(myFont = windowsFont("微软雅黑")) 
cq =ggplot(
	cq_data, 
	aes(x = long, y = lat)) +
geom_polygon(aes(group=group,fill=number),colour="white")+
scale_fill_gradientn(colours = colorRampPalette(brewer.pal(5,'GnBu'))(24))+
coord_map("polyconic") +
geom_text(aes(label=name),size =3,family="myFont",fontface="plain",data=centres)+ 
ggtitle("2020年重庆xxx数据分布")+
theme(plot.title=element_text(family="myFont",size=19,hjust = 0.5), 
,   panel.grid=element_blank(),
   panel.background=element_blank(),
   axis.text=element_blank(),
   axis.ticks=element_blank(),
   axis.title=element_blank())

library(ggsn)     ## 提供指北针
north2(cq,.82,.88)## 数字是指定指北针的位置

在这里插入图片描述

不过我的比例尺加不上去

也可以通过重庆主城区，填充颜色

主城区

cq_map <- rgdal::readOGR("重庆主城区.shp")
plot(cq_map)

在这里插入图片描述
主城新区

cq_map_12 <- rgdal::readOGR("主城新区.shp")
plot(cq_map_12)

cq_map_121 <- fortify(cq_map_12) #转化为数据框
cq_map_12@data

library(readxl)
b <- read_excel("C:/Users/Administrator/Desktop/a - 副本.xlsx")
b
c <- join(cq_map_121, b, type="full") #合并两个数据框

ggplot(
	c, 
	aes(x = long, y = lat, group = group, fill = number)) +
geom_polygon(colour="grey40") +
scale_fill_gradient(low="#99FF99",high="#FF6666") +  #指定渐变填充色，可使用RGB
theme(
   panel.grid=element_blank(),
   panel.background=element_blank(),
   axis.text=element_blank(),
   axis.ticks=element_blank(),
   axis.title=element_blank(),
   #legend.position = c(0.9,0.6)
)

在这里插入图片描述

参考：

http://blog.sina.com.cn/s/blog_6bc5205e0102vma9.html
https://zhuanlan.zhihu.com/p/29157375
https://zhuanlan.zhihu.com/p/26880814
https://zhuanlan.zhihu.com/p/187531793

卖山楂啦prss

关注

10
点赞
踩
125

收藏

觉得还不错? 一键收藏
6
评论
绘图杂记【5】R语言绘制热力地图（重庆热力地图）

安装和加载包install.packages(mapdata)install.packages(maptools)install.packages(ggplot2)install.packages(plyr)install.packages(rgdal)install.packages(mapproj)library(ggplot2)library(maps)library(mapdata)library(maptools)library(plyr)library(rgdal)
复制链接

扫一扫