<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Feb 14, 2015 at 7:29 PM, Vaclav Petras <span dir="ltr"><<a href="mailto:wenzeslaus@gmail.com" target="_blank">wenzeslaus@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Sat, Feb 14, 2015 at 11:47 AM, Paulo van Breugel <span dir="ltr"><<a href="mailto:p.vanbreugel@gmail.com" target="_blank">p.vanbreugel@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br><div class="gmail_quote"><div>For a quick solution, what about using r.tile to split the input data in tiles and compute the mahalanobis distance per tile.<span><font color="#888888"><br></font></span></div><span></span></div></blockquote><div><br></div></span><div>See PyGRASS GridModule class which will do tiling (and a lot of other things) for you.<br><br><a href="http://grass.osgeo.org/grass70/manuals/libpython/pygrass.modules.grid.html" target="_blank">http://grass.osgeo.org/grass70/manuals/libpython/pygrass.modules.grid.html</a> <br></div></div></div></div></blockquote><div><br></div>The Grid Module does indeed what I want, except that I want to run a python function (the mahalanobis distance computation) on the tiles, while the GridModule requires a raster GRASS r.* command. Any way around that? Writing an intermediate addon/function that I can then call with the GridModule?<br><br></div><br></div></div>