<div dir="ltr"><div>Hi Even,</div><div><br></div><div>thank you, it turned out that the right projection of the raster is ETRS89 Lambert Azimutal Equal Area (IGNF:ETRS89LAEA), sorry for the noise :)<br></div><div></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Mar 21, 2022 at 5:04 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<p>Alessandro,</p>
<p>the truth (and sometimes issues) is in projinfo:</p>
<p>$ projinfo -s EPSG:4258 -t EPSG:4326<br>
Candidate operations found: 1<br>
-------------------------------------<br>
Operation No. 1:<br>
<br>
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.<br>
<br>
PROJ string:<br>
+proj=noop<br>
<br>
WKT2:2019 string:<br>
COORDINATEOPERATION["ETRS89 to WGS 84 (1)",<br>
VERSION["EPSG-eur"],<br>
SOURCECRS[<br>
GEOGCRS["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>
CS[ellipsoidal,2],<br>
AXIS["geodetic latitude (Lat)",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["geodetic longitude (Lon)",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4258]]],<br>
TARGETCRS[<br>
GEOGCRS["WGS 84",<br>
ENSEMBLE["World Geodetic System 1984 ensemble",<br>
MEMBER["World Geodetic System 1984 (Transit)"],<br>
MEMBER["World Geodetic System 1984 (G730)"],<br>
MEMBER["World Geodetic System 1984 (G873)"],<br>
MEMBER["World Geodetic System 1984 (G1150)"],<br>
MEMBER["World Geodetic System 1984 (G1674)"],<br>
MEMBER["World Geodetic System 1984 (G1762)"],<br>
MEMBER["World Geodetic System 1984 (G2139)"],<br>
ELLIPSOID["WGS 84",6378137,298.257223563,<br>
LENGTHUNIT["metre",1]],<br>
ENSEMBLEACCURACY[2.0]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
CS[ellipsoidal,2],<br>
AXIS["geodetic latitude (Lat)",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["geodetic longitude (Lon)",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4326]]],<br>
METHOD["Geocentric translations (geog2D domain)",<br>
ID["EPSG",9603]],<br>
PARAMETER["X-axis translation",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8605]],<br>
PARAMETER["Y-axis translation",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8606]],<br>
PARAMETER["Z-axis translation",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8607]],<br>
OPERATIONACCURACY[1.0],<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>
<br>
</p>
<p>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.<br>
</p>
<p>Even<br>
</p>
<div>Le 21/03/2022 à 16:56, Alessandro
Pasotti via QGIS-Developer a écrit :<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>Hi,</div>
<div><br>
</div>
<div>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. <br>
</div>
<div><br>
</div>
<div>The issue can be reproduced from the console too:</div>
<div><br>
</div>
<div>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">ct = QgsCoordinateTransform(QgsCoordinateReferenceSystem(4258), QgsCoordinateReferenceSystem(4326), QgsProject.instance().transformContext())
</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">ct.transformBoundingBox(QgsRectangle(4153600, 1296000, 5000000, 2000000))</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">ct.isValid()</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">True
</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">
</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">Result:
</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap"><QgsRectangle: 4153600 1296000, 5000000.00000001396983862 1999999.99999999720603228></p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">
</p>
<p style="margin:0px;text-indent:0px;white-space:pre-wrap">Any ideas?</p>
</div>
<div><br>
</div>
<div>-- <br>
<div dir="ltr">
<div dir="ltr">
<div>Alessandro Pasotti</div>
<div>QCooperative: <a href="https://www.qcooperative.net" target="_blank">www.qcooperative.net</a><br>
</div>
ItOpen: <a href="http://www.itopen.it" target="_blank">www.itopen.it</a></div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
QGIS-Developer mailing list
<a href="mailto:QGIS-Developer@lists.osgeo.org" target="_blank">QGIS-Developer@lists.osgeo.org</a>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-developer</a>
</pre>
</blockquote>
<pre cols="72">--
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div>Alessandro Pasotti</div><div>QCooperative: <a href="https://www.qcooperative.net" target="_blank">www.qcooperative.net</a><br></div>ItOpen: <a href="http://www.itopen.it" target="_blank">www.itopen.it</a></div></div>