this bug's URL: http://intevation.de/rt/webrt?serial_num=3172
Request number 3172 was commented on by 'pkelly' (Paul Kelly).
Responding to this message will send mail to the requestor.
Request Tracker
rt@intevation.de
--------------------------------------------------------------
Cc: grass5@grass.itc.it
There definitely is a bug in there somewhere (and it may be 6.x-specific
rather than PROJ version-specific; I haven't been able to re-produce it in
5.4). d.where uses the PROJ function pj_latlong_from_proj() to convert from
the UTM projection definition to a latlong definition on the same ellipsoid.
This function is returning a value of es (eccentricity squared) with a very
low precision, only a few decimal places, and that seems to be what is
causing the error. I still haven't tracked down what is causing the low
number of decimal places.
Adding a call to pj_print_proj_params(&iproj, &oproj) before the pj_do_proj
line in display/d.where/where.c is helping to track this down.
-------------------------------------------- Managed by Request Tracker