<div dir="ltr">Hi Micha <div><br></div><div>    I understand what you mean. But it requires another step to manually identify depressions from these pre-conditioned DEM. </div><div><br></div><div>Cheers</div><div>Ming </div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> 于2021年4月3日周六 下午12:43写道：<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
On 4/2/21 5:37 PM, ming han wrote:<br>
> Maybe I am the only one who has this demand. Following is just a <br>
> recommendation to GRASS r.watershed function.<br>
> Maybe it is worth having an option to avoid r.watershed overcome <br>
> depressions.<br>
> The reasons are 1) there are many hydrologically pre condition DEM <br>
> data available globally, such as:HydroSHEDS, MERIT<br>
>                             2) the depression in these DEM are real <br>
> depressions, overcome these depressions will make the entire drainage <br>
> system<br>
<br>
<br>
Regarding Hydrosheds, the documentation[1] in section 3.4 explains how <br>
they overcame the problem of sinks. They performed a regular "fill <br>
sinks" operation on areas that were SRTM artifacts. True natural <br>
depressions were identified manually, then another manual procedure of <br>
carving rivers was done to force flow thru these depressions and produce <br>
hydrologically correct streams and basins. So pre-conditioning to <br>
overcome depressions is not a magic bullet...<br>
<br>
<br>
In my opinion, the best results are obtained when true depressions <br>
(pits, salt playas or karst regions) are identified, and set to NULL in <br>
the elevation raster. That will allow r.watershed to stop routing at <br>
those locations, and produce correct stream and basin layers.<br>
<br>
<br>
[1]<a href="https://hydrosheds.org/images/inpages/HydroSHEDS_TechDoc_v1_2.pdf" rel="noreferrer" target="_blank">https://hydrosheds.org/images/inpages/HydroSHEDS_TechDoc_v1_2.pdf</a><br>
<br>
<br>
> incorrectly.<br>
><br>
> I understand GRASS has other functions to solve this problem, but just <br>
> a user recommendation. I use GRASS a lot.<br>
><br>
> Thanks<br>
> Ming<br>
><br>
> ming han <<a href="mailto:dustming@gmail.com" target="_blank">dustming@gmail.com</a> <mailto:<a href="mailto:dustming@gmail.com" target="_blank">dustming@gmail.com</a>>> <br>
> 于2021年3月30日周二 上午8:06写道：<br>
><br>
>     Got it, thanks everyone~<br>
>     Ming<br>
><br>
>     Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">tsvibar@gmail.com</a> <mailto:<a href="mailto:tsvibar@gmail.com" target="_blank">tsvibar@gmail.com</a>>><br>
>     于2021年3月29日周一 下午2:40写道：<br>
><br>
>         Hello:<br>
><br>
>         You might try `r.param.scale`, or even better `r.geomorphons`<br>
>         modules to<br>
>         identify geomorphology features, then filter out all pixels<br>
>         identified<br>
>         as pits.<br>
><br>
><br>
>         r.watershed is purposely designed to overcome depressions, and<br>
>         find flow<br>
>         routing thru these spots. So I don't think you can use that<br>
>         module to<br>
>         identify depressions.<br>
><br>
><br>
>         On 3/27/21 8:49 PM, ming han wrote:<br>
>         > Hi  Everyone<br>
>         ><br>
>         >      When I do watershed delineation using r.watershed for<br>
>         great salt<br>
>         > lake watershed. I found r.watershed always tried to assign<br>
>         an outlet<br>
>         > for a great salt lake, which does actually not exist because<br>
>         it is an<br>
>         > inland lake and the great salt lake has no watershed outlet<br>
>         at all.<br>
>         ><br>
>         >       I noticed that there is a depression option. But is<br>
>         there any<br>
>         > way that  r.watershed can automatically identify depressions<br>
>         while<br>
>         > defining flow accumulation and stream network?<br>
>         ><br>
>         > Thanks<br>
>         > Ming<br>
>         ><br>
>         > _______________________________________________<br>
>         > grass-user mailing list<br>
>         > <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a> <mailto:<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>><br>
>         > <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
>         <<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a>><br>
><br>
>         -- <br>
>         Micha Silver<br>
>         Ben Gurion Univ.<br>
>         Sde Boker, Remote Sensing Lab<br>
>         cell: +972-523-665918<br>
><br>
-- <br>
Micha Silver<br>
Ben Gurion Univ.<br>
Sde Boker, Remote Sensing Lab<br>
cell: +972-523-665918<br>
<br>
</blockquote></div>