graphics.off()
windows(record=TRUE)

library(maptools)
library(ggmap)

precisa.baixar = !('brasilmun2010_new.shp' %in% dir())

if (precisa.baixar) {
  download.file('http://ufrn.schmert.net/dados/brasilmun2010_new.shp',
                destfile='brasilmun2010_new.shp',mode='wb')
  download.file('http://ufrn.schmert.net/dados/brasilmun2010_new.shx',
                destfile='brasilmun2010_new.shx',mode='wb')
  download.file('http://ufrn.schmert.net/dados/brasilmun2010_new.dbf',
                destfile='brasilmun2010_new.dbf',mode='wb')
}

mapa = readShapePoly('brasilmun2010_new')

##############################################################

fundo = qmap('Brazil', zoom=4, maptype='roadmap')
fundo

# lista de estados
estados = unique( mapa@data$NAME2_)

# para cada estado, acresenta uma camada ao mapa de fundo, com uma cor aleatória
# do arco-íris (rainbow)

for (nome in estados) {
  indices = which( mapa@data$NAME2_ == nome)
  est.df  = fortify( mapa[indices,])     # converter num data.frame (ggmap precisa)
  est.cor = sample( rainbow(25), 1)
  fundo = fundo + geom_polygon(data=est.df, aes(x=long,y=lat,group=id), fill=est.cor,
                      col=NA,alpha=.50)
}

fundo