Ciao a tutti, 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.
finora: - import dei cxf in postgis secondo le indicazioni del wiki [0]: FATTO - script python per importare i dati alfanumerici dell'agenzia del territorio in un db postgres: FATTO - decodifica dei tracciati record documentati qui [1] in modo: FATTO - query per collegare le particella alla rendita, anno per anno e ai proprietari: FATTO - reperire le tabelle di corrispondenza dei fogli urbano/terreno per tutta la provincia di TV: FATTO Ora ci sarebbe da affrontare il problema del calcolo della percentuale fabbricabile di un mappale che ho scelto di fare in una apposita vista: ST_area(particelle.the_geom) as area_mappale ST_area(ST_intersection(particelle.the_geom, zone.the_geom) as area_edificabile e poi: area_edificabile/area_mappale*100 as perc_edificabile Tra i dati alfanumerici del catasto c'è la superficie in ettari, are, centiare che io trasformo in mq Quindi la superificie edificabile che farà testo sarà superficie_catastale * perc_edificabile 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: ad esempio: id:1095616 area_calcolata: 2929.05078125 (postgis) derivato>area: 2.925,565 m² id:1005647 area_calcolata: 12342.7646484375 derivato>area: 1,233 ha id:1556723 area_calcolata: 46504.70703125 derivato>area:4,645 ha 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... 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? ---- [0] http://wiki.gfoss.it/index.php/Catasto [1] http://www.agenziaterritorio.it/sites/territorio/files/servizi/ServiziComuniIstituzioni/PortalePerComuni/ES-23-IS-05_100909.pdf
_______________________________________________ Gfoss@lists.gfoss.it http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. Non inviate messaggi commerciali. I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it. 605 iscritti al 10.7.2012