[R-br] União de poligonos

[Paulo Justiniano]

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
> 
>