请选择 进入手机版 | 继续访问电脑版
楼主: 1050658607
1367 1

[作图问题求助] #中国地图+ggplot2可视化 [推广有奖]

  • 0关注
  • 0粉丝

高中生

40%

还不是VIP/贵宾

-

威望
0
论坛币
38 个
通用积分
27.7227
学术水平
1 点
热心指数
1 点
信用等级
0 点
经验
1715 点
帖子
28
精华
0
在线时间
24 小时
注册时间
2016-5-4
最后登录
2017-7-19

1050658607 在职认证  发表于 2016-12-19 09:02:17 |显示全部楼层 |坛友微信交流群
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
#中国化热图
library(maptools)
library(ggplot2)
#data <- read.csv("C:\\map\\data.csv",header=TRUE,sep=",")
china_map=readShapePoly('C:\\map\\bou2_4p.shp')
plot(china_map) map.rar (1.44 MB)
theme_set(theme_bw())
ggplot(china_map,aes(x=long,y=lat,group=group)) +geom_polygon(fill="white",colour="grey")#+coord_map("orthographic")
?coord_map
x <- china_map@data           #读取行政信息
xs <- data.frame(x,id=seq(0:924)-1)          #含岛屿共925个形状
china_map1 <- fortify(china_map)           #转化为数据框
library(plyr)
china_map_data <- join(china_map1, xs, type = "full")
unique(china_map@data$NAME)
#######################
shuju<-read.csv('C:\\map\\shuju.csv',header = T)#相应的数据读入,改变数据的值即可
#######################
china_data<-join(china_map_data,shuju,type="full")
ggplot(china_data, aes(x = long, y = lat, group = group, fill = zhibiao))+
  geom_polygon(colour="grey40")+scale_fill_gradient(low="white",high="steelblue")+   #指定渐变填充色,可使用RGBcoord_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.2,0.3)
  )
midpos <- function(x) mean(range(x,na.rm=TRUE))#取形状内的平均坐标
centres <- ddply(china_data,.(NAME),colwise(midpos,.(long,lat)))
ggplot(china_data,aes(long,lat))+#此处语法与前面不同,参考ggplot2一书P85
  geom_polygon(aes(group=group,fill=zhibiao),colour="black")+
  scale_fill_gradient(low="white",high="steelblue") +
  geom_text(aes(label=NAME),data=centres)+
  theme(
    panel.grid = element_blank(),
    panel.background = element_blank(),
    axis.text = element_blank(),
    axis.ticks = element_blank(),
    axis.title = element_blank()
  )
province_city <-read.csv("C:\\map\\data.csv",header=TRUE,sep=",")
ggplot(china_data,aes(long,lat))+
  geom_polygon(aes(group=group,fill=zhibiao),colour="grey60")+ylim(c(18,55))+
  scale_fill_gradient(low="white",high="red",name="店铺数量")+
  geom_text(aes(x = Jd,y = Wd,label = 城市), data =province_city)+
  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.1,0.2)
  )+geom_point(aes(Jd,Wd),data = province_city,size=1)





二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:ggplot2 gplot 中国地图 plot 可视化 中国地图 library

已有 1 人评分论坛币 学术水平 热心指数 收起 理由
飞天玄舞6 + 40 + 1 + 1 精彩帖子

总评分: 论坛币 + 40  学术水平 + 1  热心指数 + 1   查看全部评分

1050658607 在职认证  发表于 2016-12-19 09:13:36 |显示全部楼层 |坛友微信交流群

结果

K3A5_`1I8UMYPDR~RYICD9V.png









使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-28 23:56