<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">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>