Static Layers: Meta issue

[scottstanie]

Tracking progress and ideas for adding the static layers workflow in the next few months.

Some initial thoughts:

Included products

All products are at 30 meters on the same UTM grid as DISP-S1. The grids should match to make applications possible with no resampling.

• Line of sight vectors
• NISAR DEM v1.1
• layover shadow masks
• ? Water Mask (or distance file?)
• ? Local incidence angle

Format:

• Geotiffs may be easiest, as it's nice on RTC to be able to search for a product and then in the same menu click the link.
• A 3-band geotiff (interleaved by pixel) for LOS seems nice. You would be able to pick out single pixel "E,N,U" remotely without downloading files.
• If we are asked to include the local incidence angle (which has fewer Displacement applications), perhaps we could just link/copy the RTC one (example here). It's already on 30 meter grid, and I'm not confident that there's a simple way to convert the 5x10 meter CSLC local incidence angle to 30 meters. I think Gustavo does more fancy estimation on the 30m grid by using the full res data.

Examples from RTC

For saving the COGs: https://github.com/opera-adt/RTC/blob/main/src/rtc/core.py#L59
For appending metadata: https://github.com/opera-adt/RTC/blob/main/src/rtc/rtc_s1_single_job.py#L527

An example product has the following metadata:

In [3]: ds = gdal.Open("/Users/staniewi/Downloads/OPERA_L2_RTC-S1_T064-135546-IW1_20241216T015159Z_20241216T090410Z_S1A_30_v1.0_VV.tif")

In [4]: ds.GetMetadataDomainList()
Out[4]: ['IMAGE_STRUCTURE', '', 'DERIVED_SUBDATASETS']

