<div dir="ltr"><div><br><br>On Thu, Mar 30, 2017 at 2:46 AM, Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>><br>> On Wed, 29 Mar 2017, Vaclav Petras wrote:<br>><br>>> As Soeren mentioned in the other part of the thread, use r.univar to get<br>>> the statistics. Don't rely on just looking at the resulting raster map.<br>>> There can be things like outliers which may throw the visualization off.<br>><br>><br>> Vaclav,<br>><br>>   Here are the results of r.univar on the output of r.mapcalc "output_diff =<br>> blocked_h_max - open_h_max":<br>><br>> [...]<br>> min=-1.89342188835144<br>> max=0.00093058729544282<br>> [...]<br>><br>>  [...] I see a small range and many null<br>> cells, but I expect that because 4 cells on one input DEM are set about 0.6m<br>> higher elevation than those on the other input DEM.<br><br></div>Apparently something went wrong because min should be 0 and max should be about 0.6.<br><div><br></div><div>The point locations used to edit the original DEM are probably wrong (not close to the center), and the computational region might have been changed somewhere in the process of editing the original DEM.<br><br></div><div>Markus M<br></div></div>