Estava vendo vossa discussão, tentei rodar o script e não encontra esta função readShapePolys(), e já tenho estas bibliotécas instaladas!

       
> require(maptools)
> readShapePolys()
Erro: não foi possível encontrar a função "readShapePolys"

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

---

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
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.