In [5]: ds.GetMetadata()
Out[5]:
{'ABSOLUTE_ORBIT_NUMBER': '57011',
 'ACQUISITION_MODE': 'IW',
 'BOUNDING_BOX': '[240240.0, 4240680.0, 339390.0, 4286520.0]',
 'BOUNDING_BOX_EPSG_CODE': '32611',
 'BOUNDING_BOX_PIXEL_COORDINATE_CONVENTION': 'Edges/corners. Xmin, Ymin, Xmax, Ymax.',
 'BOUNDING_POLYGON': 'POLYGON ((-119.8436808985369 38.33404367679761, -119.7909786104352 38.34196964658739, -119.7369020202512 38.35007421124919, -119.6746198797831 38.35936768361647, -119.6216243687353 38.36725428014282, -119.5644118998329 38.37573530956657, -119.5201240924068 38.38229055325848, -119.4684165352311 38.38991290528675, -119.422205919929 38.39670917089858, -119.3665939423655 38.40485206896142, -119.3160797143292 38.41222815844881, -119.2628276484003 38.41997587425311, -119.2086476082604 38.42783076279792, -119.1737822680648 38.43288877127543, -119.1275189649066 38.43956689670897, -119.0786608779893 38.446595292295, -119.0309059370584 38.4534448441843, -118.9826762391049 38.46034071843796, -118.9351950376387 38.46710933154327, -118.8833792003521 38.47446761343528, -118.8362915448103 38.48113731839624, -118.8758951998748 38.64681797641664, -118.929477048253 38.6392523242201, -118.9752431549253 38.6327616975538, -119.0233005518837 38.62592779949581, -119.0619549322919 38.62040577450281, -119.1105351090313 38.6134591746191, -119.1646338523148 38.60570341054071, -119.2139210736612 38.5986095101243, -119.2640075528816 38.59137829211034, -119.3076954702201 38.58504566665224, -119.3579810714041 38.57774231928702, -119.4189082768631 38.56887236611593, -119.4653704193609 38.56207356777865, -119.5168854096413 38.5545172747626, -119.5655959712712 38.54734713159461, -119.6179202917312 38.53962414486553, -119.6666283960955 38.53240869556407, -119.7169988927933 38.52492555270837, -119.7706217087444 38.51693645778429, -119.8275667765037 38.50842604002693, -119.8872796420559 38.49947195340746, -119.8436808985369 38.33404367679761))',
 'BOUNDING_POLYGON_EPSG_CODE': '4326',
 'BURST_ID': 't064_135546_iw1',
 'CENTER_FREQUENCY': '5405000454.33435',
 'CEOS_ANALYSIS_READY_DATA_DOCUMENT_IDENTIFIER': 'https://ceos.org/ard/files/PFS/NRB/v5.5/CARD4L-PFS_NRB_v5.5.pdf',
 'CEOS_ANALYSIS_READY_DATA_PRODUCT_TYPE': 'Normalised Radar Backscatter',
 'CONTACT_INFORMATION': '[email protected]',
 'INPUT_ANNOTATION_FILES': "['calibration-s1a-iw1-slc-vv-20241216t015159-20241216t015224-057011-070161-004.xml', 'noise-s1a-iw1-slc-vv-20241216t015159-20241216t015224-057011-070161-004.xml']",
 'INPUT_DEM_SOURCE': 'Digital Elevation Model (DEM) for the NASA OPERA project version 1.1 (v1.1) based on the Copernicus DEM 30-m and Copernicus 90-m referenced to the WGS84 ellipsoid',
 'INPUT_L1_SLC_GRANULES': "['S1A_IW_SLC__1SDV_20241216T015157_20241216T015224_057011_070161_F3FE.zip']",
 'INPUT_ORBIT_FILES': "['S1A_OPER_AUX_RESORB_OPOD_20241216T035820_V20241216T000246_20241216T032016.EOF']",
 'INSTITUTION': 'NASA JPL',
 'INSTRUMENT_NAME': 'Sentinel-1A CSAR',
 'ISCE3_VERSION': '0.15.0',
 'LAYER_DESCRIPTION': 'Radiometric terrain corrected Sentinel-1 VV backscatter coefficient normalized to gamma0',
 'LAYER_NAME': 'RTC-S1 VV Backscatter',
 'LOOK_DIRECTION': 'right',
 'ORBIT_PASS_DIRECTION': 'ascending',
 'PLATFORM': 'Sentinel-1A',
 'PROCESSING_DATETIME': '2024-12-16T09:04:15.684054Z',
 'PROCESSING_INFORMATION_BISTATIC_DELAY_CORRECTION_APPLIED': 'True',
 'PROCESSING_INFORMATION_BURST_GEOGRID_SNAP_X': '30.0',
 'PROCESSING_INFORMATION_BURST_GEOGRID_SNAP_Y': '30.0',
 'PROCESSING_INFORMATION_DEM_EGM_MODEL': 'Earth Gravitational Model 2008 (EGM2008)',
 'PROCESSING_INFORMATION_DEM_INTERPOLATION_ALGORITHM': 'biquintic',
 'PROCESSING_INFORMATION_FILTERING_APPLIED': 'False',
 'PROCESSING_INFORMATION_GEOCODING_ALGORITHM': 'Area-Based SAR Geocoding with Adaptive Multilooking (GEO-AP)',
 'PROCESSING_INFORMATION_GEOCODING_ALGORITHM_REFERENCE': 'Gustavo H. X. Shiroma, Marco Lavalle, and Sean M. Buckley, "An Area-Based Projection Algorithm for SAR Radiometric Terrain Correction and Geocoding," in IEEE Transactions on Geoscience and Remote Sensing, vol. 60, pp. 1-23, 2022, Art no. 5222723, doi: 10.1109/TGRS.2022.3147472.',
 'PROCESSING_INFORMATION_INPUT_BACKSCATTER_NORMALIZATION_CONVENTION': 'beta0',
 'PROCESSING_INFORMATION_MULTILOOKING_APPLIED': 'False',
 'PROCESSING_INFORMATION_NOISE_CORRECTION_APPLIED': 'True',
 'PROCESSING_INFORMATION_NOISE_REMOVAL_ALGORITHM_REFERENCE': 'Thermal Denoising of Products Generated by the S-1 IPF. ESA document. Accessed: May 29, 2023. [Online]. Available: https://sentinels.copernicus.eu/documents/247904/2142675/Thermal-Denoising-of-Products-Generated-by-Sentinel-1-IPF.pdf',
 'PROCESSING_INFORMATION_OUTPUT_BACKSCATTER_DECIBEL_CONVERSION_EQUATION': '10*log10(backscatter_linear)',
 'PROCESSING_INFORMATION_OUTPUT_BACKSCATTER_EXPRESSION_CONVENTION': 'linear backscatter intensity',
 'PROCESSING_INFORMATION_OUTPUT_BACKSCATTER_NORMALIZATION_CONVENTION': 'gamma0',
 'PROCESSING_INFORMATION_RADIOMETRIC_TERRAIN_CORRECTION_ALGORITHM': 'Area-Based SAR Radiometric Terrain Correction (RTC-AP)',
 'PROCESSING_INFORMATION_RADIOMETRIC_TERRAIN_CORRECTION_ALGORITHM_REFERENCE': 'Gustavo H. X. Shiroma, Marco Lavalle, and Sean M. Buckley, "An Area-Based Projection Algorithm for SAR Radiometric Terrain Correction and Geocoding," in IEEE Transactions on Geoscience and Remote Sensing, vol. 60, pp. 1-23, 2022, Art no. 5222723, doi: 10.1109/TGRS.2022.3147472.',
 'PROCESSING_INFORMATION_RADIOMETRIC_TERRAIN_CORRECTION_APPLIED': 'True',
 'PROCESSING_INFORMATION_STATIC_TROPOSPHERIC_GEOLOCATION_CORRECTION_APPLIED': 'True',
 'PROCESSING_INFORMATION_WET_TROPOSPHERIC_GEOLOCATION_CORRECTION_APPLIED': 'False',
 'PROCESSING_TYPE': 'NOMINAL',
 'PRODUCT_DATA_ACCESS': 'https://search.asf.alaska.edu/#/?dataset=OPERA-S1&productTypes=RTC',
 'PRODUCT_LEVEL': 'L2',
 'PRODUCT_SPECIFICATION_VERSION': '1.0',
 'PRODUCT_TYPE': 'RTC-S1',
 'PRODUCT_VERSION': '1.0',
 'PROJECT': 'OPERA',
 'QA_GEOMETRIC_ACCURACY_BIAS_X': '-0.72',
 'QA_GEOMETRIC_ACCURACY_BIAS_Y': '-0.67',
 'QA_GEOMETRIC_ACCURACY_STDDEV_X': '0.7',
 'QA_GEOMETRIC_ACCURACY_STDDEV_Y': '0.62',
 'QA_RFI_BURST_REPORT_AZIMUTH_TIME': '2024-12-16T01:51:59.139351Z',
 'QA_RFI_BURST_REPORT_SWATH': 'IW1',
 'QA_RFI_INFO_AVAILABLE': 'True',
 'QA_RFI_IN_BAND_OUT_BAND_POWER_RATIO': '8.31855',
 'QA_RFI_MITIGATION_DOMAIN': 'TimeAndFrequency',
 'QA_RFI_MITIGATION_PERFORMED': 'BasedOnNoiseMeas',
 'RADAR_BAND': 'C',
 'S1_READER_VERSION': '0.2.4',
 'SOFTWARE_VERSION': '1.0.2',
 'SOURCE_DATA_ACCESS': 'https://search.asf.alaska.edu/#/?dataset=SENTINEL-1&productTypes=SLC',
 'SOURCE_DATA_AVERAGE_ZERO_DOPPLER_SPACING_IN_METERS': '13.96045',
 'SOURCE_DATA_AZIMUTH_RESOLUTION_IN_METERS': '22.5',
 'SOURCE_DATA_INSTITUTION': 'ESA',
 'SOURCE_DATA_NUMBER_OF_ACQUISITIONS': '1',
 'SOURCE_DATA_NUMBER_OF_AZIMUTH_LINES': '1497',
 'SOURCE_DATA_NUMBER_OF_RANGE_SAMPLES': '21576',
 'SOURCE_DATA_PROCESSING_CENTER': 'organization: "ESA", site: "Production Service-SERCO", country: "Italy"',
 'SOURCE_DATA_PROCESSING_DATETIME': '2024-12-16T06:37:11.000000',
 'SOURCE_DATA_PRODUCT_LEVEL': 'L1',
 'SOURCE_DATA_RANGE_BANDWIDTH': '56500000.0',
 'SOURCE_DATA_SLANT_RANGE_RESOLUTION_IN_METERS': '2.7',
 'SOURCE_DATA_SLANT_RANGE_SPACING': '2.329562114715323',
 'SOURCE_DATA_SLANT_RANGE_START': '800082.3526126802',
 'SOURCE_DATA_SOFTWARE_VERSION': '3.90',
 'SOURCE_DATA_ZERO_DOPPLER_END_TIME': '2024-12-16T01:52:02.214463Z',
 'SOURCE_DATA_ZERO_DOPPLER_START_TIME': '2024-12-16T01:51:59.139351Z',
 'SOURCE_DATA_ZERO_DOPPLER_TIME_SPACING': '0.002055556299999998',
 'STATIC_LAYERS_DATA_ACCESS': 'https://search.asf.alaska.edu/#/?dataset=OPERA-S1&productTypes=RTC-STATIC&operaBurstID=T064-135546-IW1&end=2024-12-16',
 'SUB_SWATH_ID': 'IW1',
 'TRACK_NUMBER': '64',
 'ZERO_DOPPLER_END_TIME': '2024-12-16T01:52:02.214463Z',
 'ZERO_DOPPLER_START_TIME': '2024-12-16T01:51:59.139351Z',
 'AREA_OR_POINT': 'Area'}