Excuse me, I am not familiar enough with Matlab. How can I print xy or bperp_mat?

Hi,sir. I meet a new problem, when I used overview, the result looked puzzling. please look at the place I marked.

Then I made Sep17 as the master image, the congestion result have same error

I use snap6.0 under unbutu. Does anyone have this problem?

Good time!

I tried to save xy, bperp_mat, index_NaN as binary files. But I am not sure it is wery usefull.

size of files

xy - 66 693 969

bperp_mat 133 387 938

index_NaN - 61 186

Hello everyone.

dear all user,

after all processing complete, now im facing new problem. the problem happened with stamps and vistamps

in stamps :

viStaMPS

Writing output files…

GETPARM: small_baseline_flag=‘n’

55886 ref PS selected

Error using load

Unable to read file ‘mean_v’: no such file or directory.

Error in ps_output (line 92)

meanv=load(meanvname);

ps_plot(‘v’)

PS_PARMS_DEFAULT: max_topo_err = 5

PS_PARMS_DEFAULT: quick_est_gamma_flag = y

PS_PARMS_DEFAULT: filter_grid_size = 50

PS_PARMS_DEFAULT: filter_weighting = P-square

PS_PARMS_DEFAULT: gamma_change_convergence = 0.005

PS_PARMS_DEFAULT: slc_osf = 1

PS_PARMS_DEFAULT: clap_win = 32

PS_PARMS_DEFAULT: clap_low_pass_wavelength = 800

PS_PARMS_DEFAULT: clap_alpha = 1

PS_PARMS_DEFAULT: clap_beta = 0.3

PS_PARMS_DEFAULT: select_method = DENSITY

PS_PARMS_DEFAULT: density_rand = 20

PS_PARMS_DEFAULT: percent_rand = 20

PS_PARMS_DEFAULT: gamma_stdev_reject = 0

PS_PARMS_DEFAULT: weed_time_win = 730

PS_PARMS_DEFAULT: weed_max_noise = Inf

PS_PARMS_DEFAULT: weed_standard_dev = 1

PS_PARMS_DEFAULT: weed_zero_elevation = n

PS_PARMS_DEFAULT: weed_neighbours = y

PS_PARMS_DEFAULT: unwrap_method = 3D

PS_PARMS_DEFAULT: unwrap_patch_phase = n

PS_PARMS_DEFAULT: drop_ifg_index = []

PS_PARMS_DEFAULT: unwrap_la_error_flag = y

PS_PARMS_DEFAULT: unwrap_spatial_cost_func_flag = n

PS_PARMS_DEFAULT: unwrap_prefilter_flag = y

PS_PARMS_DEFAULT: unwrap_grid_size = 200

PS_PARMS_DEFAULT: unwrap_gold_n_win = 32

PS_PARMS_DEFAULT: unwrap_alpha = 8

PS_PARMS_DEFAULT: unwrap_time_win = 730

PS_PARMS_DEFAULT: unwrap_gold_alpha = 0.8

PS_PARMS_DEFAULT: scla_drop_index = []

PS_PARMS_DEFAULT: scn_wavelength = 100

PS_PARMS_DEFAULT: scn_time_win = 365

PS_PARMS_DEFAULT: scn_deramp_ifg = []

PS_PARMS_DEFAULT: ref_lon = -Inf Inf

PS_PARMS_DEFAULT: ref_lat = -Inf Inf

PS_PARMS_DEFAULT: ref_centre_lonlat = 0 0

PS_PARMS_DEFAULT: ref_radius = Inf

PS_PARMS_DEFAULT: plot_dem_posting = 90

PS_PARMS_DEFAULT: plot_scatterer_size = 120

PS_PARMS_DEFAULT: plot_pixels_scatterer = 3

PS_PARMS_DEFAULT: plot_color_scheme = inflation

PS_PARMS_DEFAULT: shade_rel_angle = 90 45

PS_PARMS_DEFAULT: lonlat_offset = 0 0

PS_PARMS_DEFAULT: merge_resample_size = 0

PS_PARMS_DEFAULT: merge_standard_dev = Inf

PS_PARMS_DEFAULT: scla_method = L2

PS_PARMS_DEFAULT: scla_deramp = n

PS_PARMS_DEFAULT: insar_processor = doris

PS_PARMS_DEFAULT: subtr_tropo = n

PS_PARMS_DEFAULT: tropo_method = a_l

115278 ref PS selected

Error using sum

Trailing string input must be ‘double’,‘native’, or ‘default’.

Error in mean (line 82)

y = sum(x,dim,flag)/size(x,dim);

Error in ps_plot (line 1078)

ph_uw=ph_uw-repmat(mean(ph_uw(ref_ps,:),1,‘omitnan’),n_ps,1);

Please help me to solved this problem.

Im using ubuntu 16.04 and matlab 2014b ubuntu.

Hi @thho

How did you manage to plot colored dots on a map? In other words, what is the command you used in StaMPS to get these dots to show up? Mine are just square pixels with displacement values.

Thanks,

Chong

Hi @hchong005,

in order to visualize the data in a manner where you can see and click single PS and have the deformation time series plotted, I exported data from matlab with this commands

```
load parms.mat;
ps_plot('v-dos', -1)
load ps_plot_v-dos.mat
lon2_str = cellstr(num2str(lon2));
lat2_str = cellstr(num2str(lat2));
lonlat2_str = strcat(lon2_str, lat2_str);
lonlat_str = strcat(cellstr(num2str(lonlat(:,1))), cellstr(num2str(lonlat(:,2))));
ind = ismember(lonlat_str, lonlat2_str);
disp = ph_disp(ind);
export_res = [lon2 lat2 disp ts];
metarow = [ref_centre_lonlat NaN transpose(day)-1];
k = 0;
export_res = [export_res(1:k,:); metarow; export_res(k+1:end,:)];
export_res = table(export_res);
writetable(export_res,'stamps_tsexport.csv')
```

