[Gfoss] Idrogramma di piena

Buon giorno a tutti,
vorrei calcolare l'idrogramma di piena di un torrente con un bacino di
piccole dimensioni (qualche decina di km2) provando a ipotizzare
differenti usi del suolo per valutare la risposta ai cambiamenti.

Ho provato a usare r.topmodel ottenendo senza fatica dei risultati a
grandi linee sovrapponibili alla registrazione dell'unica stazione
idrometrica. Ma topmodel non gestisce differenti raster per l'uso del
suolo.

Esiste qualche altro metodo per ottenere un idrogramma?

Ciao e grazie

Simone

2008/11/24 Simone Pittaluga <pitta@ibidem.it>:

Buon giorno a tutti,
vorrei calcolare l'idrogramma di piena di un torrente con un bacino di
piccole dimensioni (qualche decina di km2) provando a ipotizzare
differenti usi del suolo per valutare la risposta ai cambiamenti.

Ho provato a usare r.topmodel ottenendo senza fatica dei risultati a
grandi linee sovrapponibili alla registrazione dell'unica stazione
idrometrica. Ma topmodel non gestisce differenti raster per l'uso del
suolo.

Esiste qualche altro metodo per ottenere un idrogramma?

Sembra con JGRASS:
http://geomatica.como.polimi.it/workbooks/n5/abstract/idro_jgrass_abs_it.pdf

Markus

Ciao a tutti,

Qualcuno ha mai provato ad applicare una griglia di trasformazione di datum per l'Italia in formato NTv2 all'interno di Mapserver??

Ho fatto un programmino che trasforma i file gr2 in gsb (NTv2 binario). Il file *.gsb funziona con il programma GDaiT e con cs2cs (PROJ4) ma con mapserver mi fa sparire il layer.

Praticamente ho impostato il MAPFILE con la PROJECTION epsg:32632 ... (UTM-WGS84) ... un LAYER è però nel sistema Roma Monte-Mario (quindi PROJECTION epsg:3003). Nel file EPSG ho aggiunto alla stringa di parametri il parametro "+nadgrids=__142.gsb" e ho ovviamente messo quel file gsb nella cartella dove vengono letti i grigliati. Il file viene letto sicuramente, solo che succedono le seguenti cose: non avviene NESSUNA trasformazione (neanche quella standard) e gli elementi non si vedono più.

Grazie.
--------------------------------------
Francesco Pirotti
Dept. Te.S.A.F./CIRGEO
University of Padua
Viale dell'Università 16
I-35020 Legnaro (PD)
tel. +39 049 827 2710
fax. +39 049 827 2686
mob. +39 349 55 39 361
@mail francesco.pirotti@unipd.it
skypeID francesco197576
http://www.cirgeo.unipd.it/cirgeo/francescopirotti.htm

On Mon, Nov 24, 2008 at 5:30 PM, francesco <francesco.pirotti@unipd.it> wrote:

Ciao a tutti,

Qualcuno ha mai provato ad applicare una griglia di trasformazione di datum per l'Italia in formato NTv2 all'interno di Mapserver??

Ho fatto un programmino che trasforma i file gr2 in gsb (NTv2 binario). Il file *.gsb funziona con il programma GDaiT e con cs2cs (PROJ4) ma con mapserver mi fa sparire il layer.

Praticamente ho impostato il MAPFILE con la PROJECTION epsg:32632 ... (UTM-WGS84) ... un LAYER è però nel sistema Roma Monte-Mario (quindi PROJECTION epsg:3003). Nel file EPSG ho aggiunto alla stringa di parametri il parametro "+nadgrids=__142.gsb" e ho ovviamente messo quel file gsb nella cartella dove vengono letti i grigliati. Il file viene letto sicuramente, solo che succedono le seguenti cose: non avviene NESSUNA trasformazione (neanche quella standard) e gli elementi non si vedono più.

