t_0=90 +lon_0=0 +lat_ts=60 +a=6378137 +b=6356752" -t_srs EPSG:4326
And gives as result
0 90 0
This is really great! Let me send a message as well to the meteorological
agency so they can update their projection information.
Thanks again!
Regards,
Mattijn
On Thu, 3 Dec 2020 at 11:26, mattijn
Dear Even,
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?
Regards,
Mattijn
(forgot the Replay All, to send to the full
tly without any advanced reprojection operators.
I've prepared all this as well in this GitHub Notebook
<https://github.com/mattijn/KNMI-radar-PyGMT/blob/main/reproject%20knmi%20radar%20using%20rasterio.ipynb>
in Python, that works with GDAL2 installed, but fails similar as above with
GDAL
Dear,
I've the same error/bug appearing. Also using the GDAL python package
(GDAL‑2.2.4‑cp35‑cp35m‑win_amd64.whl) from
https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal
Regards,
Mattijn
--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093
