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

Hello,
Let me step in this discussion.

@SUN-Ravenclaw If i well understand you have separated stacks of images paires and interfergrams.
You should export all your images+interferograms to the same directory : INSAR_20190208 that you create manually

On “Stamps Export” enter : Your_Working_Directory/INSAR_20190208
in the field “Target Folder” as shown bellow :


Use the same directory each time you will export a image+interferogram pair.
It would warn you that the directory exists already and it will be overwritten. Just ignore that and make YES. It will not overwrite the files but increment theme in each export

1 Like

thank you for commenting.
Another solution is to create a stack of all image pairs (AB, AC, AD), both for the backscatter and the interferogram, and use these two stacks in the StaMPS export.

It’s been a while i have not followed this topic.
Is one Stack for multiple slaves working fine now ?

now I’m not sure anymore either :grin:
I have only been using this python package for its automation: Snap2stamps package: a free tool to automate the SNAP-StaMPS Workflow
I’m not entirely sure how it handles the different pairs.

1 Like

I checked the snap2stamps scripts again and this is how it is solved there, too. So you were right: The export has to be repeated for each pair, but the output folder is the same.

3 Likes

In a time i did it manually for a huge amount of images.
Thanks to snap2stamps authors this come to be more easier.
@SUN-Ravenclaw you can also use this tool. it can help you a lot Snap2stamps package: a free tool to automate the SNAP-StaMPS Workflow

2 Likes

Hi Abdel, thank you for the answer above, I have completed steps 1 to 8 of StaMPS in matlab. However, I encountered some warning messages in step 6 (shown below).
Warning: A’s rank is not sufficient to reach the computer precision range.

In lscov (line 200)
In uw_sb_unwrap_space_time (line 304)
In uw_3d (line 161)
In ps_unwrap (line 235)
In stamps (line 503)
Do you know why? Thanks again.

Dear @ABraun ,
I tried what you recommended to do but I still cant get anu results! would you please take a look at these:

thank you

I don’t understand why you apply Create Stack when you have already applied Back Geocoding.

I found this tutorial in forum, with these steps, I may be misunderstood:
" 1) For S1, apply TOPS split and then apply precise orbit
2) Add all products to the Backgeocoding
3) Create Stack of Backgeocoded images (master_Stack)
Create a stack of Interferograms (Follow the step one of Prepare for STACK AVERAGING) on the resulting stack. { STEP 1 : master (SAME MASTER)+ (Different Slaves)first slave -> TOPS Corregistration, if necessary Deburst (and possibly subset), Interferogram, topographic phase removal, } master_Stack_ifg­_dinsar .
Deburst both the Backgeocoded stack (master_Stack_Deb), and the Interferogram stack (master_Stack_Deb_ifg¬_dinsar .).
5) Add an elevation band to the Interferogram stack (by right-click > add elevation raster.).
6) Add the Debursted stack and Debursted Interferogram stack to the stamps export "

yes, sorry. You are right. It was a misunderstanding from my side yesterday.
Do you use the same master for all coregistered stacks?

By the way, to add elevation and lat/lon rasters can now be done in the interferogram generation step already. Unfortunately, there are many outdated instructions on StaMPS preparation in this forum. This one is the latest I know: Documentation on SNAP-StaMPS workflow

We collected them once in this topic: About the StaMPS category

Please make sure you are referring to the latest version.

Yes, I use a 201912 image as master for each of coregistered stacks, I still don’t know whats wrong
thank you so much I’m gonna use this tutorial and check if I did sth by mistake

Hi all, the following figure is my map in GMT with the results of my DInSAR PSI processing…
using the ps_mean_v file!!! @mdelgado

2 Likes

Hi again, I want to check with details the contents of folders PATCH 1… 2… of others people here! can somebody to list this please??? I am checking errors presented during the processing…

if you are free , could you please share some experiences about how to use Global Mapper and Surface to draw the pictures from data in StaMPS.
Thank you !

Oh,a beautiful graph!could you tell how you create it ?Thank you !

Hello,
It is about saving your results and drawing each of result as shaded relief map into the Surfer program. It takes a bit time according to your data. But first you need to plot your result into Global mapper and export the data as Surfer .grd data. These programs are easy to use but not free. If you find a way, I can assist you to plot your results as same as mine.
Fikret

i have got these two programs.it is need to take me some time to learn.so if you are free ,i hope to get some advices and help.Thank you!

Hello again,
First of all, did you obtain result figures ? For example, mean_LOS velocity map? If you obtain that, you can export the values using matlab and you have LAT and LON values too. Then you need to creat LAT, LON, mean_LOS values respectively into the .dat file. If you make that tell me, afterwards it is time to use global mapper.
Fikret

hello!i have got the graph you mentioned!but ,sorry ,i don’t know where and how to export the result!i have read the mannual again but couldn’t find.could you give me some hints?thank you!