中国国内航线信息的可视化

 

上图是对国内机场和航线信息进行了一个简单的可视化。圆点表示了中国163个机场的位置，线条显示了5381条航线。之前曾在这个网站上见到了作者用R语言来对全世界的航线进行可视化。正所谓见贤思齐，本图就是模仿山寨的结果。但是这个图的生成没有原文那么复杂，所用到的地理图形包和步骤也与原例略有不同，比较失败的是没有展现出原图的夜景效果。具体实施的步骤如下：

    从这个网站下载到机场数据和航线数据；
    从中挑选出中国的机场和国内航线，并加以整理；
    用ggmap包读取谷歌地图；
    将机场和航线信息绘制在地图上。

 

 

library(ggmap)
data.port <- read.csv( ' d:\\airports.dat ',F)
data.line <- read.csv( ' d:\\routes.dat ',F)
library(stringr)
# 找到中国的机场
portinchina <- str_detect(data.port[, ' V4 '],  " China")
chinaport <- data.port[portinchina,]
# 去除少数几个没有编号的机场
chinaport <-chinaport[chinaport$V5!= '',
                      c( ' V3 ', ' V5 ', ' V7 ', ' V8 ', ' V9 ')]
names(chinaport) <- c( ' city ', ' code ', ' lan ', ' lon ', ' att ')
 
# 找出国内航班
lineinchina <- (data.line[, ' V3 '] %in% chinaport$code) & (data.line[, ' V5 '] %in% chinaport$code)
chinaline <- data.line[lineinchina,c( ' V3 ', ' V5 ', ' V9 ')]
names(chinaline) <- c( ' source ', ' destination ', ' equipment ')
 
#构建一个函数，根据机场编码得到经纬度
findposition <-  function(code) {
    find <- chinaport$code==code
    x <- chinaport[find, ' lon ']
    y <- chinaport[find, ' lan ']
    return(data.frame(x,y))
}
 
# 将机场编码转为经纬度
from <- lapply(as.character(chinaline$source),findposition)
from <- do.call( ' rbind ',from)
from$group <-  1:dim(from)[ 1]
names(from) <- c( ' lon ', ' lan ', ' group ')
 
to <- lapply(as.character(chinaline$destination),findposition)
to <- do.call( ' rbind ',to)
to$group <- 1:dim(to)[ 1]
names(to) <-c( ' lon ', ' lan ', ' group ')
data.line <- rbind(from,to)
temp<- data.line[data.line$group< 100,]
# 用ggmap包从google读取地图数据，并将之前的数据标注在地图上。
ggmap(get_googlemap(center =  ' china ', zoom= 4,
                    maptype= ' roadmap '),extent= ' device ')+
    geom_point(data=chinaport,aes(x=lon,y=lan),
               colour =  ' red4 ',alpha= 0.8)+
    geom_line(data=data.line,aes(x=lon,y=lan,group=group),
              size= 0.1,alpha= 0.05,color= ' red4 ') 

 

 

有兴趣的同学还可以进一步研究，哪个航线最远，分布如何，哪条航线最多，航线飞机档次如何，哪个机场海拔最高，哪个机场最忙.

from: http://xccds1977.blogspot.com/2012/07/blog-post_26.html

转载于:https://www.cnblogs.com/vivizhyy/archive/2012/09/11/2680270.html