<div dir="ltr">Dear Even et al.,<div><br></div><div>To answer my own question. Based on your comment on the units, I was able to make it. By changing the Affine transformation as well. So the gdaltransform example command becomes like this:</div><div><div style="display:inline-block;white-space:pre;font-family:"Cascadia Mono",monospace;padding:4px"><span style="color:rgb(131,148,150);background-color:rgb(255,255,255)">echo "0 -3650000" | gdaltransform -s_srs "+proj=stere +x_0=0 +y_0=0 +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378137 +b=6356752" -t_srs EPSG:4326</span></div></div><div>And gives as result</div><div><div style="display:inline-block;white-space:pre;font-family:"Cascadia Mono",monospace;padding:4px"><span style="color:rgb(131,148,150);background-color:rgb(255,255,255)">0 90 0</span></div></div><div><div style="display:inline-block;white-space:pre;padding:4px"><span style="background-color:rgb(255,255,255)"><font face="arial, sans-serif" style="" color="#000000">This is really great! Let me send a message as well to the meteorological agency so they can update their projection information.</font></span></div></div><div><div style="display:inline-block;white-space:pre;padding:4px"><span style="background-color:rgb(255,255,255)"><font face="arial, sans-serif" style="" color="#000000">Thanks again!</font></span></div></div><div><div style="display:inline-block;white-space:pre;padding:4px"><span style="background-color:rgb(255,255,255)"><font face="arial, sans-serif" style="" color="#000000">Regards,</font></span></div></div><div><div style="display:inline-block;white-space:pre;padding:4px"><span style="background-color:rgb(255,255,255)"><font face="arial, sans-serif" style="" color="#000000">Mattijn</font></span></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 3 Dec 2020 at 11:26, mattijn <<a href="mailto:mattijn@gmail.com">mattijn@gmail.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 dir="ltr">Dear Even,<div><br></div><div>Thanks for your response! It's much appreciated. I control the creation of the source file, so I can adjust the unit of the ellipsoid axis in the source. Would that help with avoiding the need of a  -ct pipeline?</div><div><br></div><div>Regards,</div><div>Mattijn</div><div> </div><div>(forgot the Replay All, to send to the full list, instead of replying to Even only.)</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 3 Dec 2020 at 11:19, Even Rouault <<a href="mailto:even.rouault@spatialys.com" target="_blank">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">Mattijn,<br>
<br>
> Question: how do I define my source raster instance, to avoid using the<br>
> advance `-ct` operator in gdalwarp in GDAL3 during reprojection?<br>
<br>
I don't think you can avoid it here. There seems to be an issue in either your <br>
source file or in the HDF5 driver regarding the unit of the ellipsoid axis. <br>
Normally the a and b parameters should be in meters, not kilometers.<br>
<br>
> But now with GDAL3 using the same command I get:<br>
> <br>
> ERROR 1: PROJ: proj_create_operations: Source and target ellipsoid do not<br>
> belong to the same celestial body<br>
<br>
yes it detects that the source and target ellipsoids are completely different, <br>
which is in most cases involve a non sensical attempt.<br>
<br>
> And using the `-ct` operator with GDAL3 I get this:<br>
> <br>
> echo "0.000 -3650.000" | gdaltransform -s_srs '+proj=stere +x_0=0 +y_0=0<br>
> +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.137 +b=6356.752' -t_srs EPSG:4326 -ct<br>
> '+proj=axisswap +order=2,1 +step' --debug on<br>
<br>
Your -ct pipeline should include the full pipeline, not just an axisswap <br>
operation. For gdaltransform, you can remove the -s_srs and -t_srs parameters <br>
when you specify -ct. For gdalwarp / ogr2ogr, you can keep the -t_srs so that <br>
the output dataset has correct CRS information.<br>
<br>
Your pipeline should be something like<br>
<br>
+proj=pipeline +step +inv +proj=stere +x_0=0 +y_0=0<br>
+lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.137 +b=6356.752 +step +proj=axisswap <br>
+order=2,1<br>
<br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div>
</blockquote></div>