União de poligonos
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
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
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.
erro de digitacao foi com um "s" a amis require(maptools) ?readShapePoly On Mon, 12 Sep 2011, Mauro Sznelwar wrote:
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.
participantes (3)
-
Jose Vilton Costa -
Mauro Sznelwar -
Paulo Justiniano