geodesymiami / MiNoPy

MIami NOn linear phase linking in PYthon
Other
19 stars 5 forks source link

shadowMask missing issue #47

Closed stoormgeo closed 3 years ago

stoormgeo commented 3 years ago

Hi, I am processing a CSM Stripmap dataset and encountered the error of missing shadowMask in the generate_ifgram and unwrap_ifgram steps. As far as I understood that shadowMask is not a requirement? The ISCE preprocessor does not prepare such data during processing (in the geom_reference folder), unlike Sentinel output. What is the best way to solve? Thx a lot. Vit

I'm using the last version MiNoPy release version V0.2.0-7, release date 2021-10-29

minopyApp.py /home/insar/1/Fast/YLPG21/test.template --dostep generate_ifgram --dir /home/insar/1/Fast/YLPG21/proc/merged/minopy

20211107:185436 * generate_unwrap_mask.py --geometry /home/insar/1/Fast/YLPG21/proc/merged/minopy/inputs/geometryRadar.h5
input geometry file: /home/insar/1/Fast/YLPG21/proc/merged/minopy/inputs/geometryRadar.h5
read /home/insar/1/Fast/YLPG21/proc/merged/minopy/inputs/geometryRadar.h5 shadowMask
Traceback (most recent call last):
File "/insar/MiNoPy/minopy/generate_unwrap_mask.py", line 111, in
main()
File "/insar/MiNoPy/minopy/generate_unwrap_mask.py", line 39, in main
mintpy.generate_mask.main(args_shm.split())
File "/insar/MintPy/mintpy/generate_mask.py", line 318, in main
inps.outfile = create_threshold_mask(inps)
File "/insar/MintPy/mintpy/generate_mask.py", line 170, in create_threshold_mask
data, atr = readfile.read(inps.file, datasetName=inps.dset)
File "/insar/MintPy/mintpy/utils/readfile.py", line 235, in read
data = read_hdf5_file(fname,
File "/insar/MintPy/mintpy/utils/readfile.py", line 314, in read_hdf5_file
raise ValueError('input dataset {} not found in file {}'.format(datasetName, fname))
ValueError: input dataset ['shadowMask'] not found in file /home/insar/1/Fast/YLPG21/proc/merged/minopy/inputs/geometryRadar.h5

minopyApp.py /home/insar/1/Fast/YLPG21/test.template --dostep unwrap_ifgram --dir /home/insar/1/Fast/YLPG21/proc/merged/minopy

snaphu -f /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/interferograms_single_reference/20190814_20210620/config_all -d /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/interferograms_single_reference/20190814_20210620/filt_fine.int 1302 -o /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/interferograms_single_reference/20190814_20210620/filt_fine.unw
b"Can't open file /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/mask_unwrap\nAbort\n"
Time spent: 0.0018973906834920248 m
20211107:185420 * unwrap_ifgram.py --ifg /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/interferograms_single_reference/20190814_20201007/filt_fine.int --coherence /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/quality --unwrapped_ifg /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/interferograms_single_reference/20190814_20201007/filt_fine.unw --max_discontinuity 1.2 --init_method MCF --l ength 1444 --width 1302 --height 633730.1420114406 --num_tiles 1 --earth_radius 6393121.60859229 --wavelength 0.031228381041666666 -m /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/mask_unwrap
2

mirzaees commented 3 years ago

Hi @stoormgeo , Thank you for letting me know, I will correct this issue for the next release but now if you don't have the shadow mask, the easiest way is to create a mask with all values equal to 1: generate_mask.py inverted/quality -m -1 -o shadow_mask.h5 After that you need to run the first command in run_03 : generate_unwrap_mask.py ...

And then you can unwrap without problem

stoormgeo commented 3 years ago

It seems that shadow_mask should be in the geometry stack. I tried using -m flag with the same result. I will try to use add.py for combine mask and geometry files.

ValueError: input dataset ['shadowMask'] not found in file /home/insar/1/Fast/YLPG21/proc/merged/minopy/inverted/shadow_mask.h5

Hi @stoormgeo , Thank you for letting me know, I will correct this issue for the next release but now if you don't have the shadow mask, the easiest way is to create a mask with all values equal to 1: generate_mask.py inverted/quality -m -1 -o shadow_mask.h5 After that you need to run the first command in run_03 : generate_unwrap_mask.py ...

And then you can unwrap without problem

mirzaees commented 3 years ago

@stoormgeo , you may also comment this line:

https://github.com/geodesymiami/MiNoPy/blob/c8e7abbbd51c5e26c743229dfecfdd364a2b291e/minopy/generate_unwrap_mask.py#L39

and it should work with -m

stoormgeo commented 3 years ago

Yep it works! Thx

@stoormgeo , you may also comment this line:

https://github.com/geodesymiami/MiNoPy/blob/c8e7abbbd51c5e26c743229dfecfdd364a2b291e/minopy/generate_unwrap_mask.py#L39

and it should work with -m