<div dir="ltr">Ciao <br><div class="gmail_extra"><br><div class="gmail_quote">Il giorno 12 luglio 2017 20:30, SC <span dir="ltr"><<a href="mailto:elyparker1@gmail.com" target="_blank">elyparker1@gmail.com</a>></span> ha scritto:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF"><span class="gmail-">
    <div class="gmail-m_-2291815349748193989moz-cite-prefix">Il 11/07/2017 19:29, Umberto Zulian ha
      scritto:<br>
    </div>
    <blockquote type="cite">
      
      
      <div id="gmail-m_-2291815349748193989divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
        <p>Forse perche' i valori geometrici (geom) tra due layer con
          proiezioni diverse non sono sovrapponibili?</p>
        <br>
      </div>
    </blockquote></span>
    in che senso? quando li vedi nel canvas di qgis le vedi sovrapposte
    <br>
    <br>
    forse vuoi dire che da codice bisogna effettuare una riproiezione?<br></div></blockquote><div><br></div><div>Beh certo, come hai detto  giustamente tu, *nel canvas* si vedono sovrapposti, ma qui sei nella console e non gli hai mai detto di riproiettare uno dei due layer.</div><div>Io ho riscontrato un problema simile con una funzione che mi serviva per contare gli oggetit in un poligono: lo scoglio su cui mi ero arenato se non sbaglio era che non riuscivo a recuperare il crs di un layer, anche se dovrebbe essere:</div><div> layer.crs()</div><div><br></div><div>Pertanto dovresti applicare una trasformazione di cocordinate [0] per ciascuna feature prima di confrontarle, ma usando la riproiezione al volo non mi fiderei troppo delle aree risultanti. </div><div><br></div><div>E' sicuramente molto meglio lavorare in un unico geodatabase Postgis o Sqlite e creare delle viste di intersezione. Potresti scrivere del codice che importa tutti i layer i nun db sqlite, crea la vista e la esporta come layer temporaneo... mi sembra molto più semplice e gestibile</div><div><br></div><div>amefad</div><div>[0] <a href="https://qgis.org/api/classQgsCoordinateTransform.html">https://qgis.org/api/classQgsCoordinateTransform.html</a></div></div></div></div>