How to prepare Sentinel-1 images stack for PSI/SBAS in SNAP

you can check any parameter’s value with getparm and change it with setparm


https://homepages.see.leeds.ac.uk/~earahoo/stamps/StaMPS_Manual_v3.1.pdf (chapter 6)

About how to change the weeding thresholds:


https://homepages.see.leeds.ac.uk/~earahoo/stamps/StaMPS_Manual_v3.1.pdf (chapter 4)

Oh, that’s weed_time_win in the codes now.
Thanks for the tips, I try that and I let you know if any improvement.

Hi,

now i have GUI (a whole new world :smile:), and tried to plot some coordinates from the ps2.mat->lonlat file:
I have 259101 PS points, the values in the first column are 0 to 2,13(lon), and 0 to 0,068 in the second column (lat). The lower left corner is the “origo” [0,0]. Here is a screenshot:
k%C3%A9p

What should I check now?
Thank you for your answers :slight_smile:
Andrew

If you proceeded to at least step 6, you can just plot the whole map with
ps_plot('v-d')

Ok, that is what i got:

as you described, the lat/lon information wasn’t transferred correctly.
I hoped this was not the case for the GUI part.

I have no idea at which step the geocoding went lost, to be honest. If you exported lat/lon bands with at the interferogram step and the StaMPS export didn’t give any warning it should be alright.

Some thoughts:

  • Have you set the processor to ‘gamma’ before starting the matlab processing?
    You can check with getparm, sometimes it is still ‘doris’. I don’t know if this caused the error but it could be worth a try.
  • Have you replaced the matlab script and the shell script provided by catherine here: Workflow between SNAP and StaMPS
    Make sure you have a copy of the old ones and then rename them accordingly so they are addressed under the original name.

Thank you for your help, ABraun.

I had set the insar_processor to ‘gamma’ (checked it), and i replaced the matlab and shell scripts to the new one (checked it). I had to be edit the ps_load_initial_gamma file at line 180 and 183, i added the ‘-v7.3’ parameter, because of a well known file save issue.

The Stamps export didn’t drop any error or warning, and i got all off the data and dir that was described before in this topic (dem, diff, geo, rslc). I’m just copied to a new dir, which is the input of the mt_prep_gamma. The structure of the input dir:

  • dem: contains 2 files (rslc and par), from one of the stamps output dir (chosed randomly, becauseall of the same)
  • diff0: contains 24 base, 24 diff and 24 par files (i had 24 stacks).
  • geo: contains 6 files, 20160519.lat, 20160519.lon, 20160519.lat.par, 20160519.lon.par, elevation_dem.rdc, elevation_dem.rdc.par (20160519 is the master date).
  • rslc: contains 25 rslc and 25 par files (i had 25 images).

Maybe something went wrong, when i copied the Stamps outputs (geo directory)?
I ran the whole STAMPS process under terminal, without GUI, so i didn’t checked the results after the steps. Next time i will run from mt_prep_gamma to stamps(8,8), step by step, and i will check the outputs in matlab.
(if there is a better topic for my problem, please let me know :slight_smile: )

Andrew

Seems you have done nothing wrong - good luck! :slight_smile:

You should had checked the log information coming out from the mt_prep_gamma command.

What you show there is seems that the lat/lon were not considered… even if then inside MATLAB you got some coordinates. Stilll which are those values 1.0 and 0.00 ??? Are those coordinates matching with your AOI?

Just wondering…

For what you said I am not sure the mt_prep_gamma got the right values for those coordinates.

Btw, had you used the mt_prep_gamma after modifying the gamma format to read? Please use the mt_prep_snap which is provided with the latest stamps version available on github.

Let us know.
Good luck!

the image is broken.

Hi guys !

I have retreated my data by playing with weeding parameters and dropping noisy interferograms.
But I still face the same problem with an important dispersion in my time series.

May this problem come from SNAP processing?

did you maybe multi-look your data?

If you want to go sure everything is fine, I can recommend the scripts provided by @mdelgado here: Snap2stamps package: a free tool to automate the SNAP-StaMPS Workflow

You simply place all slave zip files into a slave folder, create a master folder where you apply orbit files and TOPS Split and then the scripts do the rest for you including the export. I was able to get nice results with it.

No, I have not.

Oh yeah, I have read that, it really helps to process quicker your data.

But, from what I understood, the only difference in the processing chain is the “enhance spectral diversity” step.
Else if, it seems that I have exactly done the same path.

Well, the dispersion comes from the APS of the slaves.
I have compared SNAP derived interferograms with StaMPS and other software and results were similar. With Sentinel-1 the issue comes on the APS, which should be taken carefully into account.

Had you done anything to reduce it?

No, I did not.
I assumed there was no effect because my area of study is flat, South Alsatian plane, and small, less than 50 km per 50 km.
So, yes, may be this assumption is not right.

But, I did not applied the ESD step, what is the impact of missing this step?

The Netherlands is also really flat, but atmosphere believe me affects the results.

That means, atmosphere only can be neglected if the area covers few Km2 (1km x 1km or 2km x 2km) or big deformation is expected (Earthquake events, or huge subsidence expected)
From my personal experience atmosphere matters, mainly when deformation to measure is small. If you measure deformation in the Mexico city center, atmosphere effects can be neglected.

Cheers

2 Likes

It happened to me when I used the old way (using SNAP GUI). Now I’m using the snap2stamps pack made by @mdelgado and the latlon and dem inputs are always ok!
Regards

1 Like

I am a freshman at SNAP&StaMPS. I have learned a lot from here, thank you all.
I successfully ran mt_prep_snap by cygwin64 (using a modified version provided by someone else on the web). But after running mt_prep_snap, the matlab automatically jumped out, said ""Warning: Name is nonexistent or not a directory: \StaMPS-4.1-beta\matlab:


After that, I opened the matlab from the terminal and show as the second picture.

I don’t know what to do. I hope that someone can help me,please . Thank you very much.
ps. my matlab is R2016a in windows
Jamie
2018.09.25

The link does not work for me too, please share another link.
Alternatively, can I share my email for you to send?

Your screenshots are too small to be readable.