SNAP - StaMPS Workflow Documentation

For matlab running on linux,

could be that matlab is out of memory, permission/ownership problems, or a foreign filesystem.

matlab is running under ubuntu, and the errors are :slight_smile:
STAMPS: ########################################
STAMPS: ####### StaMPS/MTI Version 4.0b6 #######
STAMPS: ####### Beta version, Jun 2018 #######
STAMPS: ########################################

STAMPS: Will process current directory only

STAMPS: ########################################
STAMPS: ################ Step 6 ################
STAMPS: ########################################
STAMPS: Directory is sf_INSAR_20200107_data

PS_UNWRAP: Starting
Phase-unwrapping…
GETPARM: small_baseline_flag=‘n’
GETPARM: unwrap_patch_phase=‘n’
GETPARM: scla_deramp=‘n’
GETPARM: subtr_tropo=‘n’
GETPARM: tropo_method=‘a_l’
Error using load
Unable to read file ‘psver’. No such file or directory.

Error in ps_unwrap (line 31)
load psver

Error in stamps (line 503)
ps_unwrap
if you can see the errors it will be nice

hallo,
when i did again stamp(1,5) i got the same errors but I just saw that when i ran matlab i got the following message, is something that i must to add to the installation '" Failed to load module “canberra-gtk-module”
thks

this is not a StaMPS related error, but to go sure you could install these graphic packages by
sudo apt install libcanberra-gtk-module libcanberra-gtk3-module

I recommend to run each step separately: step(1,1) step(e,2) and closely check the messages. Sometimes indications for potential errors are given there (too few PS candidates, empty interferograms…)

Hallo,
yesterday I installed the gtk-module and I ran stamps(1,4), he stops for patche6 the last one and i got the following errors
TAMPS: ########################################
STAMPS: ################ Step 4 ################
STAMPS: ########################################
STAMPS: Directory is PATCH_6

PS_WEED: Starting
PS_WEED: Weeding selected pixels…
GETPARM: weed_time_win=730
GETPARM: weed_standard_dev=1
GETPARM: weed_max_noise=Inf
GETPARM: weed_zero_elevation=‘n’
GETPARM: weed_neighbours=‘n’
GETPARM: small_baseline_flag=‘n’
PS_WEED: 134013 low D_A PS, 0 high D_A PS
PS_WEED: 1 PS with duplicate lon/lat dropped

PS_WEED: Estimating noise for all arcs…
PS_WEED: Estimating max noise for all pixels…
PS_WEED: 133098 PS kept after dropping noisy pixels
psver currently: 1
psver now set to: 2
PS_WEED: Finished
Error using movefile
Text file busy

Error in stamps (line 462)
movefile(‘patch.list’,‘patch.list_old’);

let’s don’t mix things up. A solution to this was suggested here yesterday:


Please give feedback on our suggestions, we cannot look into your system.

1 Like

Dear all,
At first I would like to thank you, i just finished the running of stamps after your help and the help of Sarto28, I was late but I had problem with the expiration of Matlab trial version and i found problem in the installation of the other one, anyway this is de map of ps velocity . I don’t know to comment it, must I have a reference point to calculate exactly the velocity of landsubsidence.

2 Likes

I cant select the time series in order to exporting to CSV file I used the k largest radius of 1000000m

please help me. thank you in advanced

How can we find the deformation value of Individual PS point after completing stamps(8,8)…?

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

this opens the displacement map and allows you to select an area with a defined radius (100 m are predefined but you can reduce it to 10 m).
Once you click in the image, a time-series graph is displayed for the selected area.

1 Like

Appreciated Thank you so much @ABraun

I have use ps_plot(‘v-do’, ‘ts’) for individual pixel i have use 20 meter radius and 5 pixels are in the circle… i want to check the value of these 5 pixels one by one , can we check values one by one ??

hi @atif, I built a tool for exactly this task:

2 Likes

Thankyou so much Sir @thho … i will follow this tool now…

I have used 91 interferograms of 3 years to find trend of subsidence in one area, can i plot kml of 1 year ?? i mean i want to plot kml year to year, for example 1st i want to plot kml for 2017 then 2018 and then 2019…is it possible to do so??

At this stage, it does not matter where we click? On what basis is this point selected?

In mt_prep_snap command, you just give less amplitude dispersion value. Then number of PS should decrease.

1 Like

You will find a clear explanation in the below link.

1 Like

Dear members,
Recently I’m running data set of both ascending and descending of my study area. Descending was running fine and providing good results. But while processing the Ascending stamps giving below error.
Please check it once.

STAMPS: ########################################
STAMPS: ################ Step 2 ################
STAMPS: ########################################
STAMPS: Directory is PATCH_1

