Hi all,
I have added a "-l" flag to g.region to output the region settings
in lat/long format (it performs 4 transformations for that).
Example:
g.region -p
projection: 99 (Transverse Mercator)
zone: 0
datum: potsdam
ellipsoid: bessel
north: 5776000
south: 5756000
west: 3558000
east: 3578000
nsres: 5
ewres: 5
rows: 4000
cols: 4000
g.region -l
lat: 52.11605 long: 9.84682 (north/west corner)
lat: 52.11359 long: 10.13878 (north/east corner)
lat: 51.93630 long: 9.84343 (south/west corner)
lat: 51.93386 long: 10.13422 (south/east corner)
rows: 4000
cols: 4000
This shouldn't break anything, the result looks reasonable (as it is
using PROJ4 transformation).
Why that? Because many geodata web services want to know your region
of interest in lat/long, however, the users are mostly using local
coordinate systems. The "-l" flag transforms the current settings on
the fly from current projection information into lat/long.
If you want to polish the "-l" output, go ahead (g.region/main.c).
Cheers
Markus
--
Markus Neteler * University of Hannover
Institute of Physical Geography and Landscape Ecology
Schneiderberg 50 * D-30167 Hannover * Germany
Tel: ++49-(0)511-762-4494 Fax: -3984
----------------------------------------
If you want to unsubscribe from GRASS Development Team mailing list write to:
minordomo@geog.uni-hannover.de with
subject 'unsubscribe grass5'