|

»

»

Curva ipsografica

View: New views

4 Messages — Rating Filter: Alert me

	Curva ipsografica by Patti Giuseppe :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message 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@... 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.
	Re: Curva ipsografica by ivan marchesini :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message 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@... > 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@... ivan.marchesini@... tel: +39(0)755853760 fax (university): +39(0)755853756 fax (home): +39(0)5782830887 jabber: geoivan73@... _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... 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.
	Re: Curva ipsografica by Patti Giuseppe :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message 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@... >> 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. >> >> _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... 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.
	Re: Curva ipsografica by Davide-34 :: Rate this Message: Reply to Author \| View Threaded \| Show Only this Message 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@...>: > 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@... > 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@... WMS http://geoserver1.unipv.it/cgi-bin/wms? _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione Gfoss@... 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.