R-中国地图

admin

一、调用相关包

1. install.packages(c("maps", "mapdata","sp“,"maptools"))
   
2. library(maps)
   
3. library(mapdata)
   
4. library(sp)
   
5. library(maptools)
   
6. map("china")
   

复制代码

1. x=readShapePoly('D:/工作文件夹/R语言与中国地图/chinaprovinceborderdata_tar_gz/bou2_4p.shp')
   
2. plot(a)
   
3. plot(a,col=gray(924:0/924))
   
4. 
   
5. getColor=function(mapdata,provname,provcol,othercol)
   
6. {
   
7.   f=function(x,y) ifelse(x %in% y,which(y==x),0);
   
8.   colIndex=sapply(mapdata@data$NAME,f,provname);
   
9.   col=c(othercol,provcol)[colIndex+1];
   
10.   return(col);

复制代码

1. #读取数据
   
2. b<-read.table("D:/工作文件夹/R语言与中国地图/population.txt",head=T)
   
3. provname<-as.matrix(b[,2])
   
4. pop<-as.matrix(b[,3])
   
5. provcol=rgb(red=1-pop/max(pop)/10,green=1-pop/max(pop)/2,blue=0.2)
   
6. plot(x,col=getColor(x,provname,provcol,"white"),xlab="",ylab="")

复制代码

二、利用google数据

1. library(googleVis) #类似的包还有ggmap 等
   
2. 
   
3. b<-read.table("D:/工作文件夹/R语言与中国地图/population.txt",head=T)
   
4. 
   
5. G3 <- gvisGeoChart(b, locationvar='ISO代码',
   
6. colorvar='人口数量',options=list(region='CN',
   
7. displayMode="regions",resolution="provinces",
   
8. colorAxis="{colors: ['yellow','red']}" ))
   
9. 
   
10. plot(G3)
    

复制代码

http://cos.name/2009/07/drawing-china-map-using-r/

http://www.dataguru.cn/article-3062-1.html参考网站

http://platform.eeeeee.org/info/country/省份数据

三、中国航线地图

1. library(maptools)
   
2. x=readShapePoly('D:/工作文件夹/R语言与中国地图/chinaprovinceborderdata_tar_gz/bou2_4p.shp')
   
3. plot(x)
   
4. getColor=function(mapdata,provname,provcol,othercol)
   
5. {
   
6.   f=function(x,y) ifelse(x %in% y,which(y==x),0);
   
7.   colIndex=sapply(mapdata@data$NAME,f,provname);
   
8.   col=c(othercol,provcol)[colIndex+1];
   
9.   return(col);}
   
10. #自定义函数
    
11. 
    
12. Par(mar(0,0,1,0),bg=”gray”)
    
13. plot(x,col="black",border="black")
    
14. 
    
15. 
    
16. #航线数据处理
    
17. library(stringr)
    
18. 
    
19. data.port<-read.csv('D:\\工作文件夹\\R语言与中国地图\\airports.dat',F)
    
20. data.line<-read.csv('D:\\工作文件夹\\R语言与中国地图\\ routes.dat',F)
    
21. chinaport<-str_detect(data.port[,'V4'],"China")
    
22. chinaport <- data.port[chinaport,]
    
23. chinaport<-chinaport[chinaport$V5!='',c('V3','V5','V7','V8','V9')] names(chinaport)<-c('city','code','lan','lon','att')
    
24. chinaport[2,]
    
25. 
    
26. # 找出国内航班
    
27. lineinchina <- (data.line[,'V3'] %in% chinaport$code) & (data.line[,'V5'] %in% chinaport$code)
    
28. 
    
29. chinaline <- data.line[lineinchina,c('V3','V5','V9')]
    
30. names(chinaline)<-c('f','t','a’)
    
31. 
    
32. #构建一个函数，根据机场编码得到经纬度
    
33. findposition <- function(code) {
    
34.     find <- chinaport$code==code
    
35.     x <- chinaport[find,'lon']
    
36.     y <- chinaport[find,'lan']
    
37.     return(data.frame(x,y))
    
38. }
    
39. 
    
40. # 将机场编码转为经纬度
    
41. from <- lapply(as.character(chinaline$source),findposition)
    
42. from <- do.call('rbind',from)
    
43. from$group <- 1:dim(from)[1]
    
44. names(from) <- c('lon','lan','group')
    
45. 
    
46. to <- lapply(as.character(chinaline$destination),findposition)
    
47. to <- do.call('rbind',to)
    
48. to$group <-1:dim(to)[1]
    
49. names(to) <-c('lon','lan','group')
    
50. data.line <- rbind(from,to)
    
51. temp<- data.line[data.line$group<100,]
    
52. 
    
53. 
    
54. #作图
    
55. par(bg=”gray”)
    
56. plot(x,col="black",border="black")
    
57. points(data.line$lon,data.line$lan, pch = 19, col ="white",cex=0.1)
    
58. 
    
59. for(i in 1:5832){segments(data.line[i,1], data.line[i,2], data.line[i+5832,1],data.line[i+5832,2], col= 'blue')
    
60. }
    
61. 
    
62. for(i in 1:100){ gcIntermediate (data.line[i,1], data.line[i,2], data.line[i+100],data.line[i+100,2])
    
63. }

复制代码

1. ggmap(get_googlemap(center = 'china', zoom=4,
   
2.                     maptype='roadmap'),extent='device')+
   
3.     geom_point(data=chinaport,aes(x=lon,y=lan),
   
4.                colour = 'red',alpha=0.8)+
   
5.     geom_line(data=data.line,aes(x=lon,y=lan,group=group),
   
6.               size=0.1,alpha=0.05,color='blue')
   

复制代码

1. par(bg=”gray”)
   
2. plot(x,col="black",border="black")
   
3. points(data.line$lon,data.line$lan, pch = 19, col ="white",cex=0.2)
   
4. 
   
5. for(i in 1:1000){segments(data.line[i,1], data.line[i,2], data.line[i+1000,1],data.line[i+1000,2], col= 'blue' )
   
6. }

复制代码

1. #弧线及其颜色处理
   
2. library(maps)
   
3. library(maptools)
   
4. library(geosphere)
   
5. map("china",col="gray20",bg="black",mar=c(0,0,0,0))
   
6. title("中国航线图", sub = " SimSun",
   
7.       cex.main =1, col.main= "gray",)
   
8. points(data.line$lon,data.line$lan, pch = 19, col ="white",cex=0.1)
   
9. 
   
10. from <- lapply(as.character(linecount$source),findposition)
    
11. from <- do.call('rbind',from)
    
12. to <- lapply(as.character(linecount$destination),findposition)
    
13. to <- do.call('rbind',to)
    
14. 
    
15. 
    
16. map("china",col="gray20",bg="black",mar=c(0,0,0,0))
    
17. title("中国航线图", sub = " SimSun",
    
18.       cex.main =1, col.main= "gray",)
    
19. points(data.line$lon,data.line$lan, pch = 19, col ="white",cex=0.1)
    
20. #颜色调整
    
21. library(RColorBrewer)
    
22. display.brewer.all(type = "seq")
    
23. display.brewer.all()
    
24. 
    
25. color<-brewer.pal(9,"PRGn")
    
26. color<-colors()[grep("gray",colors())]
    
27. 
    
28. map("china",col="gray20",bg="black",mar=c(0,0,0,0))
    
29. title("中国航线图", sub = " SimSun", col.main= "gray",)
    
30. points(data.line$lon,data.line$lan, pch = 19, col ="white",cex=0.1)
    
31. 
    
32. 
    
33. for(i in 1:1304)
    
34. { n= ceiling(18/linecount[i,3])
    
35. lines(gcIntermediate(from[i,], to[i,], n=100,addStartEnd=T) ,col="blue")
    
36. }

复制代码

R语言交流QQ群：99598210（欢迎加入）