StaMPS-Visualizer, SNAP-StaMPS Workflow

Hello @thho, if you want, ı can give you an idea which can be helpful for your program :slight_smile:

Adding point optiong for example ? for seing which levelling or GNSS point is close to which PS point.
This can be very helpful for users.

Also if we can add the vertical velocities or coordinates of our levelling or GNSS points, and then if we can compare our terrestirial results with image results like this, it can be so much better.

Dear @firat,

thank you for your suggestion, it is similar to the one you mentioned some weeks earlier and in fact I have started working on it but have to stop, due to the fact, that I have come to a point where I have to expand the whole project and give it more time. The GPS comparison will come sooner or later, atm I am no longer brave enough to make statements about a date, because I have to work on it in my spare time and next to other projects. Anyway, smaller features will be added gradually so feel free and thereby I would like to encourage everybody, to mention features you are missing, I will put them on the list :slight_smile:

2 Likes

hi,ABraun, I have export the kml inmy patch folder,but the kml is too large to load in google earth,what should i do?

actually, the script exports a csv file. Are you sure you entered the correct file extension?
A kml is only needed when you want to subset your data as described here: StaMPS-Visualizer, SNAP-StaMPS Workflow

i means when i finish the step “output to google earth”,the kml is too large

this is not related to the StaMPS-Visualizer then, but you find answers here:

oh , i want to use Stamps-visualizer, but dont know where to start with ,my plotform is windows can i finish the visualier under windows,can you give me some advice

Although I still believe that the whole StaMPS workflow only runs without errors under Linux, it should be possible to use the Visualizer in Windows to display the results. But note that first you will have to produce correct displacements with StaMPS (until step 7 at least)

You need to export the results from Matlab first, the visualizer itself can then be used to display the results. Please start with the list in the first post and install R and RStudio to run the scripts provided there. Once the Visualizer runs, you can find a detailled description on how to export the data as CSV from Matlab to be analyzed in this software.

1 Like

i got it ,but when i run the export_csv in matlab,:
ps_plot(‘v-do’, ‘ts’);
load parms.mat;
ps_plot(‘v-do’, -1);
load ps_plot_v-do.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’)

after the ps_plot(‘v-do’,‘ts’);
do i need to select point in the figure
image
and i select a point in the window , and go ahead
but when i got to’ export_res = [lon2 lat2 disp ts];’
the error occurs: beacaue the matrix in the bracket is not consist with each other;

for examplre:
the lon2:1041
the lat2:104
1
the disp:1041
but the ts:1
10

so when i go next:export_res = [lon2 lat2 disp ts];
an error occurs beacuae the matrix is not consist with each other
how to deal with it

1 Like

hard to guess. But it is important that you run

ps_plot(‘v-do’, ‘ts’)

first, then select a point to create the plots.
Then execute the rest of the script (beginning with load parms.mat)

Maybe @thho has an idea.

Good job creating the map by the way :+1:

1 Like

dear @zhuhaxixiong,

first of all, your steps look good to me, creating the plot, selecting a point with a search radius for the PS of interest and then continue with the load parms.mat is the right way to use the script I wrote.

Concerning your problem, my first guess is, that you have a very very high density of PS/DS. However, this sometimes leads to some unmatching issues reported earlier in this topic. To avoid this, select a very large search radius like 10000m to select all points and subset this later on, all these steps and the subsetting are documented in the app, you could look them up in the manual tab using the example version here if you have not installed R, Rstudio and run the app, which you should do anyway in order to use it.

What is a bit suspicious is, that the difference between the length of all variables to the ts variable is very large…If this does not fix the problem, could you provide a screenshot of your matlab window, showing the environment with our variables produced by StaMPS and my export script on the right hand side? This would be helpful, atm I do not have access to a matlab license and have not worked with the new StamPS 4.x maybe some details have changed in this part, but I doubt that. Anyway, a screenshot with the variables created until the export step of the csv table would be helpful.

1 Like

when i reduce the radius to smaller,the number of ps point get to four,but the same error ocurrs again,
the matrix export_res is consist of four matrix:lon2 lat2 disp and ts,and the first three matrix is a n1 matrix(n is the number of ps point ),but the ts is the time serise,in my enviroment ,there are 10 different time slc to process,so the division of ts is 101,so i think this may be the error

i have read the post in this title.two methods to solv this problem:
one is increase the radius to include all points and this will export the csv sucessfully(i dont not know why);
another is like @firat said ,use old version “staMPS_v3.3b1”
i hope some one can solve this bug

Increasing the search radius is and was the solution since the problem occurred in the first place. Using stamps 3.x is not the solution since the problem occurred first in this version.

The bug, or whatever it can be called, is not yet clear. Sometimes when I have processed with stamps, I observed doubled coordinates, which should not occur as far as I understand the algorithms… Anyway, I have managed to get a MATLAB license and now I will search for an old Insar_master date folder and check some things, I will report my insights if I do have some :slight_smile:

I have solved the problem by use stamps3.1 and repeat steps from 1 to 7,finally i export the csv file in the folder, and when i choose ts comparison two issus occurs in my picture:


(1)the line is not plot through the first time to the final time,but very random
(2)the time axix seems also random
i want wo know why?

1 Like

Could you provide the first 5 rows and maybe 6 columns (or more) of your exported csv table? Maybe a screenshot or an upload of this subset. My first guess is, that the csv was not created properly I could look it up and tell you what happened exactly.

sorry , i foget to add some lines in mannual to my matlab,it has been solved

1 Like

how to subset
Error in as.data.frame.default(x[[i]], optional = TRUE) :
cannot coerce class ‘“function”’ to a data.frame

1 Like

There seems to be a mistake in line 23 (I know you copy pasted it from the manual of the application from github or the online example, I have fixed it there too)

use:

write.table(sub.csv, file = "/home/user/stamps_visualizer/stusi/newstudysite/20170310.csv",
            row.names = F, col.names = T, sep = ",")

Explanation: sub is not an object created by the script containing a data.frame but a function called sub (which is an object too). the function write.table() awaits a data.frame like object in x as input, this explains the error message. It maybe that I first worked with an object called sub, when I wrote the script but then changed sub to sub.csv, to avoid logical conflicts with the function sub…nice style be me but bad implemented by forgetting to rename all sub objects. Thanks for the report :slight_smile:

1 Like

Nevermind…Long story short:

In the new (4.x) Version of StamPS there were some changes in the ts plot which highly effects the variables used for the csv export. I changed the Matlab script in the Manual tab, now csv export is available for the new StamPS version! Look in the github repo for changes.

However, I was just able to test this with one example, it would be great if others can test the export code and report errors or problems, @ABraun, @firat @zhuhaxixiong maybe you could help me here with your processed results?

Details:

In StamPS 3.x a ts plot showed all ts (time series) for each MP (measurement point) within the search area, the data could be found in the object ts, in the newer version 4.x ts is a no longer a matrix with all results but a matrix with one row representing some mean of all selected MP, which becomes clear when looking at the ts plot, which is just a single ts averaged for all MP within the search radius. Therefore the ts plot becomes more less messy if many MP are selected…on the other hand details in single ts are lost…a good argument to use StamPS-Visualizer!
To overcome this problem the new object ph_mm containing all ts of all measurement points is subset, resulting in a object which is similar to the old ts object. Now the matrix can be made of all these objects as usual.

Thanks to @zhuhaxixiong and @firat for pointing in the direction of this change after the StaMPS 4. release!

6 Likes