[GRASS-dev] programming raster: how to get a latitude value of the pixel to be processed, even if mapset is not in Lat/lon?

GRASS GIS GRASS GIS developer
1

Hello list,
i would need to get a latitude value of the pixel to be processed,
even if it is not in a Lat/lon projection system. is there a way to
get this done from within GRASS libraries (or using some proj lib
maybe?).

thank you,
Yann

2

Yann Chemin wrote:

i would need to get a latitude value of the pixel to be processed,
even if it is not in a Lat/lon projection system. is there a way to
get this done from within GRASS libraries (or using some proj lib
maybe?).

You need to use the PROJ library, see the r.proj source code for an
example of initialising a projection from the location's PROJ_INFO and
PROJ_UNITS files.

--
Glynn Clements <glynn@gclements.plus.com>