[STATSGRASS] Re: Problem with gmeta6() in spgrass6.

Markus Neteler neteler at itc.it
Wed Jun 21 16:06:14 EDT 2006


Dear Roger,

[to also contribute something here]

to fix this problem, simply change the call/parsing in
gmeta6():

The original code does:

> system("g.region -p3")
projection: 3 (Längen- und Breitengrad)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      46:06:39N
south:      45:57:18N
west:       11:02:42E
east:       11:19:06E
top:        1.00000000
bottom:     0.00000000
nsres:      0:00:03
nsres3:     0:00:03
ewres:      0:00:03
ewres3:     0:00:03
tbres:      1
rows:       187
rows3:      187
cols:       328
cols3:      328
depths:     1

You could change to -g:

> system("g.region -g3")
n=46.11083333
s=45.955
w=11.045
e=11.31833333
t=1
b=0
nsres=0.00083333
nsres3=0.00083333
ewres=0.00083333
ewres3=0.00083333
tbres=1
rows=187
cols=328
rows3=187
cols3=328
depths=1

which is easily parsable and also works in LatLong.

Best
 Markus




More information about the grass-stats mailing list