[Help]: Initializing SBAS post-processing - open_grids does not recognize merged grd files about pygmtsar

[AlexeyPechnikov]

You can use the function sbas.pairs() to look up all the available grids. By default, it looks for grids with the 'phasefilt' name.

sbas = SBAS(DATADIR, DEMFILE, basedir=WORKDIR, force=False)
pairs = sbas.pairs()
unwraps_ra  = sbas.open_grids(pairs, 'unwrap')
unwraps_ra 

If you call sbas.pairs(dates=True), it will not only return the interferogram pair list but also provide all the scene dates as well.

pairs, dates = sbas.pairs(dates=True)

And there's an even easier way using the dump() and restore() functions. With these functions, you only need to provide the WORKDIR path to restore your state. By the way, you can still use sbas.pairs() to detect the SBAS pairs that have already been processed.

# dump current state into WORKDIR as soon as you have the interferograms pomputed and merged
sbas.dump()
# restore the current state later from other notebook to continue
# sbas = SBAS.restore(WORKDIR)

[kinreich]

Thank you!

Running pairs, dates = sbas.pairs(dates=True) this error pops up:

Traceback (most recent call last): File "/home/name/pythonProject/venv/lib/python3.10/site-packages/IPython/core/interactiveshell.py", line 3508, in run_code exec(code_obj, self.user_global_ns, self.user_ns) File "<ipython-input-10-9f9e9cb1f67c>", line 1, in <module> pairs, dates = sbas.pairs(dates=True) File "/home/name/pythonProject/venv/lib/python3.10/site-packages/pygmtsar/SBAS_base.py", line 414, in pairs pattern = self.get_filenames(None, None, f'????????_????????_{name}') File "/home/name/pythonProject/venv/lib/python3.10/site-packages/pygmtsar/SBAS_base.py", line 166, in get_filenames subswath = self.get_subswath(subswath) File "/home/name/pythonProject/venv/lib/python3.10/site-packages/pygmtsar/SBAS_base.py", line 155, in get_subswath assert len(subswaths)==1, f'ERROR: multiple subswaths {subswaths} found, merge them first using SBAS.merge_parallel()' AssertionError: ERROR: multiple subswaths [1 2 3] found, merge them first using SBAS.merge_parallel()

[AlexeyPechnikov]

It seems that you haven't merged the subswaths, which makes the continuation of the processing a bit tricky. I recommend trying the following code for your case:

sbas = SBAS(DATADIR, DEMFILE, basedir=WORKDIR, force=False)
# get all the subswaths like [1, 2, 3]
subswaths = sbas.get_subswaths()
# emulate the merging
sbas.df = sbas.df[sbas.df.subswath == subswaths[0]]
sbas.df['subswath'] = ''.join(map(str, subswaths))
# for now, we can work as expected with the merged subswaths
pairs = sbas.pairs()
unwraps_ra  = sbas.open_grids(pairs, 'unwrap')
unwraps_ra