open them in SNAP and check if you see backscatter intensities.

the band name already says it - Intensity

If the image can be displayed, everything is fine. But sometimes, coregistration ( BackGeocoding) fails and there is no slave image. No slave image means no interferogram. No interferogram means no displacement.

I recommend getting into the basics first so you have an idea about what you expect from the analysis and wher potential errors lie.

Thank you so much ABraun.

stamps(2,2)

STAMPS: ########################################

STAMPS: ####### StaMPS/MTI Version 4.0b6 #######

STAMPS: ####### Beta version, Jun 2018 #######

STAMPS: ########################################

STAMPS: Will process current directory only

psver currently: 1

psver now set to: 1

STAMPS: ########################################

STAMPS: ################ Step 2 ################

STAMPS: ########################################

STAMPS: Directory is PATCH_1

PS_EST_GAMMA_QUICK: Starting

PS_EST_GAMMA_QUICK: Estimating gamma for candidate pixels

GETPARM: filter_grid_size=50

GETPARM: filter_weighting=‘P-square’

GETPARM: clap_win=32

GETPARM: clap_low_pass_wavelength=800

GETPARM: clap_alpha=1

GETPARM: clap_beta=0.3

GETPARM: max_topo_err=20

GETPARM: lambda=0.0554658

GETPARM: gamma_change_convergence=0.005

GETPARM: gamma_max_iterations=3

GETPARM: small_baseline_flag=‘n’

Found look angle file

PS_EST_GAMMA_QUICK: n_trial_wraps=0.243690

PS_EST_GAMMA_QUICK: Initialising random distribution…

PS_EST_GAMMA_QUICK: 1884150 PS candidates to process

PS_EST_GAMMA_QUICK: iteration #1

PS_EST_GAMMA_QUICK: Calculating patch phases…

Index in position 1 is invalid. Array indices must be positive integers or logical values.

Error in ps_est_gamma_quick (line 221)

ph_grid(grid_ij(i,1),grid_ij(i,2),:)=ph_grid(grid_ij(i,1),grid_ij(i,2),:)+shiftdim(ph_weight(i,:),-1);

Error in stamps (line 326)

ps_est_gamma_quick(est_gamma_parm);

please tell me the solution.

Thank you

this error was reported a couple of times, please use the search function.

One solution was to replace the script ps_load_initial_gamma in the stamps folder with a modified one by @mdelgado

But make sure you save a copy first, before you replace the script

I already replaced this file long back…

sorry, I forgot about it. I’m afraid I have no solution then. Maybe it is better to directly ask in the StaMPS forum

currently, this is not supported by SNAP, but you can open the InSAR Stack Overview and export the baseline information of the stack to a csv file and generate the chart outside SNAP

Thank you ABraun.

All **ps_plot(‘u’) , ps_plot(‘u-d’), ps_plot(‘u-m’) and ps_plot(‘u-dm’)** are running successfully. But while running this below command it is showing error, please check it once. Thank you.

ps_plot(‘v-dm’,‘ts’);

Error using ps_plot (line 1540)

unknown value type

ps_output

Writing output files…

GETPARM: small_baseline_flag=‘n’

GETPARM: ref_velocity=0

GETPARM: lambda=0.0554658

1241101 ref PS selected

Error using load

‘mean_v’ is not found in the current folder or on the MATLAB path, but exists in:

/mnt/e/subsidence/ahmedabad_stamps/6_psiexport

Change the MATLAB current folder or add its folder to the MATLAB path.

Error in ps_output (line 95)

meanv=load(meanvname);

v is the acutal velocity, so it should be calculated by the scripts. Did you proceed until stamps(8,8)?

Yes, I processed upto stamps(1,8).

ps_plot(‘v’) and ps_plot(‘v’,‘ts’) are also running properly. but while running ps_plot(‘v-dm’,‘ts’) it is giving error. Thank you

ps_plot(‘v-dm’,‘ts’);

Error using ps_plot (line 1540)

unknown value type

m is the atmospheric orbit error (AOE) of the master. If you did not remove it, you cannot include it in the plot.

I rerun the stamps(8,8) step once again, but still I’m getting the error. Please tell me the solution. Thank you.

ps_plot(‘v-dm’,‘ts’);

Error using ps_plot (line 1540)

unknown value type

ps_plot(‘v-dm’,1);

Error using ps_plot (line 1540)

unknown value type

stamps(6,6)

STAMPS: ########################################

STAMPS: ####### StaMPS/MTI Version 4.0b6 #######

STAMPS: ####### Beta version, Jun 2018 #######

STAMPS: ########################################

STAMPS: Will process current directory only

STAMPS: ########################################

STAMPS: ################ Step 6 ################

STAMPS: ########################################

STAMPS: Directory is PATCH_1

PS_UNWRAP: Starting

Phase-unwrapping…

GETPARM: small_baseline_flag=‘n’

GETPARM: unwrap_patch_phase=‘n’

GETPARM: scla_deramp=‘n’

GETPARM: subtr_tropo=‘n’

GETPARM: tropo_method=‘a_l’

GETPARM: drop_ifg_index=

GETPARM: unwrap_hold_good_values=‘n’

PS_UNWRAP: Code to hold good values skipped

subtracting scla and master aoe…

GETPARM: unwrap_time_win=730

GETPARM: unwrap_method=‘3D’

GETPARM: unwrap_grid_size=200

GETPARM: unwrap_gold_n_win=32

**GETPARM: unwrap_prefilter_flag='y’**

GETPARM: unwrap_gold_alpha=0.8

GETPARM: unwrap_la_error_flag=‘y’

GETPARM: unwrap_spatial_cost_func_flag=‘n’

GETPARM: max_topo_err=20

GETPARM: lambda=0.0554658

PS_UNWRAP: n_trial_wraps=0.256749

Resampling phase to grid…

Number of interferograms : 52

Number of points per ifg : 1241101

Number of resampled points: 52911

**GETPARM: unwrap_prefilter_flag='y’**

In my case it is “y” so the subtraction of SCLA and master AOE is possible. Is it correct?

I think you can proceed without the m flag then, yes.

I didn’t get your point?