我有一个安大略省边界的shapefile,内容如下:
library(sf)
library(here)
ontario <- sf::st_read(here::here("data", "messy_data", "Ontario.shp"), quiet = TRUE) %>%
st_transform(4326)
它可以绘制如下:
我只想保留ontario
的{{1}}数据。但是由于它是latitude < 51
对象,所以sf
不起作用。我知道我可以使用dplyr::filter(latitude < 51)
提取坐标,但是如何将它们重新加入数据以过滤出更高的海拔呢?
另外,st_coordinates()
中是否可以使用任何函数进行过滤?我调查了帮助,但到目前为止找不到任何相关内容。