verifica calcolo % da raster a vettoriale

Buongiorno a tutti/e

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.

Grazie mille!!

marco

Invita i tuoi amici e Tiscali ti premia! Il consiglio di un amico vale
più di uno spot in TV. Per ogni nuovo abbonato 30 € di premio per te e
per lui! Un amico al mese e parli e navighi sempre gratis
http://freelosophy.tiscali.it/


_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
666 iscritti al 22.7.2013

[pcav]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Il 24/10/2013 11:06, [hidden email] ha scritto:
> Buongiorno a tutti/e
>
> 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").

Il pugin di statistica zonale di QGIS, sul raster riclassificato?
Oppure il plugin esterno Lecos.
Saluti.
- --
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.14 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iEYEARECAAYFAlJo474ACgkQ/NedwLUzIr5sOACeMCepkWZ8kAqt0Y4oobgDE472
ZKkAn3mpodgPs8ouz/FD2EMsIcUbQvHH
=7i8B
-----END PGP SIGNATURE-----
_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
666 iscritti al 22.7.2013

[Luca Delucchi]

2013/10/24  <[hidden email]>:
> 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
_______________________________________________
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.
666 iscritti al 22.7.2013