[GRASS5] [bug #3302] (grass) cs2cs with GRASS and NTv2 grid fails
Paul Kelly via RT
grass-bugs at intevation.de
Wed Jun 8 06:46:59 EDT 2005
this bug's URL: http://intevation.de/rt/webrt?serial_num=3302
Request number 3302 was commented on by 'pkelly' (Paul Kelly).
Responding to this message will send mail to the requestor.
Request Tracker
rt at intevation.de
--------------------------------------------------------------
Cc: grass5 at grass.itc.it
Yes the reason $GISBASE/etc/nad/ doesn't need to be included in the nadgrids
definition in the PROJ_INFO file is that the functions which prepare the
GRASS co-ordinate system definitions for use by PROJ (pj_get_kv() and
pj_get_string() in lib/proj/get_proj.c) call the PROJ function pj_set_finder
() to specify a name of a "finder function", which is provided by GRASS
(actually set_proj_lib() in get_proj.c; should really be re-written using
G_asprintf() if it is staying) and returns the location of the datum
gridshift files in $GISBASE/etc/nad/. PROJ then calls this function whenever
it wants during the transformation if it needs to find the files.
At the time we were putting the datum transformation support into GRASS,
Frank indicated that this was the preferred solution.
But looks like we've hit a limitation.
My peferred solution would be, as long as specifying the full path to the
gridshift file is supposed to work and isn't an undocumented feature of
PROJ, would be to change pj_get_kv() (which is called by g.proj) to prepend
$GISBASE/etc/nad/ to the name of the gridshift file when it prepares the
PROJ string, and to do away with pj_set_finder().
-------------------------------------------- Managed by Request Tracker
More information about the grass-dev
mailing list