<div dir="auto">Wow, thank you very much Even! I think this solution is just fine! <div dir="auto"><br></div><div dir="auto">Thanks again!</div><div dir="auto"><br></div><div dir="auto">Best regards,</div><div dir="auto">Pedro</div><div dir="auto"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">A sábado, 27/11/2021, 19:38, Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<p>Petro,</p>
<p>I've queued <a href="https://github.com/OSGeo/gdal/pull/4908" target="_blank" rel="noreferrer">https://github.com/OSGeo/gdal/pull/4908</a> which will
"fix" that particular case (more precisely anything based on
ETRS89). I'm not super satisfied with this, and I believe there's
no miracle solution for all potential use cases. What GDAL returns
in both 2.4.0 and 3.4.0 is actually pretty much the same object,
despite the different representation (WKT1 vs WKT2). The root
cause of the issue is with the TOWGS84[] node. Those nodes are not
part of CRS definitions. They date back to the time where the GDAL
& PROJ stack wasn't capable of more powerful transformation,
and so a particular datum transformation to WGS 84 was picked up,
and encoded in the CRS definition. When reading such definitions
nowadays, we have a dilemna: did the user really meant to use the
particular transformation that was in TOWGS84[], or was it just
there because that was the default one. We ere on the side of
caution (most of time, there are a few exceptions there and there
in GDAL), and keep the transformation information, which makes
that this CRS is formally not the official definition, and thus
considered as ""custom""<br>
</p>
<p>However, for ETRS89, I believe we can safely remove it, since as
far as I know, there's no better transformation between ETRS89
(which is known modeled as a "datum ensemble" by EPSG, meaning any
of the 11 ETRF realizations) and WGS 84 (itself a datum ensemble)
than the TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0], and it is unlikely
there will be an alternative one in the future, given the dynamic
nature of the WGS 84 datum (maybe I'm speaking too fast. When
ETRS89 and WGS 84 coordinates will have drifted too far away in
the future, maybe there will some non-null transformation added in
the database, so that users who want to transform to a 'past' or
'future' coordinate epoch of WGS 84 can do that in a not too bad
way.)</p>
<p>(I should not that this .prj file isn't really compliant with the
Esri dialect of WKT, which is supposed to be used for shapefiles,
since normally it shouldn't contain TOWGS84 or AUHTORITY code)<br>
</p>
<p>Even<br>
</p>
<p><br>
</p>
<div>Le 27/11/2021 à 19:13, Pedro Venâncio a
écrit :<br>
</div>
<blockquote type="cite">
<div dir="ltr">Hi all,<br>
<br>
I'm seeing that QGIS 3.16 / 3.22 does not recognize the CRS of
specific layers generated by GeoServer, but the old QGIS 2.18.28
recognize.<br>
<br>
For instance, this shapefile: <a href="https://si.icnf.pt/shp/rnap" target="_blank" rel="noreferrer">https://si.icnf.pt/shp/rnap</a><br>
<br>
Checking ogrinfo in different versions:<br>
<br>
ogrinfo --version<br>
GDAL 2.4.0, released 2018/12/14 (Used by QGIS 2.18.28)<br>
<br>
ogrinfo -al rnap.shp<br>
INFO: Open of `rnap.shp'<br>
using driver `ESRI Shapefile' successful.<br>
<br>
Layer name: rnap<br>
Metadata:<br>
DBF_DATE_LAST_UPDATE=1921-11-27<br>
Geometry: Polygon<br>
Feature Count: 48<br>
Extent: (-123561.201500, -300415.400000) - (162054.118700,
268136.980100)<br>
Layer SRS WKT:<br>
PROJCS["ETRS89 / Portugal TM06",<br>
GEOGCS["ETRS89",<br>
DATUM["European Terrestrial Reference System 1989",<br>
SPHEROID["GRS 1980",6378137.0,298.257222101,<br>
AUTHORITY["EPSG","7019"]],<br>
TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],<br>
AUTHORITY["EPSG","6258"]],<br>
PRIMEM["Greenwich",0.0,<br>
AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.017453292519943295],<br>
AXIS["Geodetic longitude",EAST],<br>
AXIS["Geodetic latitude",NORTH],<br>
AUTHORITY["EPSG","4258"]],<br>
PROJECTION["Transverse_Mercator",<br>
AUTHORITY["EPSG","9807"]],<br>
PARAMETER["central_meridian",-8.133108333333334],<br>
PARAMETER["latitude_of_origin",39.66825833333334],<br>
PARAMETER["scale_factor",1.0],<br>
PARAMETER["false_easting",0.0],<br>
PARAMETER["false_northing",0.0],<br>
UNIT["m",1.0],<br>
AXIS["Easting",EAST],<br>
AXIS["Northing",NORTH],<br>
AUTHORITY["EPSG","3763"]]<br>
nome_ap: String (254.0)<br>
classifica: String (254.0)<br>
publica1: String (254.0)<br>
publica2: String (254.0)<br>
jurisdicao: String (254.0)<br>
area_ha: Real (33.15)<br>
sigla: String (254.0)<br>
<br>
<br>
<br>
ogrinfo --version<br>
GDAL 3.4.0, released 2021/11/04 (Used by QGIS 3.22)<br>
<br>
ogrinfo -al rnap.shp<br>
INFO: Open of `rnap.shp'<br>
using driver `ESRI Shapefile' successful.<br>
<br>
Layer name: rnap<br>
Metadata:<br>
DBF_DATE_LAST_UPDATE=1920-08-17<br>
Geometry: Polygon<br>
Feature Count: 48<br>
Extent: (-123561.201500, -300415.400000) - (162054.118700,
268136.980100)<br>
Layer SRS WKT:<br>
BOUNDCRS[<br>
SOURCECRS[<br>
PROJCRS["ETRS89 / Portugal TM06",<br>
BASEGEOGCRS["ETRS89",<br>
ENSEMBLE["European Terrestrial Reference System
1989 ensemble",<br>
MEMBER["European Terrestrial Reference Frame
1989"],<br>
MEMBER["European Terrestrial Reference Frame
1990"],<br>
MEMBER["European Terrestrial Reference Frame
1991"],<br>
MEMBER["European Terrestrial Reference Frame
1992"],<br>
MEMBER["European Terrestrial Reference Frame
1993"],<br>
MEMBER["European Terrestrial Reference Frame
1994"],<br>
MEMBER["European Terrestrial Reference Frame
1996"],<br>
MEMBER["European Terrestrial Reference Frame
1997"],<br>
MEMBER["European Terrestrial Reference Frame
2000"],<br>
MEMBER["European Terrestrial Reference Frame
2005"],<br>
MEMBER["European Terrestrial Reference Frame
2014"],<br>
ELLIPSOID["GRS 1980",6378137,298.257222101,<br>
LENGTHUNIT["metre",1]],<br>
ENSEMBLEACCURACY[0.1]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4258]],<br>
CONVERSION["Portugual TM06",<br>
METHOD["Transverse Mercator",<br>
ID["EPSG",9807]],<br>
PARAMETER["Latitude of natural
origin",39.6682583333333,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8801]],<br>
PARAMETER["Longitude of natural
origin",-8.13310833333333,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8802]],<br>
PARAMETER["Scale factor at natural origin",1,<br>
SCALEUNIT["unity",1],<br>
ID["EPSG",8805]],<br>
PARAMETER["False easting",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8806]],<br>
PARAMETER["False northing",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8807]]],<br>
CS[Cartesian,2],<br>
AXIS["easting (X)",east,<br>
ORDER[1],<br>
LENGTHUNIT["metre",1]],<br>
AXIS["northing (Y)",north,<br>
ORDER[2],<br>
LENGTHUNIT["metre",1]],<br>
USAGE[<br>
SCOPE["Topographic mapping (medium scale)."],<br>
AREA["Portugal - mainland - onshore."],<br>
BBOX[36.95,-9.56,42.16,-6.19]],<br>
ID["EPSG",3763]]],<br>
TARGETCRS[<br>
GEOGCRS["WGS 84",<br>
DATUM["World Geodetic System 1984",<br>
ELLIPSOID["WGS 84",6378137,298.257223563,<br>
LENGTHUNIT["metre",1]]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
CS[ellipsoidal,2],<br>
AXIS["latitude",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["longitude",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4326]]],<br>
ABRIDGEDTRANSFORMATION["ETRS89 to WGS 84 (1)",<br>
VERSION["EPSG-eur"],<br>
METHOD["Geocentric translations (geog2D domain)",<br>
ID["EPSG",9603]],<br>
PARAMETER["X-axis translation",0,<br>
ID["EPSG",8605]],<br>
PARAMETER["Y-axis translation",0,<br>
ID["EPSG",8606]],<br>
PARAMETER["Z-axis translation",0,<br>
ID["EPSG",8607]],<br>
USAGE[<br>
SCOPE["(null/copy) Approximation for medium and low
accuracy applications assuming equality between plate-fixed
static and earth-fixed dynamic CRSs, ignoring static/dynamic CRS
differences."],<br>
AREA["Europe - onshore and offshore: Albania;
Andorra; Austria; Belgium; Bosnia and Herzegovina; Bulgaria;
Croatia; Cyprus; Czechia; Denmark; Estonia; Faroe Islands;
Finland; France; Germany; Gibraltar; Greece; Hungary; Ireland;
Italy; Kosovo; Latvia; Liechtenstein; Lithuania; Luxembourg;
Malta; Moldova; Monaco; Montenegro; Netherlands; North
Macedonia; Norway including Svalbard and Jan Mayen; Poland;
Portugal; Romania; San Marino; Serbia; Slovakia; Slovenia;
Spain; Sweden; Switzerland; United Kingdom (UK) including
Channel Islands and Isle of Man; Vatican City State."],<br>
BBOX[32.88,-16.1,84.73,40.18]],<br>
ID["EPSG",1149],<br>
REMARK["ETRS89 and WGS 84 are realizations of ITRS
coincident to within 1 metre. This transformation has an
accuracy equal to the coincidence figure."]]]<br>
Data axis to CRS axis mapping: 1,2<br>
nome_ap: String (254.0)<br>
classifica: String (254.0)<br>
publica1: String (254.0)<br>
publica2: String (254.0)<br>
jurisdicao: String (254.0)<br>
area_ha: Real (33.15)<br>
sigla: String (254.0)<br>
<br>
The differences between rnap.prj (left) and <a href="https://epsg.io/3763" target="_blank" rel="noreferrer">https://epsg.io/3763</a>
(right) are few:<br>
<br>
<a href="https://cld.pt/dl/download/5d9b5324-b28a-4654-a8f7-365c7cc39c10/rnap_prj.png" target="_blank" rel="noreferrer">https://cld.pt/dl/download/5d9b5324-b28a-4654-a8f7-365c7cc39c10/rnap_prj.png</a><br>
<br>
The GeoServer version that is generating this layer is old.<br>
<br>
So, could this be a bug?<br>
<br>
Thanks!<br>
<br>
Best regards,<br>
Pedro</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
QGIS-Developer mailing list
<a href="mailto:QGIS-Developer@lists.osgeo.org" target="_blank" rel="noreferrer">QGIS-Developer@lists.osgeo.org</a>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank" rel="noreferrer">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank" rel="noreferrer">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a>
</pre>
</blockquote>
<pre cols="72">--
<a href="http://www.spatialys.com" target="_blank" rel="noreferrer">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</div>
_______________________________________________<br>
QGIS-Developer mailing list<br>
<a href="mailto:QGIS-Developer@lists.osgeo.org" target="_blank" rel="noreferrer">QGIS-Developer@lists.osgeo.org</a><br>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a><br>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a><br>
</blockquote></div>