MapInfo CRS not correctly recognized

[meyerlor]

Hello,
I m using GDAL 3.4.1 with QGIS 3.22.3

When exporting any file to a MapInfo .tab into the CRS EPSG:25832, opening that exported file is incorrectly recognized as EPSG:32632.
gdalinfo of that file reports:

Layer name: BORIS_BW_2020

Geometry: Unknown (any)
Feature Count: 296
Extent: (517700.150000, 5281321.900000) - (537137.040000, 5302830.140000)
Layer SRS WKT:
PROJCRS["unnamed",
BASEGEOGCRS["unnamed",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]]],
CONVERSION["Transverse Mercator",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",9,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",0.9996,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",500000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]
Data axis to CRS axis mapping: 1,2

It should be ETRS and not WGS84 i think. I opened an issue some time ago for QGIS, but i think here it is better placed, if not please delete (qgis/QGIS#43482) !

[jratike80]

I could not reproduce with GDAL 3.4.0dev on Windows.

I had a try both with tab and mif

ogr2ogr -f "MapInfo File" states.tab -a_srs  EPSG:25832 states.shp
ogr2ogr -f "MapInfo File" states.mif -a_srs  EPSG:25832 states.shp

and then

ogrinfo states.tab -al -so

Layer name: states
Geometry: Unknown (any)
Feature Count: 49
Extent: (-124.730000, 24.950000) - (-66.970000, 49.370000)
Layer SRS WKT:
PROJCRS["unnamed",
    BASEGEOGCRS["unnamed",
        DATUM["Euref_89",
            ELLIPSOID["GRS 80",6378137,298.257222101,
 ...

Or with gdalsrsinfo

gdalsrsinfo states.mif

PROJ.4 : +proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs

[rouault]

I've a fix pending. The issue is when using the WKT2 definition of EPSG:25832, which is what QGIS does under the hood, which lacks an explicit datum code.

[jratike80]

QGIS 3.18 finds https://epsg.org/crs_7791/RDN2008-UTM-zone-32N.html from the mif file that I created but I really made that shapefile with an odd method (attach srs UTM32 into EPSG:4326 data in degrees).

[meyerlor]

QGIS 3.18 finds https://epsg.org/crs_7791/RDN2008-UTM-zone-32N.html from the mif file that I created but I really made that shapefile with an odd method (attach srs UTM32 into EPSG:4326 data in degrees).

that changed in QGIS 3.22 from EPSG7791 to EPSG32632

@rouault thats great! looking forward to this release then :)