R语言可视化实现地图填充与散点图图层叠加
本篇内容主要讲解“R语言可视化实现地图填充与散点图图层叠加”，感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷，实用性强。下面就让小编来带大家学习“R语言可视化实现地图填充与散点图图层叠加”吧!
今天跟大家分享关于如何在地图图层上添加散点图。
散点图需要精确的经纬度信息才能在叠加的图层上进行映射，因此我们选用中国省级轮廓地图以及各省省会城市的经纬度进行案例演示。
加载包：
library(ggplot2)
library(plyr)
library(maptools)
library(sp)
导入中国省界地图：
china_map<-readShapePoly("c:/rstudy/bou2_4p.shp")
data1<- china_map@data
data2<- data.frame(id=row.names(data1),data1)
数据格式转化及业务数据合并：
china_map1 <- fortify(china_map)
china_map_data <- join(china_map1,data2, type = "full")
mydata <- read.csv("c:/rstudy/geshengzhibiao.csv")
china_data <- join(china_map_data, mydata, type="full")
各省省会城市经纬度数据：
province_city <- read.csv("c:/rstudy/chinaprovincecity.csv")
省级轮廓地图上添加散点图图层：
ggplot(china_data,aes(long,lat))+
geom_polygon(aes(group=group),fill="white",colour="grey60")+
geom_point(data=province_city,aes(x=jd,y=wd),colour="red")+
coord_map("polyconic") +
theme(
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank()
)
接下来，我们可以给各省的省会城市赋值，将散点图的大小映射给连续性数值变量，使其变成气泡图。
province_city$zhibiao<-NULL
province_city$zhibiao<-rnorm(34,100,50)
windowsFonts(myFont = windowsFont("微软雅黑"))
ggplot()+
geom_polygon(data=china_data,aes(x=long,y=lat,group=group),fill="grey95",colour="grey80")+
geom_point(data=province_city,aes(x=jd,y=wd,size=zhibiao),shape=21,fill="#8E0F2E",colour="black",alpha=0.4)+
scale_size_area(max_size=8)+
coord_map("polyconic") +
guides(size=guide_legend(reverse=TRUE,title=NULL))+
ggtitle("某公司2015~2016年度营业状况分布图")+
theme(
title=element_text(family="myFont"),
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank(),
legend.position =c(0.15,0.4),
legend.background=element_rect(colour="white",fill="white"),
legend.text.align=1
)
最后我们来处理标签问题：
因为这里使用了散点图（气泡图）作为数据展示方式，而作为底图的地图图层仅仅是作为定位信息，没有包含任何的数据信息，而且页面整体上没有太多地方放置省名称标签，所以我们有选择性的显示前五个数据所代表的省份标签，以防标签太多导致页面杂乱。
labelper<-province_city[order(province_city[,5],decreasing=T),][1:10,]
ggplot()+
geom_polygon(data=china_data,aes(x=long,y=lat,group=group),fill="grey95",colour="grey80")+
geom_point(data=province_city,aes(x=jd,y=wd,size=zhibiao),shape=21,fill="#8E0F2E",colour="black",alpha=0.4)+
scale_size_area(max_size=8)+
coord_map("polyconic") +
geom_text(aes(x=jd+2.3,y=wd,label=city),size =3,family="myFont",fontface="plain",data=labelper) +
guides(size=guide_legend(reverse=TRUE,title=NULL))+
ggtitle("某公司2015~2016年度营业状况分布图")+
theme(
title=element_text(family="myFont"),
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank(),
legend.position =c(0.15,0.4),
legend.background=element_rect(colour="white",fill="white"),
legend.text.align=1
)
如果甲你想要添加全部的标签，直接使用province_city中的city标签即可：
ggplot()+
geom_polygon(data=china_data,aes(x=long,y=lat,group=group),fill="grey95",colour="grey80")+
geom_point(data=province_city,aes(x=jd,y=wd,size=zhibiao),shape=21,fill="#8E0F2E",colour="black",alpha=0.4)+
scale_size_area(max_size=8)+
coord_map("polyconic") +
geom_text(aes(x=jd+2.3,y=wd,label=city),size =3,family="myFont",fontface="plain",data=province_city) +
guides(size=guide_legend(reverse=TRUE,title=NULL))+
ggtitle("某公司2015~2016年度营业状况分布图")+
theme(
title=element_text(family="myFont"),
panel.grid = element_blank(),
panel.background = element_blank(),
axis.text = element_blank(),
axis.ticks = element_blank(),
axis.title = element_blank(),
legend.position =c(0.15,0.4),
legend.background=element_rect(colour="white",fill="white"),
legend.text.align=1
)
到此，相信大家对“R语言可视化实现地图填充与散点图图层叠加”有了更深的了解，不妨来实际操作一番吧！这里是亿速云网站，更多相关内容可以进入相关频道进行查询，关注我们，继续学习！