Buongiorno a tutti. Sono alle prese con la determinazione di alcuni parametri di bacino e tra questi vorrei ricavare in automatico la curva ipsografica con grass. Posto che ho già un dem con cella di 20 m tagliato sul bacino (che è di circa 500 kmq), al momento l'unico modo che mi è venuto in mente per farlo è di interrogare il raster per le quote comprese in step di altitudine prefissati, ad esempio ogni 500 m noti i minimi e i massimi, e interrogare il risultato per conoscere l'area di ogni step, quindi rapportato all'area totale del bacino plottare le percentuali contro le altitudini e disegnare la curva unendo i puntini. Ovviamente non mi sembra una soluzione elegante ma credo porti al risultato. Accetto volentieri altri suggerimenti (magari anche sul calcolo automatico dell'altitudine media), ho già provato anche a usare jgrass-udig ma mi trovo in serio imbarazzo rispetto a grass, per cui preferirei quest'ultimo.
Grazie in anticipo!
Ciao ..
sul libro di Markus c'è un esempio specifico su questo...
usa la funzione ecdf di R ...
ciao
Il giorno lun, 21/07/2008 alle 09.00 +0200, Patti Giuseppe ha scritto:
Buongiorno a tutti. Sono alle prese con la determinazione di alcuni
parametri di bacino e tra questi vorrei ricavare in automatico la curva
ipsografica con grass. Posto che ho già un dem con cella di 20 m
tagliato sul bacino (che è di circa 500 kmq), al momento l'unico modo
che mi è venuto in mente per farlo è di interrogare il raster per le
quote comprese in step di altitudine prefissati, ad esempio ogni 500 m
noti i minimi e i massimi, e interrogare il risultato per conoscere
l'area di ogni step, quindi rapportato all'area totale del bacino
plottare le percentuali contro le altitudini e disegnare la curva unendo
i puntini. Ovviamente non mi sembra una soluzione elegante ma credo
porti al risultato. Accetto volentieri altri suggerimenti (magari anche
sul calcolo automatico dell'altitudine media), ho già provato anche a
usare jgrass-udig ma mi trovo in serio imbarazzo rispetto a grass, per
cui preferirei quest'ultimo.
Grazie in anticipo!_______________________________________________
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.
--
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format
Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a
06125
Perugia (Italy)
Socio fondatore GFOSS "Geospatial Free and Open Source Software" http://www.gfoss.it
e-mail: marchesini@unipg.it
ivan.marchesini@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoivan73@jabber.org
Vero! A pag. 349 della prima edizione! Ce l'avevo sotto il naso....scusate e grazie (adesso mi tocca masticare R)!
ivan marchesini ha scritto:
Ciao ..
sul libro di Markus c'è un esempio specifico su questo...
usa la funzione ecdf di R ...ciao
Il giorno lun, 21/07/2008 alle 09.00 +0200, Patti Giuseppe ha scritto:
Buongiorno a tutti. Sono alle prese con la determinazione di alcuni parametri di bacino e tra questi vorrei ricavare in automatico la curva ipsografica con grass. Posto che ho già un dem con cella di 20 m tagliato sul bacino (che è di circa 500 kmq), al momento l'unico modo che mi è venuto in mente per farlo è di interrogare il raster per le quote comprese in step di altitudine prefissati, ad esempio ogni 500 m noti i minimi e i massimi, e interrogare il risultato per conoscere l'area di ogni step, quindi rapportato all'area totale del bacino plottare le percentuali contro le altitudini e disegnare la curva unendo i puntini. Ovviamente non mi sembra una soluzione elegante ma credo porti al risultato. Accetto volentieri altri suggerimenti (magari anche sul calcolo automatico dell'altitudine media), ho già provato anche a usare jgrass-udig ma mi trovo in serio imbarazzo rispetto a grass, per cui preferirei quest'ultimo.
Grazie in anticipo!_______________________________________________
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.
Ciao Giuseppe, io avevo preso questo da un pub. di Alberto Viglione
#########################
# Curva ipsografica #
#########################
# File da dar in pasto a R
echo "heigths <- scan(\"$1\")
ipsographic <- sort(heigths, decreasing = TRUE)
n <- length(ipsographic)
percent <- c(1:n)
percent <- percent/n
intervalli <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50,
60, 70, 80, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100)
intervalli <- intervalli/100
indici <- round(intervalli*n)
quote <- ipsographic[indici]
cat(quote, file=\"file_temp_ipso3.txt\")" > file_temp_ipso2.txt
R -q --slave --vanilla < file_temp_ipso2.txt
a1percento=`awk '{print $1}' file_temp_ipso3.txt`
a2percento=`awk '{print $2}' file_temp_ipso3.txt`
a3percento=`awk '{print $3}' file_temp_ipso3.txt`
a4percento=`awk '{print $4}' file_temp_ipso3.txt`
a5percento=`awk '{print $5}' file_temp_ipso3.txt`
a6percento=`awk '{print $6}' file_temp_ipso3.txt`
a7percento=`awk '{print $7}' file_temp_ipso3.txt`
a8percento=`awk '{print $8}' file_temp_ipso3.txt`
a9percento=`awk '{print $9}' file_temp_ipso3.txt`
a10percento=`awk '{print $10}' file_temp_ipso3.txt`
a20percento=`awk '{print $11}' file_temp_ipso3.txt`
a30percento=`awk '{print $12}' file_temp_ipso3.txt`
a40percento=`awk '{print $13}' file_temp_ipso3.txt`
a50percento=`awk '{print $14}' file_temp_ipso3.txt`
a60percento=`awk '{print $15}' file_temp_ipso3.txt`
a70percento=`awk '{print $16}' file_temp_ipso3.txt`
a80percento=`awk '{print $17}' file_temp_ipso3.txt`
a90percento=`awk '{print $18}' file_temp_ipso3.txt`
a91percento=`awk '{print $19}' file_temp_ipso3.txt`
a92percento=`awk '{print $20}' file_temp_ipso3.txt`
a93percento=`awk '{print $21}' file_temp_ipso3.txt`
a94percento=`awk '{print $22}' file_temp_ipso3.txt`
a95percento=`awk '{print $23}' file_temp_ipso3.txt`
a96percento=`awk '{print $24}' file_temp_ipso3.txt`
a97percento=`awk '{print $25}' file_temp_ipso3.txt`
a98percento=`awk '{print $26}' file_temp_ipso3.txt`
a99percento=`awk '{print $27}' file_temp_ipso3.txt`
a100percento=`awk '{print $28}' file_temp_ipso3.txt`
cp file_temp_ipso3.txt $1_ipso.txt
Ciao
D.
2008/7/21 Patti Giuseppe <geognu@infinito.it>:
Buongiorno a tutti. Sono alle prese con la determinazione di alcuni
parametri di bacino e tra questi vorrei ricavare in automatico la curva
ipsografica con grass. Posto che ho già un dem con cella di 20 m tagliato
sul bacino (che è di circa 500 kmq), al momento l'unico modo che mi è venuto
in mente per farlo è di interrogare il raster per le quote comprese in step
di altitudine prefissati, ad esempio ogni 500 m noti i minimi e i massimi, e
interrogare il risultato per conoscere l'area di ogni step, quindi
rapportato all'area totale del bacino plottare le percentuali contro le
altitudini e disegnare la curva unendo i puntini. Ovviamente non mi sembra
una soluzione elegante ma credo porti al risultato. Accetto volentieri altri
suggerimenti (magari anche sul calcolo automatico dell'altitudine media), ho
già provato anche a usare jgrass-udig ma mi trovo in serio imbarazzo
rispetto a grass, per cui preferirei quest'ultimo.
Grazie in anticipo!_______________________________________________
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.
--
dott. Davide Zizioli
Dipartimento di Scienze della Terra
Università di Pavia
Via Ferrata 1, 27100 Pavia
Tel. +393498345440
Tel. Skype: Riverlab
davide.zizioli@manhattan.unipv.it
WMS http://geoserver1.unipv.it/cgi-bin/wms?