Skip to content

Instantly share code, notes, and snippets.

@erzk
Created April 1, 2016 22:08
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
Star You must be signed in to star a gist
Save erzk/b41acc653a47b5c5cc55d55a837b8d85 to your computer and use it in GitHub Desktop.
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