[GRASS5] LatLong E/W problem in southern hemisphere

Markus Neteler neteler at itc.it
Fri Oct 15 09:19:09 EDT 2004


Hi,

we have a map in LatLong, located somewhere in Argentina:

g.region rast=noaa -p
projection: 3 (Latitude-Longitude)
zone:       0
datum:      wgs84
ellipsoid:  a=6378137 es=0.00669438
north:      30:24S
south:      35:42:00.0036S
west:       58:17:59.9964W
east:       65:11:59.9892W        <- note this correct value
nsres:      0:03:47.1429
ewres:      4:42:28.800096
rows:       84
cols:       75
                                                                                                          Mapset <PERMANENT> in Location <santafell>                                                                GRASS 5.7.cvs > g.region -g
n=-30.4
s=-35.700001
w=-58.299999
e=294.800003          <- note the wrong value
nsres=0.06309525
ewres=4.70800003

The bug seems to be somewhere in libgis, because:

r.in.bin ~/pub/avhrrpf.ndvi.1ntfsa.840101.11705583 output=noaa bytes=1 north=-30.400000  south=-35.700001  east=-65.199997 west=-58.299999 rows=84 cols=75
D0/0: E: 294.800003  N:-58.299999
D0/0: E: 294.800003  N:-58.299999
D0/0: E: 294.800003  N:-58.299999
Percent Complete:  100%
CREATING SUPPORT FILES FOR noaa

I have added the debug output to
lib/gis/adj_cellhd.c
(line 103)
   G_debug(0,"E: %f  N:%F",cellhd->east,cellhd->west);

Anyone working in the southern hemisphere and observing similar
problems?

Markus




More information about the grass-dev mailing list