
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>

</head>

<body dir="ltr">

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Hei,</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Before investing time in getting r.terraflow to work on larger regions I actually tried r.fill.dir.</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

The problem is that the results are significantly different and those of r.terraflow and r.hydrodem (where no breaching is performed) are more suitable for my needs. Here, r.terraflow fills most, r.fill.dir least and r.hydrodem inbetween.</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Here a little comparison based on NC data:<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

</div>

<span>g.extension extension=r.hydrodem operation=add<br>

</span>

<div>g.region -p raster=elevation<br>

</div>

<div><br>

</div>

<div>r.fill.dir --overwrite --verbose input=elevation output=elevation_fill_dir_depressionless direction=elevation_fill_dir_direction areas=elevation_fill_dir_pas<br>

</div>

<div><br>

</div>

<div>r.hydrodem -a --overwrite input=elevation memory=5000 output=elevation_hydrodem_depressionless<br>

</div>

<div><br>

</div>

<div>r.terraflow --overwrite --verbose elevation=elevation filled=elevation_terraflow_depressionless memory=5000<br>

</div>

<div><br>

</div>

<div>for m in terraflow fill_dir hydrodem<br>

</div>

<div>do<br>

</div>

<div>r.mapcalc --o expression="${m}_effect=if((elevation_${m}_depressionless-elevation)>0,elevation_${m}_depressionless-elevation,null())"<br>

</div>

<div>done<br>

</div>

<div><br>

</div>

<div>r.univar map="terraflow_effect"<br>

</div>

<div>total null and non-null cells: 2025000<br>

</div>

<div>total null cells: 1901290<br>

</div>

<div><br>

</div>

<div>Of the non-null cells:<br>

</div>

<div>----------------------<br>

</div>

<div>n: 123710<br>

</div>

<div>minimum: 7.62939e-06<br>

</div>

<div>maximum: 12.5168<br>

</div>

<div>range: 12.5168<br>

</div>

<div>mean: 1.49711<br>

</div>

<div>mean of absolute values: 1.49711<br>

</div>

<div>standard deviation: 1.97871<br>

</div>

<div>variance: 3.91528<br>

</div>

<div>variation coefficient: 132.169 %<br>

</div>

<div>sum: 185206.976940155<br>

</div>

<div><br>

</div>

<div>r.univar map="fill_dir_effect"<br>

</div>

<div>total null and non-null cells: 2025000<br>

</div>

<div>total null cells: 1964852<br>

</div>

<div><br>

</div>

<div>Of the non-null cells:<br>

</div>

<div>----------------------<br>

</div>

<div>n: 60148<br>

</div>

<div>minimum: 7.62939e-06<br>

</div>

<div>maximum: 10.9003<br>

</div>

<div>range: 10.9003<br>

</div>

<div>mean: 0.499953<br>

</div>

<div>mean of absolute values: 0.499953<br>

</div>

<div>standard deviation: 0.939429<br>

</div>

<div>variance: 0.882526<br>

</div>

<div>variation coefficient: 187.903 %<br>

</div>

<div>sum: 30071.1971092224<br>

</div>

<div><br>

</div>

<div>r.univar map="hydrodem_effect"<br>

</div>

<div>total null and non-null cells: 2025000<br>

</div>

<div>total null cells: 1020396<br>

</div>

<div><br>

</div>

<div>Of the non-null cells:<br>

</div>

<div>----------------------<br>

</div>

<div>n: 1004604<br>

</div>

<div>minimum: 3.8147e-06<br>

</div>

<div>maximum: 4.74686<br>

</div>

<div>range: 4.74685<br>

</div>

<div>mean: 0.00587574<br>

</div>

<div>mean of absolute values: 0.00587574<br>

</div>

<div>standard deviation: 0.08715<br>

</div>

<div>variance: 0.00759513<br>

</div>

<div>variation coefficient: 1483.22 %<br>

</div>

<div>sum: 5902.79515457153<br>

</div>

<span></span>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

<br>

</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Cheers</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Stefan<br>

</div>

<div id="appendonsend"></div>

<hr style="display:inline-block;width:98%" tabindex="-1">

<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Fra:</b> Markus Metz <markus.metz.giswork@gmail.com><br>

<b>Sendt:</b> fredag 14. juni 2019 14.25<br>

<b>Til:</b> Stefan Blumentrath<br>

<b>Kopi:</b> grass-dev@lists.osgeo.org<br>

<b>Emne:</b> Re: [GRASS-dev] r.terraflow vs. r.hydrodem</font>

<div> </div>

</div>

<div>

<div dir="ltr"><br>

<br>

On Fri, Jun 14, 2019 at 1:23 PM Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@nina.no">Stefan.Blumentrath@nina.no</a>> wrote:<br>

><br>

> Hi,<br>

><br>

> While trying to make r.terraflow work with larger regions, I was wondering if it would be doable to add a flag to r.hydrodem that forces the module to perform only sink filling and no breaching / carving?<br>

><br>

<div>> Sink filling is the function I need from r.terraflow but if it would be an option to add this function to r.hydrodem, that would be fine as well of course...</div>

<div><br>

</div>

<div>what about r.fill.dir? It should work as is with larger regions.</div>

<div><br>

</div>

<div>Markus M<br>

</div>

<br>

</div>

</div>

</body>

</html>