[geotk] comparing CRS

[Martin Desruisseaux]

Hello Philipe

The first CRS:

Le 26/08/10 13:18, Philippe Kruschitz a écrit :
> PROJCS["MGI / Austria GK M31",
>   [...snip...]
> AXIS["Northing", NORTH],
> AXIS["Easting", EAST],
> AUTHORITY["EPSG","31258"]]

has (northing, easting) axes. But the CRS read from the shapefile:

>
> An here is the problematic shapefile projection (created by ArcGIS):
>
> PROJCS["MGI / Austria GK M31",
>    [...snip...]
> UNIT["m",1.0]]

does not declare any axis order. In such case, the WKT specification mandates 
that the default order is (easting, northing).

Did you tried to create the first CRS in Geotk as below?

     CoordinateReferenceSystem crs = CRS.decode("EPSG:31258", true);

It should work (the 'true' boolean argument means "force (x,y) axis order"). If 
it doesn't work, then I would need to look at the issue more closely.

    Regards,

        Martin