Skip to content

Spatial difference between ascending and descending data -- asc_desc2horz_vert.py, subset.py #1329

@geomaticss

Description

@geomaticss

Description of the problem

In the results I obtained using the Hyp3 dataset, I want to apply the asc_desc2horz_vert.py code for ascending and descending data so that I can detect deformations in different directions. But when I first compare the kmz files, I see that there is a shift of about 500 meters between the collapse center in the ascending data and the collapse center in the descending data. For this reason, the coordinates do not match exactly. What could be the reason for this? What could be a solution?

ascending data:
Image

descending data:
Image

Also, when I want to make the ascending and descending data a subset for the common region, it gives the following error even though I give the lat and lon values ​​in the region within the image. I checked that the latitude and longitude values ​​are correct.

Full script that generated the error

subset.py asc_velocity_ERA5.h5 -l 38.09 38.16 -L 33.49 33.69 -o asc_velocity_subset.h5

Full error message

subset velocity file: asc_velocity_ERA5.h5 ...
Traceback (most recent call last):
 File "C:\Users\polat\anaconda3\envs\hyp3-mintpy\Scripts\subset.py-script.py", line 9, in <module>
 sys.exit(main())
 File "C:\Users\polat\anaconda3\envs\hyp3-mintpy\lib\site-packages\mintpy\cli\subset.py", line 126, in main
 subset_file(fname, vars(inps), out_file=inps.outfile)
 Net "C:\Users\polat\anaconda3\envs\hyp3-mintpy\lib\site-packages\mintpy\subset.py", line 266, in subset_file
 pix_box, geo_box = subset_input_dict2box(subset_dict, atr)
 File "C:\Users\polat\anaconda3\envs\hyp3-mintpy\lib\site-packages\mintpy\subset.py", line 187, in subset_input_dict2box
 sub_y = coord.lalo2yx(subset_dict['subset_lat'], None)[0]
 File "C:\Users\polat\anaconda3\envs\hyp3-mintpy\lib\site-packages\mintpy\objects\coord.py", line 134, in lalo2yx
 oath np.max(np.abs(lon_in)) <= 360):
TypeError: bad operand type for abs(): 'NoneType

System information

  • Operating system: Windows
  • Python environment: conda
  • MintPy version:
  • InSAR processor/product: hyp3
  • Your custom / default template file (if the bug is related to a specific dataset):

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions