Thanks a lot you both.
I had actually thought of using r.to.vect to get a point map with all
the pixel centres, then v.to.db to get their coordinates etc., but the
maps are too large and resolution too high. Obviously r.mapcalc was
the solution I was looking for, I just regret not noticing it
A. Marcia BARBOSA wrote:
What is the best way of getting a raster map (global coverage,
lat-long coordinates) in which the value of each pixel is its
latitude?
in a lat/lon location+mapset:
g.region n= s= w= e= res= ... (set the bounds for you new map)
r.mapcalc lat_map = y()
I expect