#######################################################
Warning
The workflow described in this post (and the ones referring to it later on) is outdated and no longer recommended, because new versions of SNAP and StaMPS have been released. Please have a look at this updated list of instructions: StaMPS - Detailled instructions
#######################################################
Hi @Alex_P, here the steps I had followed,
- Apply TOPS split and then Apply orbit file.
- Add scenes (Master+Slave1 and Master+Slave2 and so on ) to do the Backgeocoding separately for each Master-Slave pair.
- Create Stack of Backgeocoded images (master_Stack1, master_stack2 and so on).
NB: I have used SRTM 1sec HGT (Auto download) DEM and Cubic-Convolution method
… - Create a stack of Interferograms (master_stack_ifg1, master_stack_ifg2 and so on for all pairs), Deburst and Topographic phase removal (master_Stack_ifg_Deb_dinsar1, master_Stack_ifg_Deb_dinsar2 and so on).
NB: If you see EAP error then you have to apply EAP for the scene before Back-gocoding (you have to apply EAP either for master or slave that you have to check, images before March 2015 need EAP correction) - Make sure to deburst both the Backgeocoded stack (master_Stack_Deb1, master_Stack_Deb2 and so on ), and the Interferogram stack already debursted**(master_Stack _ifg_Deb_dinsar1, master_Stack _ifg_Deb_dinsar2 and so on).**
- Add an elevation band to the Interferogram stack for all (by right-click > add elevation raster).
NB: I have used SRTM 1sec HGT (Auto download) DEM and Cubic-Convolution method - Add the Debursted Backgeocoded stack first and then Debursted Interferogram stack to the stamps export.
NB: If you do subset then do it with same Geo-coordinates for master_Stack_Deb1_subset and master_Stack_ifg_Deb_dinsar_subset (you can also do subset before dinsar to save some time)
…
Export to StaMPS the two files: - master_Stack_Deb AND master_Stack _ifg_Deb_dinsa
- StaMPS Export will generate four directories (/diff0, /geo, /dem, /rslc)
NOTE: For all the Stacks export each Stack and then copy results in the directories (/diff0 and /rslc), no need to copy for (/geo, /dem) is same for all pairs. - To have PS geocoded add (Lon & Lat) bands to the result of 5 (master_Stack_Deb) [do it only for one pair]
NOTE: (Right click -band maths, then edit expressions, select lon from constants (lat in the second time). It will add Lat and Lon band in your master_Stack_Deb data and then Export it to GAMMA format. Now rename the lon and lat file as 20150118.lon and 20150118.lat (20150118 is the date of your master data), copy this two file and place into /geo directory.
… - Place the new scripts mt_prep_gamma_snap to the /bin directory of StaMPS (you can rename it as mt_prep_gamma and replace the existing file).
- Place ps_load_initial_gamma_snap.m to the StaMPS /matlab directory (rename ps_load_initial_gamma.m and replace the existing file).
- Run mt_prep_gamma from the data directory (cd to the folder where you have your all exported data).
HINTS: mt_prep_gamma_snap [masterdate in format yyyymmdd] [data directory] [amp_threshold] [rg_patches] [az_patches] [az_overlap] [rg_overlap]
{for instance: mt_prep_gamma 20150118 /Insardata/Stamps_export 0.4 1 1 50 200}
NOTE: You can use mt_prep_gamma_snap (changed) and ps_load_initial_gamma_snap (changed).m provided by FeiLiu. If the script ps_load_initial_gamma_snap.m and ps_load_initial_gamma_snap.m provided by Katherine are workable then everything is also OK. - When the process is finished, type matlab in the console and in the matlab type stamps.
I hope this will help
Cheers