Skip to content

this is a set of programs that enable you to run StaMPS with output file from GMTSAR or GMT5SAR

Notifications You must be signed in to change notification settings

dedetmix/gmtsar2stamps

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

5 Commits
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

gmtsar2stamps

This is a set of programs that enable you to run StaMPS with output file from GMTSAR or GMT5SAR It has been tested on ERS and ENVISAT and ALOS-1 data. Xiaopeng Tong, University of Washington, Seattle, Feb 2016

Steps:

  1. run dispersion.csh in the SLC/ folder to make amplitude dispersion index (full resolution).
  2. modify filter.csh to generate interferograms with real.grd and imag.grd file (full resolution).
  3. make a stack/ folder parallel to intf/, raw/, topo/ and SLC/
  4. run link.bash to move all the real.grd and imag.grd into stack/
  5. modify/run cut.bash to crop the real.grd and imag.grd to appropriate size and region
  6. link the scatter.grd to stack/
  7. modify/run mt_extract_info_gmtsar
  8. modify/run mt_extract_cands_gmtsar
  9. start matlab/stamps in stack/ and run stamps(1,8)

Because of the way the phase data is extracted here, you need to modify the way Stamps read the data.
The original code has the real and imaginary data ordered such that real is in odd position and imaginary is in even position. When I used GMT to read in the SLC grids, I read all the real number first and then the imaginary number.

The phase is loaded in ps_load_initial.m, ps_load_initial_gamma.m or sb_load_initial.m. So you need to change the code:

original: ph(:,i)=complex(ph_bit(1:2:end),ph_bit(2:2:end));

modified: ph(:,i)=complex(ph_bit(1:n_ps),ph_bit(n_ps+1:end));

About

this is a set of programs that enable you to run StaMPS with output file from GMTSAR or GMT5SAR

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Shell 100.0%