Create a gist now

Instantly share code, notes, and snippets.

What would you like to do?
plot ESRI shapefiles - vanilla
library(maptools)
# read the shapefile data
shpfile1 <- "POL_adm1.shp" # 2 - voivodeship-level
sh1 <- readShapePoly(shpfile1)
shpfile2 <- "POL_adm2.shp" # 2 - powiat-level
sh2 <- readShapePoly(shpfile2)
# create the overlapping plots
plot(sh2)
plot(sh1,
border = "red",
add=TRUE)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment