[gdal-dev] MSG-CPP satellite projection to wgs84 with gdal

[Tibor Szabo]

  Hi Even,

Thank you for your quick answer.

I have tried many variants in the syntax and I got output nc file only with
the following:

+proj=geos +lon_0=0.000000 +lat_0=0 +h=35807.414063 +a=6378169 +b=6356583.8
+units=km

One error msg was also remains:

szabot at manicminer:~/gdal_knmi_test$ gdalwarp -of NetCDF -s_srs "+proj=geos
+lon_0=0.000000 +lat_0=0 +h=35807.414063 +a=6378169 +b=6356583.8 +units=km"
-t_srs EPSG:4326 -te 0 50 10 60
NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cot
TEST_MOD_COT_NL.nc
ERROR 5: OSRCalcInvFlattening(): Wrong input values
Creating output file that is 3125P x 3125L.
Processing
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc:cot
[1/1] : 0Using internal nodata values (e.g. -1) for image
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc:cot.
Copying nodata values from source
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc:cot
to destination TEST_MOD_COT_NL.nc.
ERROR 1: Too many points (529 out of 529) failed to transform, unable to
compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,3125,3125,
skipping.
...10...20...30...40...50...60...70...80...90...100 - done.

BR

Tibor


Even Rouault <even.rouault at spatialys.com> ezt írta (időpont: 2020. júl.
16., Cs, 21:18):

> On jeudi 16 juillet 2020 21:00:40 CEST Tibor Szabo wrote:
>
> > Hi All,
>
> >
>
> > I have tried to wrap a netcdf file for changing projection from satellite
>
> > view to EPSG:4326 with gdal_warp in Ubuntu and windows also with error.
>
> > Sample file:
>
> >
> https://z003.s3-us-west-2.amazonaws.com/SEVIR_OPER_R___MSGCPP__L2__20200712T
>
> > 220000_20200712T221500_0001.nc D:\meteo_model>gdalwarp -of NetCDF -s_srs
>
> > "+proj=geos +lon_0=0.000000 +lat_0=0 +h=35807.414063 +a=6378.169
>
> > +b=6356.5838" -t_srs EPSG:4326 -te 0 50 10 60
>
>
>
> +a and +b values must be expressed in metres, not kilometres.
>
>
>
> Even
>
>
>
> --
>
> Spatialys - Geospatial professional services
>
> http://www.spatialys.com
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200717/86010770/attachment.html>