Ricardo, Sem acesso aos dados é difícil entender completamente o que está acontecendo. Em questões futuras, por favor forneça (uma amostra d) os dados para que o problema seja reprodutível. Aparentemente, pelo menos um polígono dos seus shapefiles apresenta mal formação. Para resolver isso, a minha dica é usar a função checkPolygonsHoles do pacote maptools. Tipicamente, você digitaria: library(maptools) pls <- slot(wdpa, "polygons") pls1 <- lapply(pls, checkPolygonsHoles) slot(wpda, "polygons") <- pls1 Repita a operação nos outros shapefiles e então tente fazer a interseção novamente. Espero que ajude, Thiago. ________________________________ From: Ricardo Dobrovolski <rdobrovolski@gmail.com> To: r-br@listas.c3sl.ufpr.br; r-ecologia@googlegroups.com Sent: Friday, December 7, 2012 9:21 AM Subject: [R-br] Erro: Intersecção entre dois shapefiles no R Oi Pessoal, Tentei realizar a intersecção entre dois shapes por meio da função gIntersect do pacote rgeos. No entanto, ocorreu o seguinte erro, conforme mostrado abaixo:
wdpa<- readShapePoly(file.choose()) grid_05<- readShapePoly(file.choose()) inter_wdpa_grid<-gIntersection(wdpa,grid_05) Erro em createPolygonsComment(p) : rgeos_PolyCreateComment: orphaned hole, cannot find containing polygon for hole at index 2
Em 6 de dezembro de 2012 20:22, Ricardo Dobrovolski <rdobrovolski@gmail.com> escreveu: Oi Pessoal,
Eu gostaria de saber se vocês conhecem pacotes no R nos quais eu possa fazer operações entre shapefiles como intersecções, clips, etc, de maneira semelhante às operações que se podem realizar no ArcGis.
Atenciosamente,
Ricardo Dobrovolski
_______________________________________________ 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.