dbekaert / stamps

Stanford Method for Persistent Scatterers

StaMPS - Stanford Method for Persistent Scatterers

Description

StaMPS is a software package that allows to extract ground displacements from time series of synthetic aperture radar (SAR) acquisitions. The package incorporates persistent scatterer and small baseline methods plus an option to combine both approaches. It is compatible with the TRAIN software and therefore allows to incorporate various tropospheric correction methods in the processing workflow.

Installation

See the StaMPS manual for installation details and dependencies.

Supported InSAR pre-processors:

• ISCE
• GAMMA
• SNAP
• DORIS + ROI_PAC

Development

The original version was developed at Stanford University but subsequent development has taken place mainly at the University of Iceland, Delft University of Technology and the University of Leeds. In addition, there have been community contributions.

Contributors

• Andy Hooper (Lead)
• David Bekaert
• Karsten Spaans
• Ekbal Hussain
• Mahmut Arikan
• Anneleen Oyen
• Miguel Caro Cuenca
• Jose Manuel Delgado Blasco
• other commmunity members

https://github.com/dbekaert/stamps/graphs/contributors

Community contributions are welcomed.

Contact details

Andy Hooper
COMET
School of Earth and Environment
University of Leeds
Leeds, LS2 9JT

stamps's Issues

step6

when i run step6, the problem as following make me confusion, how should i solve this if i want to go on next step? please give me some suggestions about it , thank you!

ps_plot('v-d',2) problem for ISCE2STAMPS

Hi,
I think there is bug in order to plot ps_plot('v-d',2)
Best,
MMA

Warning: Usage of DATENUM with empty date character vectors is not supported.
Results may change in future versions.

In datenum (line 86)
In plot_all_ifgs (line 146)
Error using axes
Value must be numeric
Error in plot_all_ifgs (line 152)
axes('position',[imX(i),imY(i),h_x,h_y])

Error in processing using SLCs from ISCE2

I encountered an error after I ran stamps.m using the coregistered SLCs using ISCE.

STAMPS: ########################################
STAMPS: ####### StaMPS/MTI Version 4.0b6 #######
STAMPS: #######  Beta version, Jun 2018  #######
STAMPS: ########################################

STAMPS: Will process current directory only

STAMPS: ########################################
STAMPS: ################ Step 1 ################
STAMPS: ########################################
STAMPS: Directory is PATCH_1

Found inc angle file, will load the data 
Loading data into matlab...
SETPARM: heading = -11.1257
SETPARM: lambda = 0.055466
Unable to perform assignment because the left and
right sides have a different number of elements.

Error in ps_load_initial_isce (line 126)
        caldate(i)=bb;

Error in stamps (line 281)
                    ps_load_initial_isce(data_inc)

export points and their values to a list for further data analysis

Hi, I'm studying on a Dam displacement in a 6 years period. I need to fusion the displacement rate with daily Dam water level so I could find the possible coherence between these two variables. I have exported the results to StaMPS Visualizer as it has been shown below, but I need the raw displacement data for each point, I must mention It's Important That I must know each displacement data belongs to to which point on the Dam.

make_single_reference_stack_isce error

Stack path is: /mnt/d/InSAR/bukavo/merged/SLC
SLC has a suffix as: .full
Geometry files have a suffix as: .full
SM: Reference Date = 20141124

ERROR 4: 20141124.slc.full: No such file or directory
This is the Open Source version of ISCE.
Some of the workflows depend on a separate licensed package.
To obtain the licensed package, please make a request for ISCE
through the website: https://download.jpl.nasa.gov/ops/request/index.cfm.
Alternatively, if you are a member, or can become a member of WinSAR
you may be able to obtain access to a version of the licensed sofware at
https://winsar.unavco.org/software/isce
Traceback (most recent call last):
File "/home/boudhane/miniconda3/envs/InSARenv/lib/python3.11/site-packages/isce/applications/gdal2isce_xml.py", line 137, in
raise Exception('Input file is not found ....')
Exception: Input file is not found ....
xml file for the SLC does not exist, you need to have these

