Compare spatial polygons and keep or delete common boundaries in R
问题 I use to plot maps of Belgium available on GADM (select Belgium) using R. I import and plot the outside boundary of Belgium using : belgium <-readRDS("gadm36_BEL_0_sp.rds") plot(belgium) Which gives me : I import and plot the boundaries of provinces of Belgium using : provinces <-readRDS("gadm36_BEL_2_sp.rds") plot(provinces) Which gives me : What i'm trying to have is a dataframe with boundaries of provinces that are NOT outside boundaries of Belgium : I tried using over() , intersect() ,