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

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

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

期待高手指点

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

1. 4049          NA          NA
   
2. 4050   112.44239    5.316195
   
3. 4051   112.43666    5.346023
   
4. 4052   112.42864    5.342843
   
5. 4053   112.44239    5.316195
   
6. 4054          NA          NA
   
7. 4055   112.57360    5.177665
   
8. 4056   112.56042    5.218391
   
9. 4057   112.55240    5.214976
   
10. 4058   112.57360    5.177665
    
11. 4059          NA          NA
    
12. 4060   112.29858    4.250258
    
13. 4061   112.29228    4.261734
    
14. 4062   112.28712    4.255352
    
15. 4063   112.29858    4.250258

复制代码

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

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


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


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

重庆的边界也初步解决~

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

1. ######################################
   
2. require(maptools)
   
3. require(mapdata)
   
4. require(maps)
   
5. x= readShapePoly("bou2_4p.shp")# 从cos 中下载的链接
   
6. #shp的数据结构是够复杂的，下面都是在搞数据
   
7. summary(x)
   
8. class(x)
   
9. SpatialPolygonsDataFrame#从summary 中看到这个对象，查看下
   
10. #第一次遇到 @ 符号。。。此代码中有@。。。？slot
    
11. slotNames(x)
    
12. attributes(x@data)
    
13. x@data
    
14. plot(x)
    
15. 
    
16. 
    
17. #####要重上述数据中怎么提取 重庆的坐标呢？########
    
18. x@data[x@data$NAME=="重庆市",]#在第207位置，看看面积
    
19. 
    
20. 
    
21. slotNames(x@polygons[[925]])#还是为了查看数据结构
    
22. 
    
23. slot(x@polygons[[208]],"area")#尝试207,208后比较面积，得出208估计为重庆
    
24. 
    
25. chongqing=slot(x@polygons[[208]],"Polygons")[[1]]@coords##尝试再尝试怎么提取 坐标
    
26. plot(chongqing)#看看是不是重庆地图
    
27. 
    
28. ###确定是重庆坐标后。提取经纬度，生成文件。
    
29. write.table(chongqing[,1],"chongqing_jing.txt",sep=",",row.names=F,col.names=F,eol=",")
    
30. write.table(chongqing[,2],"chongqing_wei.txt",sep=",",row.names=F,col.names=F,eol=",")
    
31. 
    
32. ###把重庆的，经纬度 修改到 mapdata 中的 china 中。
    
33. ###map("china",fill=F,plot=T,interior=T)
    
34. ###ch=map("china")
    
35. ###fix(ch)
    
36. ###chinaCQ.r的 数据为structure(list(x=c(重庆东经坐标，NA,钓鱼岛东经坐标,NA,....),y=c(重庆北纬坐标，NA,钓鱼岛北纬坐标,NA,....))，。。。）
    
37. 
    
38. chi <-source("chinaCQ.r")
    
39. str(chi$value)
    
40. map(chi$value,fill=F,plot=T,interior=T,myborder = 0.001)
    
41. 
    
42. ###可以看到重庆和四川分开了， 其中重庆的地图数据精度要高的多（比R自带的），如果觉得线太重， 可以
    
43. ###稀释一些点，让它看上去和周围平缓些。
    
44. 
    
45. ###其实可以直接那个shp文件， 画图， 但为了数据结构的简单性， 觉得mapdata的数据形式更简单，直观些
    
46. ###而进行的修补。

复制代码

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

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"))

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

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

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

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

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

查看  mapdata china 数据

1. List of 4
   
2. $ x    : num [1:4508] 123 124 124 124 NA ...
   
3. $ y    : num [1:4508] 46.3 46.3 46 45.9 NA ...
   
4. $ range: num [1:4] 73.44 135.1 4.25 53.57
   
5. $ names: chr [1:712] "1" "2" "3" "4" ...
   
6. - attr(*, "class")= chr "map"

复制代码

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