[Qgis-user] wrong interpretation of .prj-file

[Even Rouault]

Uwe,

do you know how / by which software this .prj is generated ? The 
',AUTHORITY["ESRI",102328]' in it is quite suspicious as it is malformed 
WKT (102328 should be quoted in double quotes), and Esri software 
doesn't no generate such authority node.

Looking in the current official projected CRS list of Esri at 
https://raw.githubusercontent.com/Esri/projection-engine-db-doc/master/csv/pe_list_projcs.csv, 
which is incorporated in PROJ, ESRI:102328 is 7 Stellen and there's no 
trace of a 6 Stellen CRS definition:

102328,102328,"PE_PCS_ETRS_1989_UTM_32N_7STELLEN","ETRS_1989_UTM_Zone_32N_7stellen","PROJCS[""ETRS_1989_UTM_Zone_32N_7stellen"",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"",2500000.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]]","PROJCRS[""ETRS_1989_UTM_Zone_32N_7stellen"",BASEGEOGCRS[""GCS_ETRS_1989"",DATUM[""D_ETRS_1989"",ELLIPSOID[""GRS_1980"",6378137.0,298.257222101,LENGTHUNIT[""Meter"",1.0]]],PRIMEM[""Greenwich"",0.0,ANGLEUNIT[""Degree"",0.0174532925199433]],CS[ellipsoidal,2],AXIS[""Latitude 
(lat)"",north,ORDER[1]],AXIS[""Longitude 
(lon)"",east,ORDER[2]],ANGLEUNIT[""Degree"",0.0174532925199433]],CONVERSION[""Transverse_Mercator"",METHOD[""Transverse_Mercator""],PARAMETER[""False_Easting"",2500000.0,LENGTHUNIT[""Meter"",1.0]],PARAMETER[""False_Northing"",0.0,LENGTHUNIT[""Meter"",1.0]],PARAMETER[""Central_Meridian"",9.0,ANGLEUNIT[""Degree"",0.0174532925199433]],PARAMETER[""Scale_Factor"",0.9996,SCALEUNIT[""Unity"",1.0]],PARAMETER[""Latitude_Of_Origin"",0.0,ANGLEUNIT[""Degree"",0.0174532925199433]]],CS[Cartesian,2],AXIS[""Easting 
(X)"",east,ORDER[1]],AXIS[""Northing 
(Y)"",north,ORDER[2]],LENGTHUNIT[""Meter"",1.0]]","ETRS 1989 UTM Zone 
32N 7Stellen","Esri","10.0.1","no","Germany - west of 
12~E",47.27,55.9,5.5,12.0

However googling for ETRS_1989_UTM_Zone_32N_6Stellen, I can find a few 
references to it using the same definition as the one you point to, but 
none of them seem to reference it as ESRI:102328.

If you replace the content of the .prj file with the following one 
(which uses equivalent EPSG:25832 definition), that should work:

PROJCS["ETRS_1989_UTM_Zone_32N",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",500000.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]]

I've also submitted a fix/improvement in PROJ per 
https://github.com/OSGeo/PROJ/pull/3532, so that it trusts less 
authority codes. With that change, QGIS will use EPSG:25832 for that file.

Even

Le 02/01/2023 à 12:00, Uwe via QGIS-User a écrit :
>
> Hello list,
>
> first: best wishes for a peaceful 2023 to everyone!
>
> I often have to process shapefiles that come with a prj file including 
> the „Authority: ESRI 102328“ statement (complete content of that file 
> see below).
>
> This .prj is processed wrong by QGIS 3.28, because QGIS sets the false 
> easting to 2.500.000 (which has to be 500.000) and it interpretes the 
> projection as „7 Stellen = decimal places“ which is „6 Stellen“ in fact.
>
> Can I change the way QGIS interpretes the .prj-File? It is annoying to 
> change the projection setting for every single shapefile I load in the 
> legend for each new shape delivery. In ESRI-Software, the .prj is 
> being read correct.
>
> Here is the complete .prj which is being interpreted wrong:
>
> PROJCS["ETRS_1989_UTM_Zone_32N_6Stellen",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",500000.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],AUTHORITY["ESRI",102328]]
>
> Mit freundlichen Grüßen,
>
> Uwe Fischer
>
> Ingenieurbüro Fischer
>
> Esbecker Str. 8
>
> *31036 Eime*
>
> /T. 05182/8325/
>
> /T. 0172/8876934/
>
>
> _______________________________________________
> QGIS-User mailing list
> QGIS-User at lists.osgeo.org
> List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
> Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user

-- 
http://www.spatialys.com
My software is free, but my time generally not.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-user/attachments/20230102/8652ef8a/attachment.htm>