Я не знаю, какую географическую область c или водный объект вы выбрали на этом сайте. Я загрузил Цифровой Файл Границы Провинций / территорий. Я прочитал шейп-файл Briti sh Columbia с использованием следующей команды
library(rgdal)
df <- readOGR(dsn = "C:\\Users\\User\\Desktop\\gpr_000a11a_e", layer = "gpr_000a11a_e")
#To see the data head
head(df@data)
#To see a particular column of the shapefile
df@data$PRENAME
#Subsetting British Columbia from the shapefile
BC = subset(df, PRENAME=="British Columbia")
#Plotting it
plot(BC)