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