Skip to content

GUNW product - gridmapping/z-values missing as meta-data for radarGrid fields #7

@dbekaert

Description

@dbekaert

Hi,

The GUNW product works well with the netcdf reader for the data layers, for plotting in QGIS, panopoly etc.
Confirm this both for the public product and a product created internally (as of 20 Feb). Below is consistent between both the public and internal version of GUNW created:

For the radarGrid:

  • there seems to be a meta-data variable missing that identifies the projection variable in the same folder to be the grid.
  • I also cdid not see a mapping to the heightAboveEllipsoid for the z-values.

Screenshots for radarGrid illustration

Details

see the difference between unwrapped layer meta-data (has projection and x-y coordinates referenced):
Screenshot 2024-02-21 at 3 40 03 PM

parallelBaseline meta-date layer (does not have projection, also z-coordinates aka heightAboveEllipsoid are not included):
Screenshot 2024-02-21 at 3 39 35 PM

Following could be used to replicate expected behavior

Details

For Unwrapped phase this is the output as expected using the netcdf drives (given CF compliance)
gdalinfo NETCDF:"GUNW_0000424561_001001_ALOS2355610360-201223_0000426271_001001_ALOS2392870360-210901.h5"://science/LSAR/GUNW/grids/frequencyA/unwrappedInterferogram/HH/unwrappedPhase

PROJCRS["WGS 84 / UTM zone 11N",
    BASEGEOGCRS["WGS 84",
        DATUM["World Geodetic System 1984",
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4326]],
    CONVERSION["UTM zone 11N",
        METHOD["Transverse Mercator",
            ID["EPSG",9807]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",-117,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",0.9996,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",500000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    ID["EPSG",32611]]
Corner Coordinates:
Upper Left  (  365480.000, 3913615.000) (118d28'50.04"W, 35d21'24.39"N)
Lower Left  (  365480.000, 3789215.000) (118d27'38.33"W, 34d14' 7.21"N)
Upper Right (  456440.000, 3913615.000) (117d28'46.19"W, 35d21'53.61"N)
Lower Right (  456440.000, 3789215.000) (117d28'22.96"W, 34d14'35.25"N)
Center      (  410960.000, 3851415.000) (117d58'24.23"W, 34d48' 3.87"N)

For the cubes we would expect something alike recognizing projection but also incorporating the Z-dimension
Below is the current produced output:

Corner Coordinates:
Upper Left ( 362317.178, 3915829.957) 
Lower Left ( 362317.178, 3788829.957) 
Upper Right ( 472317.178, 3915829.957) 
Lower Right ( 472317.178, 3788829.957) 
Center   ( 417317.178, 3852329.957)

Separate note for GUNW product
Suspect this is already fixed (only applied to the public sample and not internally created product).

Science --> LSAR --> GUNW --> identification --> trackNumer* (trackNumber)

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions