[QGIS-Developer] issue transforming coords from EPSG:4258 'ETRS89' to 4326
Alessandro Pasotti
apasotti at gmail.com
Mon Mar 21 09:21:13 PDT 2022
Hi Even,
thank you, it turned out that the right projection of the raster is ETRS89
Lambert Azimutal Equal Area (IGNF:ETRS89LAEA), sorry for the noise :)
On Mon, Mar 21, 2022 at 5:04 PM Even Rouault <even.rouault at spatialys.com>
wrote:
> Alessandro,
>
> the truth (and sometimes issues) is in projinfo:
>
> $ projinfo -s EPSG:4258 -t EPSG:4326
> Candidate operations found: 1
> -------------------------------------
> Operation No. 1:
>
> EPSG:1149, ETRS89 to WGS 84 (1), 1.0 m, 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.
>
> PROJ string:
> +proj=noop
>
> WKT2:2019 string:
> COORDINATEOPERATION["ETRS89 to WGS 84 (1)",
> VERSION["EPSG-eur"],
> SOURCECRS[
> GEOGCRS["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]],
> CS[ellipsoidal,2],
> AXIS["geodetic latitude (Lat)",north,
> ORDER[1],
> ANGLEUNIT["degree",0.0174532925199433]],
> AXIS["geodetic longitude (Lon)",east,
> ORDER[2],
> ANGLEUNIT["degree",0.0174532925199433]],
> ID["EPSG",4258]]],
> TARGETCRS[
> GEOGCRS["WGS 84",
> ENSEMBLE["World Geodetic System 1984 ensemble",
> MEMBER["World Geodetic System 1984 (Transit)"],
> MEMBER["World Geodetic System 1984 (G730)"],
> MEMBER["World Geodetic System 1984 (G873)"],
> MEMBER["World Geodetic System 1984 (G1150)"],
> MEMBER["World Geodetic System 1984 (G1674)"],
> MEMBER["World Geodetic System 1984 (G1762)"],
> MEMBER["World Geodetic System 1984 (G2139)"],
> ELLIPSOID["WGS 84",6378137,298.257223563,
> LENGTHUNIT["metre",1]],
> ENSEMBLEACCURACY[2.0]],
> PRIMEM["Greenwich",0,
> ANGLEUNIT["degree",0.0174532925199433]],
> CS[ellipsoidal,2],
> AXIS["geodetic latitude (Lat)",north,
> ORDER[1],
> ANGLEUNIT["degree",0.0174532925199433]],
> AXIS["geodetic longitude (Lon)",east,
> ORDER[2],
> ANGLEUNIT["degree",0.0174532925199433]],
> ID["EPSG",4326]]],
> METHOD["Geocentric translations (geog2D domain)",
> ID["EPSG",9603]],
> PARAMETER["X-axis translation",0,
> LENGTHUNIT["metre",1],
> ID["EPSG",8605]],
> PARAMETER["Y-axis translation",0,
> LENGTHUNIT["metre",1],
> ID["EPSG",8606]],
> PARAMETER["Z-axis translation",0,
> LENGTHUNIT["metre",1],
> ID["EPSG",8607]],
> OPERATIONACCURACY[1.0],
> 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."]]
>
> So basically ETRS89 and WGS 84 are considered the same with a 1 metre
> accuracy when considering the generic/datum ensemble versions of EPSG:4258
> and EPSG:4326. You might possibly get time-dependent transformations if
> considering one of the ETRFxxxx realization and a WGS 84 (Gyyyy) / ITRFzzzz
> realization instead.
>
> Even
> Le 21/03/2022 à 16:56, Alessandro Pasotti via QGIS-Developer a écrit :
>
> Hi,
>
> I'm not sure if this is a local issue on my dev env or if it is a real
> issue, I have a raster layer in EPSG:4258 ('ETRS89') from copernicus and
> the project's CRS is 4326, the transformation silently fails and the raster
> is not reprojected.
>
> The issue can be reproduced from the console too:
>
> ct = QgsCoordinateTransform(QgsCoordinateReferenceSystem(4258),
> QgsCoordinateReferenceSystem(4326),
> QgsProject.instance().transformContext())
>
> ct.transformBoundingBox(QgsRectangle(4153600, 1296000, 5000000, 2000000))
>
> ct.isValid()
>
> True
>
> Result:
>
> <QgsRectangle: 4153600 1296000, 5000000.00000001396983862
> 1999999.99999999720603228>
>
> Any ideas?
>
> --
> Alessandro Pasotti
> QCooperative: www.qcooperative.net
> ItOpen: www.itopen.it
>
> _______________________________________________
> QGIS-Developer mailing listQGIS-Developer at lists.osgeo.org
> List info: https://lists.osgeo.org/mailman/listinfo/qgis-developer
> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-developer
>
> -- http://www.spatialys.com
> My software is free, but my time generally not.
>
>
--
Alessandro Pasotti
QCooperative: www.qcooperative.net
ItOpen: www.itopen.it
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-developer/attachments/20220321/0a862518/attachment.html>
More information about the QGIS-Developer
mailing list