Обрезка полигона поверхности довольно неприятным большим пространственным фреймом данных полигона

Question

Я столкнулся с хитрой проблемой.

Я пытаюсь выполнить простой многоугольный клип, используя пакет sp, используя либо функцию fuction st_difference(st_union(x),st_union(y)) (или ее вариант), либо st_intersection, в зависимости от того, что лучше всего работает.

Хотя это легко с двумя поверхностными полигонами, мне нужно обрезать его в ужасно большой скачанный Large SpatialPolygonsDataFrame, это просто шейп-файл для Великобритании, скачанный из: https://gadm.org/download_country_v3.html

Шейп-файлы имеют следующий вид (нанесено на leaflet):

    > 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, причина этого в том, что слишком много времени занимает рендеринг шейп-файла из Великобритании (и Франции в этом отношении) в листовку.

Answer 1

Это, вероятно, подойдет?

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) )

обновление

Если вы хотите вырезать коробку с шейп-файлом из Великобритании, используйте st_difference()

mapview::mapview( st_difference (box, UK) )