[GRASS-dev] Extracting coordinates of RasterNumpy objects (pygrass)

[Pierre Roudier]

Hi,

I am using (with enthusiasm!) Pietro's pygrass library to develop a raster
module. I am using Numpy/Scipy as my working horse, so I manipulate a lot
of the RasterNumpy objects that have been introduced with pygrass.

In a specific step, I am identifying pixels using a test (this would be
similar to my_array <  100). I would like to extract the points that
satisfy the test, and access not only their values and index in the array,
but I would also go back to their coordinates to extract them as (x, y, z).

Is there a way to do this using RasterNumpy and pygrass?

Cheers,

Pierre
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20121217/5a0be636/attachment.html>