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/
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Il 24/10/2013 11:06, marco.donnini@tiscali.it 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-----
2013/10/24 <marco.donnini@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