用R语言绘制地理热图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

参考来源：http://blog.sina.com.cn/s/blog_6bc5205e0102vma9.html
数据来源：中国地图数据文件——可到微信公众号 iamExcelPro 发送 shapefile 获得下载地址
          各省、直辖市经纬度坐标——可到微信公众号 iamExcelPro 发送 shapefile 获得下载地址
          业务数据为2014年各省GDP——国家统计局
R语言绘制地理热图.rar (1.44 MB, 需要: 2 个论坛币) 本附件包括：

• bou2_4p.dbf
• bou2_4p.shp
• bou2_4p.shx
• china-province-border-data.tar.gz
• chinaprovinceborderdata_tar_gz.zip
• chinaprovincecity.csv
• geshengzhibiao.csv
• 用R绘制地图.R

2016-8-22 08:44:17 上传

程序附件
需要: 2 个论坛币  [购买]

1. ###############设置工作目录#####################
   
2. setwd("D:/RScript")
   
3. ##################################################
   
4. 
   
5. ##############加载程序包，读取地图数据##############
   
6. library("maptools")
   
7. library(ggplot2)
   
8. library(plyr)
   
9. china_map = readShapePoly("bou2_4p.shp")       # 读取地图空间数据
   
10. 
    
11. 
    
12. ###################数据预处理#######################
    
13. x <- china_map@data          #读取行政信息
    
14. xs <- data.frame(x,id=seq(0:924)-1)          #含岛屿共925个形状
    
15. china_map1 <- fortify(china_map)           #转化为数据框
    
16. china_map_data <- join(china_map1, xs, type = "full")       #合并两个数据框
    
17. 
    
18. ##读取指标数据，csv格式
    
19. mydata <- read.csv("geshengzhibiao.csv")         
    
20. china_data <- join(china_map_data, mydata, type="full")          #合并两个数据框
    
21. 
    
22. 
    
23. ####################绘制地图，并清除背景色、坐标轴、经纬线,图例放左下角##########
    
24. province_city <- read.csv("chinaprovincecity.csv")  #读取省会城市坐标
    
25. china_data <- join(china_data, province_city, type="full")          #合并两个数据框
    
26. 
    
27. ##图1：用数据填充，GDP量越大颜色越深
    
28. pp <- ggplot(china_data,aes(long,lat))+
    
29.   labs(title="2014年全国各省GDP（无香港和台湾数据）单位：亿元")+
    
30.   geom_polygon(aes(group=group,fill=X2014年),colour="grey60")+
    
31.   scale_fill_gradient(name="2014年GDP",low="white",high="steelblue") +
    
32.   coord_map("polyconic") +
    
33.   #coord_map(projection = "azequidistant")+
    
34.   geom_text(aes(x = jd,y = wd,label = province), data =province_city)+
    
35.   theme(
    
36.     #panel.grid = element_blank(),
    
37.     #panel.background = element_blank(),
    
38.     #axis.text = element_blank(),
    
39.     #axis.ticks = element_blank(),
    
40.     axis.title = element_blank(),
    
41.     legend.position = c(0.2,0.3)
    
42.   )
    
43. pp
    
44. 
    
45. ##图2：将数据映射为散点绘制在地图上
    
46. pp <- ggplot(china_data,aes(long,lat,group=group))+geom_path()
    
47. pp <- pp+labs(title="2014年全国各省GDP（无香港和台湾数据）单位：亿元")
    
48. pp <- pp+coord_map("polyconic")
    
49. pp <- pp+geom_point(data=china_data,aes(x=jd,y=wd,size=X2014年),color="steelblue")  #aes() mapping the data to the graph
    
50. pp <- pp+scale_size_continuous(name="2014年GDP")
    
51. #pp <- pp+scale_colour_continuous(low="white",high="black")
    
52. pp <- pp+  theme(
    
53.   axis.title = element_blank(),
    
54.   legend.position = c(0.2,0.3)
    
55. )
    
56. pp

复制代码

扫码加我 拉你入群

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

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

分享0 收藏14 回帖

关键词：R语言 Shapefile Library ggplot2 script 中国地图 下载地址 library 统计局 直辖市

支持分享。

不错，收藏了

求问X2014年表示什么

gpnrnvyy 发表于 2016-9-7 09:50
求问X2014年表示什么

是2014年的GDP，原数据是存入CSV文件中的，读入数据后R自动在前面添加了一个X

Jack00878 发表于 2016-8-22 08:38
参考来源：http://blog.sina.com.cn/s/blog_6bc5205e0102vma9.html
数据来源：中国地图数据文件——可到微 ...

谢谢你

这个看起来不错哦很实用

请问为什么地图还是灰色的没有显示渐变颜色呀

谢谢分享