示例:
bbox <- c(-0.1178, 51.4232, -0.0185, 51.5147) # I know it needs to be sf df object
# we have
df
#> Geometry set for 300 features
#> geometry type: LINESTRING
#> dimension: XY
#> bbox: xmin: -0.113894 ymin: 51.49739 xmax: -0.0764779 ymax: 51.59839
#> epsg (SRID): 4326
#> proj4string: +proj=longlat +datum=WGS84 +no_defs
#> LINESTRING (-0.113894 51.50631, -0.1135137 51.5...
#> LINESTRING (-0.0767875 51.59837, -0.0764779 51....
#> ....
我该怎么做
df[bbox]
并保留bbox中的线串。谢谢。
答案 0 :(得分:0)
在此示例中,使用sf
中的tigris
对象只是为了提高可重复性。我正在使用康涅狄格州纽黑文县的城镇,将其绘制成图。然后使用st_crop
将其裁剪到我组成的边界框内,我相信它是最近才添加到{{1 }}。如果我将bbox作为形状而不是坐标向量,则可以使用sf
。
我没有方便的线串对象,但我认为它的工作方式相同。
st_intersection
library(tidyverse)
library(sf)
# selecting just to limit the amount of data in my sf
ct_sf <- tigris::county_subdivisions(state = "09", county = "09", cb = T, class = "sf") %>%
select(NAME, geometry)
plot(ct_sf)