4. L2a Products Format Specification

4.1 Main Product Header

The Main Product Header (MPH) is an XML file present in all the BIOMASS L2 products.

BIOMASS L2 products are folder based and the MPH is top level leaf of the product tree.

The MPH main scope is to provide all the information related to the:

• Genesis of the product;

• Identification of the product for example for cataloguing purposes;

• Product consistency (filesize).

It must be noted that, while the footprint is inherited from lower level products (thus vertices represent the extent of the acquired data), DGG tiles (tileId list) are related to the full extent of L2a product (that is geocoded).

For the detailed description of the MPH content and format refer to [AD3]. Here just a

schematic representation is reported, , where this notation is applicable:

This is a container of other elements

This element or block of elements can occur more than one time (i.e. it is a list)

This is a container of actual metadata

This is a mandatory element for the standard which is left empty by this specification

This element or block of elements is optional

bio:EarthObservation Root element of metadata representation tailored for BIOMASS SAR products

om:phenomenonTime Time coverage of the product

gml:TimePeriod Time coverage of the product as a period

gml:beginPosition Start date and time of the product

gml:endPosition Stop date and time of the product

om:resultTime Time when the sensing data has become available

gml:TimeInstant Identifiable position in time

gml:timePosition Date and time of data availability

om:validTime Time validity of the product

gml:TimePeriod Time validity of the product as a period

gml:beginPosition Validity start date and time of the product

gml:endPosition Validity stop date and time of the product

om:procedure Procedure used to sense the data

eop:EarthObservationEquipment Equipment used to sense the data

eop:platform Platform description

eop:Platform Nested element for platform description

eop:shortName Satellite name

eop:instrument Instrument description

eop:Instrument Nested element for instrument description

eop:shortName Sensor name

eop:sensor Sensor description

eop:Sensor Nested element for sensor description

eop:sensorType Sensor type

eop:operationalMode Sensor operational mode

eop:swathIdentifier Swath identifier

eop:acquisitionParameters Acquisition parameters

bio:Acquisition Acquisition description

eop:orbitDirection Orbit direction

eop:wrsLongitudeGrid Track number

eop:wrsLatitudeGrid Frame number

sar:antennaLookDirection SAR antenna look direction

bio:missionPhase BIOMASS mission phase

bio:globalCoverageID BIOMASS global coverage identifier

bio:majorCycleID BIOMASS major cycle identifier

bio:tileID List of DGG geographic cell (tile) identifiers

bio:basinID Retrieval basin identifier

om:observedProperty Observed property (empty)

om:featureOfInterest Observed area

eop:Footprint Footprint description

eop:multiExtentOf Footprint representation structure, coordinates in posList

gml:MultiSurface

gml:surfaceMember

gml:Polygon

gml:exterior

gml:LinearRing

gml:posList Footprint points in NE,SE,SW,NW,NE order

eop:centerOf Acquisition centre representation structure

gml:Point

gml:pos Coordinates of the centre of the acquisition

om:result Observation result

eop:EarthObservationResult

eop:product

bio:ProductInformation

eop:fileName

ows:ServiceReference Product filename

ows:RequestMessage download request (empty)

eop:version Product baseline

eop:product (list)

bio:ProductInformation

eop:fileName

ows:ServiceReference Product filename (internal files)

ows:RequestMessage download request (empty)

eop:size Product size

bio:rds product schema

eop:metaDataProperty Observation metadata

bio:EarthObservationMetaData

eop:identifier Product identifier

eop:doi Digital Object Identifier

eop:acquisitionType Type of acquisition

eop:productType Product type

eop:status Status of the product

eop:processing Data processing information

bio:ProcessingInformation

eop:processingCenter Processing centre code

eop:processingDate Processing date and time

eop:processorName Processor name

eop:processorVersion Processor version

eop:processingLevel Processing level

eop:auxiliaryDataSetFileName Auxiliary DS filename

eop:processingMode File Class

bio:sourceProduct BIOMASS source product identifier

bio:refDoc Reference document

Main product header L2a (metadata schematic representation)

4.2 Measurement Data Sets (MDS)

The MDS of the L2a product is constituted by two images:

• the biophysical measurement consisting of the Forest Height binary raster, representing Top Canopy Height (TCH);

• a layer of quality, consisting of a percentage bias value for each of the Forest Height image pixel, indicating the inversion performance.

The two quantities are organised in separated Cloud Optimized GeoTIFF (COG) files [RD1] and coded using a 4-byte float per pixel (float32).

A Cloud Optimized GeoTIFF (COG) is a regular GeoTIFF file, aimed at being hosted on a

HTTP file server, with an internal organization that enables more efficient workflows on the

cloud. It does this by leveraging the ability of clients issuing HTTP GET range requests to ask for just the parts of a file they need.

The COG format is natively supported by GDAL library [RD2] [RD3] and by libraries and tools based on GDAL library (Dask and Rasterio packages for Python, …).

