Hi all,
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.
For instance, this shapefile: https://si.icnf.pt/shp/rnap
<https://si.icnf.pt/shp/rnap>
Checking ogrinfo in different versions:
ogrinfo --version
GDAL 2.4.0, released 2018/12/14 (Used by QGIS 2.18.28)
ogrinfo -al rnap.shp
INFO: Open of `rnap.shp'
using driver `ESRI Shapefile' successful.
Layer name: rnap
Metadata:
DBF_DATE_LAST_UPDATE=1921-11-27
Geometry: Polygon
Feature Count: 48
Extent: (-123561.201500, -300415.400000) - (162054.118700, 268136.980100)
Layer SRS WKT:
PROJCS["ETRS89 / Portugal TM06",
GEOGCS["ETRS89",
DATUM["European Terrestrial Reference System 1989",
SPHEROID["GRS 1980",6378137.0,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0.0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.017453292519943295],
AXIS["Geodetic longitude",EAST],
AXIS["Geodetic latitude",NORTH],
AUTHORITY["EPSG","4258"]],
PROJECTION["Transverse_Mercator",
AUTHORITY["EPSG","9807"]],
PARAMETER["central_meridian",-8.133108333333334],
PARAMETER["latitude_of_origin",39.66825833333334],
PARAMETER["scale_factor",1.0],
PARAMETER["false_easting",0.0],
PARAMETER["false_northing",0.0],
UNIT["m",1.0],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","3763"]]
nome_ap: String (254.0)
classifica: String (254.0)
publica1: String (254.0)
publica2: String (254.0)
jurisdicao: String (254.0)
area_ha: Real (33.15)
sigla: String (254.0)
ogrinfo --version
GDAL 3.4.0, released 2021/11/04 (Used by QGIS 3.22)
ogrinfo -al rnap.shp
INFO: Open of `rnap.shp'
using driver `ESRI Shapefile' successful.
Layer name: rnap
Metadata:
DBF_DATE_LAST_UPDATE=1920-08-17
Geometry: Polygon
Feature Count: 48
Extent: (-123561.201500, -300415.400000) - (162054.118700, 268136.980100)
Layer SRS WKT:
BOUNDCRS[
SOURCECRS[
PROJCRS["ETRS89 / Portugal TM06",
BASEGEOGCRS["ETRS89",
ENSEMBLE["European Terrestrial Reference System 1989
ensemble",
MEMBER["European Terrestrial Reference Frame 1989"],
MEMBER["European Terrestrial Reference Frame 1990"],
MEMBER["European Terrestrial Reference Frame 1991"],
MEMBER["European Terrestrial Reference Frame 1992"],
MEMBER["European Terrestrial Reference Frame 1993"],
MEMBER["European Terrestrial Reference Frame 1994"],
MEMBER["European Terrestrial Reference Frame 1996"],
MEMBER["European Terrestrial Reference Frame 1997"],
MEMBER["European Terrestrial Reference Frame 2000"],
MEMBER["European Terrestrial Reference Frame 2005"],
MEMBER["European Terrestrial Reference Frame 2014"],
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]],
ENSEMBLEACCURACY[0.1]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4258]],
CONVERSION["Portugual TM06",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",39.6682583333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",-8.13310833333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",1,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",0,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["easting (X)",east,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["northing (Y)",north,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["Topographic mapping (medium scale)."],
AREA["Portugal - mainland - onshore."],
BBOX[36.95,-9.56,42.16,-6.19]],
ID["EPSG",3763]]],
TARGETCRS[
GEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["latitude",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["longitude",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4326]]],
ABRIDGEDTRANSFORMATION["ETRS89 to WGS 84 (1)",
VERSION["EPSG-eur"],
METHOD["Geocentric translations (geog2D domain)",
ID["EPSG",9603]],
PARAMETER["X-axis translation",0,
ID["EPSG",8605]],
PARAMETER["Y-axis translation",0,
ID["EPSG",8606]],
PARAMETER["Z-axis translation",0,
ID["EPSG",8607]],
USAGE[
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."],
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."],
BBOX[32.88,-16.1,84.73,40.18]],
ID["EPSG",1149],
REMARK["ETRS89 and WGS 84 are realizations of ITRS coincident
to within 1 metre. This transformation has an accuracy equal to the
coincidence figure."]]]
Data axis to CRS axis mapping: 1,2
nome_ap: String (254.0)
classifica: String (254.0)
publica1: String (254.0)
publica2: String (254.0)
jurisdicao: String (254.0)
area_ha: Real (33.15)
sigla: String (254.0)
The differences between rnap.prj (left) and https://epsg.io/3763
<https://epsg.io/3763> (right) are few:
https://cld.pt/dl/download/5d9b5324-b28a-4654-a8f7-365c7cc39c10/rnap_prj.png
<https://cld.pt/dl/download/5d9b5324-b28a-4654-a8f7-365c7cc39c10/rnap_prj.png>
The GeoServer version that is generating this layer is old.
So, could this be a bug?
Thanks!
Best regards,
Pedro
_______________________________________________
QGIS-Developer mailing list
[email protected]
List info: https://lists.osgeo.org/mailman/listinfo/qgis-developer
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-developer
