[R-br] Recorte de mapa de MG

[Paulo Justiniano]

Só vejo duas alternativas

1. na linha do que voce está tentando
entratanto se tiver como importar de algum lugar uma coluna
adicional que já tenha a subregião do estado isto ficaria mais facil

2. A outra   seria verificar se consegue
definir a zona da mata a partir de coordenadas de um poligono
(pode ser aproximado)

Com isto veria quais os centroides de menuicipios que estariam dentro
deste poligono
coordinates(OBJETO)
te dá as goordenadas e basta um poligono aproximado
qua contenha os centroides dos municipios
e com iuso seleciona subset com over()




Paulo Justiniano Ribeiro Jr
LEG (Laboratorio de Estatistica e Geoinformacao)
Universidade Federal do Parana
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR  -  Brasil
Tel: (+55) 41 3361 3573
VOIP: (+55) (41) (3361 3600) 1053 1066
Fax: (+55) 41 3361 3141
e-mail: paulojus AT  ufpr  br
http://www.leg.ufpr.br/~paulojus


On Mon, 14 Dec 2015, Andre Oliveira wrote:

> Pessoal plotei o mapa de MG, gostaria de  selecionar  e colocar em cima só as cidades da zona da mata de MG! São em torno de 75 cidades das 853.
> Estou usando o comando subset. Tem algo mais simples de selecionar estas 75 cidades?
> 
> require(rgdal)
> mg=readOGR("31MUE250GC_SIR.shp",layer="31MUE250GC_SIR",encoding = "latin1")
> plot(mg,axes=T)
> zmata=subset(mg, NM_MUNICIP%in%c("ABRE CAMPO","ACAIACA","ALVINOPOLIS","AMPARO DO SERRA"))
> plot(zmata,col="blue",add=TRUE)
> 
> 
> André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito
> Santo.  IFES
> 
>