[Gfoss] unione DEM quota decimale e conseguente problema con r.statistic

salve, ho un problema dovuto al fatto che ho un dem dove le quote non
sono numeri interi ma decimali...

partiamo dall'antefatto:
dal sito
della nasa ho scaricato i dem, con qgis li ho uniti con il comando
"raster - > unisci" ottenendo un dem "grande"
poi con grass l'ho
importato nel mio sistema di coordinate (utm wgs84 fuso 32N)
ho cos'
ottenuto un dem dove le quote non sono espresse come numeri interi, ma
come decimali...

fino a qui non mi è mai sembrato un problema, il
problema si è proposto quando ho creato un'area vettoriale che mi
rappresenta un bacino idrologico.
vorrei porre fra gli attributi di
tale vettoriale la quota media.

per prima cosa ho trasformato il
vettoriale in raster con il comando "v.to.rast"
poi con il comando "r.
mapcalc" ho "ritagliato" il dem sulla base del mio bacino, creando - in
sostanza - un raster della stessa area del bacino idrografico, che ha
per ciascun pxl la quota come da dem (non so se mi sono spiegato).

questo raster l'ho chiamato "elevation"

ho usato il comando "r.
statistic" per creare un raster con la quota media.
ho impostato come
"base raster map" il raster del bacino
e come "cover map raster" la
mappa "elevation" da me creata con "r.mapcalc"
come metodo ho settato
"average" ed ho settato un nome all'output

facendo così sulla shell mi
compare il seguente msg:
ERRORE: This module currently only works for
integer (CELL) maps

riprendo un altro dem (creato tempo fa con esri
che però non copre tutta la mia area di studio con tutti i bacini che
ho creato) e vedo che le quote sono numeri interi.
ritaglio con "r.
mapcalc" e "r.statistic" funziona correttamente

come posso risolvere
il problema??
grazie mille!!
marco

Risparmia con Tutto Incluso Light: telefono + adsl 8 mega a soli 19,95 € al mese per un anno!

SCONTO DI 120 EURO!

http://abbonati.tiscali.it/telefono-adsl/prodotti/tc/tuttoincluso_light/?WT.mc_id=01fw

Ciao Marco,
potresti esportare il dem con r.out.ascii con l'opzione (di default) che il numero di cifre significative dp sia = 0;
dopo di che, lo riimporti con r.in.ascii.

Ciao,
     Bianca

Il 27/09/2010 19.01, marco.donnini@tiscali.it ha scritto:

salve, ho un problema dovuto al fatto che ho un dem dove le quote non
sono numeri interi ma decimali...

partiamo dall'antefatto:
dal sito
della nasa ho scaricato i dem, con qgis li ho uniti con il comando
"raster -> unisci" ottenendo un dem "grande"
poi con grass l'ho
importato nel mio sistema di coordinate (utm wgs84 fuso 32N)
ho cos'
ottenuto un dem dove le quote non sono espresse come numeri interi, ma
come decimali...

fino a qui non mi è mai sembrato un problema, il
problema si è proposto quando ho creato un'area vettoriale che mi
rappresenta un bacino idrologico.
vorrei porre fra gli attributi di
tale vettoriale la quota media.

per prima cosa ho trasformato il
vettoriale in raster con il comando "v.to.rast"
poi con il comando "r.
mapcalc" ho "ritagliato" il dem sulla base del mio bacino, creando - in
sostanza - un raster della stessa area del bacino idrografico, che ha
per ciascun pxl la quota come da dem (non so se mi sono spiegato).

questo raster l'ho chiamato "elevation"

ho usato il comando "r.
statistic" per creare un raster con la quota media.
ho impostato come
"base raster map" il raster del bacino
e come "cover map raster" la
mappa "elevation" da me creata con "r.mapcalc"
come metodo ho settato
"average" ed ho settato un nome all'output

facendo così sulla shell mi
compare il seguente msg:
ERRORE: This module currently only works for
integer (CELL) maps

riprendo un altro dem (creato tempo fa con esri
che però non copre tutta la mia area di studio con tutti i bacini che
ho creato) e vedo che le quote sono numeri interi.
ritaglio con "r.
mapcalc" e "r.statistic" funziona correttamente

come posso risolvere
il problema??
grazie mille!!
marco

Risparmia con Tutto Incluso Light: telefono + adsl 8 mega a soli 19,95 € al mese per un anno!

SCONTO DI 120 EURO!

http://abbonati.tiscali.it/telefono-adsl/prodotti/tc/tuttoincluso_light/?WT.mc_id=01fw
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
Gfoss@lists.gfoss.it
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 rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
474 iscritti al 18.9.2010

--
Bianca Federici
DICAT - Dipartimento di Ingegneria delle Costruzioni, dell'Ambiente e del Territorio
Laboratorio di Geodesia, Geomatica e GIS
Università degli Studi di Genova

- Via Montallegro, 1 - 16145 Genova
   Tel: +39 010 353 2421 - Fax: +39 010 353 2527
- Polo accademico di Savona - Palazzina Delfino
   Via A. Magliotto, 2 - 17100 Savona
   Tel: +39 019 219 45381 - Fax: +39 019 219 45382
- E-mail: bianca.federici@unige.it - www.dicat.unige.it

---------------------------------------------------------------------------------
L'unica razza che conosco è quella umana - A. Einstein
Inizio e fine della ricerca è la conoscenza - A.M. Manzino
---------------------------------------------------------------------------------

Please consider the environment before printing this email!