<div dir="ltr">Hi,<br><br><div class="gmail_quote"><div dir="ltr">El mié., 5 dic. 2018 a las 11:04, Moritz Lennert (<<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On 5/12/18 10:47, Margherita Di Leo wrote:<br>
> Hi,<br>
> <br>
> I have a strds and I try to get statistics for an area within a mask. A <br>
> couple of images have no data within the masked area, and t.rast.univar <br>
> gives the warning message:  "Unable to get statistics for raster map". <br>
> In the resulting csv, the missed dates are skipped, and this is sort of <br>
> unexpected and creates hurdles when for example I compare it with the <br>
> behavior of t.rast.what, that keeps the date and associates a no-data to <br>
> it. Wouldn't it be better if also t.rast.univar would yield a no-data <br>
> rather than skip the date?<br></blockquote><div><br></div><div>+1 <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
This can probably somehow be solved within t.rast.univar, but IIUC the <br>
actual issue comes from the underlying r.univar call which has the same <br>
behaviour, for example when call with the 'zones' parameter and the '-t' <br>
flag: data is just absent if there are only null values in a given zone. <br>
I've been struggling with that, for example for addon modules such as <br>
i.segment.stats. I don't know how easy it would be to change that within <br>
r.univar, but it would be nice.<br>
<br>
IIUC, the issue comes from stats.c [1]:<br>
<br>
123             if (stats[z].n == 0)<br>
124                 continue;<br>
<br>
Maybe a flag to fill the stats with NULL values, instead of ignoring <br>
them, would be appropriate ? </blockquote><div><br></div><div>+1!!</div><div> <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">But I don't know what an good NULL value would be here.<br></blockquote><div> </div><div>IIRC, <a href="http://r.out.xyz">r.out.xyz</a> and <a href="http://t.rast.out.xyz">t.rast.out.xyz</a> use the '*' character. <br></div><div>It would be good to create a ticket for this, no? As to not forget afterwards...<br></div></div><div class="gmail_quote"><div><br></div>my 0.01 cents</div><div class="gmail_quote">Vero<br></div></div>