JdP*_*JdP 6 r transform geospatial r-sf
有没有一种简单的方法可以将简单特征(对象)的边界转换sf为另一个 CRS?
该函数的结果st_bbox()属于 类bbox。无法将st_transform()其转换为另一个 CRS。
我正在使用基于计算的边界框EPSG: 28992:
sf::st_bbox(xmin = 187470, xmax =194587,
ymin = 409753, ymax = 412715,
crs = st_crs(28992))
Run Code Online (Sandbox Code Playgroud)
现在我想把这个盒子改造成EPSG:4326
对象有一个st_as_sfc方法bbox,所以我们可以像这样转换 a bbox:
library(sf)
bb = sf::st_bbox(c(xmin = 187470, xmax =194587,
ymin = 409753, ymax = 412715),
crs = st_crs(28992))
bb_ll = st_bbox(
st_transform(
st_as_sfc(bb),
4326
)
)
# or pipey
library(magrittr)
bb_ll = bb %>%
st_as_sfc() %>%
st_transform(crs = 4326) %>%
st_bbox()
bb_ll
xmin ymin xmax ymax
5.856639 51.675176 5.959866 51.702257
Run Code Online (Sandbox Code Playgroud)