<div dir="ltr">Dear all, <div>after several tests and a bit of inverse-engineering, Longzhu and I have understood the reason of such as a small difference in the delineation of the polygons basins.</div><div><br></div><div>All the process beyond r.watershed are computed in 3x3 moving window (e.g. slope, flow accumulation, etc.), therefore if we change the extent of study area we change the starting point of the mowing window and of course the values of the flow accumulation change a bit. These changes in the flow accumulation cause variation on the delineations of the basin borders, of adjacent tiles. </div><div><br></div><div>The reason, it is a bit obvious but did not come into my mind, immediately. </div><div><br></div><div>Now I'm trying to figure out a clever way to merge these flow accumulations that come from different tiles and that are slightly different. </div><div><br></div><div>Ciao <br>Giuseppe </div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div> </div><div><br></div><div><br></div><div><br></div><div> </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 5 Jun 2020 at 02:21, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Giuseppe,<br>
<br>
On Wed, May 13, 2020 at 2:16 AM Giuseppe Amatulli<br>
<<a href="mailto:giuseppe.amatulli@gmail.com" target="_blank">giuseppe.amatulli@gmail.com</a>> wrote:<br>
><br>
> Hi Markus,<br>
><br>
> I use GRASS 7.6.0<br>
> under a HPC running<br>
> Distributor ID: RedHatEnterpriseServer<br>
> Description: Red Hat Enterprise Linux Server release 7.7 (Maipo)<br>
> Release: 7.7<br>
> Codename: Maipo<br>
><br>
> The compression is:<br>
><br>
> GRASS 7.6.0 (nc_spm_08_grass7):~ > r.compress -p map=flow<br>
> <flow> is compressed (method 5: ZSTD). Data type: DCELL<br>
> <flow> has a compressed NULL file<br>
> [Raster MASK present].<br>
><br>
> My option in selecting r.watershed with out the -m is manly due to be able to use the -b option and also be able to run the tile-computation in 6 hours vs several days (probably weeks) of the full globe (or continents) run in segmentation mode - made quit hard to do debugging and re-running.<br>
><br>
> Do you have any suggestions for the mismatching of the borders?<br>
<br>
I have no idea but it is recommended to update to GRASS GIS 7.8.3.<br>
<br>
Best,<br>
markusN<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Giuseppe Amatulli, Ph.D.<br><br>Research scientist at<br><span>School of Forestry & Environmental Studies<br></span><span>Center for Research Computing</span><br>Yale University</div><div>New Haven, CT, USA</div><div>06511<br><div>
Teaching: <a href="http://spatial-ecology.net" target="_blank">http://spatial-ecology.net</a></div>
Work: <a href="https://environment.yale.edu/profile/giuseppe-amatulli/" target="_blank">https://environment.yale.edu/profile/giuseppe-amatulli/</a> <br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>