[GRASS-dev] Compute mahalanobis distance using Scipy

[Paulo van Breugel]

Hi Javier,

This looks really useful, thanks! My python skills are very limited, so I
will have to look at it better to see what I can understand. But to start,
the parallel function, does that also solve the handling of potentially
very large raster layers?

Btw, this does look like a very useful set of functions, any plan to
develop this into a GRASS addon or function?

Paulo


On Sun, Feb 15, 2015 at 9:22 PM, Javier Martínez-López <
javi.martinez.lopez at gmail.com> wrote:

> Hi Paulo,
>
> you can see an implementation of the mahalanobis distance computed in
> parallel (using all computer processors) for a single image here:
>
> https://github.com/javimarlop/eHabpy/blob/master/ehab.py
>
> see lines 32-86 and 597. I hope it helps! Cheers, Javier
>
> On Sat, Feb 14, 2015 at 7:29 PM, Vaclav Petras <wenzeslaus at gmail.com>
> wrote:
> >
> > On Sat, Feb 14, 2015 at 11:47 AM, Paulo van Breugel <
> p.vanbreugel at gmail.com>
> > wrote:
> >>
> >>
> >> For a quick solution, what about using r.tile to split the input data in
> >> tiles and compute the mahalanobis distance per tile.
> >
> >
> > See PyGRASS GridModule class which will do tiling (and a lot of other
> > things) for you.
> >
> >
> http://grass.osgeo.org/grass70/manuals/libpython/pygrass.modules.grid.html
> >
> > _______________________________________________
> > grass-dev mailing list
> > grass-dev at lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/grass-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20150215/13b0a59c/attachment.html>