[GRASS-dev] PIxel-wise regression contribution

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.

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

(attachments)

ogrady.tar.gz (5.11 KB)

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.

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.

could you add your script to grass-addons? so it wouldn't be maybe lost in
the mailing list.

grass-addons:

   http://grass.osgeo.org/wiki/GRASS_AddOns

you need a osgeo-userid to add things to grass addons:

   http://www.osgeo.org/osgeo_userid

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/PIxel-wise-regression-contribution-tp4975696p4975706.html
Sent from the Grass - Dev mailing list archive at Nabble.com.

I am not sure from the documentation how to do this. I have an
osgeo-userid, but no write access to the addons list (I assume). I
would be grateful if you would point me in the right direction.

Regards

Damien

On Sat, 2012-05-19 at 02:51 -0700, Helmut Kudrnovsky wrote:

>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.
>
>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.

could you add your script to grass-addons? so it wouldn't be maybe lost in
the mailing list.

grass-addons:

   http://grass.osgeo.org/wiki/GRASS_AddOns

you need a osgeo-userid to add things to grass addons:

   http://www.osgeo.org/osgeo_userid

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/PIxel-wise-regression-contribution-tp4975696p4975706.html
Sent from the Grass - Dev mailing list archive at Nabble.com.
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

On 21 May 2012 02:18, “Damien O'Grady” <damien.ogrady@my.jcu.edu.au> wrote:

I am not sure from the documentation how to do this. I have an
osgeo-userid, but no write access to the addons list (I assume). I
would be grateful if you would point me in the right direction.

Please read this

http://trac.osgeo.org/grass/wiki/HowToContribute#WriteaccesstotheGRASS-Addons-SVNrepository

Regards

Damien


Ciao
Luca

http://www.lucadelu.org
http://gis.cri.fmach.it/delucchi

I am not sure from the documentation how to do this. I have an
osgeo-userid, but no write access to the addons list (I assume). I
would be grateful if you would point me in the right direction.

just ask for addons-write-access in the dev-ML (and giving your
osgeo-userid), that's all.

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/PIxel-wise-regression-contribution-tp4975696p4975923.html
Sent from the Grass - Dev mailing list archive at Nabble.com.

On Mon, May 21, 2012 at 3:33 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

I am not sure from the documentation how to do this. I have an
osgeo-userid, but no write access to the addons list (I assume). I
would be grateful if you would point me in the right direction.

just ask for addons-write-access in the dev-ML (and giving your
osgeo-userid), that's all.

Well, please to the grass-psc list... and we also need the statement
about RFC2 as explained in

http://trac.osgeo.org/grass/wiki/HowToContribute#WriteaccesstotheGRASS-Addons-SVNrepository

Best
Markus

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@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Yes, you are quite right, and r.regression.series is far more elegant.
I wish I'd seen it three years ago, when I really needed it!

Thanks

Damien

On Tue, 2012-05-22 at 12:04 +0200, Markus Metz wrote:

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@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev

Markus Metz wrote:

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

Damien wrote:

Yes, you are quite right,

let's race 'em!

I wish I'd seen it three years ago, when I really needed it!

it's only 8 months old :slight_smile:

Hamish