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

Different result between ascending and descending #334

Open
NabilAttar opened this issue Mar 1, 2024 · 26 comments
Open

Different result between ascending and descending #334

NabilAttar opened this issue Mar 1, 2024 · 26 comments

Comments

@NabilAttar
Copy link

NabilAttar commented Mar 1, 2024

Greeting Mr. Yumorishita

I have done to process frame 171A_09587_041313 and 120D_09443_131313 in 105.16/105.36/-5.53/-5.31 geo clip area, but it produces different velocity. In ascending we see that there is an area with up lift and subsidence. But in descending we see there're no area with subsidence, almost all areas experienced uplift. Is my processing wrong? Why does it happen?

Thank you.....

Ascending
ascending2
descending

@yumorishita
Copy link
Owner

Maybe the time period is different?

@NabilAttar
Copy link
Author

Maybe the time period is different?

Yes sir, ascending in 2015 - 2023. Descending 2017 - 2023. does it have any effect?

@yumorishita
Copy link
Owner

I think so. Also, you should check the network and gaps.

@NabilAttar
Copy link
Author

I think so. Also, you should check the network and gaps.

in the normal case, should they both have the same velocity? Thank you

@yumorishita
Copy link
Owner

They are never same because the LOS geometry is different, but they should be similar if there is no horizontal component.

@NabilAttar
Copy link
Author

NabilAttar commented Mar 3, 2024

greeting Mr. Yumorishita

this is network between ascending (up) and descending (down). i think that descending is bad because infrenquently IFG, Do you have any suggestions regarding processing results that are not similar? thank you

ascending network13_nobad
descending network13_nobad

@yumorishita
Copy link
Owner

I suggest using the common time period (Aug 2019 - May 2022).

@NabilAttar
Copy link
Author

I have tried time period (Aug 2019 - May 2022) in descending. However, there is still only uplift in the surrounding area. Maybe it's already from the IFG data

@NabilAttar
Copy link
Author

excuse me sir, sorry for out of topic. I have done ascending and descending processing. but I think the results of descending are quite bad. i have use LiCSBAS_decompose but i think the result of UD is bad caused descending frame. Is it possible to get a UD velocity value just by using ascending? I only need UD and not EW. thank you sir

@ditafaith
Copy link

Hi @NabilAttar
If you want to use only 1 orbit, and ensure that horizontal movement very small contribution, then you may extract time series at certain point and divided by cosinus incidence angle. If you want to get an image, may be you can operate gdal_calc for ascending.tiff file.

@NabilAttar
Copy link
Author

Hi @NabilAttar If you want to use only 1 orbit, and ensure that horizontal movement very small contribution, then you may extract time series at certain point and divided by cosinus incidence angle. If you want to get an image, may be you can operate gdal_calc for ascending.tiff file.

thank you sir, is gdal_calc same like raster calculator in QGIS? so, formula for get velocity upward/vertical direction tiff is only [velocity LOS.tiff / cos(incidence angle] ?? thank you

@NabilAttar
Copy link
Author

Hi @NabilAttar
If you want to use only 1 orbit, and ensure that horizontal movement very small contribution, then you may extract time series at certain point and divided by cosinus incidence angle. If you want to get an image, may be you can operate gdal_calc for ascending.tiff file.

excuse me sir, is my statement correct? Thank you

@yumorishita
Copy link
Owner

I think your statement is correct.

@NabilAttar
Copy link
Author

Thank you sir

@NabilAttar
Copy link
Author

ascending network13_nobad descending network13_nobad

greeting Mr. Yumorishita

I want to ask about ascending and descending baselines which I have uploaded. In ascending use dummy unit (look in left side picture), and descending use Bperp unit (baseline perpendicular). why the results are different? in your paper the network result in Bperp unit, what is dummy? thank you

@yumorishita
Copy link
Owner

If baselines file does not contain the values for all dates, the dummy baseline is used because the correct figure cannot be drawn from the incomplete data. Please check the log in step02 to find what date is lacking in baselines.

print('ERROR: bperp for {} not found!'.format(imd), file=sys.stderr)

print(' baselines file found, but not complete. Make dummy', flush=True)

@NabilAttar
Copy link
Author

image

excuse sir, here is "not baseline found" do you have any suggestions? is it a good idea if i delete 20170814 in GEOC and start batch licsbas from step 02? thank you

@NabilAttar
Copy link
Author

image
in metadata network png on GEOC, baseline 20170814 is avaible

@yumorishita
Copy link
Owner

What is the frame ID?

is it a good idea if i delete 20170814 in GEOC and start batch licsbas from step 02?

Yes, you will get the correct network figure with non-dummy y-axis.

@NabilAttar
Copy link
Author

frame ID is 171A_09587_041313

image

i have tried by deleting 20170814 in GEOC and delete GEOCml1. but another IFG without baseline is found. do you have any suggestions? thank you

@yumorishita
Copy link
Owner

This is the baseline file used for the network figure:
https://gws-access.jasmin.ac.uk/public/nceo_geohazards/LiCSAR_products/171/171A_09587_041313/metadata/baselines

I found some dates have no data in the 3rd column, e.g.,

20190512 20170802 -65.0 -648
20190512 20170814  -636
20190512 20170826 -52.0 -624

You can manually edit the file and fill the 3rd column, or ask COMET to fix this issue.

@NabilAttar
Copy link
Author

how to determine the value of the 3rd column? are there any provisions for fill the 3rd column?

@ditafaith
Copy link

how to determine the value of the 3rd column? are there any provisions for fill the 3rd column?

You might try to download the images and check it manually using SNAP software. Also, you can check it through Alaska Satellite Portal. It seems you are working on subsidence in Lampung.

@NabilAttar
Copy link
Author

how to determine the value of the 3rd column? are there any provisions for fill the 3rd column?

You might try to download the images and check it manually using SNAP software. Also, you can check it through Alaska Satellite Portal. It seems you are working on subsidence in Lampung.

yes sir, i have research about this area. also i have visited alaska and check baseline manually. I equated it with baseline data from licsar, the perpendicular baseline difference from the missing third column is only about 3 meters. Thank you for the advice

@luvkesh5511
Copy link

Dear Prof.,
I had a doubt regarding how t o use a combination of ascending and descending pass interferograms to analyse ground deformation? Does the platform support combining two different passes (Asc and Desc) ? If so, how can we integrate them, should we add two frame ids in line 72: (p01_frame = "129A_xxxxx" "135D_xxxxx") ? Kindly guide.

Thank you

@yumorishita
Copy link
Owner

  1. Run LiCSBAS step01-16 for the ascending frame (129A_xxxxx) using batch_LiCSBAS.sh
  2. Run LiCSBAS step01-16 for the descending frame (135D_xxxxx) using batch_LiCSBAS.sh
  3. Use LiCSBAS_decomposeLOS.py with the input from both ascending and descending

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants