search

opera-adt / disp-s1

OPERA Displacement workflows
Other
5 stars 7 forks source link

Copy orbit, wavelength, and times to Compressed slc HDF5 #103

Closed scottstanie closed 1 month ago

scottstanie commented 1 month ago

Before:

$ h5ls -r output/compressed_slcs/compressed_t027_056721_iw2_20180107_20180107_20180425.h5
/                        Group
/complex64               Type
/data                    Group
/data/VV                 Dataset {4737, 20628}
/data/amplitude_dispersion Dataset {4737, 20628}
/data/spatial_ref        Dataset {SCALAR}
/data/x                  Dataset {20628}
/data/y                  Dataset {4737}

after:

In [13]: product.copy_opera_cslc_metadata( 'input_slcs/OPERA_L2_CSLC-S1_T027-056721-IW2_20180107T132747Z_20240701T182533Z_S1A_VV_v1.1.h5', 'output/compressed_slcs/compressed_t027_056721_iw2_20180107_20180107_20180425.h5')

(mapping-311) staniewi:calval$ h5ls -r output/compressed_slcs/compressed_t027_056721_iw2_20180107_20180107_20180425.h5
/                        Group
/complex64               Type
/data                    Group
/data/VV                 Dataset {4737, 20628}
/data/amplitude_dispersion Dataset {4737, 20628}
/data/spatial_ref        Dataset {SCALAR}
/data/x                  Dataset {20628}
/data/y                  Dataset {4737}
/identification          Group
/identification/zero_doppler_end_time Dataset {SCALAR}
/identification/zero_doppler_start_time Dataset {SCALAR}
/metadata                Group
/metadata/orbit          Group
/metadata/orbit/orbit_direction Dataset {SCALAR}
/metadata/orbit/orbit_type Dataset {SCALAR}
/metadata/orbit/position_x Dataset {12}
/metadata/orbit/position_y Dataset {12}
/metadata/orbit/position_z Dataset {12}
/metadata/orbit/reference_epoch Dataset {SCALAR}
/metadata/orbit/time     Dataset {12}
/metadata/orbit/velocity_x Dataset {12}
/metadata/orbit/velocity_y Dataset {12}
/metadata/orbit/velocity_z Dataset {12}
/metadata/processing_information Group
/metadata/processing_information/input_burst_metadata Group
/metadata/processing_information/input_burst_metadata/wavelength Dataset {SCALAR}

and the opera_utils functions work


In [9]: opera_utils.get_radar_wavelength('output/compressed_slcs/compressed_t027_056721_iw1_20180107_20180107_20180425.h5')
Out[9]: 0.05546576

In [10]: opera_utils.get_zero_doppler_time('output/compressed_slcs/compressed_t027_056721_iw1_20180107_20180107_20180425.h5')
Out[10]: datetime.datetime(2018, 1, 7, 13, 27, 46, 767654)

Note that this is using the sample data created from this script, where we zero out all the big rasters to make he HDF5 file ~3 MB