The exported .csv was than imported to a R Shiny application I am currently developing, that is what you see on my screenshots. The application allows you to click single dots to look up the ts for the dot and compare up to 5 ts in one graph…the status of this app is somewhat beta and I am changing and adding features all the way, give me a few days and I can share a beta version. To prepare yourself, you need R, Rstudio, and the packages shiny, colorRamps, lubridate and leaflet

some impressions

Hi @FeiLiu,

I changed some parameters mainly in step 6 in stamps and the PS data looks now like that (my “best” result) even that the trend becomes more clear and the selected points look good, in the end, the last point is unwrapped wrong or disturbed.

Until now I do not use any other software to subtract atmosphere, I know about TRAIN but I have not implemented it yet…maybe this would help. As you know I am trying hard on SBAS approach for this case study and hope I get a more satisfying result for the last and most interesting points in the TS plot. I think this could work because in PS the btemp is the largest for the last point, which is clear and a problem of the approach and in SBAS this problem may be solved…so lets see what happens next week

Hi @hchong005 ,

in this thread you will find the beta version of the application I use to visualize my data

Hi

My name is Daniel

So I’ve been working based on workflow described by amazing people such as @katherine , @thho, @FeiLi

And I get this result using ps_plot(‘w’) command

so using 7 SLC data to create 6 interferograms

my question is since I’m very new at this and still struggling and learning on inSAR maybe you can take a look at my result and give some interpretation.

And I have question for @thho in order to export the point to CSV using script in your post How to prepare Sentinel-1 images stack for PSI/SBAS in SNAP

is it run on Matlab command window or write an m file and then run it on matlab?

or is there a default command in StamPS to export PS to CSV?

good to see that the workflow works

I can not really help you to interpret the images, since I am new to this topic too, I am struggling myself what is deformation, what atmosphere etc. Maybe someone with more knowledge can help you.

Regarding your question, I would suggest to use the lines step by step copy pasting them to the matlab console, because when you do ps_plot(‘v-do’, ‘ts’) you have to do some manual input in defining the plot area and radius. I really recommend to do this just when a reference point is set! Your results will always be heavily affected when relying on the mean of all points of the image…creating and using the csv output should be discussed here I suggest. No there is no default command to export from StaMPS as far as I know, I used a combination of objects created during the StaMPS process and by some ps_plot() commands to build table which can then be exported as csv, it is somehow similar with the google earth kml export which you can use in StaMPS…

hi all

so what is the different between v-ds and v-d because when I plot using ps_plot. I got result as follows:

ps_plot(‘v-ds’)

ps_plot(‘v-d’)

I’ve read the manual and there’s no description about v-d and v-ds variable

hello guys

i was able to run the psinsar test thx to members from the group and special thx to @thho. but trying to run a new stack i am getting the following error. anyone have any idea why this is coming and how to correct it

Dear FeiLiu, I looked you question about stack averaging but I didn’t found an answer. Do you know why it is necessary put dinsar averaging stack and not dinsar stack?

All the best,

Antonella

Dear Sir

I am also working on Landslide deformation detection. Wish to do SBAS or PSI analysis. I do not worked with Stamps.

Now SNAP is also being made capable to handle PSI/Stamps.

What would you suggest ?

Will you please help me in Stamps processing?

Regards

ASVirk

Hi @virk_rana,

SNAP offers you to perform a preprocessing and export to StaMPS for PSI, SBAS is not possible yet. If you want to join a group of people testing the SBAS workflow with SNAP-StaMPS please have a look here.

I I recommend to read this thread, the papers published by Andy Hooper in which he presents StaMPS and the literature cited in this papers. There are some slides which helped me to become more familiar with the idea, here some more for SBAS. There are other sources that helped me to understand SBAS/PS, this thesis was very helpful for me, but here I will stop recommend literature.

when you feel comfortable with the theory, I suggest following one of the summaries for SNAP-StaMPS workflow given in this Forum, like this one.

The processing in StaMPS somehow depends on your dataset, how many images, how large is your studysite, which topographic and atmospheric conditions, how large is the area where you think a landslide occurs etc. I can not help you answering this questions, it is to specific (for my level of knowledge) and I have to dig deep into each dataset which I process, to set parameters right in StaMPS. But do not feel discouraged, StaMPS is a great tool, you have to invest some time, but you will become more and more familiar with your data and the parameters. So reading some recommended theory and starting preprocessing in SNAP and processing in StaMPS step by step is my advice. Many questions are already answered here or in MAINSAR google group. Gathering the information and hints is another good way of learning the process

Good luck and I am curious about some of your results!

Hello! I coming back to my quastion. In Some cases (most of all) mix index in Index_NaN is more than lendth of bperp_mat.

For example max index in Index_NaN is 24454944 and length of bperp_mat is 24454943.

Maybe you know the some reason of it and maybe you have some decision?

PS: some addition information. If the pixel size of product is odd - processing works. If it is even - send error in ps_load_initial_gamma in string 191 - index Index NaN is bigger then size of bperp_mat. Or in 205 string - diferent size cann’t concatinate.

No, my results are related to co-seismic and post-seismic results means after the big earthquake causes big load deformation.

Hi, @thho. Inoted that you have mentioned that Intrieri et al. 2017 used SBAS(SqueeSAR) technology. Actualy, SqueeSAR is one of the PSInSAR. see the picture below: