[Gfoss] catasto e aree fabbricabili: calcolo area differente in Postgis e Qgis

[Amedeo Fadini]

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
-------------- parte successiva --------------
Un allegato HTML ? stato rimosso...
URL: <http://lists.gfoss.it/pipermail/gfoss/attachments/20121021/f13fd49c/attachment.html>