使用ggplot2 / geom_sf进行绘图时，修复osm道路中的间隙

Question

使用osmdata绘制从osm下载的道路数据时，如果在size中使用较大的geom_sf值，则所得的图会有空白（参见图像）。

以下是使用SW伦敦的一段路的可复制示例。绘图时如何去除线条中的白色间隙？

library(tidyverse)
library(sf)
library(osmdata)

# define bounding box for osm data
my_bbox <- 
  matrix(c(-0.2605616, -0.2605616,
           -0.2004485, -0.2004485,
           -0.2605616, 51.4689943,
           51.4288980, 51.4288980,
           51.4689943, 51.4689943),
         ncol = 2)
bbox_sf <- st_geometry(st_polygon(x = list(my_bbox)))
st_crs(bbox_sf) <- 4326

#get osm road data for bounding box
osm_roads_secondary_sf <- 
  opq(bbox = st_bbox(bbox_sf)) %>%
  add_osm_feature(key = 'highway', value = 'secondary') %>%
  osmdata_sf()

ggplot() + 
  geom_sf(data=osm_roads_secondary_sf$osm_lines,size=4)

会话信息：

R version 3.5.0 (2018-04-23)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS High Sierra 10.13.6

other attached packages:
 [1] osmdata_0.0.7        sf_0.6-3             forcats_0.3.0        
stringr_1.3.1       
 [5] dplyr_0.7.5          purrr_0.2.5          readr_1.1.1          
tidyr_0.8.1         
 [9] tibble_1.4.2         ggplot2_3.0.0        tidyverse_1.2.1.9000

Answer 1

理想的解决方案是将lineend = "round"传递给geom_sf，然后将其传递给geom_path（实际上是其下一层）以使行尾四舍五入，从而导致重叠且外观光滑。遗憾的是，这行不通：

ggplot(osm_roads_secondary_sf$osm_lines) + 
    geom_sf(size = 4, lineend = "round")
#> Warning: Ignoring unknown parameters: lineend

我已经在GitHub上提交了an issue，但是由于ggplot刚刚发布，因此任何修复都不会在一段时间内发布到CRAN。

同时，解决方法包括使用st_coordinates从“几何”列中提取路径。可以使用geom_path绘制强制转换为数据帧的结果矩阵，该矩阵愉快地接受lineend参数：

osm_roads_secondary_sf$osm_lines %>% 
    st_coordinates() %>% 
    as.data.frame() %>% 
    ggplot(aes(X, Y, group = L1)) + 
    geom_path(size = 4, lineend = "round") + 
    coord_sf(crs = 4326)

将颜色更改为合适的灰色阴影，以使外观更像geom_sf。

一种更简单的方法是将线段合并为连续的线，这些线当然没有间隙。 st_line_merge的工作量很大，但是您需要事先将它们汇总到多行中，以便获得必要的数据：

osm_roads_secondary_sf$osm_lines %>% 
    st_union() %>% 
    st_line_merge() %>% 
    ggplot() + 
    geom_sf(size = 4)

请注意，这主要是但并非完全更好。线之间的间隙消失了，但是st_line_join不知道如何固定三向交叉点，因此那里仍然有一个很小的间隙。如果您的真实数据有很多这样的交集（这很有可能），那么这种方法将不会产生良好的结果。

最后一种方法是简单地使用基本sf绘图，默认为圆线末端：

plot(osm_roads_secondary_sf$osm_lines$geometry, lwd = 10)

这种方法是否可行，取决于该图还有哪些其他工作要做，以及您对基本图的适应程度如何。