<div dir="ltr"><div>Don't have time right now, but running it with valgrind might show the problem.</div><div><br></div><div>Anna</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 9, 2020 at 8:27 AM Vincent Bain <<a href="mailto:bain@toraval.fr">bain@toraval.fr</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">Hello Vincent,<br>
<br>
in case it could be related to my recent issue with r.param.scale:<br>
<br>
<a href="https://lists.osgeo.org/pipermail/grass-user/2020-October/081788.html" rel="noreferrer" target="_blank">https://lists.osgeo.org/pipermail/grass-user/2020-October/081788.html</a><br>
<br>
In my case there was most probably something wrong with my OpenMP<br>
configuration.<br>
<br>
HTH,<br>
Vincent.<br>
<br>
Le mercredi 09 décembre 2020 à 13:59 +0100, Vincent Godard a écrit :<br>
> Dear Grass users,<br>
> <br>
> I've been using r.param.scale to compute various type of curvatures<br>
> over <br>
> a 1 m LiDAR DEM (FCELL) (from <a href="http://opentopography.org" rel="noreferrer" target="_blank">opentopography.org</a>) :<br>
> <br>
> r.param.scale --o -c input=dem output=curvature size=15 method=profc<br>
> <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r curvature<br>
> <br>
> I usually get reasonable results like (checked with independent <br>
> calculations) :<br>
> <br>
> min=-0.203124601340989<br>
> max=0.0685868278408049<br>
> <br>
> But when running this command repeatedly, I randomly (usually every<br>
> 3rd <br>
> or 4th time) get very different results like :<br>
> <br>
> min=-127112309.998249<br>
> max=131128979.492382<br>
> <br>
> I am running GRASS 7.8.4 on Ubuntu 20.04. The problem is also<br>
> present <br>
> when running r.param.scale from QGIS or rgrass7 (in a R session).<br>
> <br>
> I have replicated the problem with another high resolution DEM of <br>
> different source (still FCELL), and also with a 30 m SRTM DEM<br>
> (CELL), <br>
> but in this case the occurrence of these odd results is less<br>
> frequent <br>
> and they are not contrasting as much with the expected outcome<br>
> ("only" <br>
> by a 1 or 2 order of magnitude). The problem is present with or<br>
> without <br>
> the -c option, and seems only to affect the computation of<br>
> curvatures <br>
> (profc,planc,longc,crosc,minic,maxic).<br>
> <br>
> Has anyone encountered something like that, or have any idea about<br>
> what <br>
> might be going on?<br>
> <br>
> Thanks in advance,<br>
> <br>
> Vincent<br>
> <br>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<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>
</blockquote></div></div>