这种好看的航线图怎么绘制？超简单~~

大家好，欢迎来到IT知识分享网。

今天这篇推文小编给大家介绍一个我一直想绘制的可视化作品-航线地图，网上有很多类似的教程，但在小编看来，都或多或少有些不完美，小编这次就对其进行完善和美化，主要内容如下：

• 航线地图所需数据处理
• 航线地图绘制完善美化
• 更多详细的数据可视化教程，可订阅我们的店铺课程：
• 

航线地图所需数据处理

这一部分主要是对数据进行处理，我们直接给出代码步骤，文中用到的数据，文末都有免费获取方式。

• 导入所需包

library(tidyverse) library(data.table) library(geosphere) library(grid) library(jpeg) library(plyr) library(ggtext) 

• 数据导入和处理

# 导入数据 airlines <- fread("airlines.dat", sep = ",", skip = 1) airports <- fread("airports.dat", sep = ",") routes <- fread("routes.dat", sep = ",") # 使用谷歌夜间影像图片作为背景图 earth <- readJPEG("BlackMarble_2016_01deg.jpg", native = TRUE) earth <- rasterGrob(earth, interpolate = TRUE) # 数据转换和整理 #添加列名 colnames(airlines) <- c("airline_id", "name", "alias", "iata", "icao", "callisign", "country", "active") colnames(airports) <- c("airport_id", "name", "city", "country","iata", "icao", "latitude", "longitude","altitude", "timezone","dst","tz_database_time_zone","type", "source") colnames(routes) <- c("airline", "airline_id", "source_airport", "source_airport_id","destination_airport","destination_airport_id","codeshare", "stops","equipment") #类型转换 routes$airline_id <- as.numeric(routes$airline_id) # airlines与routes数据融合 flights <- left_join(routes, airlines, by="airline_id") # flights与airports数据融合 airports_orig <- airports[,c(5,7,8)] colnames(airports_orig) <- c("source_airport","source_airport_lat", "source_airport_long") airports_dest <- airports[, c(5, 7, 8)] colnames(airports_dest) <- c("destination_airport", "destination_airport_lat", "destination_airport_long") flights <- left_join(flights, airports_orig, by = "source_airport") flights <- left_join(flights, airports_dest, by = "destination_airport") #剔除缺失值 flights <- na.omit(flights, cols = c("source_airport_long", "source_airport_lat", "destination_airport_long", "destination_airport_lat")) # 按航空公司拆分数据集 flights_split <- split(flights, flights$name) # 计算两点之间的距离(geosphere包计算) flights_all <- lapply(flights_split, function(x) gcIntermediate(x[, c("source_airport_long", "source_airport_lat")], x[, c("destination_airport_long", "destination_airport_lat")], n=100, breakAtDateLine = FALSE, addStartEnd = TRUE, sp = TRUE)) # 转换为数据框 flights_fortified <- lapply(flights_all, function(x) ldply(x@lines, fortify)) # Unsplit lists flights_fortified <- do.call("rbind", flights_fortified) flights_fortified$name <- rownames(flights_fortified) flights_fortified$name <- gsub("\\..*", "", flights_fortified$name) flights_points <- flights_fortified %>%   group_by(group) %>%   filter(row_number() == 1 | row_number() == n()) 

由于绘制全部机场的数据较大(小编的本子太菜，耗时太长)，这里我们选择两个机场线路进行绘制，首先我们进行站点筛选和对应的数据选择：

flights_subset <- c("Delta Air Lines","Lufthansa") flights_subset <- flights_fortified[flights_fortified$name %in% flights_subset, ] flights_subset_points <- flights_subset%>%   group_by(group)%>%   filter(row_number()==1|row_number()==n()) 

接下来，我们就将选择的数据进行可视化展示。

航线地图绘制完善美化

在已有的教程中，可视化结果都不够完美(或颜色单一、或没有明显图例)，这里小编进行了必要绘图元素的添加，代码如下：

library(extrafont) # 添加字体 windowsFonts(Cinzel = windowsFont("Cinzel")) fonts() # 设置颜色 color_set <- c("#9b5de5", "#f15bb5") # 可视化绘制 light_map2 <- ggplot() +   annotation_custom(earth, xmin = -180, xmax = 180, ymin = -90, ymax = 90) +   geom_path(data = flights_subset,aes(long, lat, group = id, colour = name),              alpha = 0.1, lineend="round",lwd=0.1) +    geom_point(data = flights_subset_points, aes(long, lat), alpha = 0.8, size = 0.1, colour = "white") +   scale_color_manual(values =color_set,name="Airlines")+   labs(     title = "Example of charts about <span style='color:#D20F26'>Flights at Night</span>",     subtitle = "processed charts with <span style='color:#1A73E8'>geosphere library</span>",     caption = "Visualization by <span style='color:#DD6449'>DataCharm</span>") +   theme_minimal(base_family = "Cinzel")+   theme(         plot.title = element_markdown(hjust = 0.5,vjust = .5,color = "grey60",                                       face = 'bold',size = 22, margin = margin(t = 5, b = 5)),         plot.subtitle = element_markdown(color = "grey60",hjust = 0,vjust = .5,size=12,                                         margin = margin(t = 3, b = 3)),         plot.caption = element_markdown(color = "grey60",face = 'bold',size = 8),                text = element_text(family = "Cinzel"),         panel.background = element_rect(fill='#05050F',colour='#05050F'),          plot.background = element_rect(fill='#05050F',colour='#05050F'),         panel.grid.major = element_blank(),          panel.grid.minor = element_blank(),         axis.text.x = element_blank(),         axis.text.y = element_blank(),         legend.position = c(0.1,0.15),         legend.box.margin = margin(0.2, 0, 0, 0, "cm"),          legend.box.spacing = unit(0.1, "cm"),          legend.key.size = unit(1.5, "lines"),         legend.background = element_rect(colour = NA, fill = NA),         legend.key = element_rect(colour = NA, fill = NA, size = 10),         legend.title = element_text(family = "Cinzel",                                      color = "grey60",                                     size = 10,                                     face = "bold"),         legend.text = element_text(family = "Cinzel",                                     color = "grey60",                                    size = 9,                                    face = "bold")         ) +      guides(color = guide_legend(title.position = "top",                               title.hjust = 0.5, nrow =2 ,                              label.position = "right",                              override.aes = list(size = 0.5,                                                  alpha = 0.8)))  

大家可以重点看theme() 关于图例 (legend) 的参数设置和guides() 中对color图例映射的定制化设置，注意，这里用到了override.aes(),小伙伴们需要掌握哦~~ 可视化结果如下：

Two airline charts example

后续：后来小编试着将结果保存成pdf格式，保存时间较短，但存在保存时字体不存在错误，可通过如下样啦进行解决：

library(tidyverse) library(extrafont) fonts() df <- data.frame(x = c(1:10), y = c(1:10)) # Dummy data  windowsFonts(Calibri = windowsFont("Calibri")) plot <- ggplot(df, aes(x, y)) + # Dummy plot   geom_point() +   labs(title = "Correct font in R, NO fonts at all in pdf :-(") +   theme(text = element_text(size=15, family= "Tw Cen MT Condensed Extra Bold")) ggsave("df_plot.pdf",         plot = plot,         device = cairo_pdf,         dpi = 320) 

注意ggsave()部分，小伙伴们可尝试绘制多个机场的航线图，方法都是一样的哈~~

总结

今天的这篇推文算是满足小编一直想绘制的一幅可视化作品，当然，也有读者私信小编要求绘制类似图表，感兴趣的小伙伴可在获取数据后自行进行练习再现哦~~