[GRASS-dev] spline interpolation over time series

Markus Metz markus.metz.giswork at gmail.com
Sun Jun 16 13:33:14 PDT 2013


On Sun, Jun 16, 2013 at 10:20 PM, Margherita Di Leo
<dileomargherita at gmail.com> wrote:
> Hi,
>
>
> On Sun, Jun 16, 2013 at 9:19 PM, Markus Metz <markus.metz.giswork at gmail.com>
> wrote:
>>
>> On Sun, Jun 16, 2013 at 8:50 PM, Sören Gebbert
>> <soerengebbert at googlemail.com> wrote:
>> > Hi,
>> > Nikos is absolutely right, the ideal candidate for spline interpolation
>> > in
>> > time for raster maps is r.series.interp. There is no other module that
>> > performs temporal interpolation.
>>
>> r.hants does temporal interpolation, to a degree. Assuming you have
>> e.g the raster maps tseries.1, tseries.3, tseries.5, tseries.7,
>> tseries.9, you could create dummy tseries maps tseries.2, tseries.4,
>> tseries.6, tseries.7 with
>> r.mapcalc "tseries.2 = null()"
>> etc, then
>>
>> r.hants
>> in=tseries.1,tseries.2,tseries.3,tseries.4,tseries.5,tseries.6,tseries.7,tseries.8,tseries.9
>>
>> with the appropriate filtering options.
>>
>>
> Thank you for all your inputs. Actually what I should do is not a proper
> interpolation but a regression

That sounds like r.series (one time series) or r.regression.series
(regression of two time series).

Markus M


>
> --
> Best regards,
>
> Margherita DI LEO
> Postdoctoral Researcher
>
> European Commission - DG JRC
> Institute for Environment and Sustainability (IES)
> Via Fermi, 2749
> I-21027 Ispra (VA) - Italy - TP 261
>
> Tel. +39 0332 78 3600
> margherita.di-leo at jrc.ec.europa.eu
>
> Disclaimer: The views expressed are purely those of the writer and may not
> in any circumstance be regarded as stating an official position of the
> European Commission.


More information about the grass-dev mailing list