Hi Menno,
I have also tested one of your files on Windows here, with MS4W 4.0.3
(GDAL 2.4.4, NetCDF 4.7.3, HDF5 1.10.5), and it also gives me errors:
gdalinfo
--------
gdalinfo --format netcdf
GDAL_HAS_HDF4=YES
GDAL_HAS_HDF5=YES
NETCDF_CONVENTIONS=CF-1.5
NETCDF_HAS_NC2=YES
NETCDF_VERSION=4.7.3 of Jan 23 2020 15:16:48 $
gdalwarp
-------
gdalwarp -geoloc -r cubic -t_srs "+proj=merc +a=6378137 +b=6378137
+lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null
+wktext +no_defs"
NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200305T000000_20200305T001500_0001.nc":precip
"output.tiff"
output
------
Warning 1: NetCDF driver detected file type=5, but libnetcdf detected type=3
ERROR 5: OSRCalcInvFlattening(): Wrong input values
Warning 1: NetCDF driver detected file type=5, but libnetcdf detected type=3
ERROR 5: OSRCalcInvFlattening(): Wrong input values
Warning 1: NetCDF driver detected file type=5, but libnetcdf detected type=3
ERROR 5: OSRCalcInvFlattening(): Wrong input values
Creating output file that is 2790P x 4447L.
Processing
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200305T000000_20200305T001500_0001.nc:precip
[1/1] : 0Using internal nodata values (e.g. -1) for image
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200305T000000_20200305T001500_0001.nc:precip.
Copying nodata values from source
NETCDF:SEVIR_OPER_R___MSGCPP__L2__20200305T000000_20200305T001500_0001.nc:precip
to destination output.tiff.
...10...20...30...40...50...60...70...80...90...100 - done.
thoughts
--------
MS4W does include C Sharp GDAL bindings but I did not test them with
your data.
-jeff
--
Jeff McKenna
MapServer Consulting and Training Services
https://gatewaygeomatics.com/
On 2020-03-12 9:44 a.m., Menno van Scheers - HUSS wrote:
Dear list,
While in the process of upgrading from 2.2.4 to 2.4.4 we encountered an
issue with a specific NETCDF file which can no longer be warped.
Sample files can be found at : ftp://msgcpp-ogc-realtime.knmi.nl/
The file is a multiband dataset from which we want to extract the
*precip* data, it contains 2 bands for georeferencing Lat/Lon.
So in 2.2.4 we used the -geoloc warp parameter to correctly warp it to a
Webmercator or Mercator projection, this was done with the c-sharp swig
bindings.
With 2.4.4 the program just halts and does not continue, the same can be
reproduced with the commandline gdalwarp tool :
gdalwarp -geoloc -r cubic -t_srs "+proj=merc +a=6378137 +b=6378137
+lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null
+wktext +no_defs" NETCDF:"samplefile ":precip "output.tiff"
When interrupting this process HDF5 gives an error :
HDF5: infinite loop closing library
D,T,FD,P,FD,P,FD,P,E,E,SL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL,FL
Unsure if this is the actual problem or a result of the interrupt.
Best Regards,
Menno Van Scheers
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
