[Gfoss] Curva ipsografica
Davide
riverlab a gmail.com
Lun 21 Lug 2008 16:07:17 CEST
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 a 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 a 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 a manhattan.unipv.it
WMS http://geoserver1.unipv.it/cgi-bin/wms?
Maggiori informazioni sulla lista
Gfoss