[geotk] Polar sterographic projection axes reversed in Geotk 3.10 CRS database

[Jonathan Blower]

Hi,

Up until recently, I have been using Geotk 3.04 - I downloaded the "referencing" bundle and put an epsg-wkt "database" on my classpath.  This all worked fine.  I've upgraded to 3.10, which seems to come with its own CRS database (I'm not sure what form this takes), but I'm having problems converting from polar stereographic to lon-lat coordinates.  Essentially the coordinates are reversed: when I convert the coordinates [x,y] from polar stereo to lon-lat, I get the result I would have expected from converting the coordinates [y,x].

The same happens for south polar stereographic projection.

Here's the result of CRS.decode("EPSG.32661").toString() in Geotk 3.10 (referencing bundle):

PROJCS["WGS 84 / UPS North", 
  GEOGCS["WGS 84", 
    DATUM["World Geodetic System 1984", 
      SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], 
      AUTHORITY["EPSG","6326"]], 
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
    UNIT["degree", 0.017453292519943295], 
    AXIS["Geodetic longitude", EAST], 
    AXIS["Geodetic latitude", NORTH], 
    AUTHORITY["EPSG","4326"]], 
  PROJECTION["Polar Stereographic (variant A)", AUTHORITY["EPSG","9810"]], 
  PARAMETER["central_meridian", 0.0], 
  PARAMETER["latitude_of_origin", 90.0], 
  PARAMETER["false_easting", 2000000.0], 
  PARAMETER["false_northing", 2000000.0], 
  PARAMETER["scale_factor", 0.994], 
  UNIT["metre", 1.0], 
  AXIS["Northing", "South along 180°"], 
  AXIS["Easting", "South along 90°E"], 
  AUTHORITY["EPSG","32661"]]


And here's the result from Geotk 3.04 with the embedded epsg-wkt database:

PROJCS["WGS 84 / UPS North", 
  GEOGCS["WGS 84", 
    DATUM["WGS_1984", 
      SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], 
      AUTHORITY["EPSG","6326"]], 
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
    UNIT["degree", 0.017453292519943295], 
    AXIS["Longitude", EAST], 
    AXIS["Latitude", NORTH], 
    AUTHORITY["EPSG","4326"]], 
  PROJECTION["Polar_Stereographic"], 
  PARAMETER["central_meridian", 0.0], 
  PARAMETER["latitude_of_origin", 90.0], 
  PARAMETER["false_easting", 2000000.0], 
  PARAMETER["false_northing", 2000000.0], 
  PARAMETER["scale_factor", 0.994], 
  UNIT["metre", 1.0], 
  AXIS["x", EAST], 
  AXIS["y", NORTH], 
  AUTHORITY["EPSG","32661"]]

As you can see, the axis definitions are different, and in reverse order.  I think the ordering from 3.04/epsg-wkt is correct.

Is there a problem with the database used in the 3.10 referencing bundle?  If so, is there a way to force it to use the epsg-wkt database until a fix is found?

Thanks, Jon

--
Dr Jon Blower
Technical Director, Reading e-Science Centre
Environmental Systems Science Centre
University of Reading
Harry Pitt Building, 3 Earley Gate
Reading RG6 6AL. UK
Tel: +44 (0)118 378 5213
Fax: +44 (0)118 378 6413
j.d.blower at reading.ac.uk
http://www.nerc-essc.ac.uk/People/Staff/Blower_J.htm