[GRASS-dev] more proj4 inconsistencies in g.proj

[Michael Barton]

Paul and Glyn,

I've been tearing my hair out for 2 days over what I thought was a Python
problem only to discover that it's another inconsistency with the projection
information files--with ellipse.table again.

For example, the line for wgs84 is...

wgs84        "World Geodetic System 1984"    a=6378137.000
f=1/298.257223563

But for it to be read into g.proj in acceptable format, it needs to be...

wgs84        "World Geodetic System 1984"    a=6378137.000
rf=298.257223563

That is, 1/298.257223563 needs to be 298.257223563

Why is it 1/298.257223563? Is this correct in some kind of format acceptable
to g.proj?

Anyway, with this solved, the location wizard works for all except xy
locations (suggestions on how to make?) and extents.

Michael



__________________________________________
Michael Barton, Professor of Anthropology
Director of Graduate Studies
School of Human Evolution & Social Change
Center for Social Dynamics & Complexity
Arizona State University

phone: 480-965-6213
fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-dev/attachments/20070705/5c5d3027/attachment.html