No subject


Wed Nov 14 13:37:45 EST 2007


a AVHRR data set which I want to store in a Goode Homolosine LOCATION.

The lat/long boundaries are (stored into file 'latlong'):
-9.1 47.3
-9.1 -45.1
12.9 47.3
12.9 -45.1

RESOLUTION: 8000m

To get the boundaries in Goode projection, I used m.proj2 (enter in one
line):

m.proj2 inproj="proj=ll" outproj="proj=goode,name=goode,
ellps=wgs84,lat_0=0.0,lon_0=20.0,unfact=1.0"
input=latlong output=goode

Result:
-2289443.208536 5251097.248241
-2347834.155174 -5014324.412190
-558592.672873 5251097.248241
-572839.261228 -5014324.412190

Accordingly I selected:
N: 5251097.248   -> 5251100
S: -5014324.412  -> -5014300
W: -2347834.155  -> -2347800
E: -558592.672   -> -558600
East-West: 8000
North-South: 8000

  north:       5251100
  south:       -5014300
  east:        -558600
  west:        -2347800

  e-w res:     7987.5  (Changed to conform to grid)
  n-s res:     8001.09119252  (Changed to conform to grid)

  total rows:             1283
  total cols:              224
  total cells:         287,392


Now the test:
g.region -l
long: -9.84175 lat: 47.30003 (north/west corner)
long: 12.89990 lat: 47.30003 (north/east corner)
long: 13.07650 lat: -45.09977 (south/east corner)
long: -9.09950 lat: -45.09977 (south/west corner)
rows:       1283
cols:       224
Center Longitude: 1:45:31.635386E [1.75879]
Center latitude:  1:06:00.451679N [1.10013]

This is very strange, as the result must be within 
N: 47.3, S: -45.1, E: 12.9, W: -9.1

Perhaps someone has a hint for me what's wrong above.

Thanks

 Markus



More information about the grass-dev mailing list