Petascope_missing geo-reference when encode multibands in NetCDF

[Bang Pham Huu]

When use crsTranform() in WCPS it return a coverage in NetCDF without georeference.

Update: The problem only happens if use a coverage with multibands (ranges) and encode all of these ranges in NetCDF. In this example: test_projection has 3 bands (Red, Green, Blue).

Describe Coverage metadata: https://www.diffchecker.com/9uk36tkw

If only select 1 band (Red) then it will have the full geo-referenced metadata.

WCPS query with metadata (select 1 band)

for c in (test_projection),  d in (eobstest) 
return encode(crsTransform((c.Red),
 { Lat:"http://www.opengis.net/def/crs/EPSG/0/3542", Long:"http://www.opengis.net/def/crs/EPSG/0/3542"}, {}), "netcdf")

WCPS query without metadata (multibands)

for c in (test_projection),  d in (test_projection) 
return encode(crsTransform((c + d),
 { Lat:"http://www.opengis.net/def/crs/EPSG/0/3542", Long:"http://www.opengis.net/def/crs/EPSG/0/3542"}, {}), "netcdf")

or not use crsTransform() with sample data: irr_cube_2

for c in (irr_cube_2)
return encode(c[ansi(148654)], "netcdf")

Driver: netCDF/Network Common Data Format
Files: /home/rasdaman/Downloads/ows(126)
Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#Conventions=CF-1.5
  NC_GLOBAL#GDAL=GDAL 1.11.2, released 2015/02/10
  NC_GLOBAL#history=Wed May 11 11:00:33 2016: GDAL CreateCopy( /tmp/rasdaman-uZVa17, ... )
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"/home/rasdaman/Downloads/ows(126)":Band1
  SUBDATASET_1_DESC=[260x261] Band1 (8-bit integer)
  SUBDATASET_2_NAME=NETCDF:"/home/rasdaman/Downloads/ows(126)":Band2
  SUBDATASET_2_DESC=[260x261] Band2 (8-bit integer)
  SUBDATASET_3_NAME=NETCDF:"/home/rasdaman/Downloads/ows(126)":Band3
  SUBDATASET_3_DESC=[260x261] Band3 (8-bit integer)
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)