【版块活动】一月一包主题讨论

6851

收藏 2014-04-14

本月主题：地图主题包

提议人：jmpamao

讨论框架：包括但不限于以下几个方面(欢迎补充)
(1)R中地图的绘制方法(空间数据可视化)，包括利用R本身的作图功能实现或是利用接口实现等
(2)空间数据的建模与分析，包括数据的读取、操作、处理等，有实例则更佳

参与方法(以跟帖形式)：
(1)分享资料
(2)分享使用和学习心得
(3)直接参与讨论(也欢迎提问)

奖励方法
(1)将有价值的回帖设为精彩回复或直接编辑到主题帖中，供坛友学习
(2)只要是具有学习价值或者讨论价值的回帖，基于论坛币的奖励
(3)请不要进行无意义的回帖，诸如顶、学习一下等，表达支持可以通过给回帖评分的形式
(4)对参与活动特别热心且有突出贡献的坛友可推荐为本版热心会员或后备版主

温馨提示
(1)分享资料包括好的书籍、博文以及之前论坛的精彩回帖等，考虑到R的更新速度，请在分享前确认资料的有效性。
(2)鼓励免费分享，不建议收取论坛币，对于好的资料会有丰厚的论坛币补偿。
(3)R毕竟只是一个工具，在分享使用经验欢迎附上使用过程中的背景知识。
(4)建议对于回帖正文，除了部分难以翻译的英文关键词，请不要使用英文。

感谢jmpamao对活动的倡议和努力，欢迎大家提出更多好的意见建议。

扫码加我拉你入群

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

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

全部回复

求证1加1

2014-4-14 15:59:32

先抛砖引玉一下

推荐几篇不错的博文
(1)关于R for map的一个系列：
http://site.douban.com/182577/widget/notes/10568279/note/257898418/
有代码有图有说明，值得一读，不过部分代码已经过时，但在看了帮助文档之后很容易知道怎么改过来，并且内容基本上已经可以满足一般的可视化需求了

(2)COS上的几篇老文(不止这些)
http://cos.name/2009/07/drawing-china-map-using-r/
http://cos.name/2013/01/drawing-map-in-r-era/
http://cos.name/2013/06/weibo-fans-map-visualization/

(3)这篇是介绍ggmap的，介绍的很仔细且易懂
http://www.maenchi.com/?p=199

另外这是我前几天看到的一篇文章，这个蛮不错的，反正也有例子，有兴趣的朋友可以去玩玩
http://blog.revolutionanalytics.com/2014/04/animated-choropleths-in-r.html

大家接着分享和补充吧~奖励多多哈~

扫码加我拉你入群

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

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

yywan0913

2014-4-14 16:04:24

我能先顶一下么？

扫码加我拉你入群

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

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

yywan0913

2014-4-14 16:17:51

以上整理就这些吧？欢迎添加
maps #很基础
mapdata  maptools #利用已知数据
ggplot2
ggmap  mapproj    #借助谷歌，有限制。不过挺好用
rMaps  #借助软件

扫码加我拉你入群

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

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

hlw

2014-4-14 17:12:09

不错
很好

扫码加我拉你入群

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

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

dxystata

2014-4-14 21:32:46

画中国地图应包括钓鱼岛、台湾、南海等！

扫码加我拉你入群

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

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

点击查看更多内容…

dxystata

2014-4-14 21:36:11

建议最后由jmpamao大神整理并通过实例展示！

扫码加我拉你入群

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

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

jmpamao

2014-4-14 23:09:27

maps 包中 num.polygons 代码怎么看？

扫码加我拉你入群

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

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

ywh19860616

2014-4-15 08:26:21

求证1加1 发表于 2014-4-14 15:59
先抛砖引玉一下

推荐几篇不错的博文

哇，非常不错
版主，有没有看到在地图中画柱形图的？就是现在有多个国家或者地区，
然后根据某个指标，把相应的柱形图标注上去，来直观比较大小。

扫码加我拉你入群

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

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

jmpamao

2014-4-15 21:09:02

dxystata 发表于 2014-4-14 21:32
画中国地图应包括钓鱼岛、台湾、南海等！

复制代码

中国地图.rar
大小:(53.8 KB)

马上下载

本附件包括：

中国地理位置_经纬度.doc
china.r

扫码加我拉你入群

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

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

jmpamao

2014-4-15 21:22:18

