ShapeFile Paraná

require(maps)
require(maptools)
 
map <- readShapePoly("41mu500gc.shp") #lendo o Shapefile
dbf <- read.dbf("41mu500gc.dbf") #lendo o arquivo dbf
plot(map, col=colors()[1:30])