This is a draft version. Please visit the CEOS-ARD website for the latest endorsed version of this document.
Product Family Specification, Synthetic Aperture Radar, Normalised Radar Backscatter
Proposed revisions may be provided to: ard-contact@lists.ceos.org
Note: This document is the successor of the former CEOS-ARD for SAR PFS v1.3.1 for product type Normalised Radar Backscatter (NRB).
Justification: Migration to building blocks.
Editor: Matthias Mohr
CEOS Analysis Ready Data (CEOS-ARD) are satellite data that have been processed to a minimum set of requirements and organized into a form that allows immediate analysis with a minimum of additional user effort and interoperability both through time and with other datasets.
Product Family Specification: Synthetic Aperture Radar, Normalised Radar Backscatter (NRB)
Version: 5.6.0-draft
Applies to: Data collected by Synthetic Aperture Radar sensors
This PFS is specifically aimed at users interested in exploring the potential of SAR but who may lack the expertise or facilities for SAR processing.
The CEOS-ARD Normalised Radar Backscatter (NRB) specification describes products that have been subject to Radiometric Terrain Correction (RTC) and are provided in the Gamma-Nought () backscatter convention (Small 2011), which mitigates the variations from diverse observation geometries and is recommended for most land applications. An additional metadata layer can be optionally provided for conversion of to Sigma-Nought () backscatter layer for compatibility with legacy software or numerical models.
As the NRB product contains backscatter values only, it cannot be directly used for SAR polarimetry or interferometric applications that require relative polarization phase or local phase estimates respectively. However, as an option, a “flattened” phase data layer can be provided with an NRB product for enabling InSAR analysis. The flattened phase is the interferometric phase, with respect to a reference orbit and to a digital elevation model (DEM), for which the topographic phase contribution is removed.
The NRB specification accommodates basic mosaic products generated from more than one input data source, where each pixel value in the NRB product uniquely corresponds to the pixel value of one of its input data sources, indicated in the Acquisition ID per-pixel metadata. For products where the pixel value is dependent on multiple input sources, see the Composite Backscatter (CB) PFS.
WARNING: The section numbers in front of the title (e.g. 1.1) are not stable and may change or may be removed at any time. Do not use the numbers to refer back to specific requirements! Instead, use the textual identifier that is provided below the title.
1.
General MetadataThese are metadata records describing a distributed collection of pixels. The collection of pixels referred to must be contiguous in space and time. General metadata should allow the user to assess the overall suitability of the dataset, and must meet the requirements listed below.
1.1.
TraceabilityIdentifier: meta-trace-sar
Not required.
Data must be traceable to SI reference standard.
Notes:
1.2.
Metadata Machine ReadabilityIdentifier: meta-memare-sar
Metadata is provided in a structure that enables a computer algorithm to be used to consistently and automatically identify and extract each component/variable for further use.
As threshold, but metadata is formatted in accordance with the latest corresponding CEOS-ARD SAR Metadata Specifications, or in a community endorsed standard that facilitates machine-readability, such as ISO 19115-2, Climate and Forecast (CF) convention and the Attribute Convention for Data Discovery (ACDD), etc.
1.3.
Product TypeIdentifier: meta-protype
CEOS-ARD product type name – or names in case of compliance with more than one product type – and, if required by the data provider, copyright.
As threshold.
1.4.
Document IdentifierIdentifier: meta-pfsurl
Reference to CEOS-ARD PFS document as URL.
As threshold.
1.5.
Data Collection TimeIdentifier: meta-time-sar
Number of source data acquisitions of the data collection is identified. The start and stop UTC time of data collection is identified in the metadata, expressed in date/time. In the case of composite or mosaic products, the dates/times of the first and last data takes is provided with the product.
As threshold, but using ISO 8601 time format.
2.
Source MetadataMetadata describing (detailing) each acquisition used to generate the ARD product.
Source data attribute information can refer to other products for higher level ARD derived from those, under the condition of their availability (Source Metadata: Source Data Access).
2.1.
Acquisition IDIdentifier: src-macqid
Source data attribute information are described for each acquisition and sequentially identified e.g. as acqID = 1, 2, 3, …
As threshold.
2.2.
Source Data AccessIdentifier: src-daccess-src
The metadata identifies the location from where the source data can be retrieved, expressed as a URL or DOI.
The metadata identifies an online location from where the data can be consistently and reliably retrieved by a computer algorithm without any manual intervention being required.
2.3.
InstrumentIdentifier: src-instru-sar
The instrument used to collect the data is identified in the metadata:
As threshold, but using CEOS Mission-Instruments-Measurements (MIM) database as reference.
2.4.
Source Data Acquisition TimeIdentifier: src-time-src
The start date and time of source data is identified in the metadata, expressed in UTC in date and time, at least to the second.
As threshold.
2.5.
Source Data Acquisition ParametersIdentifier: src-acqpar
Acquisition parameters related to the SAR antenna:
As threshold.
2.6.
Source Data Orbit InformationIdentifier: src-orbit
Information related to the platform orbit used for data processing:
Note:
As threshold, including also:
2.7.
Source Data Processing ParametersIdentifier: src-propar
Processing parameters details of the source data:
Note:
As threshold, plus additional relevant processing parameters, e.g., range- and azimuth look bandwidth and LUT applied.
2.8.
Source Data Image AttributesIdentifier: src-imgatt-sar
Image attributes related to the source data:
Note:
Geometry of the image footprint expressed in WGS84 in a standardised format (e.g., WKT).
2.9.
Sensor CalibrationIdentifier: src-sencal-sar
Not required.
Sensor calibration parameters are identified in the metadata or can be accessed using details included in the metadata. Ideally this would support machine-to-machine access.
2.10.
Performance IndicatorsIdentifier: src-perfind
Provide performance indicators on data intensity noise level ( and/or and/or , i.e., noise equivalent Sigma- and/or Beta- and/or Gamma-Nought). Provided for each polarization channel when available.
Parameter may be expressed as the mean and/or minimum and maximum noise equivalent values of the source data.
Values do not need to be estimated individually for each product, but may be estimated once for each acquisition mode, and annotated on all products.
Provide additional relevant performance indicators (e.g., ENL, PSLR, ISLR, and performance reference DOI or URL).
2.11.
Polarimetric Calibration MatricesIdentifier: src-polcalm
Not required.
The complex-valued polarimetric distortion matrices with the channel imbalance and the cross-talk applied for the polarimetric calibration.
2.12.
Mean Faraday Rotation AngleIdentifier: src-farotan
Not required.
The mean Faraday rotation angle estimated from the polarimetric data and/or from models with reference to the method or paper used to derive the estimate.
2.13.
Ionosphere IndicatorIdentifier: src-ionind
Not required.
Flag indicating whether the backscatter imagery is “significantly impacted” by the ionosphere (0 – false, 1 – true). Significant impact would imply that the ionospheric impact on the backscatter exceeds the radiometric calibration requirement or goal for the imagery.
3.
Product MetadataInformation related to the CEOS-ARD product generation procedure and geographic parameters.
3.1.
Product Data AccessIdentifier: prd-daccess-prod
Processing parameters details of the CEOS-ARD product:
The metadata identifies an online location from where the data can be consistently and reliably retrieved by a computer algorithm without any manual intervention being required.
3.2.
Auxiliary DataIdentifier: prd-auxdat-sar
Not required.
The metadata identifies the sources of auxiliary data used in the generation process, ideally expressed as DOIs.
Note:
3.3.
Product Sample SpacingIdentifier: prd-samspa
CEOS-ARD product processing parameters details:
As threshold.
3.4.
Product Equivalent Number of LooksIdentifier: prd-enl
Not required.
Equivalent Number of Looks (ENL)
3.5.
Product ResolutionIdentifier: prd-resol
Not required.
Average spatial resolution of the CEOS-ARD product along:
3.6.
Product FilteringIdentifier: prd-spekfil
Flag if speckle filter has been applied (true/false).
Metadata should include:
As threshold.
3.7.
Product Bounding BoxIdentifier: prd-geobbox
Two opposite corners of the product file (bounding box, including any zero-fill values) are identified, expressed in the coordinate reference system defined in Product Metadata: Product Coordinate Reference System.
Four corners of the product file are recommended for scenes crossing the Antemeridian, or the North or the South Pole.
As threshold.
3.8.
Product Geographical ExtentIdentifier: prd-geoarea-sar
The geometry of the SAR image footprint expressed in longitude/latitude based on WGS84 (EPSG 4326), in a standardised format (e.g., WKT Polygon).
As threshold.
3.9.
Product Image SizeIdentifier: prd-imgsize
Image attributes of the CEOS-ARD product:
As threshold.
3.10.
Product Pixel Coordinate ConventionIdentifier: prd-pixcoco
Coordinate referring to the centre, the upper left corner, or the lower left corner of a pixel. Values are pixel centre, pixel ULC or pixel LLC.
As threshold.
3.11.
Product Coordinate Reference SystemIdentifier: prd-crs-sar
The metadata lists the map projection (or geographical coordinates, if applicable) that was used and any relevant parameters required to geolocate data in that map projection, expressed in a standardised format (e.g., WKT).
Indicate EPSG code, if defined for the CRS.
As threshold.
3.12.
Reference OrbitIdentifier: prd-reorbit-nrb-pol
Usage: Only when Flattened phase per-pixel metadata (see Radiometrically Corrected Measurements: Flattened Phase) is provided.
Not required.
Provide the absolute orbit number used as reference for topographic phase flattening. In case a virtual orbit has been used, provide orbit parameters or orbit state vectors as DOI or URL.
4.
Per-Pixel MetadataThe following minimum metadata specifications apply to each pixel. Whether the metadata is provided in a single record relevant to all pixels or separately for each pixel is at the discretion of the data provider. Per-pixel metadata should allow users to discriminate between (choose) observations on the basis of their individual suitability for application.
4.1.
Metadata Machine ReadabilityIdentifier: pxl-memare-sar
Metadata is provided in a structure that enables a computer algorithm to be used to consistently and automatically identify and extract each component/variable for further use.
As threshold, but metadata is formatted in accordance with the latest corresponding CEOS-ARD SAR Metadata Specifications, or in a community endorsed standard that facilitates machine-readability, such as ISO 19115-2, Climate and Forecast (CF) convention and the Attribute Convention for Data Discovery (ACDD), etc.
4.2.
Data Mask ImageIdentifier: pxl-damaski
Mask image indicating:
File format specifications/contents provided in metadata:
Notes:
As threshold, including additional bit value representations, e.g.:
4.3.
Scattering Area ImageIdentifier: pxl-piscata
Usage: Recommended for scenes that include land areas.
Not required.
DEM-based scattering area image used for Gamma-Nought terrain normalisation is provided. This quantifies the local scattering area used to normalise for radiometric distortions induced by terrain to the measured backscatter. The terrain-flattened is best understood as divided by the local scattering area.
File format specifications/contents provided in metadata:
Notes:
4.4.
Local Incident Angle ImageIdentifier: pxl-ploinca
DEM-based Local Incident angle image is provided.
File format specifications/contents provided in metadata:
Note:
As threshold.
4.5.
Ellipsoidal Incident Angle ImageIdentifier: pxl-pelinca
Not required.
Ellipsoidal incident angle is provided.
File format specifications/contents provided in metadata:
Note:
4.6.
Noise Power ImageIdentifier: pxl-pinopow
Not required.
Estimated Noise Equivalent (or or , as applicable) used for noise removal, if applied, for each channel. and are both based on either an ellipsoid Earth model or the local topography.
File format specifications/contents provided in metadata:
4.7.
Gamma-to-Sigma Ratio ImageIdentifier: pxl-gasiri
Not required.
Ratio of the integrated area in the Gamma projection over the integrated area in the Sigma projection (ground). Multiplying RTC by this ratio results in an estimate of RTC .
File format specifications/contents provided in metadata:
Note:
4.8.
Acquisition ID ImageIdentifier: pxl-pacqidm
Usage: Required for mosaic products only.
Acquisition ID, or acquisition date, for each pixel is identified.
In case of multi-temporal image stacks, use source acquisition ID (i.e., Source Metadata: Acquisition ID) to list contributing images.
In case of date, data represent (integer or fractional) day offset to reference observation date (in UTC). Date used as reference (“Day 0”) is provided in the metadata.
Pixels not representing a unique date or ID (e.g., pixels averaged in image overlap zones) are flagged with a pixel value referencing a date range that is provided in the metadata.
File format specifications/contents provided in metadata:
As threshold.
4.9.
Per-Pixel DEMIdentifier: pxl-pidem
Not required.
Provide DEM or DSM as used during the geometric and radiometric processing of the SAR data, resampled to an exact geometric match in extent and resolution with the CEOS-ARD SAR image product.
File format specifications/contents provided in metadata:
Note:
5.
Radiometrically Corrected MeasurementsThe requirements indicate the necessary outcomes and, to some degree, the minimum steps necessary to be deemed to have achieved those outcomes. Radiometric corrections must lead to normalised measurement(s) of backscatter intensity and/or decomposed polarimetric parameters. As for the per-pixel metadata, information regarding data format specification needs to be provided for each record. The requirements below must be met for all pixels/samples/observations in a collection.
5.1.
Backscatter Measurements (NRB)Identifier: rcm-backsca-nrb
“Terrain-flattened” Radiometrically Terrain Corrected (RTC) Gamma-Nought backscatter coefficient () is provided for each polarization.
File format specifications/contents provided in metadata:
Note:
As threshold.
5.2.
Scaling ConversionIdentifier: rcm-scaconv
If applicable, indicate the equation to convert pixel linear amplitude/power to logarithmic decibel scale, including, if applicable, the associated calibration (dB offset) factor, and/or the equation used to convert compressed data (int8/int16/float16) to float32.
As threshold, but use of float32.
5.3.
Noise RemovalIdentifier: rcm-noiser
Flag if noise removal (see note) has been applied (Y/N). Metadata should include the noise removal algorithm and reference to the algorithm as URL or DOI.
Note:
As threshold.
5.4.
Radiometric Terrain Correction AlgorithmIdentifier: rcm-radtalg-appl
Adjustments were made for terrain by modelling the local contributing scattering area using the preferred choice of a published peer-reviewed algorithm to produce radiometrically terrain corrected (RTC) backscatter estimates.
Metadata references, e.g.
Note:
As threshold.
5.5.
Radiometric AccuracyIdentifier: rcm-radacc-sar
Not required.
Uncertainty (e.g., bounds on or ) information is provided as document referenced as URL or DOI. SI traceability is achieved.
5.6.
Flattened PhaseIdentifier: rcm-flapha
Usage: Alternative to GSLC product for NRB and POL products
Not required.
The Flattened Phase is the interferometric phase for which the topographic phase contribution is removed. It is derived from the range-Doppler SLC product using a DEM and the orbital state vectors with respect to a reference orbit (see annex “Topographic phase removal” in the applicable PFS). The use of the Flattened Phase with the NRB or POL intensity (Radiometrically Corrected Measurements) provides the GSLC equivalent, as follows:
File format specifications/contents provided in metadata:
In case of polarimetric data, indicate the reference polarization.
6.
Geometric CorrectionsGeometric corrections are steps that are taken to place the measurement accurately on the surface of the Earth (that is, to geolocate the measurement) allowing measurements taken through time to be compared. This section specifies any geometric correction requirements that must be met in order for the data to be analysis ready.
6.1.
Geometric Correction AlgorithmIdentifier: gcor-geocalg
Not required.
Metadata references, e.g.:
Note:
6.2.
Digital Elevation ModelIdentifier: gcor-cdem
Usage: For products including land areas.
6.3.
Geometric AccuracyIdentifier: gcor-geomacc-sar
Accurate geolocation is a prerequisite to radar processing to correct for terrain and to enable interoperability between radar sensors.
The absolute geolocation error (ALE) for a sensor is typically assessed through analysis of Single Look Complex (SLC) imagery and measured along the slant range and azimuth directions (case A: SLC ALE). The end-to-end “ARD” ALE of the final CEOS-ARD product could be measured directly in the final image product in the chosen map projection, i.e., in the map coordinate directions: e.g., Northing and Easting (case B: ARD ALE). Providing accuracy estimates based on measurements following at least one scheme (A or B or both) meets the threshold requirement.
Estimates of the ALE is provided as a bias and a standard deviation, with (Case A) SLC ALE expressed in slant range and azimuth, and (Case B) ARD ALE expressed in map projection dimensions.
For composite products, when sources come from different SAR platforms or different beam modes, provide averaged ALE or averaged ARD ALE.
Notes:
Output product sub-sample accuracy should be less than or equal to 0.1 (slant range) pixel radial root mean square error (rRMSE).
Provide documentation of estimates of ALE as DOI or URL.
6.4.
Geometric Refined AccuracyIdentifier: gcor-georacc
Not required.
Values provided under Geometric Corrections: Geometric Accuracy are provided by the SAR mission Cal/Val team.
CEOS-ARD processing steps could include method refining the geometric accuracy, such as cross-correlation of the SAR data in slant range with a SAR scene simulated from a DSM or DEM.
Methodology used (name and reference), quality flag, geometric standard deviation values should be provided.
For composite products, provide averaged ALE or averaged ARD ALE estimated from all sources.
6.5.
Gridding ConventionIdentifier: gcor-gridconv
A consistent gridding/sampling frame is used. The origin is chosen to minimise any need for subsequent resampling between multiple products (be they from the same or different providers). This is typically accomplished via a “snap to grid” in relation to the most proximate grid tile in a global system.
Note:
Provide DOI or URL to gridding convention used.
When multiple providers share a common map projection, providers are encouraged to standardise the origins of their products among each other.
In the case of UTM/UPS coordinates, the upper left corner coordinates should be set to an integer multiple of sample intervals from a 100 km by 100 km grid tile of the Military Grid Reference System’s 100k coordinates (“snap to grid”).
For products presented in geographic coordinates (latitude and longitude), the origin should be set to an integer multiple of samples in relation to the closest integer degree.
This section aims to provide background and specific information on the processing steps that can be used to achieve analysis ready data for a specific and well-developed Product Family Specification. This Guidance material does not replace or override the specifications.
CEOS-ARD are products that have been processed to a minimum set of requirements and organized into a form that allows immediate analysis with a minimum of additional user effort. In general, these products would be resampled onto a common geometric grid (for a given product) and would provide baseline data for further interoperability both through time and with other datasets.
CEOS-ARD products are intended to be flexible and accessible products suitable for a wide range of users for a wide variety of applications, including particularly time series analysis and multi-sensor application development. They are also intended to support rapid ingestion and exploitation via high-performance computing, cloud computing and other future data architectures. They may not be suitable for all purposes and are not intended as a replacement for other types of satellite products.
The CEOS-ARD branding is applied to a particular product once:
Agencies or other entities considering undertaking an assessment process should consult the CEOS-ARD Governance Framework.
A product can continue to use CEOS-ARD branding as long as its generation and distribution remain consistent with the peer-reviewed assessment.
Threshold (Minimum) requirements are the minimum that is needed for the data to be analysis ready. This must be practical and accepted by the data producers.
Goal (Desired) requirements (previously referred to as “Target”) are the ideal; where we would like to be. Some providers may already meet these.
Products that meet all threshold requirements should be immediately useful for scientific analysis or decision-making.
Products that meet goal requirements will reduce the overall product uncertainties and enhance broad-scale applications. For example, the products may enhance interoperability or provide increased accuracy through additional corrections that are not reasonable at the threshold level.
Goal requirements anticipate continuous improvement of methods and evolution of community expectations, which are both normal and inevitable in a developing field. Over time, goal specifications may (and subject to due process) become accepted as threshold requirements.
As can be seen from the individual PFS descriptions, only a few minor details in terms of generated parameters and/or the addition of supplemental data distinguish these CEOS-ARD products. In part, they are to a large extent all backward-compatible. For example, POL products implicitly include NRB products, while a coastal NRB or POL product can simply be made compatible with other ORB products by applying gamma-to-sigma conversion. Just as GSLC can be converted to NRB (given that terrain-flattening was applied, a goal-requirement for GSLC, the inverse conversion can be made true by including the optional topographically flattened phase. In this way a NRB or POL product can be used like a GSLC for InSAR applications. Consequently, it becomes obvious that they all can follow a common approach, in terms of content and structure, in order to optimize their interoperability.
The radiometric interoperability of CEOS-ARD SAR products is ensured by a common processing chain during production. The recommended processing roadmap involves the following steps:
Table 1 lists possible sequential steps and existing software tools (e.g., Gamma software (GAMMA, 2018)) and scripting tasks that can be used to form the CEOS-ARD SAR processing roadmap.
| Step | Implementation option |
|---|---|
| 1. Orbital data refinement | Check xml date and delivered format. RADARSAT-2, pre EDOT (July 2015) replace. Post July 2015, check if ‘DEF’, otherwise replace. (Gamma - RSAT2_vec) |
| 2. Apply radiometric scaling Look-Up Table (LUT) to Beta-Nought | Specification of LUT on ingest. (Gamma - par_RSAT2_SLC/SG) |
| 3. Generate covariance matrix elements | Gamma – COV_MATRIX |
| 4. Radiometric terrain normalisation | Gamma - geo_radcal2 |
| 5. Speckle filtering (Boxcar or Sigma Lee) | Custom scripting |
| 6. Geometric terrain correction/Geocoding | Gamma – gc_map and geocode_back |
| 7. Create metadata | Custom scripting |