To check for issues, I need a look at the raw csv as presented in a text editor. I have to check for the col delimiter if it is a , or ; . in Libreoffice this is not possible to see.
Hi dear @thho
As you mentioned I checked the manual of the visualizer (Create a coregistered stack product (InSAR stack) in SNAP. Be sure that the prime image is the first in the stack!) but when I try to create a coregistered stack I got this error
1- Topsar Split
2- Apply the orbit file
3- Back Geocoding
4- ESD
The pic is just an example, and for the baseline plot, I got 34 SLC images. Do I need to perform more steps?
hi @masoodhp
in the first script of your earlier post, you have to delete the . in line 13 and 14 when referring to the cols
in the second script, change ; to , in line 16 and 27
That should do the trick
Hi @masoodhp ,
hmm havenāt done this for a very long time, @ABraun or @mdelgado can you have a quick look at it?
Thanks for your response.
I changed the lines 16 and 27 now I got this error.
and about the first one (subset ts plot export).
@masoodhp this is just a warning, the shapefile was most probably written to disk. ESRI shapefile has its drawbacks, one of it is, that you have a quite limited number of characters which you can use as column headers in the attribute tableā¦the warning is about, that they were shortened automatically during the write process
The problem has been solved following your instructions and installing the latest version of STAMPS-Visualizer. Also, I re-run the whole process, dividing the dataset to further Patches and all works fine. Thank you very much!
1 Like
Dear @thho ,
I have 5 GPS stations in my Study area. I want to plot these station locations in Stamps_visualizer plot using latitude and longitude ?
Is it possible to plot those stations in the PS_InSAR plot generated by using StaMPS_Visualizer?
hi @suribabu,
the current version of the visualizer supports the import of custom gejson geometries, thus adding a point layer should work out of the box.
In order to make these point also clickable, in order to plot their time series is more complicated and you would have to add some options in the code I guess.
1 Like
Hey, @Gijs
1- Can you explain how did you manage to combine Asc and Des data, please? I saw your codes on GitHub but I donāt know how to use them.
2- Actually here (StaMPS-Visualizer, SNAP-StaMPS Workflow - #199 by Gijs) you said that āAfter the reference point you can use the formula: LOS displacement / (cos(radians(incidence_angle)) to convert the LOS displacement to vertical displacementā. So If I want to calculate vertical displacements, no combination is needed, right?
If yes, should I use the velocity column from the CSV file or the displacements of the points?
@thho @ABraun Can you help me with these?
Greetings
Masood
Hello, dear @ABraun . How can I explain this graph? In this case, interpret it for me a little, please.
Sincerely.
This is a typical sime series of deformation extracted from multi-temporal DInSAR techniques (PSI or SBAS). Each blue dot represents the displacement at one time relative to a reference date. The red line is optional and indiaces the overall trend which is negative (decrease of the surface), but the RĀ² of 0.82 indicates that the observed trend is not a perfect line, that means the displacement has phases of acceleration and deceleration.
Thank you, dear @ABraun , but what factors can affect the distribution and trend of blue dots (apart from displacement), can rainfall and groundwater feeding be among these factors?
Two kinds of factors can have an impact:
- actual surface changes (groundwater extraction, tectonics, soil compaction)
- phase errors because of atmospheric disturbance or unwrapping errors.
Itās not always easy to decide which of them are observed in the time-series, but long time-series tend to be more robust towards errors.
After assign the longitude, latitude and radius successfully, while importing the CSV file with respect to that particular lat long and radius is giving error, please check it once.
K>> setparm(āref_cenā,[72.012 21.630])
SETPARM: ref_centre_lonlat = 72.012 21.63
K>> setparm(āref_radā,1000)
SETPARM: ref_radius = 1000
K>> ps_plot(āv-doā)
Deramping computed on the fly.
**** z = ax + by+ c
861 ref PS selected
Color Range: -6.20838 to 8.16557 mm/yr
ans =
Figure (1: v-do) with properties:
Number: 1
Name: 'v-do'
Color: [0.9400 0.9400 0.9400]
Position: [680 558 560 420]
Units: 'pixels'
Show all properties
K>> a new window will open in the new window select a radius and location of the radius center to select the PS to export
load parms.mat;
the 'v-doa' 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-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);
disp_ts = ph_mm(ind,:);
export_res = [lonlat(ind,1) lonlat(ind,2) disp 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);
% you can specify the location and name of the .csv export by renaming the second parameter
writetable(export_res,āreference_deva.csvā)
Deramping computed on the fly.
**** z = ax + by+ c
861 ref PS selected
savename =
'ps_plot_v-do'
Color Range: to 109 m/yr
Error using keyboard
Unrecognized function or variable ālon2ā.
Did you mean:
K>> lon2_str = cellstr(num2str(log2));
Error using keyboard
Not enough input arguments.
Dear @thho,
Itās an excellent explanation please update it. If you have already updated the manual please forward the link.
Thank you.
this error appeared after attempting to install the leaflet as required after running the ui.R script
And this error appeared after clicking the run app button