PS_EST_GAMMA_QUICK: Starting
PS_EST_GAMMA_QUICK: Estimating gamma for candidate pixels
GETPARM: filter_grid_size=50
GETPARM: filter_weighting=‘P-square’
GETPARM: clap_win=32
GETPARM: clap_low_pass_wavelength=800
GETPARM: clap_alpha=1
GETPARM: clap_beta=0.3
GETPARM: max_topo_err=20
GETPARM: lambda=0.0554658
GETPARM: gamma_change_convergence=0.005
GETPARM: gamma_max_iterations=3
GETPARM: small_baseline_flag=‘n’
Found look angle file
PS_EST_GAMMA_QUICK: n_trial_wraps=0.355834
PS_EST_GAMMA_QUICK: Initialising random distribution…
PS_EST_GAMMA_QUICK: 2563970 PS candidates to process
PS_EST_GAMMA_QUICK: iteration #1
PS_EST_GAMMA_QUICK: Calculating patch phases…
PS_EST_GAMMA_QUICK: Estimating topo error…
PS_EST_GAMMA_QUICK: 100000 PS processed
PS_EST_GAMMA_QUICK: 200000 PS processed
PS_EST_GAMMA_QUICK: 300000 PS processed
PS_EST_GAMMA_QUICK: 400000 PS processed
PS_EST_GAMMA_QUICK: 500000 PS processed
PS_EST_GAMMA_QUICK: 600000 PS processed
PS_EST_GAMMA_QUICK: 700000 PS processed
PS_EST_GAMMA_QUICK: 800000 PS processed
PS_EST_GAMMA_QUICK: 900000 PS processed
PS_EST_GAMMA_QUICK: 1000000 PS processed
PS_EST_GAMMA_QUICK: 1100000 PS processed
PS_EST_GAMMA_QUICK: 1200000 PS processed
PS_EST_GAMMA_QUICK: 1300000 PS processed
PS_EST_GAMMA_QUICK: 1400000 PS processed
PS_EST_GAMMA_QUICK: 1500000 PS processed
PS_EST_GAMMA_QUICK: 1600000 PS processed
PS_EST_GAMMA_QUICK: 1700000 PS processed
PS_EST_GAMMA_QUICK: 1800000 PS processed
PS_EST_GAMMA_QUICK: 1900000 PS processed
PS_EST_GAMMA_QUICK: 2000000 PS processed
PS_EST_GAMMA_QUICK: 2100000 PS processed
PS_EST_GAMMA_QUICK: 2200000 PS processed
PS_EST_GAMMA_QUICK: 2300000 PS processed
PS_EST_GAMMA_QUICK: 2400000 PS processed
PS_EST_GAMMA_QUICK: 2500000 PS processed
PS_EST_GAMMA_QUICK: gamma_change_change=0.000000
GETPARM: gamma_change_convergence=0.005
GETPARM: gamma_max_iterations=3
PS_EST_GAMMA_QUICK: Finished

STAMPS: ########################################
STAMPS: ################ Step 3 ################
STAMPS: ########################################
STAMPS: Directory is PATCH_1

PS_SELECT: Starting
PS_SELECT: Selecting stable-phase pixels…
GETPARM: slc_osf=1
GETPARM: clap_alpha=1
GETPARM: clap_beta=0.3
GETPARM: clap_win=32
GETPARM: select_method=‘DENSITY’
GETPARM: density_rand=20
GETPARM: gamma_stdev_reject=0
GETPARM: small_baseline_flag=‘n’
GETPARM: drop_ifg_index=
Warning: Not enough random phase pixels to set gamma threshold - using default threshold of 0.3

In ps_select (line 184)
In stamps (line 357)
PS_SELECT: Initial gamma threshold: 0.300 at D_A=0.08 to 0.300 at D_A=0.50
PS_SELECT: 0 PS selected initially
PS_SELECT: Reestimation gamma threshold: 0.300 at D_A=0.08 to 0.300 at D_A=0.50
PS_SELECT: 0 ps selected after re-estimation of coherence
No PS points left. Updating the stamps log for this*
PS_SELECT: Finished

Good afternoon to everyone,
I’m asking for help to understand how the displacement time series export works. The displacement in mm that can be exported after selecting the center and radius on the ps_plot(‘v-dao’, ‘ts’) figure, I know they are referred to the master date, meaning that its displacement is zero. I usually have a master that is in the middle of the time series. But my question is: do I have to change the sign of all the displacements AFTER the master date or they are already correct?
This is because as far as I have understood, when computing the phase difference (and as a consequence also the displacement) the equation is: master - slave, therefore once we are after the master date, the sign of the result changes.
The code that I use to export the time series is the one suggested in the StaMPS-Visualizing discussion. My question is arising because of data that I’m actually processing, that I’m trying to validate with in-situ GNSS receivers, and if I apply this correction to the displacements after the master date the two time series are very similar, which is not true in the other case.
From the StaMPS manual I can say that positive values mean uplift while negative subsidence (“For plotting of velocities, the units are mm/year with positive values being towards the satellite”).
Thank you in advance!