[GRASS-dev] PIxel-wise regression contribution

[Markus Metz]

Damien O'Grady wrote:
> Hi there
>
> I would like to contribute a script that I have used in my research to
> the broader field.  The script carries out linear regression, except
> that rather than using two rasters for the x and y variables, it uses
> any number of pairs of rasters, allowing the regression to be carried
> out across a time series, resulting in a slope, intercept and R^2 values
> for each pixel.

It seems that your script is doing the same like the C module
r.regression.series in grass-addons [0] ?

Markus M


[0] https://trac.osgeo.org/grass/browser/grass-addons/grass7/raster/r.regression.series

>
> The crunch is that it is written in perl, as this is what I know best
> and it serves me well.  I have done my best to conform to the stated
> standards, and I have produced a fairly detailed manual page.  This will
> give you a better idea of what the script does, with examples.  I attach
> the script and html here.  I hope that you will consider making it
> available to others.
>
> Here at James Cook University, by the way, we are encouraging a switch
> from problematic and costly proprietary software to GRASS GIS.  I am a
> major fan of GRASS, particularly on a Linux platform, and I am pleased
> that more and more students are listening.  You should see a sharp rise
> in its use here in Queensland in the near future.
>
> Thanks
>
> Damien O'Grady
>
>
> _______________________________________________
> grass-dev mailing list
> grass-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev