我遇到了一个棘手的问题。
我正在尝试使用sp
包,使用功能st_difference(st_union(x),st_union(y))
(或其中的变体)或st_intersection
函数来填充一个简单的多边形剪辑,以哪种效果最好。
虽然使用两个表面多边形很容易,但是我需要将其裁剪到一个可怕的大型下载的Large SpatialPolygonsDataFrame,它只是英国的shapefile,可从以下位置下载: https://gadm.org/download_country_v3.html
> str(uk)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
..@ data :'data.frame': 1 obs. of 70 variables:
.. ..$ ID_0 : Factor w/ 1 level "242": 1
.. ..$ ISO : Factor w/ 1 level "GBR": 1
.. ..$ NAME_0 : Factor w/ 1 level "United Kingdom": 1
# .....etc.
#
> str(box)
sfc_POLYGON of length 1; first list element: List of 1
$ : num [1:5, 1:2] -7.237 0.126 0.126 -7.237 -7.237 ...
- attr(*, "class")= chr [1:3] "XY" "POLYGON" "sfg"
我想将box
(蓝色)剪切到uk
,原因是在英国(以及法国)的小叶中渲染shapefile花费的时间太长。 。
答案 0 :(得分:2)
这可能会做什么?
sf::st_intersection(UK, box)
完整代码
library(sf)
UK <- readRDS("./gadm36_GBR_0_sf.rds")
#create box since it was not provided in question
box <- c("POLYGON((-7.237 48,0 48,0 52,-7.237 52, -7.237 48))") %>%
st_as_sfc(crs = "+proj=longlat +datum=WGS84")
mapview::mapview(list(UK,box))
mapview::mapview( st_intersection(UK, box) )
如果要使用英国的shapefile剪切框,请使用st_difference()
mapview::mapview( st_difference (box, UK) )