Bom, seu problema é mais ou menos isso.
#abrindo pacote pra plotar mapas, assumo que vc ta usando esse, o maptools
library(maptools)
#lendo um mapa pra exemplo
mapa<- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1],
IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))
#mapa
summary(mapa)
#informações aqui
mapa@data
#vamo pega o nome das cidades
mapa@data$NAME
#o problema que vc ta fazendo isso
plot(mapa, axes=TRUE, col=mapa@data$NAME=="Alamance")
plot(mapa, col=mapa@data$NAME=="Avery", add=TRUE)
plot(mapa, col=mapa@data$NAME=="Buncombe", add=TRUE)
#no col vc ta pondo essa informação
mapa@data$NAME=="Alamance"
#que é sempre 1, logo vc colore de preto, 1 na palleta default é
preto, true é igual 1, tem 0 0 0 que é branco e 1 que é preto
sum(mapa@data$NAME=="Alamance")
#vc precisa de numeros diferentes pra cores
#ou o nome da cor pra cidade
#ou faze uma palleta
#mas assim vc tem numeros diferentes
as.numeric(mapa@data$NAME)
#e vai sair cores diferentes, olha esse plot, note o que to usando no col.
plot(mapa, axes=TRUE, col=as.numeric(mapa@data$NAME))
#nesse exemplo vc ta usando essas cores (numeros) pra esses lugares
data.frame(cores=as.numeric(mapa@data$NAME),localidade=mapa@data$NAME)
#então faz um vetor, com o numero de cidades, e escrever "black",
"blue", as cores que vc quer por cidade, ou usa numeros, ou usa uma
palleta pronta, tipo as que vem no pacote Rcolorbrewer.
Bem espero ter ajudado.
Em 25 de junho de 2012 13:00, Wirton Macedo Coutinho
<[email protected]> escreveu:
> Boa tarde pessoal,
>
> Gostaria de saber se alguém tem alguma dica de como colorir alguns
> municipios em um mapa estadual (extensão shp) com cores distintas.
>
> Preciso colorir 7 municipios no mapa do Estado de Goiás com cores distintas
> (cada município com uma cor diferente), mas não estou conseguindo.
>
> Tenho tentado, sem sucesso, da seguinte forma:
>
> plot(mapa, axes=TRUE, col=mapa$Nome_Munic=="Santa Helena de Goiás")
> plot(mapa, col=mapa$Nome_Munic=="Chapadão do Céu", add=TRUE)
> plot(mapa, col=mapa$Nome_Munic=="Montividiu", add=TRUE)
> plot(mapa, col=mapa$Nome_Munic=="Acreúna", add=TRUE)
> plot(mapa, col=mapa$Nome_Munic=="Perolândia", add=TRUE)
> plot(mapa, col=mapa$Nome_Munic=="Palmeiras de Goiás", add=TRUE)
> plot(mapa, col=mapa$Nome_Munic=="Cristalina", add=TRUE)
>
> Apesar de conseguir colorir todos os 7 municipios, com os comando acima,
> todos os 7 ficam com a mesma cor.
>
> Certo de contar com a boa vontade de todos, agradeço antecipadamente.
>
> Att.,
>
> --
> Wirton Macedo Coutinho
> Pesquisador Fitopatologia
> Embrapa Algodão
> Rua Oswaldo Cruz, 1143, Centenário
> Campina Grande PB
> CEP 28428-095
>
> _______________________________________________
> R-br mailing list
> [email protected]
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
> mínimo reproduzível.
--
Grato
Augusto C. A. Ribas
Site Pessoal: http://augustoribas.heliohost.org
Lattes: http://lattes.cnpq.br/7355685961127056
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
