Workflow between SNAP and StaMPS

Does someone have an idea what goes wrong,
All my 146120 PS are in the red box, and I think that they should be distributed throughout the whole processing area?

interferogram represents whole processing area in StaMPS-SNAP,
and I have got PermanentScatterers only in the red box (white box is my AOI),
Thanks in advance,


Hi all,

I am trying to export two files: (Stack with ESD) and (Stack with ESD, Interferogram, DInSAR). But, I am having following error.

When I try StaMPS export for individual stacks, I get error that says “Input requires a coregistered stack (1st product) and an interferogram of at least 4 interferograms”. (Figure attached for reference).

Can anyone please tell me how to get rid of this error? Do I need to generate 4 interferograms again or am I doing something wrong?

Thank you for help!

both files must be inside the File list: the intensity stack and the interferogram stack

Can you please tell me how to do it? I moved both stacks in same folder. What else do I need to do?

I mean you have to add both to the ProductSet-Reader grafik at once

1 Like

Ok, it seems to be working. Thank you for help!

I will let you know if I have any other question. I believe I am going to bother you and other people here a lot! :sweat_smile:

this is what this forum is for!

The more precise you ask, the more likely are people which can help you :+1:

1 Like

Hello again!

I am having trouble understand the procedure for using stamps. I have used Stamps export and I have 4 directories. I have downloaded the folder with all graphs from Github. I also read the readme file and some other publications but I don’t understand what should I do next now. I also have terrain corrected stacks for ESD and ESD with Interferogram and DInSAR. Can you please tell me how should I process now?

Thank you!

the stamps export (or the python scripts snap2stamps) creates four directories:

  • dem
  • geo
  • rslc
  • dif0

After you downloaded the stamps binaries, you have to make sure that you modify the StaMPS_CONFIG file to match all required folders (MATLABPATH, TRIANGLE, SNAPHU) and source it as described in the manual (chapter 2.1) to make all files available to the Linux shell.

Then you create a new folder and run the script mt_prep_snap as described in the manual (chapter 4). This creates interferograms as they are required by matlab. Once this script is executed, you open matlab and proceed with the actual stamps processing, basically just entering the commands stamps(1,1), stamps(2,2), and so on.

1 Like

Is this process same for doing Pre-processing in SNAP 7.0?

no, as stated in the overview, this description is outdated.

Please follow one of them:


Hello again!

I am processing stamps(2,2). I got the error in line 221 of ps_est_gamma_quick.m that says “Subscript indices must either be real positive integers or logicals.” Can you please tell me what is wrong and how to solve it?

Thank you!

hoch many patches did you create?
according to this answer, it happens for very small patches

I used 1 patch for range and azimuth, which was set by default. What value should I set in mt_prep_snap?

then the size of the patch is not the problem.

How did you create the latitude and longitude bands?

I created the lat long bands in SNAP in topo phase removal step while pre-processing.

hmm, also no source of error then. I have no idea to be honest.

the variable grid(i,j) is having some NaN values. So, we have currently used if ~isnan condition to bypass the error. Stamps(2,2) seems to work for now.

had you used the ps_load_initial_gamma.m that I have put in the forum in March 2019? This solved not only for NaN.

Good that you solved.

1 Like

could you please give a link to this post? I will include it in the overview post.