<div dir="ltr"><div><div><div><div>Bom dia Milú,<br><br></div>Se os dados são da CAOP, muito provavelmente o resultado correto é o da ferramenta do Processing, porque existem error topológicos (self-intersection) na CAOP (pelo menos em algumas versões que tenho testado). Assim como na COS.<br><br></div>A questão é que esses erros, no modelo de dados OGC, não são erros no modelo de dados da ESRI. Ou seja, no primeiro, os "buracos" devem ser ilhas; já no segundo são admitidas as self-intersections. Pessoalmente considero a abordagem topológica do modelo OGC mais correta e mais lógica, mas a verdade é que existem os dois modelos. Vê aqui graficamente a diferença:<br><br><a href="http://www.h2gis.org/docs/1.3/ST_IsValidReason/">http://www.h2gis.org/docs/1.3/ST_IsValidReason/</a></div><div><br></div><div>No Postgis corrige-se facilmente esse problema com o ST_MakeValid.</div><div><br></div><div><br></div>Cumprimentos,<br></div>Pedro Venâncio<br><div><div><br><br><br><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">No dia 7 de dezembro de 2017 às 10:24, MiluCaridadeGaspar <span dir="ltr"><<a href="mailto:milugaspar@gmail.com" target="_blank">milugaspar@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Bom dia a todos,<br>
<br>
No desenvolvimento de um plugin preciso de usar a ferramenta de<br>
geoprocessamento intersecção.<br>
Em ambiente desktop, consigo aceder a duas ferramentas diferentes, uma no<br>
menu Vector/Ferramentas de Geoprocessamento/Intersectar que me devolve<br>
resultados satisfatórios e outra que aparece na Caixa de Ferramentas -<br>
Intersecção - e que me devolve o erro /"Erro no geoprocessamento GEOS: Um ou<br>
mais elementos de entrada têm geometria inválida"/.<br>
<br>
Em python a instrução que uso é processing.runalg("qgis:<wbr>intersection",<br>
.......)<br>
<br>
/processing.runalg("qgis:<wbr>intersection",LayerCAOP,<wbr>Temp01,saida)/<br>
<br>
E obtenho o mesmo erro de geometria iniválida.<br>
<br>
Nenhum dos layers terá geometria inválida, inclusivamente já corri um<br>
"reparador de geometrias". Para além disso parece-me que  o erro está<br>
relacionado com a complexidade, dimensão ou numero de vértices  dos<br>
shapefiles envolvidos na intersecção...<br>
<br>
Penso que conseguiria contornar o problema se conseguisse aceder via código<br>
python à ferramenta "intersectar" no menu Vector. Sabem como o posso fazer?<br>
Ou qualquer dica/sugestão que me possa ajudar a resolver esta questão?<br>
<br>
Muito obrigada!<br>
<br>
Com os melhores cumprimentos,<br>
Milú<br>
<br>
<br>
<br>
<br>
<br>
<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/QGIS-pt-f5128248.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.<wbr>nabble.com/QGIS-pt-f5128248.<wbr>html</a><br>
______________________________<wbr>_________________<br>
QGIS-pt mailing list<br>
<a href="mailto:QGIS-pt@lists.osgeo.org">QGIS-pt@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/qgis-pt" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/qgis-pt</a><br>
</blockquote></div><br></div>