Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

hyp3_insar_stack_for_ts_analysis.ipynb permutes corners #469

Open
feigl opened this issue Sep 21, 2024 · 2 comments
Open

hyp3_insar_stack_for_ts_analysis.ipynb permutes corners #469

feigl opened this issue Sep 21, 2024 · 2 comments
Labels
bug Something isn't working

Comments

@feigl
Copy link

feigl commented Sep 21, 2024

In the example notebook https://github.com/ASFHyP3/hyp3-docs/blob/main/docs/tutorials/hyp3_insar_stack_for_ts_analysis.ipynb
The function named get_common_overlap permutes the corners.
This issue appeared in an ascending swath where the area of interest was slightly off the western edge, i.e. with longitude between the longitudes of the NW and SW corners. This case is rare.

Here is the corrected code:

    # 2024/09/21 This looks strange 
    # ulx = max(corner['upperLeft'][0] for corner in corners)
    # uly = min(corner['upperLeft'][1] for corner in corners)
    # lrx = min(corner['lowerRight'][0] for corner in corners)
    # lry = max(corner['lowerRight'][1] for corner in corners)
    
    # 2024/09/24 x coordinate is UTM easting, mininum is to our left, i.e. west edge
    ulx = min(corner['upperLeft'][0] for corner in corners)
    uly = max(corner['upperLeft'][1] for corner in corners)
    lrx = max(corner['lowerRight'][0] for corner in corners)
    lry = min(corner['lowerRight'][1] for corner in corners)

search3utm

@feigl feigl added the bug Something isn't working label Sep 21, 2024
@asjohnston-asf
Copy link
Member

The get_common_overlap() function returns the intersection of the input geometries. Your revisions would instead return the union of the input geometries. Either is a reasonable approach, and I agree the behavior of the function could be made more clear in the docstring.

Using the intersection of the input interferograms restricts the time series to the region where every input interferogram contributes and guarantees a fully-connected network for the entire region.

Using the union of the input interferograms expands the time series to the region where any input interfergoram contributes, but outlying areas not covered by all interferograms (like your example) will have a less-connected network through time. You could get similar results by removing interferograms not covering your area of interest from your list prior to section #3 of the notebook.

@asjohnston-asf
Copy link
Member

If your area of interest is small, burst-based InSAR may be a better fit than SLC-based InSAR. Sentinel-1 bursts are consistently geolocated through time, eliminating the case where your area of interest is covered by some interferograms in your stack but not by others.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants