Uma possibilidade é o uso das funcionalidades do pacote maptools que voce já está usando As funcoes relevantes seriam require(maptools) readShapePolys() gpclibPermit() # pode ser evitado se voce tiver a biblioteca rgeos instalada unionSpatialPolygons() O aRT (www.leg.ufpr.br/aRT) tb tem funcoes para isto A ideia básica é que na funcao de uniao voce passa os ID's que quer unir, ele vai gerar um novo cunjunto de IDs no objeto resultante On Fri, 9 Sep 2011, Jose Vilton Costa wrote:
Pessoal, Gostaria de saber como faço para unir poligonos de um mapa e mante-los agregado no mapa. Explicando melhor, o meu mapa está subdividido em 49 regiões (id=1...49) e preciso unir os poligonos "24" e "36" em um único poligono, bem como, os poligonos "26" e "46". Desta forma o meu mapa deverá ter 47 regiões. No final de tudo, gostaria de ter um novo id para os poligonos que foram unidos". Para ler o meu mapa estou usando o seguinte código: mapa=readShapePoly("CS-region.shp",ID="ID", proj4string=CRS("+proj=utm +zone=23 +datum=SAD69")) Att, José Vilton