Ciao a tutti!
Nell'ultimo tempo ho fatto nuovi provi, ma la situazione è la stessa. Perché è
già passato tanto tempo, riprendo il problema: dovevo calcolare il percentuale
del bosco (oppure solo la superficie boscata, perché il percentuale si può
ottenere tramite il calcolatore dei campi) per alcuni comuni/bacini o
suddividendo un comune/bacino in maglie quadrate (1 kmq ciascuna). Sono partita
da due vettoriali (.shp): lo strato dei comuni (o dei bacini o una griglia) e
lo strato dei boschi, ambedue poligonali.
In lista Paolo e Markus mi avevano suggerito di usare ZonalStats, Statistiche
zonali e v.rast.stats. Dopo le prime prove mi sono resa conto che Statistiche
zonali e v.rast.stats erano proprio gli strumenti che cercavo, essendo molto
facile da usarli, però, paragonando con le aree calcolate sia in QGIS (sui
vettoriali) sia in gvSIG (per controllo), i risultati ottenuti erano sbagliati.
Devo precisare che per usare Statistiche zonali e v.rast.stats ho dovuto
trasformare uno degli strati da vettore in raster. Dunque, ho usato: lo strato
poligonale (.shp) dei comuni/o una griglia (1 kmq ciascun quadrato, realizzata
con fTools) e lo strato raster dei boschi.
Risultati:
1. Quando ho usato lo strato raster dei boschi ad una risoluzione di 10
(stabilita in GRASS) le aree calcolate sia con Statistiche zonali (SZ) sia con
v.rast.stats erano 100 volte più piccole rispetto all'area reale (per esempio
un quadrato di 1 kmq = 1 000 000 mq era calcolato come 10 000 mq).
2. Cambiando la risoluzione del raster si ottengono valori differenti: con la
risoluzione 1 (stabilita in GRASS) si ottengono i valori quasi reali, ma non
per tutte le aree, per esempio:
- 1 000 000 mq (valore reale, QGIS, gvSIG) - 1 000 000 mq (calcolato con SZ e
con v.rast.stats - evviva!!!)
- 999 990,76 mq (valore reale, QGIS, gvSIG) - 999 992 mq (calcolato con SZ e
con v.rast.stats)
- 762 933,02 mq (valore reale, QGIS, gvSIG)- 763 336 mq (calcolato con SZ e con
v.rast.stats)
- 335 316,97 mq (valore reale, QGIS, gvSIG)- 335 581 mq (calcolato con SZ e con
v.rast.stats)
- 241 915,87 mq (valore reale, QGIS, gvSIG)- 242 047 mq (calcolato con SZ e con
v.rast.stats)
3. Con una risoluzione del raster di 0,5 si ottengono valori quattro volte più
grandi (!!!) rispetto al valore reale (comunque, per grandi regioni è difficile
lavorare con questa risoluzione e anche con quella di 1,0).
4. Anche per la lunghezza totale dei fiumi (.shp di tipo linea, rasterizzato,
risol. 1,0) per ogni comune e per ogni maglia quadrata di 1 kmq si ottengono
valori sbagliati (i più vicini alla realtà con risol. di 1,0): 12 m invece di
15 m ecc.
5. Ho provato a rasterizzare anche lo strato dei comuni/e della griglia alla
risoluzione del raster e trasformarla di nuovo in vettore ma niente è cambiato.
6. Osservo che sia ZS che v.rast.stats e ZonalStatscalcolano solo valori INTERI
e non decimali (come r.report) per l'area. Ho provato con un DEM (risol. 10,
risoluzione che per il calcolo delle area risulta completamente sbagliata) e
per l'altitudine massima, minima e la differenza tra loro (range) può calcolare
e inserire nella tabella i valori decimali corretti (per la carta dell'energia
del rilievo).
Ho fatto centinaia di provi e non capisco dove sbaglio...qui casca l'asino.
Mi potete aiutare? Qualcuno conosce come devo fare? Mi sembra una cosa così
semplice per un programma GIS...
Grazie,
Gabriela
________________________________
Da: Paolo Cavallini <[email protected]>
A: [email protected]
Inviato: Giovedì 26 Aprile 2012 13:25
Oggetto: Re: [Gfoss] calcolo e rappresentazione del percentuale del bosco_QGIS
o GRASS
Il 20/03/2012 16:15, Paolo Cavallini ha scritto:
>
> i fiumi, i valori della loro lunghezza per ogni kmq sono completamente
> sbagliati.
> Per esempio, per una lunghezza reale di 492,27 m è stato calcolato 97 (?), per
> 257,106 m - 49. Ho calcolato i valori reali dei fiumi usando linea di misura
> esistente in QGIS. Con zonalstatist e statistiche zonali si ottengono gli
> stessi
> valori. Dove sbaglio?
>> Puo' darsi che sia un bug (lo strumento e' nuovo): qualcuno puo' per favore
>> confermare, ed aprire un ticket, meglio se con un piccolissimo campione di
>> esempio?
Salve.
Se ci mandate un piccolissimo esempio di dati di prova, con i risultati
giusti e quelli sbagliati, verifichiamo ed eventualmente correggiamo.
Saluti, e grazie.
--
Paolo Cavallini
See: http://www.faunalia.it/pc
_______________________________________________
[email protected]
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 rispecchiano necessariamente
le posizioni dell'AssociazioneGFOSS.it.
584 iscritti al 7.4.2012
_______________________________________________
[email protected]
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 rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
584 iscritti al 7.4.2012
