Hi,

Hope someone can help me... I did get a shapefile, I believe produced by QGIS, 
with .prj file:

PROJCS["ETRS_1989_UTM_Zone_N32",
    GEOGCS["GCS_ETRS_1989",
       DATUM["D_ETRS_1989",
          SPHEROID["GRS_1980",6378137.0,298.257222101]
       ],
       PRIMEM["Greenwich",0.0],
       UNIT["Degree",0.0174532925199433]
    ],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",32500000.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",9.0],
    PARAMETER["Scale_Factor",0.9996],
    PARAMETER["Latitude_Of_Origin",0.0],
    UNIT["Meter",1.0]
]

If I feed this WKT to the code:

CoordinateReferenceSystem crs = CRS.parseWKT(wkt);
Integer lookupResult = CRS.lookupEpsgCode(crs, true);

.. I do get number 5652, which is correct, but has northing-easting axis order 
(https://epsg.io/5652). I'd rather prefer 4647, which is the same, but with 
easting-northing. Is there a way to get it? I'm using gt-epsg-hsql plugin, here 
is the gradle classpath:

implementation 'org.geotools:gt-api'
implementation 'org.geotools:gt-referencing'
implementation 'org.geotools:gt-main'
implementation 'org.geotools:gt-epsg-hsql'
implementation 'org.geotools:gt-coverage'

Geotools version is 34.1.

// Strange thing, if I copy-paste Implementation of CRS class and set 
FORCE_LONGITUDE_FIRST_AXIS_ORDER to false, I still don't get 4647, but null. 
Might it be a bug, maybe 4647 is not in the hsql database?

With best regards, Slava


_______________________________________________
GeoTools-GT2-Users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users

Reply via email to