Ciao a tutti, <br>in tempi di start-up, visto che non riesco a stare dietro a lavoro e università, ho pensato anch'io di partire con un nuovo progetto e sto cercando di mettere a punto un metodo efficace per il calcolo delle aree fabbricabili... spero possa diventare presto un plugin di qgis, se qualcuno vuole dare una mano, specie i dipendenti pubblici che conoscono bene la materia è il benvenuto.<br>

<br>finora:<br><br>- import dei cxf in postgis secondo le indicazioni del wiki [0]: FATTO <br>- script python per importare i dati alfanumerici dell'agenzia del territorio in un db postgres: FATTO<br>- decodifica dei tracciati record documentati qui [1] in modo: FATTO<br>

- query per collegare le particella alla rendita, anno per anno e ai proprietari: FATTO<br>- reperire le tabelle di corrispondenza dei fogli urbano/terreno per tutta la provincia di TV: FATTO<br><br><br>Ora ci sarebbe da affrontare il problema del calcolo della percentuale fabbricabile di un mappale che ho scelto di fare in una apposita vista:<br>

<br>ST_area(particelle.the_geom) as area_mappale<br>ST_area(ST_intersection(particelle.the_geom, zone.the_geom) as area_edificabile<br><br>e poi:<br><br>area_edificabile/area_mappale*100 as perc_edificabile<br><br>Tra i dati alfanumerici del catasto c'è la superficie in ettari, are, centiare che io trasformo in mq<br>

Quindi la superificie edificabile che farà testo sarà <br><br>superficie_catastale * perc_edificabile<br><br><br>Il problema è che le aree calcolate in POSTGIS sono differenti da quelle che Qgis calcola nel campo derivato, pensavo dipendesse dal sistema di riferimento (qgis usa 3003 con +twgs84 e postgis senza) ma anche andando a impostare la stessa stringa proj4 in qgis  trovo risultati discordanti:<br>

<br>ad esempio:<br><br>id:1095616<br>area_calcolata: 2929.05078125 (postgis)<br>derivato>area: 2.925,565 m²<br><br><br>id:1005647<br>area_calcolata: 12342.7646484375<br>derivato>area: 1,233 ha<br><br>id:1556723<br>
area_calcolata: 46504.70703125<br>derivato>area:4,645 ha<br><br><br>So bene che parliamo circa di 10 mq ogni ettaro, cioè lo 0,1% ma vorrei sapere se qualcuno mi spiega perché il calcolo è differente...<br>Disabilitando la priezione al volo sembra che la differenza svanisca, ma non è chiaro perché il campo in ettari  ha solo tre decimali.. c'è modo di mostrare sempre in metri quadri?<br>
<br>----<br>[0] <a href="http://wiki.gfoss.it/index.php/Catasto" target="_blank">http://wiki.gfoss.it/index.php/Catasto</a><br>[1] <a href="http://www.agenziaterritorio.it/sites/territorio/files/servizi/ServiziComuniIstituzioni/PortalePerComuni/ES-23-IS-05_100909.pdf" target="_blank">http://www.agenziaterritorio.it/sites/territorio/files/servizi/ServiziComuniIstituzioni/PortalePerComuni/ES-23-IS-05_100909.pdf</a><br>