R 使用ggplot2/geom_sf打印时固定osm道路中的间隙

R 使用ggplot2/geom_sf打印时固定osm道路中的间隙,r,ggplot2,openstreetmap,sf,overpass-api,R,Ggplot2,Openstreetmap,Sf,Overpass Api,当使用osmdata绘制从osm下载的道路数据时,当使用geom\u sf中较大的size值时,生成的图中有间隙(见图) 这是一个使用伦敦西南部一段道路的可复制示例。打印时如何删除线中的白色间隙 library(tidyverse) library(sf) library(osmdata) # define bounding box for osm data my_bbox <- matrix(c(-0.2605616, -0.2605616, -0.2004

当使用
osmdata
绘制从osm下载的道路数据时,当使用
geom\u sf
中较大的
size
值时,生成的图中有间隙(见图)

这是一个使用伦敦西南部一段道路的可复制示例。打印时如何删除线中的白色间隙

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)

理想的解决方案是将
lineend=“round”
传递到
geom_sf
,然后将其传递到
geom_路径
(实际上是它的底层)以使线的末端变圆,从而产生轻微的重叠和平滑的外观。但遗憾的是,这不起作用:

ggplot(osm_道路、辅助道路、sf$osm_线路)+
几何尺寸(尺寸=4,lineend=“圆形”)
#>警告:忽略未知参数:lineend

我已经在GitHub上提交了文件,但由于ggplot刚刚发布,任何修复都不会在一段时间内进入CRAN

同时,解决方法包括使用
st_坐标
从几何体列中提取路径。强制为数据帧的结果矩阵可以用
geom_path
绘制,它愉快地接受
lineend
参数:

osm\u道路\u次要\u sf$osm\u线路%>%
st_坐标()%>%
as.data.frame()%>%
ggplot(aes(X,Y,组=L1))+
几何路径(大小=4,lineend=“round”)+
合作基金(crs=4326)

将颜色更改为合适的灰色,以获得更像geom_sf的外观

一种更自然的方法是将线段合并成连续的直线,当然,这些直线没有间隙
st_line_merge
完成了繁重的工作,但您需要事先将它们聚合为多行,以便它具有必要的数据:

osm\u道路\u次要\u sf$osm\u线路%>%
st_union()%>%
st_line_merge()%>%
ggplot()+
geom_sf(尺寸=4)

请注意,这在很大程度上是好的,但不是完全好的。线内的间隙消失了,但是
st_line_join
不知道如何修复三向交叉点,因此仍然有一个很小的间隙。如果你的真实数据有很多这样的交叉点(这是很有可能的),这种方法将不会产生好的结果

最后一种方法是简单使用基本sf打印,默认为圆线端点:

绘图(osm_道路_次要道路_sf$osm_线$geometry,lwd=10)


这种方法是否可行取决于绘图还有哪些工作要做,以及基本绘图的方便程度。

提供
lineend
参数
“round”
可能会给您提供最好的连接。@alistaire谢谢,
geom_sf
看起来没有
lineend
参数,虽然是的,它应该传递到
geom_line
,但它看起来工作不正常。如果您
GG以高分辨率保存图像,间隙是否仍然存在?@camille turning
dpi
up没有区别。应该可以将
geom_路径
解决方案与
coord_sf()
相结合,以获得具有适当分划和轴刻度标签的结果。
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