<div dir="ltr"><div dir="auto">What exactly do you mean by median-filtering in time only, Nikos? Get the median and then filtering out all pixels above/below that value? If that's the case, then maybe some of the examples here <a href="https://grasswiki.osgeo.org/wiki/Temporal_data_processing#Spatio-temporal_algebra_with_STRDS">https://grasswiki.osgeo.org/wiki/Temporal_data_processing#Spatio-temporal_algebra_with_STRDS</a> might be of help<div dir="auto"><br></div><div>Basically, you first estimate the median of the series (or the period you are interested in) either with r.series or t.rast.series and then use t.rast.mapcalc or t.rast.algebra to get only pixels below/above that value represented pixelwise in the median map. It can also be a single value of course.</div><div><br></div><div>Does that make sense?</div><div><br></div><div>best, <br></div><div>Vero<br></div><div dir="auto"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mié., 27 ene. 2021 08:44,  <<a href="mailto:nik@nikosalexandris.net" target="_blank">nik@nikosalexandris.net</a>> escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I am trying to median-filter time series in time and not in space.<br>
<br>
Not sure `r.series.filter` is the right tool.<br>
I skimmed through the manual, the paper and the source code. Yet I am <br>
asking for a confirmation.<br>
Its (odd integer) parameter `winsize=` implies a moving window (of size <br>
winsize^2).<br>
Does this rather filter images, of a time series, in space?<br>
Or does the `winsize` here refer to the length of a moving temporal <br>
window?<br>
<br>
Then, there is `r.series.lwr`.<br>
Can I use the tool with `order=0` and `delta=1` (or `order=3` and <br>
`delta=1`)<br>
and perform a moving average filtering in time?<br>
It's not median. Still I'd like to see a filtering in time-only effect.<br>
<br>
And, now I recall `r.hants`.<br>
Is there a way to use it for filtering-out only lows and highs (`-l` and <br>
`-h` respectively) pixel-wise in time, without interpolating nor <br>
extrapolating?<br>
<br>
Else, I guess there should be a way to build a 'custom' filter using <br>
`t.rast.algebra` using its temporal neighborhood modifier.<br>
It's convenient if one the tools above does already median-filtering in <br>
time-only.<br>
<br>
Kind regards, Nikos<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" rel="noreferrer" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div>