search for: st_area

Scotland is the second feature in the UK data, so get it and split this one MULTIPOLYGON feature into individual POLYGONS scot = st_cast(the_uk$geometry[2],"POLYGON") # which is the largest polygon? which.max(st_area(scot)) [1] 1 # the first one. ok... plot(scot[[1]]) # mainland # add the rest of the islands for context, in grey, maybe to show they're outside our study area: for(i in 2:length(scot)){plot(scot[[i]], col="grey", add=TRUE)} There are 2794 polygons in Scotland (according to this...