The MDS images ground coverage is corresponding to the geocoded spatial stack coverage of inputs L1c used to compute L2a product: the L2a MDS images are geocoded to the minimum required latitude-longitude coverage to accommodate for the product on ground, following the same sampling requirements specified in [RD9] ensuring a ground sampling always better than the resolution.

Compression strategy

The MDS images are composed by real values: it is possible to introduce with benefit a strategy for data compression, making use of native LibTIFF libraries capabilities [RD4] adopting a GeoTIFF LERC+ZSTD compression, where:

• LERC compression algorithm can be configured to decide how lossy the compression will be, through the “MAX_Z_ERROR” option (loss-less only if set to zero);

• ZSTD is a loss-less compression algorithm with a configurable compression level from 0 to 9, where zero means no compression.

Both set by configuration. More details in sec. 4.3.1.4.

Metadata

The TIFF and GeoTIFF tags used to annotate the MDS files are described in the following table. For the complete description of each tag refer to official TIFF format documentation [RD5] and official GeoTIFF documentation [RD11].

Table 2 — MDS – TIFF and GeoTIFF tags

Tag Name	Code	Value	Description
ImageWidth	256	Variable	Number of columns in the image, i.e. the number of pixels per row.
ImageLength	257	Variable	Number of rows of pixels in the image.
BitsPerSample	258	32	Number of bits per component.
Compression	259	34887, 34926	Compression scheme used on the image data (LERC and ZSTD)
PhotometricInterpretation	262	1	Colour space of the image. Always set to 1 (minimum value is black).
ImageDescription	270	Variable	A string that describes the image. Valid Strings are composed by prefix: “BIOMASS L2a FP_FH__L2A: ” or “BIOMASS L2b FP_FH__L2B: ” and: “Forest Height [m]” “Forest height quality” “Forest Mask” (L2b only) “Heat Map” (L2b only) “acquisition id image” (L2b only)
SamplesPerPixel	277	1	The number of components per pixel. Always set to 1, since each product is stored in a different file.
PlanarConfiguration	284	1	How the components of each pixel are stored. Always set to 1 (contiguous).
Software	305	Variable	Software name and version that created the image.
DateTime	306	Variable	Date and time of image creation in the format: YYYY:MM:DD HH:MM:SS
TileWidth	322	Variable	The tile width in pixels. This is the number of columns in each tile.
TileLength	323	Variable	The tile length (height) in pixels. This is the number of rows in each tile.
TileOffsets	324	Variable	For each tile, the byte offset of that tile, as compressed and stored on disk.
TileByteCounts	325	Variable	For each tile, the number of (compressed) bytes in that tile.
SampleFormat	339	3 or 1	Specifies how to interpret each data sample in a pixel. Set to 3 (IEEE floating point data)
ModelPixelScaleTag	33550	Variable	List of doubles, (ScaleX, ScaleY, ScaleZ) where ScaleX and ScaleY give the horizontal and vertical spacing of raster pixels. The ScaleZ is zero.
ModelTiePointTag	33922	N * (pixel, line, 0, lat, long, 0)	Maps image pixels to geographic coordinates, where N is the number of tie points in the scene.
GeoKeyDirectoryTag	34735	Variable	Contains the value of all GeoTIFF keys of type SHORT and the location of the values for GeoTIFF keys of other types.
GeoDoubleParamsTag	34736	Variable	Contains the value of all GeoTIFF keys of type DOUBLE.
GeoAsciiParamsTag	34737	Variable	Contains the value of all GeoTIFF keys of type ASCII.
GDAL_METADATA	42112	Variable	Contains GDAL metadata (see table below)
GDAL_NODATA	42113	“-9999” or “255”	Pixel value used in case of invalid data. ASCII string encoded value. (See 3)
LercParameters	50674	Variable	LERC compression parameters: version, additional compression. e.g. 4,2 for LERC v4 with ZSTD

The GDAL_METADATA tag content is listed in the following table

Table 3 — MDS – GDAL_METADATA tiff tag content

Tag Name	Value	Description
basinID	list of strings	List of DGG basin identifiers (for querying purposes). More than one ID in case of tile covering more than one geographic location. Private Tag code.
MAX_Z_ERROR	float	Absolute maximum error present in the image data, due to LERC compression. Private Tag code.
Swath	string	Image swath. Private Tag code (L2A only)
tileID	list of strings	List of DGG geographic cell (tile) identifiers (for querying purposes). Private Tag code.

Overviews

The exploitation of GDAL library allows to easily add to COG file one or more overviews of the L2 image.

Overviews are sub-sampled versions of an L2 image. Adding one or more overviews at different resolution levels allows to use pyramidal viewers directly on L2 image file, without the need of transcoding operations or of further processing.

The following (minimum) number of overviews is foreseen:

• multi-look factors (latitude x longitude) of 2x2.

• multi-look factors (latitude x longitude) of 4x4.

