GDAL_missing geo-reference when encode multibands in NetCDF from TIFF file (Know Fails)

[Bang Pham Huu]

Update 2: The real problem is from GDAL does not support to convert a multibands (> 2 bands) GeoTIFF to NetCDF. But if subset to create 1 band from this TIFF file then can convert to NetCDF with full metadata.

Example with 1 band NetCDF from 1 band GeoTIFF.

gdal_translate '/home/rasdaman/tmp/test_projection/mean_summer_airtemp
/irr_cub2_b1.tif' 
-of NetCDF output.nc

Metasdata

Driver: netCDF/Network Common Data Format
Files: output.nc
Size is 63, 36
Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],

if use 2 bands GeoTIFF to transform to NetCDF it will lost these information

gdal_translate '/home/rasdaman/test_server/NEW/rasdaman/systemtest
/testcases_services/test_wcs/testdata/irr_cube_2/AUT-PM10_2008-01-01.geo.tif'
-of NetCDF output.nc

Coordinate System is `'
Metadata:
  NC_GLOBAL#Conventions=CF-1.5
  NC_GLOBAL#GDAL=GDAL 1.11.2, released 2015/02/10
  NC_GLOBAL#GDAL_AREA_OR_POINT=Area
  NC_GLOBAL#history=Wed May 18 15:47:27 2016: GDAL CreateCopy( output.nc, ... )
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"/home/rasdaman/tmp/test_projection/mean_summer_airtemp/output.nc":Band1
  SUBDATASET_1_DESC=[36x63] Band1 (32-bit floating-point)
  SUBDATASET_2_NAME=NETCDF:"/home/rasdaman/tmp/test_projection/mean_summer_airtemp/output.nc":Band2
  SUBDATASET_2_DESC=[36x63] Band2 (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)