Unble to import an AQUA image data (Level 2)in netcdf file format

[laurent001]

GRASS versions : 7.2.0 and 7.3.0 (current).

Failure to import a netcdf file (extension : .nc)(raster data) of an AQUA (satellite) image in the L2 level.

Seems to fail with HDF5 and GDAL according to error the log file.

Not tested L3 level satellite image (raster) data, nor L0 and L1.

GDAL, netCDF and HDF5 libraries were rebuild. Everything was OK but GRASS.

Consequence : unable to process satellite image adta with GRASS.

Encl. Log file generated by GRASS 7.2.0 as a text file.

[neteler]

From the log...:

r.import input=/Volumes/Data_1/_Image-Data/_NASA-Ocean_colour/_AQUA/_AQUA-2016-11-30-Bay_of_Biscay/A2016335134000.L2_LAC_OC.nc output=A2016335134000 -o
3	HDF5-DIAG: Error detected in HDF5 (1.8.18) thread 0:
4	  #000: /Volumes/External_1/_Other-Apps-
5	GNOME/_HDF/_Source/_V5/hdf5-1.8.18/src/H5G.c line 812 in
6	H5Gclose(): not a group
7	    major: Invalid arguments to routine
8	    minor: Inappropriate type
9	HDF5-DIAG: Error detected in HDF5 (1.8.18) thread 0:
...

r.import is a wrapper script for r.in.gdal which used GDAL.

Does "gdalinfo" accept the file?

[veroandreo]

I just checked. It seems it is not supported (gdal 2.0.2)

gdalinfo A2017012160000.L2_LAC_OC.nc
Warning 1: Variable has 0 dimension(s) - not supported.
gdalinfo failed - unable to open 'A2017012160000.L2_LAC_OC.nc'.

Global Level 3 products are recognized and imported normally, though:

r.in.gdal -o -l input=NETCDF:A20030012003008.L3m_8D_CHL_chl_ocx_4km.nc:chl_ocx output=prueba_clor_new

Maybe, SeaDAS, the software provided by NASA for oceancolor, is needed to read L2 products and export into a different format. Just guessing, no idea.

[laurent001]

Yes gdalinfo accepts and reads an AQUA image data. Encl. 1 log file (text format)

[laurent001]

Yes gdalinfo accepts and reads an AQUA image data. Encl. 1 log file (text format)

Encl. 1 log file

[mmetz]

Replying to laurent001:

Yes gdalinfo accepts and reads an AQUA image data. Encl. 1 log file (text format)

You need to select a subdataset with HDF5:file_name:subdataset, e.g. with

r.in.gdal in=HDF5:"/Volumes/Data_1/_Image-Data/_NASA-Ocean_colour/_AQUA/_2017-02-17-AQUA-Bay_of_Biscay/A2017048130000.L2_LAC_OC.nc"://geophysical_data/Kd_490

See also GDAL HDF5 and NETCDF documentation.

[neteler]

Ticket retargeted after milestone closed

[martinl]

Ticket retargeted after milestone closed

[neteler]

Closing as it will work when specifying the subdataset.