insarlab / MintPy

Miami InSAR time-series software in Python
https://mintpy.readthedocs.io
Other
608 stars 259 forks source link

Bug report when processing hyp3 with bridging method in correct_unwrap_error #832

Closed Arctic-Ambrun closed 2 years ago

Arctic-Ambrun commented 2 years ago

Description of the problem

Seems there`s no [ connectComponent ] when using bridging method with data build with hyp3 products. https://github.com/insarlab/MintPy/blob/36a083521b39d9fb557741f9c5d703b7e42448c3/mintpy/unwrap_error_bridging.py#L243

the ccName need to be transfered as 'unwrapPhase' above the Line 243

My suggestion is to add the judgment of hyp3 product, or maybe the issue due to the generation of ifgram.h5

Code:

correct unwrapping error in /mnt/Disk2/0_ASF_SBAS_WOLF/asc_whole_mintpy/inputs/ifgramStack.h5 with bridging ... read water mask from file: waterMask.h5 open /mnt/Disk2/0_ASF_SBAS_WOLF/asc_whole_mintpy/inputs/ifgramStack.h5 with r+ mode input dataset: unwrapPhase output dataset: unwrapPhase_bridging create /unwrapPhase_bridging of np.float32 in size of (1026, 2080, 1941) Traceback (most recent call last): File "/home/wjg/anaconda3/envs/AM_manual/bin/smallbaselineApp.py", line 33, in sys.exit(load_entry_point('mintpy==1.3.3', 'console_scripts', 'smallbaselineApp.py')()) File "/home/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/mintpy-1.3.3-py3.8.egg/mintpy/smallbaselineApp.py", line 1275, in main app.run(steps=inps.runSteps) File "/home/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/mintpy-1.3.3-py3.8.egg/mintpy/smallbaselineApp.py", line 1052, in run self.run_unwrap_error_correction(sname) File "/home/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/mintpy-1.3.3-py3.8.egg/mintpy/smallbaselineApp.py", line 528, in run_unwrap_error_correction mintpy.unwrap_error_bridging.main(iargs_bridge) File "/home/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/mintpy-1.3.3-py3.8.egg/mintpy/unwrap_error_bridging.py", line 300, in main run_unwrap_error_bridge(inps.ifgram_file, File "/home/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/mintpy-1.3.3-py3.8.egg/mintpy/unwrap_error_bridging.py", line 244, in run_unwrap_error_bridge 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/wjg/anaconda3/envs/AM_manual/lib/python3.8/site-packages/h5py/_hl/group.py", line 328, 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)"

PASTE CODE HERE

Full error message

PASTE ERROR MESSAGE HERE

System information

welcome[bot] commented 2 years ago

đź‘‹ Thanks for opening your first issue here! Please filled out the template with as much details as possible. We appreciate that you took the time to contribute! Make sure you read our contributing guidelines.

Arctic-Ambrun commented 2 years ago

Hello After changing the dataset name from 'connectComponent' to 'unwrapPhase' , another issue comes out. input reference point is NOT included in the connectComponent. seems the input reference point is not in the connectComponent, I wonder how to deal with it, thanks

Arctic-Ambrun commented 2 years ago

@yunjunz Dear Yunjun, May I ask a question? Does hyp3 product contains conn comp file? It is werid that some of the bridging file is written but some were failed. An image was uploaded named unwrapPhase_bridging_1.pn unwrapPhase_bridging_1 g

Arctic-Ambrun commented 2 years ago

maybe I`m wrong, the hyp3 product does not contains connectComponent data, so for HYP3 The process is stoped??

yunjunz commented 2 years ago

HyP3 InSAR products are generated using Gamma, which does not provide connected components info. Thus, the phase unwrapping error correction methods can not be used.

The same constraint applies to Gamma, GMTSAR also, as far as I am aware. https://mintpy.readthedocs.io/en/latest/dir_structure/#asf_hyp3

Arctic-Ambrun commented 2 years ago

Thank you for your kindhearted help, Mr.Yunjun