<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p><br>
</p>
<div class="moz-cite-prefix">Le 22/01/2024 à 17:33, Wilco K via
gdal-dev a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:PRAPR10MB5275537E57ED10BBEBD2829ADF752@PRAPR10MB5275.EURPRD10.PROD.OUTLOOK.COM">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style type="text/css" style="display:none;">P {margin-top:0;margin-bottom:0;}</style>
<div class="elementToProof"
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi,</div>
<div class="elementToProof"
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof"
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
these 2 lines work with GDAL 2.2.3, but not with GDAL 3.5.0
anymore.</div>
<div class="elementToProof"
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">gdal_translate
-of GTiff -a_nodata 65535 -a_srs "+proj=stere +lat_0=90
+lon_0=0 +lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0"
-a_ullr 0.0, -3649.9792, 700.000, -4414.9792
"HDF5:\"RAD_NL25_RAC_03H_202401220800.h5\"://image1/image_data"
"test-translate.tif"</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">gdalwarp
-t_srs EPSG:4326 -dstnodata 65535 -of GTiff
"test-translate.tif" "test-warp.tif"</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">Error:</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">ERROR
1: PROJ: proj_create_operations: Source and target ellipsoid
do not belong to the same celestial body<br>
ERROR 6: Cannot find coordinate operations from
`PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6378.14,298.183263207102,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Reference
meridian",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]],CONVERSION["Polar
Stereographic (variant B)",METHOD["Polar Stereographic
(variant B)",ID["EPSG",9829]],PARAMETER["Latitude of standard
parallel",60,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8832]],PARAMETER["Longitude
of
origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8833]],PARAMETER["False
easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False
northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",south,MERIDIAN[90,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],ORDER[1],LENGTHUNIT["metre",1]],AXIS["(N)",south,MERIDIAN[180,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],ORDER[2],LENGTHUNIT["metre",1]]]'
to `EPSG:4326'</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">But
when the +a and +b in the gdal_translate call are changed,
gdalwarp does work: </span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">gdal_translate
-of GTiff -a_nodata 65535 -a_srs "+proj=stere +lat_0=90
+lon_0=0 +lat_ts=60 +a=6378140 +b=6356750 +x_0=0 +y_0=0"
-a_ullr 0 -3649999.11191775 700000.903671186 -4415003.88199764
"HDF5:\"RAD_NL25_RAC_03H_202401220800.h5\"://image1/image_data"
"test-translate.tif"<br>
<br>
</span></div>
<div class="elementToProof"><span
style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">What
is the problem?</span></div>
</blockquote>
<p>The value of the +a and +b parameters must be in metres, not in
km. So you add to cheat also on the geotransform. But if datum
transformations had to be done the fact that you transform between
an ellipsoid of ~ 6000 meters to one of ~ 6000 km wouldn't work
well... It probably sort of worked because datum transformation
was skipped, but the new behaviour which checks that the shape of
the source and target ellipsoid isn't too different is definitely
saner and will avoid potential reprojection errors</p>
<span style="white-space: pre-wrap">
</span>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>