Il test del centroide è un test del pirla, scusate, perché il centroide può cadere anche fuori dalla geometria...<br>Posterò una domanda nell ml di sviluppo ma intanto lo chiedo anche qua: come si può sapere l'orientamento di una geometria tramite le api di Qgis?<br>
<br>giovanni<br><br><div class="gmail_quote">Il giorno 15 giugno 2011 12:36, G. Allegri <span dir="ltr"><<a href="mailto:giohappy@gmail.com">giohappy@gmail.com</a>></span> ha scritto:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dall'API di QGis non trovo un metodo per determinare l'orientamento della sequenza di coordinate... o se si tratta di un'isola. Un modo grezzo ma, credo, efficace è lanciare questo script nella console di Python dentro QGis, nel quale faccio un intersect tra la geometria e il suo centroide. Se fossse un'isola dovrebbe tornarmi False... o sbaglio?<br>

<br>iface = qgis.utils.iface<br>lyr = iface.activeLayer()<br>prov = lyr.dataProvider()<br>attrlist = prov.attributeIndexes()<br>prov.select(attrlist)<br>feat = QgsFeature()<br>for i in range(lyr.featureCount()):<br>    prov.nextFeature(feat)<br>

    geom = feat.geometry()<br>    cent = geom.centroid()<br>    dentro = geom.intersects(cent)<br>    if (dentro):<br>        attributes = feat.attributeMap()<br>        print 'La feature %s contiene una geometria piena' % attributes[0].toInt()[0]<br>

    else:<br>        print '->>> La feature %s sembra contenere un'isola' % attributes[0].toInt()[0]<br><br><br>Giovanni<br>        <br><br><br><br><br><div class="gmail_quote">Il giorno 15 giugno 2011 11:00, marco zanieri <span dir="ltr"><<a href="mailto:marcozanieri@gmail.com" target="_blank">marcozanieri@gmail.com</a>></span> ha scritto:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div class="h5">Salve,<br>ho un problema con alcune geometrie areali, avrei la necessità di verificare il senso di digitalizzazione (orario:interno poligono pieno; antiorario: interno poligono vuoto);<br>

esiste quest possibilità in Qgis?<br>
<br>Grazie mille,<br>marco<br clear="all"><font color="#888888"><br>-- <br>            dott. Marco Zanieri<br>   e-mail: <a href="mailto:marcozanieri@gmail.com" target="_blank">marcozanieri@gmail.com</a><br><br>           cartografia tematica<br>

          banche dati territoriali<br>
     sistemi informativi geografici <br>      applicazioni GIS e webGIS<br><br><br><br><br>
</font><br></div></div>_______________________________________________<br>
Iscriviti all'associazione GFOSS.it: <a href="http://www.gfoss.it/drupal/iscrizione" target="_blank">http://www.gfoss.it/drupal/iscrizione</a><br>
<a href="mailto:Gfoss@lists.gfoss.it" target="_blank">Gfoss@lists.gfoss.it</a><br>
<a href="http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss" target="_blank">http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss</a><br>
Questa e' una lista di discussione pubblica aperta a tutti.<br>
Non inviate messaggi commerciali.<br>
I messaggi di questa lista non rispecchiano necessariamente<br>
le posizioni dell'Associazione GFOSS.it.<br>
518 iscritti al 3.6.2011<br></blockquote></div><br><div style="padding:0px;margin-left:0px;margin-top:0px;overflow:hidden;word-wrap:break-word;color:black;font-size:10px;text-align:left;line-height:130%">
</div>
</blockquote></div><br><div style="visibility: hidden; left: -5000px; position: absolute; z-index: 9999; padding: 0px; margin-left: 0px; margin-top: 0px; overflow: hidden; word-wrap: break-word; color: black; font-size: 10px; text-align: left; line-height: 130%;" id="avg_ls_inline_popup">
</div>