我现在就想怎么把钓鱼岛的坐标放进那个 structure（list（x=...y=......), 大家帮忙看看，估计要看代码
map了，并不长，但理解不易。

扫码加我拉你入群

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

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

xingzhaoh

2014-4-16 07:54:34

jmpamao 发表于 2014-4-15 21:22
我现在就想怎么把钓鱼岛的坐标放进那个 structure（list（x=...y=......), 大家帮忙看看，估计要看代 ...

期待高手指点

扫码加我拉你入群

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

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

jmpamao

2014-4-16 09:36:33

知道怎么添加钓鱼岛坐标了，今天早上灵光一想，为什么坐标中， structure 中有那么多的 NA,这没事做吗？
所以把 structure 变为 data.frame 看。

复制代码

觉得这些数据有可能是通过 NA 来组成一个闭环，这也是为什么之前想当然的以为，只要把钓鱼岛的坐标放进 x, y 就ok，变为线段，而没有断截。所以添加 NA 后，发现 wow， works~
#红圈内的小黑点为钓鱼岛坐标，添加NA 组成闭环，看来是MAP 的处理了，还得看代码

附件列表

钓鱼岛.jpg

原图尺寸 19.38 KB

不加NA为线段.jpg

原图尺寸 26.13 KB

扫码加我拉你入群

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

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

jmpamao

2014-4-16 10:38:28

下一步，我们看看能否把四川和重庆边界的坐标找出来，捣鼓一下，看看能否把 mapdata 中 china 的数据进行完善~~~ structure(......)

cos 中也有大牛，提到四川和重庆边界的问题，是提出来了，但没有去把数据进行完善。
而我们可以趁这个讨论地图主题的机会，集合大家的力量，看看能否完善其数据。
再进行分享，觉得是有意义~

##晚上要消化下 1+1版主给的资料。。。。

扫码加我拉你入群

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

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

jmpamao

2014-4-16 19:21:37

重庆的边界也初步解决~

1+1版主提供的http://cos.name/2009/07/drawing-china-map-using-r/ 中》有从这里下载中国地图的GIS数据，这是一个压缩包，完全解压后包含三个文件（bou2_4p.dbf、bou2_4p.shp和bou2_4p.shx）

复制代码

把 chinaCQ.txt 也上传下，改chinaCQ.r
大家都来参与讨论啊~~~

附件列表

重庆边界.jpg

原图尺寸 38.68 KB

chinaCQ.txt

大小:154.7 KB

马上下载

扫码加我拉你入群

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

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

xingzhaoh

2014-4-17 08:39:35

jmpamao 发表于 2014-4-16 19:21
重庆的边界也初步解决~

1+1版主提供的http://cos.name/2009/07/drawing-china-map-using-r/ 中》有从这 ...

library(maptools)

axx<-readShapePoly("bou2_4p.shp", proj4string=CRS("+proj=longlat +ellps=clrk66"))

x=axx

getColor=function(mapdata,provname,provcol,othercol)

{

f=function(x,y) ifelse(x %in% y,which(y==x),0);

colIndex=sapply(mapdata$NAME,f,provname);

col=c(othercol,provcol)[colIndex+1];

return(col);

}

provname=c("湖北省","甘肃省","河北省","辽宁省","黑龙江省","内蒙古自治区");

provcol=c("green","green","green","green","green","green","green");

plot(x,col=getColor(x,provname,provcol,"white"))

怎么在各个省加上省的名字上去

扫码加我拉你入群

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

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

jmpamao

2014-4-17 09:28:58

xingzhaoh 发表于 2014-4-17 08:39
library(maptools)axx

我之前画的图，标注省会城市的
是从百度下的坐标。
你去找下有没有省的标注坐标（就是在每个区域内，摆放的那个点比较居中就ok）

在 shp 的数据也可以调出一个省的边界坐标，如果你愿意的话，自己查看range，然后，大体可以定一个位置（上面有怎么找重庆市，就其边界坐标，但还是麻烦些）
。。。
你有了这些数据后，就可以 text 了，请查看 mapproj 包，
####
其他，哪个包，可以直接通过参数来调用，我也想知道

扫码加我拉你入群

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

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

jmpamao

2014-4-17 09:45:41

jmpamao 发表于 2014-4-16 19:21
重庆的边界也初步解决~

1+1版主提供的http://cos.name/2009/07/drawing-china-map-using-r/ 中》有从这 ...

大家帮忙看看着个重庆边界的问题

