<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;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>
<br>
<div style="color: rgb(0, 0, 0);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Da:</b> QGIS-it-user <qgis-it-user-bounces@lists.osgeo.org> per conto di SC <elyparker1@gmail.com><br>
<b>Inviato:</b> martedì 11 luglio 2017 18:21<br>
<b>A:</b> qgis-it-user@lists.osgeo.org<br>
<b>Oggetto:</b> [QGIS-it-user] Fwd: intersezione layer e calcolo area e crs differenti</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Salve ho scritto un piccolo script per la console di QGIS per calcolare<br>
<br>
le aree intersecate dei vari geometrie di un livello rispetto l'altro, ai piu non sfuggirà la possibilità di usarlo per per calcolare certificati di destinazione urbanistica rispetto un prg<br>
<br>
c'è un problema quando i layer sono due sistemi di riferimenti<br>
differenti ,il calcolo non funziona, mentre va bene se i liveli sono<br>
sullo stesso crs<br>
<br>
sapete dirmi se è<br>
<br>
1)un bug di qgis<br>
<br>
2) scelta progettuale (perchè?) da ovviare tramite codice<br>
<br>
3)api sbagliate che ho usato si poteva scegliere altra soluzione<br>
<br>
saluti<br>
<br>
<br>
<br>
# /* copyright 2017 by Salvo Caligiore caligiore@gmail.com<br>
<br>
livelli=iface.mapCanvas().layers()<br>
<br>
<br>
def dataGeom(geom,lv):<br>
#calcola le aree delle intersezioni delle features di un determinato<br>
livello rispetto ad una data geometria e restituisce una lista di<br>
feature con il corrispondente valore di area intersecato<br>
#la feauture in lista serve per altri calcoli nella procedura principale<br>
rect = geom.boundingBox()<br>
<br>
request=QgsFeatureRequest(rect)<br>
lista=[]<br>
feat=lv.getFeatures(request)<br>
da = QgsDistanceArea()<br>
for f in feat:<br>
<br>
intersezione=f.geometry().intersection(geom)<br>
<br>
if intersezione==None : continue<br>
else :<br>
a=round(da.convertAreaMeasurement(da.measureArea(intersezione),1),2)<br>
<br>
<br>
<br>
if a>0 : lista.append((f,round(a,2)))<br>
if lv.geometryType()==0: lista.append((f," "))<br>
return lista<br>
<br>
<br>
da = QgsDistanceArea()<br>
#incrocia i livelli fra di loro e interseca ogni feature con le altre<br>
del livello incrociato stampando le informazioni<br>
for x in livelli :<br>
for y in livelli:<br>
for f in y.getFeatures():<br>
<br>
lst=dataGeom(f.geometry(),x)<br>
<br>
for z,a in lst:<br>
areafeature=round(da.convertAreaMeasurement(f.geometry().area(),1),2)<br>
print x.name(), y.name(),"ID",f.id(),"Area<br>
Totale",areafeature , "Area Intersecata", a , "ZID",z.id(),"rap % fra<br>
inter e area f" ,str(round((a/areafeature)*100,2))+'%'<br>
<br>
_______________________________________________<br>
QGIS-it-user mailing list<br>
QGIS-it-user@lists.osgeo.org<br>
<a href="https://lists.osgeo.org/mailman/listinfo/qgis-it-user" id="LPlnk708005" previewremoved="true">https://lists.osgeo.org/mailman/listinfo/qgis-it-user</a>
<div id="LPBorder_GT_14997935580870.6687223979816441" style="margin-bottom: 20px; overflow: auto; width: 100%; text-indent: 0px;">
<table id="LPContainer_14997935580830.949263158208784" role="presentation" cellspacing="0" style="width: 90%; background-color: rgb(255, 255, 255); position: relative; overflow: auto; padding-top: 20px; padding-bottom: 20px; margin-top: 20px; border-top: 1px dotted rgb(200, 200, 200); border-bottom: 1px dotted rgb(200, 200, 200);">
<tbody>
<tr valign="top" style="border-spacing: 0px;">
<td id="TextCell_14997935580850.131520166163267" colspan="2" style="vertical-align: top; position: relative; padding: 0px; display: table-cell;">
<div id="LPRemovePreviewContainer_14997935580850.1292371831367749"></div>
<div id="LPTitle_14997935580850.14031948283608564" style="top: 0px; color: rgb(0, 120, 215); font-weight: normal; font-size: 21px; font-family: wf_segoe-ui_light, "Segoe UI Light", "Segoe WP Light", "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; line-height: 21px;">
<a id="LPUrlAnchor_14997935580860.9761047411196069" href="https://lists.osgeo.org/mailman/listinfo/qgis-it-user" target="_blank" style="text-decoration: none;">Pagina di informazioni della lista QGIS-it-user</a></div>
<div id="LPMetadata_14997935580860.2910241595571892" style="margin: 10px 0px 16px; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 14px;">
lists.osgeo.org</div>
<div id="LPDescription_14997935580870.030209622450099616" style="display: block; color: rgb(102, 102, 102); font-weight: normal; font-family: wf_segoe-ui_normal, "Segoe UI", "Segoe WP", Tahoma, Arial, sans-serif; font-size: 14px; line-height: 20px; max-height: 100px; overflow: hidden;">
Lista di discussione e supporto per gli utenti italiani di QGIS. Per consultare la raccolta dei messaggi precedentemente inviati alla lista, visita gli Archivi ...</div>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<br>
</div>
</span></font></div>
</div>
</body>
</html>