4.3 Annotation data sets (ADS)

4.3.1 Main ADS

This section describes in details the content of main ADS for L2a products.

An overview is provided by figure and table below.

Fig. 3 — Main L2a ADS

Table 4 — Main L2a ADS

Name	Description	Data Type	Cardinality
product	Product L2a DSR. This DSR contains the L2a product information.	productL2aType	1
rasterImage	Raster image DSR. This DSR contains all the necessary information to exploit the L2a product raster images.	rasterImageType	1
InputInformation	Input Information DSR. This DSR contains the necessary information to identify the input data set to the L2a processing.	InputInformationL2aType	1
processingParameters	Processing parameters DSR. This DSR contains the description of L2a processing parameters.	processingParametersL2aType	1
annotationLUT	Annotation LUT DSR. This DSR contains the list of Look-Up Tables (LUTs) complementing product main annotations.	layerListType	1

The content of each DSR is described making use of tables as the one here above. All the tables share the same structure and are filled following the same conventions, i.e.:

• Name: defines the name of the element or attribute;

• Description: describes the purpose of the element or attribute including the units (if applicable). The absence of a unit designation for string data types implies a string literal and for numerical data types (integers, floats, etc.…) the absence of a unit designation implies an absolute number;

• Data Type: defines the data type of the element or attribute;

• Cardinality: for elements, this column defines the number of occurrences of the element in the form minOccurs..maxOccurs; for attributes, “optional” means the attribute may or may not be present and “required” means the presence of the attribute is mandatory.

In general, elements are displayed in normal font, attributes are displayed in italics.

Data types requiring a dedicated definition table are underlined, for all the others the common reference is sec. 7.1.1.

4.3.1.1 Product L2a DSR

This DSR contains the parameters to describe the L2a products.

An overview is provided by figure and table below.

Fig. 4 — Product L2a DSR (productL2aType)

Table 5 — Product L2a DSR (productL2aType)

Name	Description	Data Type	Cardinality
mission	Mission (BIOMASS).	missionType	1
tileID	List of DGG geographic cell (tile) identifiers.	stringListType	1
basinID	List of DGG Basin identifiers.	stringListType	1
productType	Product type (“FH_L2A”).	productType	1
startTime	Start time of the product image [UTC].	timeType	1
stopTime	Stop time of the product image [UTC].	timeType	1
radarCarrierFrequency	Radar carrier frequency [Hz].	doubleWithUnit	1
missionPhaseID	Mission phase identifier.	missionPhaseIDType	1
sensorMode	Sensor mode (always Measurement).	sensorModeType	1
globalCoverageID	Global coverage identifier.	xsd:unsignedShort	1
swath	Swath (S1, S2, S3).	swathType	1
majorCycleID	Major cycle identifier.	xsd:unsignedShort	1
absoluteOrbitNumber	List of absolute orbit numbers at start time (one for each input L1c).	integerListType	1
relativeOrbitNumber	Relative orbit number (track) at start time	xsd:unsignedShort	1
orbitPass	Orbit pass (Ascending, Descending).	orbitPassType	1
dataTakeID	List of data take identifiers (one for each input L1c).	integerListType	1
frame	Frame identifier.	xsd:unsignedShort	1
platformHeading	Platform heading relative to North [deg].	doubleWithUnit	1
forestCoveragePercentage	Forest coverage of the input STA products, compute as a percentage of all the pixels in the STA footprint.	double	1

Table 6 — stringListType

Name	Description	Data Type	Cardinality
ID	Tile identifier string	string	1 … unbounded

Table 7 — integerListType

Name	Description	Data Type	Cardinality
val	Integer value	uint32	1 … unbounded

4.3.1.2 Raster Image DSR

This DSR contains all the necessary information to exploit the raster images of the product. An overview is provided by figure and table below.

Fig. 5 — Raster image DSR (rasterImageType)

