使用ggplot2/geom_sf绘制时,修复渗透道路中的间隙

Question

在绘制从osm下载的道路数据osmdata时,当使用较大的sizein 值时,得到的绘图中存在间隙geom_sf(见图).

这是一个可重复的例子,使用伦敦西南部的一段道路.绘图时如何删除线条中的白色间隙？

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上提交了一个问题,但是由于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)

这种方法是否切实可行取决于该情节还有哪些工作要做,以及你对基础绘图的敏感程度.