Error using movefile Cannot write from source

Error using movefile
Cannot write from source: H:\ss\patch.list. Use the 'f' option to override.

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

When i use

movefile('patch.list','patch.list_old','f')

it shows the error

Error using movefile
The system cannot find the file specified.

Expand gamma description of manual

Expand gamma description with more information on:

• expected precision of files double, float etc
• inform about precision flag in the stamps_gamma scripts
• big versus Small Endian
• relative files sizes between ifg, lonlat etc.

SMALL BASELINE

Dear
I’m trying to run with STAMPS the small baseline but I have the problem in matlab

    stamps

STAMPS: ########################################
STAMPS: ####### StaMPS/MTI Version 4.0b6 #######
STAMPS: ####### Beta version, Jun 2018 #######
STAMPS: ########################################

Reference to non-existent field ‘small_baseline_flag’.

Error in ps_parms_default (line 111)
if strcmpi(parms.small_baseline_flag,‘y’)

Error in getparm (line 14)
ps_parms_default

Error in stamps (line 69)
quick_est_gamma_flag=getparm(‘quick_est_gamma_flag’);

Someone could help me?? I followed the STAMP manual

Thanks

Pedro

About matlab .m file

When running stamps, matlab prompts that there is no “strsplit” function.

Empty pscands files --- Error in llh2local (line 39) z=llh(2,:)~=0; Error in ps_load_initial_isce (line 161) xy=llh2local(lonlat',ll0)*1000; Error in stamps (line 281) ps_load_initial_isce(data_inc)

Hello,

I am using ISCE ( stackStripMap ) to create StaMPS - PS Processing with TerraSAR-X datasets. With the following steps:

1. !stackStripMap.py -s /home/toni/WCF_TerraSARX_processing/scenes/HS_044/SLC/ -d /home/toni/WCF_TerraSARX_processing/scenes/HS_044/DEM/demLat*.dem.wgs84 -t 120 -b 400 -a 1 -r 2 -m 20160627 -u snaphu -W slc --nofocus –zero

2. Ran all files run_files 1-7 using: !run_files/run_0*

3. In terminal StaMPS processing dir:

     $ make_single_reference_stack_isce  input_file
   
                             Input_file:
                             source_data slc_stack
                             slc_stack_path /home/toni/WCF_TerraSARX_processing/scenes/HS_044/merged/SLC/
                             slc_stack_reference 20170627
                             slc_stack_geom_path /home/toni/WCF_TerraSARX_processing/scenes/HS_044/merged/geom_reference/
                             slc_stack_baseline_path /home/toni/WCF_TerraSARX_processing/scenes/HS_044/merged/baselines/
                             maskfile /home/toni/WCF_TerraSARX_processing/scenes/HS_044/merged/geom_reference/shadowMask.rdr
                             azimuth_looks 1
                             range_looks 2
                             aspect_ratio 2
                             lambda 0.056
                             slc_suffix
                             geom_suffix
   

4. $ cd INSAR_20170627/

5. $ mt_prep_isce 0.4 2 1 50 50

6. $ matlab
   a. >> stamps
   b. Error message at STEP 1 of StaMPS
   STAMPS: ########################################
   STAMPS: ####### StaMPS/MTI Version 4.0b6 #######
   STAMPS: ####### Beta version, Jun 2018 #######
   STAMPS: ########################################

                                                 ans =
                                                 
                                                      0
                                                 
                                                 STAMPS: Will process patch subdirectories
                                                 
                                                 STAMPS: ########################################
                                                 STAMPS: ################ Step 1 ################
                                                 STAMPS: ########################################
                                                 STAMPS: Directory is PATCH_1
                                                 
                                                 Found inc angle file, will load the data 
                                                 Loading data into matlab...
                                                 SETPARM: heading = -10.8843
                                                 SETPARM: lambda = 0.056
                                                 Index in position 1 exceeds array bounds.
                                                 
                                                 Error in llh2local (line 39)
                                                    z=llh(2,:)~=0;
                                                 
                                                 Error in ps_load_initial_isce (line 161)
                                                 xy=llh2local(lonlat',ll0)*1000;
                                                 
                                                 Error in stamps (line 281)
                                                                     ps_load_initial_isce(data_inc)
   

My selpsc_patch command produces empty pscands.1.ij and pscands.1.da files which in return makes the pscands.1.hgt, pscands.1.ll, and pscands.1.ph files empty (literally 0-bytes). The input files (selpsc.in, patch.in, maskfile) for the selpsc_patch are not empty. Any ideas why?

The ph1.mat of one PATCH is corrupt while ph1.mat file of other PATCHes are normal.

Dear all,
After running stamps(1, 1) for all patches, there is one patch whose ph1.mat file is corrupt and can not be loaded by Matlab while ph1.mat file of other patches are normal.
Please help me if you know it. I am using StaMPS-master (Version 4.1b).
Thank you!

###############################
The error message are:

Error using load
Unable to read MAT-file
./INSAR_20160419/PATCH_13/ph1.mat. File might be corrupt.

Error in ps_plot (line 612)
            rc=load(phname);

PS_PLOT('vs') error

Hello everyone, after run the SBAS by stamps, the problem occurred when i want to plot the vs figure. Looking forward to your replies.

Issue in Step1(stamps(1,1))

wrong use ps_load_initial_gamma
./pscands.1.ll does not exit
stamps(line 263)
ps_load_initial_gamma

how to solve this problem?
before this, i have done the step :mt_prep_gamma ******** 0.6 1 1 50 50
and open matlab in the patch1

'Error opening file', while using mt_prep_isce for smallbaseline processing more than 1020 files

I have a small baseline network of 160 SLCs. This network contains 622 interferograms image pairs, which are processed and generated by ISCE software. I first used "make_small_baselines_isce" to create the small baseline interferograms listed, everything is fine. But when I use "mt_prep_isce" for preprocessing, the following error occurs like below.

`opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210605/reference.slc [file 1010]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210605/secondary.slc [file 1011]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210611/reference.slc [file 1012]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210611/secondary.slc [file 1013]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210617/reference.slc [file 1014]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210524_20210617/secondary.slc [file 1015]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210605/reference.slc [file 1016]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210605/secondary.slc [file 1017]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210611/reference.slc [file 1018]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210611/secondary.slc [file 1019]...
opening /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210617/reference.slc [file 1020]...
Error opening file /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/20210530_20210617/reference.slc
terminate called after throwing an instance of 'char const*'
Aborted (core dumped)

psclonlat /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/psclonlat.in pscands.1.ij pscands.1.ll
opening pscands.1.ij...
Error opening file pscands.1.ij

pscdem /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/pscdem.in pscands.1.ij pscands.1.hgt
opening pscands.1.ij...
pscdem: Error opening file pscands.1.ij

pscphase /mnt/f/ma_landslide/stack1/merged/INSAR_20200803/SMALL_BASELINES/pscphase.in pscands.1.ij pscands.1.ph
opening pscands.1.ij...
Error opening file pscands.1.ij`

I am quite sure that the memory and hard disk space are sufficient. And the file that need to be opened is not broken. And the "selsbc_patch" program is compiled correctly. Because when I shrink the size of small baseline network to reduce the number of interferograms pairs, the entire "mt_prep_isce" script can run smoothly when the number of opened files is less than 1020.
I was able to locate the error at the line 198 of the src/selsbc_patch.c file, but I am very unfamiliar with the C language, and I don’t know how to fix it. Does anyone have good suggestions? Thanks so much!

Issue in Stamps Step(1,1)

While running the stamps step(1,1) I've encountered this error.
Found inc angle file, will load the data
Loading data into matlab...
Index exceeds the number of array elements. Index must not exceed 0.

Error in ps_load_initial_isce (line 82)
bperp=bperp(day_ix);

Error in stamps (line 281)
ps_load_initial_isce(data_inc)****

Kindly help me with this.

isse in stamps (step 6)

Dear all, I find the error"Error using / ,Matrix dimensions must agree." when I run step 6 in StaMPS. I have no idea with it and look forward your answer.

step2 error about ph_grid is stemmed from the bug of psclonlat

I have found some comments about the ph_grid error in step 2. Actually, it is due to the nan value or error value ( lon and lat) in the pscands.1.ll, which is produced by the psclonlat code. This bug of nan value or error value will also cause the hgt error in the pscands.1.hgt, which is producted by the pscdem code. I have tested the question with sentinel-1 data, and found the data with line number < 34000 usually not cause the bug. But if the data with line number > 34000 or 35000 usually cause the bug in the PATCH_, whose line input number > 34000 or 35000 (refer to patch.in in the PATCH_). So, it seems there are some bugs on the psclonlat.c and pscdem.c code in the /src folder of StaMPS. I hope someone can help to check and resolve this question. I don't know much about programming.

K2q.m didn't include Sentinel-1 parameters

K2q.m is used for dem error estimation but it's sensor depended ans only supprot ERS/Envisat, ALOS1 and Radarsat-2 U18. And it will assume ERS/Envisat if sensor not specified.

else
%ERS/ENVISAT
slant_range=7.9048902811596;
near_range=831265;
h=794141;
lambda=0.0565647;
end

Can someone confirm if this will affect the accuracy of dem error estimation? I am willing to add some lines for Sentinel-1.

mt_prep_gamma line 143

Hi
I see that mt_prep_gamma has in line 143 an "end" and so It didn't ran completely,
I commented this and the program ran ok, but I am concerned that I should have removed that "end"

Processing time on mt_extract_cands

Hi, I am not sure if this is a bug or not. And I can't really ask the question in a more easy way. So sorry in advance.

I have prepared the files in the directory using snap2stamps and I execute the command

./mt_prep_snap 20200227 /mnt/d/INSAR_20200227 DA

I use the WSL (ubuntu under windows 10) and my files are on the d drive, while the scripts are under the linux machine. I do not see any errors so far, but the last bit of the scipr (mt_exract_cands) takes too long and I am not sure if it crashes or not (computer mentoring shows high cpu usage). How long does this process takes and how can I know if something is not going well

link_slcs execution error

I executed the following command according to the StaMPS manual:
link_slcs / home / sar / ENVISAT
But the following error message appears:info:
/home/sar/ENVISAT/ASA_IMS_1PNESA20030801_024605_000000182018_00347_07415_0000.N1
awk: Fatal error: Cannot open file "/home/sar/ENVISAT/ASA_IMS_1PNESA20030801_024605_000000182018_00347_07415_0000.N1/L" in read mode (not a directory)
mkdir: missing operand*

What should I do to fix it? Thank you.

HELP,I can't run stamps(1,1),did anyone met like that

i run stamps after finished isce, but i can't solve this problem.

or like this, in this directory

Sign Convention for Plotting Velocity and Phase

I am confused with the sign convention of velocity and phase retrieved from StaMPS processing using ascending and descending passes. What should be the sign conventions for the following situations-

Ascending Data

1. If reference image predates the secondary image- what is the meaning of positive and negative velocities and phase
2. If reference image postdates the secondary image-what is the meaning of positive and negative velocities and phase
3. If reference image is in the middle of all the secondary images - what is the meaning of positive and negative velocities and phase

Descending Data

1. If reference image predates the secondary image- what is the meaning of positive and negative velocities and phase
2. If reference image postdates the secondary image-what is the meaning of positive and negative velocities and phase
3. If reference image is in the middle of all the secondary images - what is the meaning of positive and negative velocities and phase

Please help me if you know it. I am using StaMPS-master (Version 4.1b).

Thanks in advance
Hari Shankar

Step 2 in StaMPS

Hi,
I'm working on StaMPS in step-2 I got the following error:
Error in ps_est_gamma_quick (line 221)
ph_grid(grid_ij(i,1),grid_ij(i,2),:)=ph_grid(grid_ij(i,1),grid_ij(i,2),:)+shiftdim(ph_weight(i,:),-1);

I need your help to solve it.

mt_prep_snap on Mac

I'm running the mt_prep_snap on Mac for the first time. It runs now for over 8 hours and there is no results or any error messages or even any change in the Terminal window?
I do not know if the process still running or how long it takes to be completed?
any idea?

plot_all_ifgs error

Dear all,

Please how can I solve this problem with plot_all_ifgs, I'm using snap2stamps

Warning: Usage of DATENUM with empty date character vectors is not supported.
Results may change in future versions.

In datenum (line 94)
In plot_all_ifgs (line 146)
In plot_all_ifgs (line 18)
In plot_all_ifgs (line 18)
Error using axes
Value must be numeric and finite
Error in plot_all_ifgs (line 152)
axes('position',[imX(i),imY(i),h_x,h_y])

Cheers

Unable to write files when using any StaMPS function

Hi everybody, I'm having an issue.
While running any function like setparm() or stamps(), this error happens:

Error using save Unable to write file no_ps_info.mat: Read-only file system.

I have also tried to change permissions in the main directory, but nothing changes.

mt_prep_snap bug

Dear All, When I execute mt_prep_snap 20170107 /home/... 0.4 5 5 50 50 The programme reports an error as shown below, how to solve the problem?

scihub_poeorb_download.sh > mv: cannot stat [...] File name too long

Please delete this issue, thanks

Incorrect? PS hgt (height): StaMPS PSI Processing

I plotted an histogram of the persistent scattering hgt file from StaMPS and realized that the majority of the PS heights were -28, -27, and -26. I did the same, but using the SRTM dem that was used in the ISCE2 stack processing (clipping it using each PS location to get the WGSS84 elevation). The SRTM dem seems like a more realistic distribution of height in my study area. Is there a way to fix this problem? I'm confused as to why my PS heights are mainly -28, -27, and -26. (see attached figures)

ISCE2STAMPS plot_all_ifgs

Hi David,
there is a problem for plot_all_ifgs for isce2stamps output
Warning: Usage of DATENUM with empty date character vectors is not supported.
Results may change in future versions.

In datenum (line 86)
In plot_all_ifgs (line 146)
Error using axes
Value must be numeric
Error in plot_all_ifgs (line 152)
axes('position',[imX(i),imY(i),h_x,h_y])
I know the get_quickview_isce can create what I want but I was wondering if it is possible in Matlab to see the interferograms?

mt_prep_snap issue

Hello dear all,when I run mt_prep_snap 20200722 /home/sun/export 0.4 5 5 50 50,it shows error opening file pscand.1.ij,it output some patch files without pscand.1.ij,and I have 2g
memory,Is it because the memory is too small？

How to show the azimuth and range direction

After using ps_plot('v'). is it possible to show azimuth and range direction based on information from metadata file with in plot using arrows to get idea of where the satellite is heading and how it is looking at the ground ?

NaN valyes brake stage 2

I am using STAMPS on win64. Prepare SENTINEL-1 data in SNAP. There is NaN values when export from SNAP 6 to STAMPS. For the STAMPS 3.3 i download ps_load_initial_gamma from https://forum.step.esa.int/t/how-to-prepare-sentinel-1-images-stack-for-psi-sbas-in-snap-5/4981/68. There is script with corection NaN values. On pictures is error report


Could you help me to repair this error?

Move file error

When I run stamps(1,1) on win64, get error in line 462 - somesing like "file is busy, cannot get access".
Repair problem by addin in line 129 string - fclose(fid)

GEOLOCATION CORRECTION - STAMPS

Dear,

Any ideas on how to update the lon/lat coordinates using the residual topographic information from the estimated PS during the StaMPS-steps?
I found a good approach here: https://www.mdpi.com/2072-4292/11/11/1258
I am even using that same article to validate the methodology of my work. The first image below is the result obtained in StaMPS without the lon/lat coordinate correction approach. The second is with the approach, how can I improve?

Image 1

Image 2

The abortion adopted in the article uses the equations below:

The parameters I used were:
theta - angle of view stored in la2.mat as la
phi - heading stored in params.mat
deltaH - DEM error that is provided by in rad/m, being stored in scla2.mat as K_ps_uw, to convert to meters I used the K2q.m function. However, due to the non-existence of parameters for Sentinel-1, I adopted the following values:

For the range_pixel parameter present in the K2q.m function, I used the values stored in ps2.mat as ij(:,3).

It is worth mentioning that for the InSAR processing, the DEM SRTM 1Sec HGT was not used, as it presented problems in the download, as shown in the image below. Thus, the Copernicus 30m Global DEM was used.

Thanks in advance if you can help me.

Error in PS_SB_MERGE Operation for MTI Processing

Dear All

I have processed first five steps of PSI i.e. stamps(1,5) under INSAR_Master_date and then first five steps of SBAS under /SMALL_BASELINES. And then as per the requirement of MTI processing, I have to merge/combine the pixels processed with both PSI and SBAS methods, and for this I have fired this command in matlab workspace: ps_sb_merge. While running this command, the following error message is received.
Note:- In PSI processing: Total IFGs 44, and in SBAS Processing: Total IFGs 153.

Run stamps (4, 4), the terminal reports an error

Run stamps (4, 4), the terminal reports an error

how to solve this problem？
thanks

uw_sb_unwrap_space_time.m (misunderstanding)

Hello everyone,
I am a bit confused regarding the following issue:

At uw_sb_unwrap_space_time.m in Line 229 we have: phase_residual=cpxphase.exp(-1i(K(i)bperp_sub)); This makes sense to me.
However, at uw_sb_unwrap_space_time.m in Line 242 we have: dph_space=dph_space.exp(-1iKbperp'); This is not clear to me.

According to my understanding, we calculated K values based on bperp_sub (temporal differences of bperp of each slc) and resphase (edge phase temporal differences). Consequently, I would expect to have bperp_sub instead of bperp in line 242.

Any help/suggestion is welcome :)

processor.txt: permission denied

mt_prep_snap Andy Hooper, August 2017

Amplitude Dispersion Threshold: 0.4
Processing 1 patch(es) in range and 1 in azimuth

processor.txt: Permission denied.

How to solved this problem?

step_2

Subscript indices must either be real positive integers or logicals.

Error in clap_filt (line 71)
ph_bit(1:n_win,1:n_win)=ph(i1:i2,j1:j2);

Error in ps_est_gamma_quick (line 225)
ph_filt(:,:,i)=clap_filt(ph_grid(:,:,i),clap_alpha,clap_beta,n_win0.75,n_win0.25,low_pass);

Error in stamps (line 327)
ps_est_gamma_quick(est_gamma_parm);

Mapping

Hello,
May I ask how to import the results of subsidence rate and cumulative subsidence of StaMPS into ArcGIS for mapping?
Looking forward to your reply.
Thank you.

No mt_prep_isce binary

The mt_prep_isce command is listed in the manual, but is not included in the bin directory

Content of 4.8 in the manual is missing

In the p25 of the manual, chapter 4.8 does not have a content.

4.8 Step 8: Atmospheric filtering