<p dir="ltr">Con r.watershed devi creare la mappa di drainage e quella degli stream. Il punto di chiusura del bacino deve essere poi scelto sopra uno dei pixel degli stream, altrimenti r.water.outlet risale per poche o magari per nessuna cella ed il bacino risulta minuscolo o nullo. Forse il problema era quello?  Ciao  ivan</p>

<br><div class="gmail_quote"><div dir="ltr">On Fri, 24 Jul 2015 21:09 Ignazio Giuntoli <<a href="mailto:ixg282@bham.ac.uk" target="_blank">ixg282@bham.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Salve a tutti,<br>
<br>
sto cercando di delineare dei bacini idrografici a partire da un DEM su<br>
Grass 7<br>
utilizzando la funzione r.watershed e poi r.water.outlet con cui<br>
calcolare il bacino sotteso a un punto specifico (outlet di coordinate<br>
Lon, Lat).<br>
<br>
Sto seguendo le istruzioni trovate qui<br>
(<a href="http://www.slideshare.net/madi468/lez-3-corso-di-modelli-e-gis-per-lambiente?related=2" rel="noreferrer" target="_blank">http://www.slideshare.net/madi468/lez-3-corso-di-modelli-e-gis-per-lambiente?related=2</a><br>
- pagine 24 e 44).<br>
<br>
Questa operazione la devo fare su tanti punti (stazioni idrometriche)<br>
situati su vari DEM contigui.<br>
<br>
Partendo dal caso più semplice di un solo punto su cui delineare il<br>
bacino, il raster calcolato da Grass risulta vuoto.<br>
<br>
Mi chiedevo se qualcuno nel forum ha già effettuato operazioni simili.<br>
> Per info, i dati che sto usando sono:<br>
><br>
> a. il dem: srtm1605 (da<br>
> <a href="http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp" rel="noreferrer" target="_blank">http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp</a>)<br>
> b. il punto su cui calcolare il bacino: -100.2112531, 37.032528<br>
Fino a r.watershed tutto bene, ma, come accennato: r.water.outlet non<br>
crea nessun bacino.<br>
<br>
Usando un altro punto a caso più a ovest mi crea un micro bacino<br>
[r.water.outlet input=drenaggio output=bacino coordinates=-102.5,37.5 --o],<br>
ma niente alle coordinate del punto b che dovrebbe dare il bacino<br>
arancione (da shape-file del USGS) visibile in figura<br>
<a href="https://www.dropbox.com/s/hke3d9s4bpanrdg/bacini-01.png?dl=0" rel="noreferrer" target="_blank">https://www.dropbox.com/s/hke3d9s4bpanrdg/bacini-01.png?dl=0</a>.<br>
<br>
Grazie in anticipo per i vostri consigli.<br>
<br>
Ignazio<br>
<br>
<br>
ps. stavo provando a fare queste operazioni su QGIS 2.4.0 usando<br>
funzioni GRASS dal Processing Toolbox, ma pare ci sia un bug in<br>
r.water.outlet su QGIS per cui inserendo per esempio -100.2 per easting<br>
coord of outlet point, quando si va al campo successivo diventa -1002<br>
(cosa che non avviene su GRASS).<br>
<br>
_______________________________________________<br>
<a href="mailto:Gfoss@lists.gfoss.it" target="_blank">Gfoss@lists.gfoss.it</a><br>
<a href="http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss" rel="noreferrer" target="_blank">http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss</a><br>
Questa e' una lista di discussione pubblica aperta a tutti.<br>
I messaggi di questa lista non hanno relazione diretta con le posizioni dell'Associazione GFOSS.it.<br>
750 iscritti al 18.3.2015</blockquote></div>