Sono diventato scemo con un problema simile a quello che citi. Avevo
chiesto se qualcuno in lista capisce perché ma senza risposta (e, per
mancanza di tempo, non ho approfondito). Secondo me c'è qualcosa che
non va in qualche parte di gdal-ogr-proj4 oppure ho capito male io
(ipotesi da tenere sempre presente). Se specifichi la trasformazione
utilizzando per il mapdatum roma40 l'EPSG che trovi in wikipedia *non*
succede niente. Ma se fai:

    gb = osgeo.osr.SpatialReference()
    gb.ImportFromProj4("+proj=tmerc +lat_0=0 +lon_0=15 +k=0.999600
+x_0=2520000 +y_0=0 +ellps=intl +units=m
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68")
    wgs84 = osgeo.osr.SpatialReference()
    wgs84.ImportFromEPSG(4326)
    tr = osgeo.osr.CoordinateTransformation(gb, wgs84)

allora va. Ovviamente dovresti dare la definizione in formato PROJ.4 a
mapserver per il datum del layer ma mi pare che la radice del problema
sia simile (pure a me accadeva la stessa cosa che citi).

--
Christian Pellegrin, see http://www.evolware.org/chri/
"Real Programmers don't play tennis, or any other sport which requires
you to change clothes. Mountain climbing is OK, and Real Programmers
wear their climbing boots to work in case a mountain should suddenly
spring up in the middle of the computer room."

francesco ha scritto:

Ciao a tutti,

Qualcuno ha mai provato ad applicare una griglia di trasformazione di datum per l'Italia in formato NTv2 all'interno di Mapserver??

all'interno di mapserver ancora no... se posso, preferisco trasformare a
monte i dati

Ho fatto un programmino che trasforma i file gr2 in gsb (NTv2 binario). Il file *.gsb funziona con il programma GDaiT e con cs2cs (PROJ4) ma con mapserver mi fa sparire il layer.

anche Bud Bruegger aveva fatto un python script in tal senso:
http://n2.nabble.com/generazione-di-grigliati-NTV2-personalizzati-td759181.html

Praticamente ho impostato il MAPFILE con la PROJECTION epsg:32632 ... (UTM-WGS84) ... un LAYER è però nel sistema Roma Monte-Mario (quindi PROJECTION epsg:3003). Nel file EPSG ho aggiunto alla stringa di parametri il parametro "+nadgrids=__142.gsb" e ho ovviamente messo quel file gsb nella cartella dove vengono letti i grigliati. Il file viene letto sicuramente, solo che succedono le seguenti cose: non avviene NESSUNA trasformazione (neanche quella standard) e gli elementi non si vedono più.

Sembrerebbe che le PROJ_LIB (proj\nad) siano ignorate. Prova a
verificare se sono nella path. Hai provato con il cs2cs presente nelle
proj di mapserver? Al posto tuo trasformerei prima i dati...
risparmieresti un ulteriore carico di lavoro al tuo mapserver, anche se
a quanto pare il formato NTv2 non è particolarmente oneroso a livello
computazionale.

ciao
Antonio

Riguardo il mio post vi tengo aggiornati, visto che sembra essere di interesse per tutti (e ci credo!!) .

Sembra che il problema sia che se si specifica un nadgrid allora lo utilizza ma proj ritorna un errore se i punti non ricadono dentro la finestra del file grigliato (come specificato nel FAQ). Bisogna quindi aggiungere un grid null , quindi nella stringa 3003 tolgo i parametri +towgs e inserisco +nadgrids=<ilmiofile.gsb>,null ... in questo caso i miei punti vengono disegnati, anche se non trasformati... appena ho i file dei grigliati giusti provo a fare un super file NTv2 con tanti subgrids quanti sono i grigliati e vedo se lo utilizza... vi farò sapere se con successo o no.

--------------------------------------
Francesco Pirotti
Dept. Te.S.A.F./CIRGEO
University of Padua
Viale dell'Università 16
I-35020 Legnaro (PD)
tel. +39 049 827 2710
fax. +39 049 827 2686
mob. +39 349 55 39 361
@mail francesco.pirotti@unipd.it
skypeID francesco197576
http://www.cirgeo.unipd.it/cirgeo/francescopirotti.htm

---------- Original Message ----------
To: (gfoss@faunalia.com)
From: francesco (francesco.pirotti@unipd.it)
Subject: [Gfoss] trasformazione tra datum WGS84 e Roma40
Date: 24/11/2008 17:30:40

Ciao a tutti,

Qualcuno ha mai provato ad applicare una griglia di trasformazione di datum per l'Italia in formato NTv2 all'interno di Mapserver??

Ho fatto un programmino che trasforma i file gr2 in gsb (NTv2 binario). Il file *.gsb funziona con il programma GDaiT e con cs2cs (PROJ4) ma con mapserver mi fa sparire il layer.

Praticamente ho impostato il MAPFILE con la PROJECTION epsg:32632 ... (UTM-WGS84) ... un LAYER è però nel sistema Roma Monte-Mario (quindi PROJECTION epsg:3003). Nel file EPSG ho aggiunto alla stringa di parametri il parametro "+nadgrids=__142.gsb" e ho ovviamente messo quel file gsb nella cartella dove vengono letti i grigliati. Il file viene letto sicuramente, solo che succedono le seguenti cose: non avviene NESSUNA trasformazione (neanche quella standard) e gli elementi non si vedono più.

Grazie.
--------------------------------------
Francesco Pirotti
Dept. Te.S.A.F./CIRGEO
University of Padua
Viale dell'Università 16
I-35020 Legnaro (PD)
tel. +39 049 827 2710
fax. +39 049 827 2686
mob. +39 349 55 39 361
@mail francesco.pirotti@unipd.it
skypeID francesco197576
http://www.cirgeo.unipd.it/cirgeo/francescopirotti.htm

Ciao Simone (grazie Markus),
in JGrass c'e' un modello, Peakflow, che permette il calcolo
dell'idrogramma di piena.

vorrei calcolare l'idrogramma di piena di un torrente con un bacino di
piccole dimensioni (qualche decina di km2) provando a ipotizzare
differenti usi del suolo per valutare la risposta ai cambiamenti.

Peakflow ha un approcio diverso, basato sulla saturazione del terreno,
trovi una presentazione aggiornata qui:
http://www.slideshare.net/silli/corso-pat-jgrasspeakflow-2008-10-23-presentation?type=powerpoint

Ciao
Andrea

Ho provato a usare r.topmodel ottenendo senza fatica dei risultati a
grandi linee sovrapponibili alla registrazione dell'unica stazione
idrometrica. Ma topmodel non gestisce differenti raster per l'uso del
suolo.

Esiste qualche altro metodo per ottenere un idrogramma?

Ciao e grazie

Simone

_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
Gfoss@faunalia.com
http://www.faunalia.com/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.

Grazie per i suggerimenti.
La presentazione è molto utile proverò a seguire le istruzioni.

Volevo scusarmi e ringraziarlo ancora con Markus perchè ieri gli ho
risposto in forma privata (per errore).

Ciao
Simone

------- Messaggio originale -------
Da: andrea antonello <andrea.antonello@gmail.com>
CC: gfoss <gfoss@faunalia.com>
Oggetto: Re: [Gfoss] Idrogramma di piena
Data: Tue, 25 Nov 2008 09:05:31 +0100

Peakflow ha un approcio diverso, basato sulla saturazione del terreno,
trovi una presentazione aggiornata qui:
http://www.slideshare.net/silli/corso-pat-jgrasspeakflow-2008-10-23-presentation?type=powerpoint

Ciao
Andrea

_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
Gfoss@faunalia.com
http://www.faunalia.com/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.