Quick PR to add the listOfCovarianceTerms to GCOV products (#855)

* add HDF5 dataset listOfCovarianceTerms to GCOV products

* fix function _save_list_cov_terms()

* fix function _save_list_cov_terms() (2)

* fix function _save_list_cov_terms() (3)

* extend _save_list_cov_terms() to off-diagonal terms

* extend _save_list_cov_terms() to off-diagonal terms (2)

* extend _save_list_cov_terms() to off-diagonal terms (3)

Author: gshiroma

python/packages/nisar/workflows/gcov.py

@@ -346,13 +346,19 @@ def run(cfg):
             yds = hdf5_obj[os.path.join(root_ds, 'yCoordinates')]
             cov_elements_list = [p.upper()+p.upper() for p in pol_list]
 
+            # save GCOV imagery
             _save_hdf5_dataset(temp_output.name, hdf5_obj, root_ds,
                                yds, xds, cov_elements_list,
                                long_name=output_radiometry_str,
                                units='',
                                valid_min=clip_min,
                                valid_max=clip_max)
 
+            # save listOfCovarianceTerms
+            freq_group = hdf5_obj[root_ds]
+            if not flag_fullcovariance:
+                _save_list_cov_terms(cov_elements_list, freq_group)
+
             # save nlooks
             if flag_save_nlooks:
                 _save_hdf5_dataset(temp_nlooks.name, hdf5_obj, root_ds, 
@@ -406,7 +412,8 @@ def run(cfg):
                         if (b2 <= b1):
                             continue
                         off_diag_terms_list.append(p1.upper()+p2.upper())
-    
+                _save_list_cov_terms(cov_elements_list + off_diag_terms_list, 
+                                     freq_group)
                 _save_hdf5_dataset(temp_off_diag.name, hdf5_obj, root_ds,
                                    yds, xds, off_diag_terms_list,
                                    long_name = output_radiometry_str, 
@@ -444,6 +451,17 @@ def run(cfg):
     t_all_elapsed = time.time() - t_all
     info_channel.log(f"successfully ran geocode COV in {t_all_elapsed:.3f} seconds")
 
+
+def _save_list_cov_terms(cov_elements_list, dataset_group):
+
+    name = "listOfCovarianceTerms"
+    cov_elements_list.sort()
+    cov_elements_array = np.array(cov_elements_list, dtype="S4")
+    dset = dataset_group.create_dataset(name, data=cov_elements_array)
+    desc = f"List of processed covariance terms"
+    dset.attrs["description"] = np.string_(desc)
+
+
 def _save_hdf5_dataset(ds_filename, h5py_obj, root_path,
                        yds, xds, ds_name, standard_name=None,
                        long_name=None, units=None, fill_value=None,