Dear STEP community,

I’ve found some intrigging results in the STAMPS Standard Deviation export, for what I can’t find any explation yet.

I’m using an algorithm adapted from @thho StaMPS-Visualizer manual.

It follows:

%%% the ‘v-doa’ parameter is an example you can change it to your needs

ps_plot(‘v-dao’,‘a_linear’,‘ts’);

%%% a new window will open

%%% in the new window select a radius and location of the radius center to select the PS to export

%%% the ‘v-dao’ parameter is an example you can change it to your needs

%%% but be sure that you use the same paramters as above in the ps_plot()!ps_plot(‘v-dao’,‘a_linear’,-1);

%%% loading phase displacement velocity from v-dao plot

load ps_plot_v-dao.mat;

%%% concatenation of phase velocity with time series usin lat lon

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); %% mean displacement velocity

disp_ts = ph_mm(ind,:); %% displacement time-series%%% leaving one column with zeros to fill with STD results (next)

%%% 1#id 2-3#lon-lat 4#vel 5#stdev 6-n# time-seriesexport_res_vel = [linspace(1,length(ind),length(ind))’ lonlat(ind,1) lonlat(ind,2) disp zeros(length(ind),1) disp_ts];

%%% To see time table (uncomment next line)

%%% datetime(day,‘ConvertFrom’,‘datenum’)%%% now lets load Time-series STANDAR DEVIATION for previous selected region

ps_plot(‘vs-dao’,‘a_linear’,-1);

load ps_plot_vs-dao.mat;%%% variable with STD values

export_res_vel(:,5)=ph_disp(ind);

export_res = [[zeros(1,5) m2xdate(day)'];export_res_vel];

export_res = table(export_res);

writetable(export_res,‘ts_peneda_Asc_20210111_20230218.csv’)

Every time I repeat this procedure for the following PS Time Series:

Using a 50 000 m radius and selecting the same reference lat lon (given by the small balck asterisk (*) in next image):

I always get different STD values, although the TS Velocities are always the same:

FIRST RUN results:

Export4 (is TS mean Velocity) and Export 5 (is Standard Deviation)

SECOND RUN results:

THIRD RUN results:

Does anyone have any guess about what could be trouble here? I really don’t have any idea whats is the issue