每R一点:各种画地图,全是知识点,90%人不知道!(转)
2015-01-22 09:42
274 查看
R语言绘制地图,在数据分析中经常能够用到,并且会达到非常好的展示效果,本节以例子形式,介绍如何使用R语言工具,画出理想的地图。
本节例子在 R version 2.15.3版本下运行顺畅,其他版本待定。
代码如下:第一个小例子
下面使用ggplot来画。
画一个更好看的地图 ,如何!
再画一个城市区域型的 地图
关注点:
⊙ 上面只是介绍几种画图的办法。除此之外,还有更多地图函数,各位朋友可以好好找找看。对于地图中具体的参数和相应指标,可以用R自带帮助文档查看。
转自:http://mp.weixin.qq.com/s?__biz=MzA3NDUxMjYzMA==&mid=209747952&idx=2&sn=72725a5801c0b1c14a7da6cc7e68719c#rd
本节例子在 R version 2.15.3版本下运行顺畅,其他版本待定。
代码如下:第一个小例子
#加载相应的包,并读入数据,最后画图。 library(maptools); library(ggplot2); china_map<-readShapePoly("d://map//bou2_4p.shp", proj4string=CRS("+proj=longlat +ellps=clrk66")); plot(china_map);
下面使用ggplot来画。
ggplot(china_map, aes(x = long, y = lat, group = group)) + geom_polygon(fill = "beige") + geom_path(colour = "grey40");
画一个更好看的地图 ,如何!
x<-china_map@data; # S4类数据 用@ 取子集 xs<-data.frame(x,id=seq(0:924)-1); # library(ggplot2) china_map1<-fortify(china_map); # 转化为数据框 # 没有行政区域的信息啊,这就需要前面的那个数据框xs了。 library(plyr); china_mapdata<-join(china_map1, xs, type = "full"); # 合并两个数据框 # 下面用这个新的数据框来画中国地图 ggplot(china_mapdata, aes(x = long, y = lat, group = group,fill=NAME))+geom_polygon( )+geom_path(colour = "grey40")+ scale_fill_manual(values=colours(),guide=FALSE);
再画一个城市区域型的 地图
zhejiang<-subset(china_mapdata,NAME=="浙江省"); ggplot(zhejiang, aes(x = long, y = lat, group = group,fill=NAME))+ geom_polygon(fill="beige" )+ geom_path(colour = "grey40")+ ggtitle("中华人民共和国浙江省")+ geom_point(x=120.12,y=30.16,fill=FALSE)+ annotate("text",x=118.3,y=30,label="杭州市");
关注点:
⊙ 上面只是介绍几种画图的办法。除此之外,还有更多地图函数,各位朋友可以好好找找看。对于地图中具体的参数和相应指标,可以用R自带帮助文档查看。
转自:http://mp.weixin.qq.com/s?__biz=MzA3NDUxMjYzMA==&mid=209747952&idx=2&sn=72725a5801c0b1c14a7da6cc7e68719c#rd
相关文章推荐
- Vim的各种小知识点
- Qt开发中的实用笔记三--关于各种类的零碎知识点:
- 90%的用户都不知道手机内部功能
- iOS--点击地图上某一点获取该点的经纬度
- ios 备忘录---你可能不知道的小知识点
- JavaScript 优雅的实现方式包含你可能不知道的知识点
- Android Context完全解析,你所不知道的Context的各种细节
- 给还在不断纠结各种linux的同学们分享一点体会。
- 一点思考,战棋类游戏地图核心应该有什么?
- 最好用的止咳药,就在厨房里!90%的人不知道!
- 什么是用户?估计90%人不知道
- Android Context完全解析,你所不知道的Context的各种细节
- 原来IC是这样设计的,竟然90%的人都不知道!
- 我想翻译一点关于编写控件的文章,不知道有没有人喜欢?
- Service的各种知识点
- 开通了博客,以后写一点或者转一点自己觉得有用也对大家有用的各种东西
- Android Context完全解析,你所不知道的Context的各种细节
- [导入]项目开发过程中遇到的一点小知识点
- ArcGIS地图文档MXD效率慢的一点建议
- javascript你不知道的知识点