[GRASS5] [bug #3302] (grass) cs2cs with GRASS and NTv2 grid fails

Hamish hamish_nospam at yahoo.com
Wed Jun 8 08:46:16 EDT 2005


> this bug's URL: http://intevation.de/rt/webrt?serial_num=3302
..
> 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().



temporary work-around to fix v.in.garmin's WGS84->GRASS location bug:


IN_PROJ="+proj=longlat +towgs84=0.000,0.000,0.000"
OUT_PROJ="`g.proj -jf | sed -e
's/+nadgrids=/+nadgrids=$GISBASE\/etc\/nad\//'`"

...

cs2cs -f %.7f $IN_PROJ +to $OUT_PROJ < "${TMP}.base" > "${TMP}.P_base"



I'd like to work towards a general purpose m.proj2 module again with -i
and -o flags to very simply convert a list of points to/from WGS84 & the
current location settings. As above, it's easy to do with cs2cs in a
script, but I still haven't figured out how to pipe into a script well 
'cause "read" is kind of crappy.. so perhaps C is better albeit
reinventing the wheel & something else that only a smaller part of the
crowd can maintain, etc.



Hamish




More information about the grass-dev mailing list