[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