[GRASS-user] Retrieve pixel position

Nikos Alexandris nikos.alexandris at felis.uni-freiburg.de
Tue Aug 3 17:56:55 EDT 2010


On Tuesday 03 of August 2010 19:11:27 Nikos Alexandris wrote:
> (forgot to cc to the list)
> 
> António:
> > Thanls for the tip. But isn't there any other more "interactive" way
> > (like r.what?)
> 
> I thought you have a points vector map and want to retrieve row,col for
> each point. Maybe "d.what.rast -c" is what you are looking for then?

In case I did think correctly that you want to script something, I think 
"r.stats -x" is the only way to get specific row,cols. I don't know any other 
module that will report rows and cols.

So, given that your points of interest are a vector map, then something like 
"v.to.rast" + "r.mapcalc ( to create a MASK)" + "r.stats -xg over the MASK", 
would give you coordinates and rowcol for _all_ non-NULL pixels in the MASK.

If you loop this (using "for") over your points of interest, will return for 
each point the requested values.

Hope this helps, Nikos


More information about the grass-user mailing list