Table 8 — Raster image DSR (rasterImageType

Name	Description	Data Type	Cardinality
footprint	Image footprint, expressed as a list of latitude, longitude tuples in NE, SE, SW, NW order [deg].	floatArrayWithUnits	1
firstLatitudeValue	First latitude value of the image [deg].	floatWithUnit	1
firstLongitudeValue	First longitude value of the image [deg].	floatWithUnit	1
latitudeSpacing	Latitude spacing between samples [deg].	floatWithUnit	1
longitudeSpacing	Longitude spacing between samples [deg].	floatWithUnit	1
numberOfSamples	Total number of samples in the image.	xsd:unsignedInt	1
numberOfLines	Total number of lines in the image.	xsd:unsignedInt	1
projection	Projection of the image: “Latitude longitude based on DGG”.	projectionType	1
datum	Datum used during processing.	datumType	1
pixelRepresentation	Representation of the image pixels within the image MDS.	pixelRepresentationChoiceType	1
pixelType	Data type of output pixels within the image MDS.	pixelTypeChoiceType	1
noDataValue	Pixel value in case of invalid data. The value is always set to -9999.0 float32 or 255 ubyte, depending on data type, see sec. 3, 5.	noDataValueChoiceType	1

Table 9 — pixelRepresentationChoiceType(L2a/L2b FH described here)

Name	Description	Data Type	Cardinality
FH	Representation of the image pixels within the image MDS: Forest Height [m]	pixelRepresentationType	1
Quality	Representation of the image pixels within the image MDS: Forest height quality	pixelRepresentationType	1
FH_heatMap	Representation of the image pixels within the image MDS: Heat Map (Only L2b)	pixelRepresentationType	0,1
BPS_FNF	Representation of the image pixels within the image MDS: Forest Mask (Only L2b)	pixelRepresentationType	0,1
Acquisition_Id_Image	Representation of the image pixels within the image MDS: acquisition_id_image (Only L2b)	pixelRepresentationType	0,1

Table 10 — pixelTypeChoiceType

Name	Description	Data Type	Cardinality
floatPixelType	Data type of output pixels within the image MDS: 32 bit Float	pixelTypeType	1
intPixelType	Integer data type of output pixels within the image MDS: 8 bit Unsigned Integer (L2b only)	pixelTypeType	0,1

Table 11 — noDataValueChoiceType

Name	Description	Data Type	Cardinality
floatNoDataValue	Pixel value in case of invalid data. The value is always set to -9999.0	xsd:float	1
intNoDataValue	Pixel value in case of invalid data. The value is always set to +255 (L2b only)	xsd:unsignedByte	0,1

4.3.1.3 Input Information DSR

This DSR contains information about the L1c acquisitions in input to the L2a processor.

An overview is provided by figures and tables below.

Fig. 6 — Input Information DSR (InputInformationL2aType)

Fig. 7 — Input Information DSR, acquisition list detail (acquisitionListType)

Fig. 8 — Input Information DSR, acquisition detail (acquisitionType)

Fig. 9 — Input information DSR, staQualityType

Fig. 10 — Input information DSR, staQualityParametersListType

Fig. 11 — Input information DSR, staQualityParametersType

Table 12 — Input Information DSR (InputInformationL2aType)

Name	Description	Data Type	Cardinality
productType	Input product type: always “STA” (L1c STA)	productType	1
overallProductsQualityIndex	Quality index based on the “overallProductQualityIndex” of the single STA input data, giving a qualitative overall idea of the input stack. Equal to “0” if all STA feature “overallProductQualityIndex”=False and to “1” if at least one is degraded.	xsd:unsignedByte	1
nominalStack	True if it is a nominal stack (3 STA products for INT phase or 7 STA products for TOM phase), False otherwise	bool	1
polarisationList	List of polarisations composing input L1c data: for each polarisation it is also specified, with an attribute, the decorrelation factors due to RFI strength.	polarisationListType	1
projection	Projection of the input acquisitions, always “Slant Range”.	projectionType	1
footprint	Image footprint of all the acquisitions composing the stack, expressed as a list of latitude, longitude tuples in NE, SE, SW, NW order [deg].	floatArrayWithUnits	1
verticalWavenumbers	Minimum and maximum values of the vertical wavenumbers [rad/m] in the read and used Input STA products.	minMaxTypeWithUnit	1
heightOfAmbiguity	Minimum and maximum values of the height of ambiguity (HoA) [m], computed from the vertical wavenumbers in the read and used Input STA products.	minMaxTypeWithUnit	1
acquisitionsList	Info of each acquisition composing the stack, in input to L2a.	acquisitionListType	1

The acquisitionListType is composed by 2 to 8 entries, where each entry is an L1c acquisition folder name string, as defined in [RD10]

Table 13 — Input information DSR (acquisitionListType)

Name	Description	Data Type	Cardinality
count	Number of acquisitions.	unsignedInt	required
acquisition	Info about acquisition of the stack.	acquisitionType	2…8

Table 14 — Input information DSR (acquisitionType)

Name	Description	Data Type	Cardinality
referenceImage	Flag to indicate if this is the reference image or not.	bool	1
averageWavenumber	Optional Average Wavenumber, it is not used in FH (only in L2A FD)	float	0
FolderName	Folder name which univocally identifies an acquisition of the stack.	xsd:string	1
staQuality	Folder name which univocally identifies an acquisition of the stack, as defined in [RD10].	staQualityType	1

Table 15 — Input information DSR (staQualityType)

Name	Description	Data Type	Cardinality
overallProductQualityIndex	Overall product quality index. This annotation is calculated based on specific quality parameters and gives an overall quality value to the product. Equal to 0 for valid products and to 1 for invalid ones.	xsd:unsignedInt	1
staQualityParametersList	Quality parameters list. The list contains an entry for each output product’s polarisation.	staQualityParametersListType	1

Table 16 — STA quality DSR (staQualityParametersListType)

Name	Description	Data Type	Cardinality
count		xsd:unsignedInt	required
staQualityParameters	Quality parameters for the current polarisation.	staQualityParametersType	1 … 4

Table 17 — STA quality DSR (staQualityParametersType)

Name	Description	Data Type	Cardinality
polarisation		polarisationType	required
invalidL1aDataSamples	Number of invalid L1a data samples (according to their mean and standard deviation), expressed as a normalized fraction of the total number of samples. The value is between 0 and 1.	xsd:float	1
rfiDecorrelation	Decorrelation factor due to RFI strength. The value is between 0 and 1.	xsd:float	1
rfiDecorrelationThreshold	Maximum decorrelation factor due to RFI strength admissible to use the image as the primary image. The value is between 0 and 1 and copied from the AUX-PPS.	xsd:float	1
faradayDecorrelation	Decorrelation factor due to Faraday residual. The value is between 0 and 1.	xsd:float	1
faradayDecorrelationThreshold	Maximum decorrelation factor due to Faraday residual admissible to use the image as the primary image. The value is between 0 and 1 and copied from AUX-PPS.	xsd:float	1
invalidResidualShiftsFraction	Number of invalid residual coregistration shifts, expressed as a normalized fraction of the total number of residual coregistration shifts. The value is between 0 and 1.	xsd:float	1
residualShiftsQualityThreshold	Threshold on residual coregistration shifts quality. The value is between 0 and 1 and copied from AUX-PPS. Residual shifts with a quality lower than this threshold are not used for processing.	xsd:float	1
invalidGroundPhasesScreenEstimatesFraction	Number of invalid ground phase screen estimates, expressed as a normalized fraction of the total number of ground phase screen estimates. The value is between 0 and 1.	xsd:float	1
groundPhasesScreenQualityThreshold	Threshold on ground phase screen estimation quality. The value is between 0 and 1 and copied from the AUX-PPS. Estimates with a quality lower than this threshold are not corrected. Used only in case skpPhaseCorrectionFlag is set to True.	xsd:float	1
skpDecompositionIndex	Error code for the SKP decomposition. Positive if the SKP decomposition hit a nonblocking contingency case (e.g. SVD decomposition failure etc.) and the exported SKP-related LUTs contain values obtained from a contingency handling procedure. 0 otherwise.	xsd:unsignedInt	1

Table 18 — Input information DSR (polarisationListType)

Name	Description	Data Type	Cardinality
count	Number of polarization images.	xsd:unsignedInt	required
polarisation	Polarisation (HH, VH, VV).	polarisationType	3

4.3.1.4 Processing parameters L2a DSR

This DSR contains the exhaustive list of L2a Forest Height inversion processing parameters applied. An overview is provided by figure and table below.

Fig. 12 — Processing parameters L2a DSR (processingparametersL2aType)

Table 19 — Processing parameters L2a DSR (processingparametersL2aType)

Name	Description	Data Type	Cardinality
processorVersion	Version of the processor used to generate the product.	xsd:string	1
productGenerationTime	Product generation time [UTC].	timeType	1
generalConfigurationParameters	General processing flags (not specific to the Forest Height processing).	generalConfigurationParametersType	1
verticalReflectivityOption	Specify which vertical reflectivity profile has been used: default modelled profile or tomographic profile (out of BPS scope). See also verticalReflectivityDefaultProfile.	verticalProfileOptionType	1
verticalReflectivityDefaultProfile	Used if verticalReflectivityOption is set to “default profile”, in this case the profile stored to xml is used: a modeled tomographic volume-only reflectivity profile normalized to unit height; deterministic and fixed for all the globe [AD9].	verticalReflectivityProfileType	1
modelInversion	Model inversion algorithm used between single or dual (baseline).	modelInversionType	1
spectralDecorrelationCompensationFlag	True if range spectral decorrelation compensation has been applied, False otherwise.	bool	1
SNRDecorrelationCompensation	True if SNR decorrelation compensation has been applied, False otherwise.	bool	1
correctTerrainSlopesFlag	True if terrain slope correction has been applied, False otherwise.	bool	1
normalisedHeightEstimationRange	Validity range of heights used for the canopy height estimation process, normalized.	minMaxType	1
normalisedWavenumberEstimationRange	Validity range of wavenumbers used for the canopy height estimation process, normalized.	minMaxNumType	1
groundToVolumeRatioRange	Validity ground to volume ratio range used for the canopy height estimation process.	minMaxNumType	1
temporalDecorrelationEstimationRange	Validity temporal decorrelation range used for the canopy height estimation process.	minMaxNumType	1
temporalDecorrelationGroundToVolumeRatio	Ratio of temporal decorrelation between ground and volume (0.0 means no temporal decorrelation for ground, while 1.0 means ground and volume are equally impacted by temporal decorrelation)	xsd:float	1
residualDecorrelation	Residual non-volumetric decorrelation used in error model computation.	xsd:float	1
productResolution	Value in [m] used as the resolution on ground range map and also to perform the coherence averaging in radar coordinates.	floatWithUnit	1
uncertaintyValidvaluesLimits	Uncertainty valid values limits applied, as a percentage: estimates out of this limits have been discarded and set to no data value.	minMaxType	1
verticalWavenumberValidvaluesLimits	Vertical wavenumber valid values limits in [rad/m] applied: estimates out of this limits have been discarded and set to no data value.	minMaxType	1
lowerHeightLimit	FH estimates lower this limit [m] were discarded and set to no data value	floatWithUnit	1
upsamplingFactor	Upsampling factor used for coherence.	xsd:unsignedByte	1
compressionOptions	Configurable compression options for all the L2 MDS COG and ADS NetCDF LUT variables (see 4.3.2 for LUT “ZLIB” compression details).	compressionOptionsL2Type	1

Table 20 — Processing parameters L2a DSR (verticalReflectivityProfileType)

Name	Description	Data Type	Cardinality
count	Number of values composing the function.	xsd:unsignedInt	required
val	N floating point values which compose the default modelled tomographic volume-only reflectivity profile function.	xsd:float	1…N

Fig. 13 — Processing parameters L2a DSR (generalConfigurationParametersType)

Table 21 — Processing parameters L2a DSR (generalConfigurationParametersType)

Name	Description	Data Type	Cardinality
applyCalibrationScreen	Phase calibration performed: “none”: no phase screen has been applied “geometry”: only flattening phase screen has been applied (i.e., as computed from acquisition geometry) “skp”: complete phase screen has been applied (default)	calibrationScreenType	1
forestCoverageThreshold	Minimum percentage forest coverage in L2a product footprint, used to trigger L2a processing.	xsd:float	1
forestMaskInterpolationThreshold	This parameter is a threshold which has been used to fix rounding of pixels with decimal values originated from binary FNF interpolation onto L2a grid. This created a safety buffer around forest border.	xsd:float	1
subsettingRule	In case of more than three acquisitions the user selected subsetting rule is here annotated (always “geometry” for FH L2a processor).	subsettingRuleType	1
polarisationCombinationMethod	Polarisations combination method: “HV” or “VH”, if just one of the two is selected (in addition to HH and VV ones); “Average”, if the average of HV and VH is computed and used; “None” if no combination is performed (all the four polarizations are used).	polarisationCombinationMethodType	1

Fig. 14 — Processing parameters L2a DSR (compressionOptionsL2Type)

Name	Description	Data Type	Cardinality
MDS àFH à compressionFactor MDS àQualityà compressionFactor	ZSTD algorithm compression factor for both FH and quality images MDS. From 0 to 9, where zero means no compression.	xsd:unsignedByte	1
MDS à FH à MAX_Z_ERROR MDS à Qualityà MAX_Z_ERROR	For both FH and quality images MDS, define exactly how lossy the LERC compression algorithm is allowed to be, specifying the absolute maximum error admitted. Zero means loss-less compression.	xsd:float	1
ADS àFNF à compressionFactor	ZLIB algorithm compression factor for the FNF ADS. From 0 to 9, where zero means no compression.	xsd:unsignedByte	1
MDS_blockSize	MDS COG blocking algorithm size. Same value is used for both data array dimensions.	xsd:unsignedByte	1
ADS_blockSize	ADS NetCDF chunking algorithm size. Same value is used for both data array dimensions.	xsd:unsignedByte	1

Processing parameters L2a DSR (compressionOptionsL2aType)

4.3.1.5 Annotation LUT DSR

This DSR contains the list of Look-Up Tables (LUTs) complementing product main annotations. An overview is provided by figure and table below.

Fig. 15 — Annotation LUT DSR (layerListType)

Table 22 — Annotation LUT DSR (layerListType)

Name	Description	Data Type	Cardinality
count	Number of layers.	xsd:unsignedInt	required
layer	LUT layer content.	layerType	1 … unbounded

4.3.2 Look-Up Tables (LUT) ADS

This section describes in details the content of Look-Up Tables (LUT) ADS.

This ADS complements the Main ADS one, including those annotations that can be represented as 2D maps of values. This approach has been chosen in order to avoid worsening the Main ADS readability and handling, and also to provide an easier access to such information.

The format adopted for LUT ADS is NetCDF one [RD6] and in particular the used implementation is the Python NetCDF4 [RD7]. Thanks to the flexibility provided by this format, all the components of this ADS are stored in the same file.

The LUT ADS provides the Forest/Non-Forest mask (generated outside of BPS [AD6] [AD7] [AD8] [AD9]).

Fig. 16 provides a schematic representation of the LUT ADS structure. Note that an XSD diagram has been used for this purpose, even if there is no link between the two formats. The ADS is organised as follows:

• a series of root level Attributes, containing all those global parameters referring to the product as a whole;

• a Dimensions group, where names and lengths of the dimensions used when creating data arrays are here defined;

• a series of Coordinate variables, special variables to store the axis data values, for each dimension in the Dimensions group;

• a series of components (groups of variables), one for each LUT stored in the ADS. Each component is composed by a group of variables which can be divided in:

  • the data arrays, storing the annotation Look-Up Tables themselves. A variable represents an array of values of the same data type and it is defined over two or more dimensions from the Dimensions group;

  • ancillary scalar variables containing other useful metadata to complete the LUT definition.

Fig. 16 — Look-Up Tables (LUT) ADS

Compression strategy

The NetCDF data arrays variable is composed by binary values, it is possible to introduce with benefit a strategy for data compression making use of native Python NetCDF4 capabilities [RD7], using the ZLIB (loss-less) compression algorithm, with a compression level from 0 to 9, where zero means no compression, set by configuration, more details in sec. 4.3.1.4.

4.3.2.1 Look-Up tables detailed content

The content of L2a LUT ADS is reported hereafter, where:

• the series of Attributes contain the same list of parameters described in sec. 4.3.1.1 for the Product L2a DSR (for all the details about these parameters refer to this section);

• the Dimensions group is composed by Latitude and Longitude, as the data array variable is defined in the Latitude Longitude based on DGG projection;

• there is only one data array variable, which is the Forest/Non-Forest mask;

• this notation is applicable:

  • this is the attribute or variable name

  • this is the actual content of constant values over the mission

  • this is a description of values content which depends on the processing

LUT ADS Attributes:

mission: “BIOMASS”

tileID: [string] json encoded list of DGG geographical tile identifiers (NxxEyyy)

basinID: [string] json encoded list of DGG Basin identifiers

productType: FP_FH__L2A

startTime: [str: YYYY-MM-DD hh:mm:ss.ppp] Start time of the product image, UTC

stopTime: [str: YYYY-MM-DD hh:mm:ss.ppp] Stop time of the product image, UTC

radarCarrierFrequency: [float32] Radar carrier frequency [Hz]

missionPhaseID: [str: COM, TOM or INT] Mission phase identifier

sensorMode: “Measurement”

globalCoverageID: [uint16] Global coverage identifier

swath: [str: S1, S2 or S3] BIOMASS swath identifier

majorCycleID: [uint16] Major cycle identifier

absoluteOrbitNumber: [list of uint32] List of absolute orbit numbers at start time (one for each input L1c)

relativeOrbitNumber: [uint16] Relative orbit number (track) at start time

orbitPass: [str: Ascending or Descending] Orbit pass

dataTakeID: [string] json encoded list of data take identifiers (one for each input L1c)

frame: [uint16] Frame identifier

platform_heading: [float32] Platform heading relative to North, deg

forest_coverage_percentage: [float32] Forest coverage of the input STA products, compute as a percentage of all the pixels in the STA footprint

LUT ADS Dimensions / Coordinate variables:

Name: Latitude

size: [int] size

type: dtype(float32)

units: deg

description: FNF LUT latitude axis

Name: Longitude

size: [int] size

type: dtype(float32)

units: deg

description: FNF LUT longitude axis

LUT ADS group of variables:

Group Name: FNF

FNF group contents:

Name: FNF

dimensions: (Longitude, Latitude)

shape: ([int] x shape, [int] y shape)

type: dtype(ubyte)

description: Forest/Non-Forest mask : this is an external image, defined outside of BPS from aggregated C3S LCM, representing the initial two-dimensional mask, here defined in the same latitude-longitude based on DGG projection used for L2a images and cropped to have a coverage containing the L2a product image boundaries. It defines, for each geographical pixel, the initial presence or absence of forest.

Name: firstSample

shape: (1, 1)

type: dtype(uint32)

description: Product image MDS sample the first LUT sample corresponds to

Name: firstLine

shape: (1, 1)

type: dtype(uint32)

description: Product image MDS line the first LUT line corresponds to

Name: samplesInterval

shape: (1, 1)

type: dtype(uint32)

description: Number of product image MDS samples one LUT sample corresponds to

Name: linesInterval

shape: (1, 1)

type: dtype(uint32)

description: Number of product image MDS lines one LUT line corresponds to

Name: pixelType

shape: (1, 1)

type: <class str>

description: Data type of output LUT pixels

Name: noDataValue

shape: (1, 1)

type: dtype(ubyte)

description: Pixel value in case of invalid data. The value is always set to 255. (See 3)

Name: projection

shape: (1, 1)

type: <class str>

description: LUT projection, DGG, as described in [RD8]

Name: coordinateReferenceSystem

shape: (1, 1)

type: <class str>

description: Coordinate reference systems, as per WKT representation [RD7]

Name: datum, geodeticReferenceFrame

shape: (1, 1)

type: <class str>

description: Geodetic reference frame

4.3.3 Quick-look ADS

The Quick-look ADS is composed by lower resolution versions of the MDS Forest Height and quality layer images; they are stored in PNG file format and they are always included with any BIOMASS L2 product.

The quantities are composed by real values, so the images are generated without any other processing than averaging and decimation.

4.3.4 Overlay ADS

The Overlay ADS is a Keyhole Markup Language (KML) file that describes the product coverage area and is suitable for viewing the Quick-look ADS in any application that supports KML format (like, e.g., Google Earth or GIS applications).

Its structure is described in the following figure and tables.

Fig. 17 — Overlay ADS

Table 23 — Overlay ADS

Name	Description	Data Type	Cardinality
Document	Document container for KML components.	DocumentType	1

Table 24 — Overlay ADS (DocumentType)

Name	Description	Data Type	Cardinality
name	Name of the document.	xsd:string	1
description	Description of the document.	xsd:string	1
GroundOverlay	Contains the parameters required to specify the footprint of the image and overlay the quicklook image on a map.	GroundOverlayType	1
Placemark	Contains the metadata associated to the overlay.	PlacemarkType	1

Table 25 — Overlay ADS (GroundOverlayType)

Name	Description	Data Type	Cardinality
name	Name of the ground overlay.	xsd:string	1
visibility	Specifies whether the feature is drawn in the 3D viewer when it is initially loaded (0: not visible, 1: visible).	xsd:unsignedInt	1
TimeStamp	Structure containing the time stamp of the overlay that can be used for dynamic visualizations over time.	TimeStampType	1
Icon	Structure describing the image file used on the map overlay.	IconType	1
LatLonQuad	Structure containing the latitude and longitude coordinates used to position the image overlay on the map.	LatLonQuadType	1

Table 26 — Overlay ADS (TimeStampType)

Name	Description	Data Type	Cardinality
when	Overlay time stamp in ISO 8601 format (YYYY-MM-DDThh:mm:ss).	xsd:string	1

Table 27 — Overlay ADS (IconType)

Name	Description	Data Type	Cardinality
href	A local file specification or URL used to load the desired image and overlay it on the map.	xsd:string	1

Table 28 — Overlay ADS (LatLonQuadType)

Name	Description	Data Type	Cardinality
coordinates	A string of 5 lon,lat coordinate pairs which describe the corners of the image. The string is of the form: lon,lat lon,lat lon,lat lon,lat lon,lat The coordinates must appear in the following order: last line first pixel, last line last pixel, first line last pixel, first line first pixel, last line first pixel	xsd:string	1

Table 29 — Overlay ADS (PlacemarkType)

Name	Description	Data Type	Cardinality
Name	Name of the placemark.	xsd:string	1
visibility	Specifies whether the feature is drawn in the 3D viewer when it is initially loaded (0: not visible, 1: visible).	xsd:unsignedInt	1
TimeStamp	Structure containing the time stamp of the placemark that can be used for dynamic visualizations over time	TimeStampType	1
Style	Structure describing the stype used to display the placemark.	StyleType	1
ExtendedData	Structure containing the metadata associated to the overlay.	ExtendedDataType	1
Polygon	Structure containing the latitude and longitude coordinates used to position the image overlay on the map.	PolygonType	1

Table 30 — Overlay ADS (StyleType)

Name	Description	Data Type	Cardinality
LyneStyle	Line style.	LyneStyleType	1
PolyStyle	Polygon style.	PolyStyleType	1

Table 31 — Overlay ADS (LyneStyleType)

Name	Description	Data Type	Cardinality
color	Line color in 32 bit hex format (AABBGGRR order).	xsd:string	1
width	Line width.	xsd:unsignedInt	1

Table 32 — Overlay ADS (PolyStyleType)

Name	Description	Data Type	Cardinality
color	Polygon color in 32 bit hex format (AABBGGRR order).	xsd:string	1

Table 33 — Overlay ADS (ExtendedDataType)

Name	Description	Data Type	Cardinality
Data	Metadata.	DataType	1…inf

Table 34 — Overlay ADS (DataType)

Name	Description	Data Type	Cardinality
value	Metadata value.	xsd:string	1

Table 35 — Overlay ADS ()

Name	Description	Data Type	Cardinality
tessellate	Tessellate option.	xsd:unsignedInt	1
altitudeMode	Altitude mode option.	xsd:string	1
outerBoundaryIs	Polygon outer boundary description.	outerBoundaryIsType	1

Table 36 — Overlay ADS (outerBoundaryIsType)

Name	Description	Data Type	Cardinality
LinearRing	Linear ring description.	LinearRingType	1

Table 37 — Overlay ADS (LinearRingType)

Name	Description	Data Type	Cardinality
coordinates	A string of 5 lon,lat,height coordinate tuples which describe the corners of the image. The string is of the form: lon,lat,height lon,lat,height lon,lat,height lon,lat,height lon,lat,height The coordinates must appear in the following order: last line first pixel, last line last pixel, first line last pixel, first line first pixel, last line first pixel	xsd:string	1