
14 Fev
2013
14 Fev
'13
13:02
Alexandro, seu codigo nao e' reproduzivel a partir da linha ###sp = readShapePoly("C:/r/shapefile-sp-ibge/35mu500gc.shp") mas suponha 'sp' um 'SpatialPolygonsDataFrame' com 645 areas Vai ficar realmente muito ruim visualizar esses 645 nomes + um valor! Se quiser tentar use coo <- coordinates(sp) ### centroids text(coo[,1], coo[,2], paste(cidades, variavel)) ##ou text(coo[,1], coo[,2], paste(cidades, variavel, sep='\n')) Uma alternativa e' usar plotGoogleMaps() do pacote 'plotGoogleMaps' para plotar um subconjunto de areas, suponha que 'j' seja um vetor identificando esse subconjunto plot(sp[j, ]) ### plotar area 'i' e suas vizinhas viz <- poly2nb(sp) plot(sp[c(i, viz[[i]]),])