r-下面多边形的边界线

Question

我是{tmap}软件包的新手，很难创建正确的路线图。这是我的代码(不能立即运行，因为我在磁盘上有一个导出的OpenStreetMap.osm文件(。

# Load OpenStreetMap (OSM) polygon data of the Holden area from disk
holden_bbox_shapes = sf::read_sf('./localdata/holden.osm', 'multipolygons')
# Load OSM line data from disk. These are almost all roads.
holden_lines = sf::read_sf('holden.osm', 'lines')
# Get just the administrative town of holden
holden_town_polygon = dplyr::filter(holden_bbox_shapes, name == 'Holden')
tm_shape(holden_town_polygon) +
tm_fill() +
tm_shape(holden_lines) +
tm_lines()

图像结果在这里，其中道路超出了holden_town_polygon的可见区域，但在我想象的holden_town_polygon对象的边界框内。我不知道为什么会发生这种情况，也不知道如何配置它。我尝试在所有tm_*函数中使用bbox参数进行finaling，但没有什么能改变结果。

我只想看到阴影多边形内的道路部分。这样做的正确方法是什么？

Answer 1

您可以从sf包中获得使用st_intersection。

它看起来像这样：

library(sf)
inside_lines <- st_intersection(holden_lines, holden_town_polygon)
tm_shape(holden_town_polygon) +
tm_fill() +
tm_shape(inside_lines) +
tm_lines()