1、四川可能要重新画，现在的方式还是包含重庆。
2、想请教，边界的点，比如四川和重庆的，这些点到底是属于四川还是重庆，还是两者公用，来组成一个块。请指教

扫码加我拉你入群

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

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

wfldragon

2014-4-17 11:04:58

其实核心的问题还是分组。只是边界还比较容易，如果涉及面积就还需要把重庆从四川里面减去，否则面积不对

扫码加我拉你入群

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

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

jmpamao

2014-4-17 11:57:40

wfldragon 发表于 2014-4-17 11:04
其实核心的问题还是分组。只是边界还比较容易，如果涉及面积就还需要把重庆从四川里面减去，否则面积不对

查看 mapdata china 数据

复制代码

是利用了 711个NA 来分 712个块，现在的names 是用 “1”，“2”，“3,”来表示。
也跟 1+1版主提了，看看大家有没有兴趣，添加一项names_cn : "重庆“，”北京“。。。。。这类的。以后数据调用的时候，估计有帮助

扫码加我拉你入群

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

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

xingzhaoh

2014-4-17 15:48:26

jmpamao 发表于 2014-4-17 09:28
我之前画的图，标注省会城市的
是从百度下的坐标。
你去找下有没有省的标注坐标（就是在每个区域内， ...

能说的详细一点嘛，着急用呢

扫码加我拉你入群

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

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

Livreanalystes

2014-4-17 15:54:17

dingyiding

扫码加我拉你入群

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

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

jmpamao

2014-4-17 16:29:33

xingzhaoh 发表于 2014-4-17 08:39
library(maptools)axx

其数据是目测的，

复制代码

扫码加我拉你入群

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

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

jmpamao

2014-4-19 20:44:28

今天再次搞四川和重庆的问题。
1、搞清楚， plot(c(1,1,1,1,1,1,1,1,2),c(1,1,1,1,1,1,1,1,2))和 plot(c(1,2),c(1,2))的制图效果一样，没有因为点的重复，而加深，这和text 不同。所以在重庆和四川的边界点，制图上并没有因此加深。（还纠结过很久，⊙﹏⊙b汗）

复制代码

2、计算面积
尝试通过经纬度计算面积。需要大家的帮助。从网上查阅，有些如何计算坐标 (X1，Y1)，(X2，Y2),...(X1，Y1) 组成的闭环的多边形面的例子，下面的这个网站，需要科学上网，。
http://local.wasp.uwa.edu.au/~pbourke/geometry/polyarea/

复制代码

[1] 7.715986

复制代码

尝试了四川，上海后都得到此公式的结果，和国家地图shp 文件中保存的“AREA"项，相同。
因此，请教大家，这个公式的单位是什么经纬度的平方？是平面面积非球面？这个单位怎么转换平方公里等。貌似维度的1度长度是变化的。上述公式如何解释？

bou2_4p.rar
大小:(403.63 KB)

马上下载

本附件包括：

bou2_4p.shp

扫码加我拉你入群

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

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

小洋仔

2014-4-20 18:52:11

jmpamao 发表于 2014-4-15 21:09

赞。。

扫码加我拉你入群

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

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

hyu9910

2014-4-24 00:16:03

牛叉哈

扫码加我拉你入群

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

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

统计R浪人

2014-4-24 07:53:27

来晚了，好热闹啊，以后要常来学习

扫码加我拉你入群

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

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

yywan0913

2014-4-24 08:34:56

统计R浪人发表于 2014-4-24 07:53
来晚了，好热闹啊，以后要常来学习

都等着你压轴表演呢。

扫码加我拉你入群

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

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

dxystata

2014-4-24 09:05:33

ZF网站应该有相应数据可以下载吧！

扫码加我拉你入群

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

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

统计R浪人

2014-4-24 22:26:45

jmpamao 发表于 2014-4-19 20:44
今天再次搞四川和重庆的问题。
1、搞清楚， plot(c(1,1,1,1,1,1,1,1,2),c(1,1,1,1,1,1,1,1,2))和 plot(c(1 ...

1、如果应用精度要求不高的话可以用一个比例直接转换经纬度面积与平方公里面积
2、如果精度需要很高的话可以参考该方法，看看这个如何经纬度面积转平方米
http://blog.csdn.net/lockepeak/article/details/60792063、不规则图形计算如图示：

附件列表

Snap1.jpg

原图尺寸 32.79 KB

示意图

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群