connectComponents doesn't exist

234 views
Skip to first unread message

Asmit Roy Burman

unread,
Nov 7, 2023, 12:16:55 PM11/7/23
to MintPy
Hey Yunjun,

Continuing to the conversation subject: False Signal for Subsidence, I tried execution the correct_unwrap_error step with flags set as mintpy.unwrapError.method = phase_closure.

The following error message was received:

MintPy version 1.5.2, date 2023-08-09 --RUN-at-2023-11-07 17:13:34.418656-- Current directory: /home/jovyan/Kolkata/Descending Run routine processing with smallbaselineApp.py on steps: ['correct_unwrap_error'] Remaining steps: ['invert_network', 'correct_LOD', 'correct_SET', 'correct_troposphere', 'deramp', 'correct_topography', 'residual_RMS', 'reference_date', 'velocity', 'geocode', 'google_earth', 'hdfeos5'] -------------------------------------------------- Project name: Kolkata Go to work directory: /home/jovyan/Kolkata/Descending/output read custom template file: /home/jovyan/Kolkata/Descending/output/Kolkata.txt update default template based on input custom template No new option value found, skip updating /home/jovyan/Kolkata/Descending/output/smallbaselineApp.cfg read default template file: /home/jovyan/Kolkata/Descending/output/smallbaselineApp.cfg ******************** step - correct_unwrap_error ******************** Input data seems to be geocoded. Lookup file not needed. unwrap_error_phase_closure.py /home/jovyan/Kolkata/Descending/output/inputs/ifgramStack.h5 --template /home/jovyan/Kolkata/Descending/output/smallbaselineApp.cfg --update --cc-mask /home/jovyan/Kolkata/Descending/output/maskConnComp.h5 read options from tempalte file: smallbaselineApp.cfg -------------------------------------------------- update mode: ON 1) output dataset: unwrapPhase_phaseClosure NOT found. run or skip: run. -------------------------------------------------- calculating the integer ambiguity for the common regions defined in /home/jovyan/Kolkata/Descending/output/maskConnComp.h5 reference pixel in y/x: (571, 692) from dataset: unwrapPhase number of interferograms: 404 number of triplets: 354 read common mask from /home/jovyan/Kolkata/Descending/output/maskConnComp.h5 refine common mask based on water mask file waterMask.h5 remove regions with area < 2500 number of common regions: 7 number of samples per region: 100 solving the phase-unwrapping integer ambiguity for unwrapPhase based on the closure phase of interferograms triplets (Yunjun et al., 2019) using the L1-norm regularzed least squares approximation (LASSO) ... 1/7[==================================================] 5s / 0s 2/7 skip calculation for the reference region 3/7[==================================================] 5s / 0s 4/7[==================================================] 5s / 0s 5/7[==================================================] 5s / 0s 6/7[==================================================] 5s / 0s 7/7[==================================================] 5s / 0s figure size : [6.0, 7.8] saved common regions and sample pixels to file common_region_sample.png -------------------------------------------------- correct unwrapping error in /home/jovyan/Kolkata/Descending/output/inputs/ifgramStack.h5 with phase closure ... open ifgramStack file: ifgramStack.h5 read water mask from file: waterMask.h5 open /home/jovyan/Kolkata/Descending/output/inputs/ifgramStack.h5 with r+ mode input dataset: unwrapPhase output dataset: unwrapPhase_phaseClosure create /unwrapPhase_phaseClosure of np.float32 in size of (404, 1273, 979) Traceback (most recent call last): File "/home/jovyan/.local/envs/earthscope_insar/bin/smallbaselineApp.py", line 10, in <module> sys.exit(main()) ^^^^^^ File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/cli/smallbaselineApp.py", line 208, in main run_smallbaselineApp(inps) File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 1117, in run_smallbaselineApp app.run(steps=inps.runSteps) File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 886, in run self.run_unwrap_error_correction(sname) File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 355, in run_unwrap_error_correction mintpy.cli.unwrap_error_phase_closure.main(iargs_closure) File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/cli/unwrap_error_phase_closure.py", line 153, in main run_unwrap_error_phase_closure(inps) File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/unwrap_error_phase_closure.py", line 418, in run_unwrap_error_phase_closure correct_unwrap_error_phase_closure( File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/mintpy/unwrap_error_phase_closure.py", line 358, in correct_unwrap_error_phase_closure cc = np.squeeze(f[ccName][i, :, :]) ~^^^^^^^^ File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper File "/home/jovyan/.local/envs/earthscope_insar/lib/python3.11/site-packages/h5py/_hl/group.py", line 357, in __getitem__ oid = h5o.open(self.id, self._e(name), lapl=self._lapl) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper File "h5py/h5o.pyx", line 190, in h5py.h5o.open KeyError: "Unable to open object (object 'connectComponent' doesn't exist)"

Asmit Roy Burman

unread,
Nov 7, 2023, 12:18:45 PM11/7/23
to MintPy

Posting the error message again for better understanding

Zhang Yunjun

unread,
Nov 8, 2023, 8:10:32 PM11/8/23
to MintPy
The closure_phase_bias.py does not require connected components to run, because it calls isce2 within mintpy to unwrap the generated closure phase, which can be time-consuming.

The regular smallbaselineApp.py, which does not include closure_phase_bias.py within it yet, does require the connected component for the phase unwrapping error correction step. Although advanced users could re-unwrap using isce2/snaphu for hyp3/gamma products, this is not built-in for smallbaselineApp.py.

I hope this clears things up.
Yunjun

Asmit Roy Burman

unread,
Nov 9, 2023, 12:30:43 PM11/9/23
to MintPy
Okay thanks for clearing it out. It was helpful.

Asmit

Reply all
Reply to author
Forward
0 new messages