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.