Merge pull request #1016 from jungkyoj/debug_h5prep_multipol

Debug h5prep multipol

Author: jungkyoJung

python/packages/nisar/workflows/h5_prep.py

@@ -153,8 +153,10 @@ def cp_geocode_meta(cfg, output_hdf5, dst):
                                     f'{ident_path}/{dst_data}')
 
         # Delete flag isGeocoded if exist, and assign it again
+        ident = dst_h5[ident_path]
+        if 'isGeocoded' in ident:
+            del ident['isGeocoded']
         is_geocoded = dst in ['GCOV', 'GSLC', 'GUNW', 'GOFF']
-
         dst_h5[ident_path].create_dataset('isGeocoded',
                                           data=np.string_(str(is_geocoded)))
         desc = "Flag to indicate radar geometry or geocoded product"
@@ -476,7 +478,7 @@ def prep_ds_insar(pcfg, dst, dst_h5):
                                              'RUNW'] else 'grids'
             product_path = f'{common_path}/{dst}'
             freq_path = f'{product_path}/{grid_swath}/frequency{freq}'
-            dst_h5[product_path].create_group(grid_swath)
+            dst_h5[product_path].require_group(grid_swath)
             dst_h5[f'{product_path}/{grid_swath}'].create_group(
                 f'frequency{freq}')
 
@@ -939,17 +941,17 @@ def prep_ds_insar(pcfg, dst, dst_h5):
                                          'crossCorrelationMethod',
                                          descr=descr, units=None, data=lay_cfg.get('cross_correlation_method'),
                                          long_name='cross correlation method')
-            # Add perpendicular and parallel baseline
-            descr = "Perpendicular component of the InSAR baseline"
-            _create_datasets(dst_h5[grid_path], igram_shape, np.float64,
-                             "perpendicularBaseline",
-                             descr=descr, units="meters",
-                             long_name='perpendicular baseline')
-            _create_datasets(dst_h5[grid_path], igram_shape, np.float64,
-                             "parallelBaseline",
-                             descr=descr.replace('Perpendicular', 'Parallel'),
-                             units="meters",
-                             long_name='parallel baseline')
+        # Add perpendicular and parallel baseline
+        descr = "Perpendicular component of the InSAR baseline"
+        _create_datasets(dst_h5[grid_path], igram_shape, np.float64,
+                            "perpendicularBaseline",
+                            descr=descr, units="meters",
+                            long_name='perpendicular baseline')
+        _create_datasets(dst_h5[grid_path], igram_shape, np.float64,
+                            "parallelBaseline",
+                            descr=descr.replace('Perpendicular', 'Parallel'),
+                            units="meters",
+                            long_name='parallel baseline')
 
 
 def get_off_params(pcfg, param_name, is_roff=False, pattern=None,