[Gfoss] calcolo e rappresentazione del percentuale del bosco_QGIS o GRASS

Lun 14 Maggio 2012 18:04:04 CEST

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 <cavallini a faunalia.it>
A: gfoss a lists.gfoss.it 
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

_______________________________________________
Gfoss a 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 rispecchiano necessariamente
le posizioni dell'AssociazioneGFOSS.it.
584 iscritti al 7.4.2012
-------------- parte successiva --------------
Un allegato HTML ? stato rimosso...
URL: <http://lists.gfoss.it/pipermail/gfoss/attachments/20120514/1ae4c7c8/attachment.html>