[gdal-dev] Infinite loop Netcdf/HDF5 -geoloc when warping (version 2.4.4)

[Jeff McKenna]

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