Export ps_plot ('vs-do')

I would like to export a .csv file relative to the standard deviation of velocities (and maybe also of the displacement).

It is easy to export just the velocity STD, but I would like to add the lat and long.

Do you know how to do it?

My script to export just the STD velocities is:

load parms.mat;
ps_plot(‘vs-do’, -1,0,0,[], 0, 0, 0,[], [10.93 10.94], [45.06 45.08]);
load ps_plot_vs-do.mat;
export_res = [ph_disp];
export_res = table(export_res);


I moved your topic to the #stamps category

1 Like

Does anyone know how to find the coordinates of the PS velocity standard deviation?


maybe you can adapt on the export code of the StaMPS-Visualizer, SNAP-StaMPS Workflow

Maybe I have found a way to adapt that script to ps_plot (vs’).

hi @Sarto28, would you mind telling how did you solve the problem of the exporting of vs-plot?

I have created a new script similar to that used to export v-do. But first, you must run the ps_plot (‘v’), select the area you are interested in and save three files: “lonlat”, “lon2” and “lat2”.
Then you can run the new script as follow:

load parms.mat;
ps_plot(‘vs-do’, -1,0,0,[], 0, 0, 0,[], [10.92 10.95], [45.06 45.09]);
load ps_plot_vs-do.mat;
load (‘lonlat’);
load (‘lon2’);
load (‘lat2’)
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))));
lonlat_str_sec = strcat(cellstr(num2str(lonlat(:,1))),’,’,cellstr(num2str(lonlat(:,2))));
ind = ismember(lonlat_str, lonlat2_str);
vel = ph_disp(ind);
lonlat_str_aux = lonlat_str_sec (ind);
for i=1:row
Auxvar = strsplit(lonlat_str_aux{i,1},’,’);
lon2exp(i,1)= str2num(char(Auxvar(1,1)));
lat2exp(i,1)= str2num(char(Auxvar(1,2)));

export_res = table([lon2exp lat2exp vel]);

It worked for me.


1 Like

Thank you Luca for your time and help. I tried to do it with my worksheet but it didn’t go well. I had this problem:

Probably becouse i made an error with the first plot. To find and save the coordintes, i used the same plot as for mean velocity deformation: ps_plot(‘v-do’, ‘ts’); and then i run you script. I tried to use ps_plot (‘v’) but it doesn’t work at all.