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:
- run dispersion.csh in the SLC/ folder to make amplitude dispersion index (full resolution).
- modify filter.csh to generate interferograms with real.grd and imag.grd file (full resolution).
- make a stack/ folder parallel to intf/, raw/, topo/ and SLC/
- run link.bash to move all the real.grd and imag.grd into stack/
- modify/run cut.bash to crop the real.grd and imag.grd to appropriate size and region
- link the scatter.grd to stack/
- modify/run mt_extract_info_gmtsar
- modify/run mt_extract_cands_gmtsar
- 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));