[R-br] Erro: Intersecção entre dois shapefiles no R

[Thiago V. dos Santos]

  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 em gmail.com>
To: r-br em listas.c3sl.ufpr.br; r-ecologia em 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 em 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 em 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.