[Gfoss] verifica calcolo % da raster a vettoriale

[Luca Delucchi]

2013/10/24  <marco.donnini a tiscali.it>:
> Buongiorno a tutti/e
>

ciao

> Ho una mappa raster: è il Corine Land Cover 2006 (mappa dell'uso del suolo)
> da me semplificata in modo da avere solo due valori numerici (valore
> 100: uso "antropico", valore 200: uso "non antropico").
>
> Ho poi una mappa vettoriale con alcuni bacini idrografici.
>
> Io vorrei sapere per ciascun bacino la % di uso antropico e la % di uso non
> antropico.
>
> Io ho in mente un metodo (forse un po' macchinoso), ma non so se è il metodo
> più rapido...
>
> Per prima cosa trasformare il raster in vettoriale (r.to.vect).
>
> Per prima cosa bisogna fare l'overlay della mappa nei bacini.
> (v.overlay/and).
> Poi creare il database con “createdb nome_database”
> Connettersi a postgres → db.connect → selezionare postrgres
> (nome driver: pg, nome database: db creato)
> Fare la copia dell'overlay che si connette in automatico con postgres.
>
> Aggiungere al file “overlay” una colonna (double precision) dove calcolare
> l'area (v.db.addcol)
>
> v.to.db → aggiungere l'area nella nuova colonna area
>
> db.select sql="select campo_bacino,campo_usosuolo,sum(area) as area_tot from
> layer_overlay group by campo_bacino,campo_usosuolo order by campo_bacino" >
> stat_lito_bas
>
> info:
> stat_lito_bas = file di testo creato da comando
> con  “campo_bacino” e “campo_usosuolo” si intende i nomi delle tabelle
> attributi dell'overlay inerenti i nomi del bacino e del tipo di uso del
> suolo.
> con layer_overlay si intende il vettore “overlay” creato e connesso con pg.
>

non so se ho capito bene, ma forse v.rast.stats può fare al caso tuo

http://grass.osgeo.org/grass64/manuals/v.rast.stats.html

> Grazie mille!!
>
> marco
>


-- 
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org