[GRASS-dev] latlong problem

Markus Neteler neteler at itc.it
Wed Mar 7 08:53:28 EST 2007


On Wed, Mar 07, 2007 at 02:03:58PM +0100, Martin Landa wrote:
> Hi,
> 
> there is a significant problem in LatLong location. Try the following:
> 
> g.region n=80 s=-80 w=-170 e=170 res=1;
> v.in.region out=reg;
> g.region n=90 s=-90 w=-180 e=180 res=1;
> d.erase;d.vect reg
> [1] should be [2]
> 
> I prepared a simple patch to fix this bug.
> 
> I am not sure whether it is OK. Anyway testing "(east1 - east2) > 180"
> does not make sense to me... or am I wrong?
> 
> Martin
> 
> [1] http://gama.fsv.cvut.cz/~landa/.../grass/reg.png
> [2] http://gama.fsv.cvut.cz/~landa/.../grass/reg2.png

Great work, Martin!

This patch cures the problem of
 http://mpa.itc.it/markus/tmp/admin98_oceans_buggy.png
to
 http://mpa.itc.it/markus/tmp/admin98_oceans_fixed.png

Unrelated bug (?): The boundary box at the East end isn't
visible - is the display driver cutting off by 1 pixel from
East? The region is 90N, 90S, 180W, 180E, Latlong WGS84.
See admin98_oceans_fixed.png for illustration, generated
with d.out.png.

Markus




More information about the grass-dev mailing list