How to start with StaMPS using Sentinel-1 data

[mehdidar...@gmail.com]

Hi everybody,

I am a newbie on the StaMPS/MTI and going to use StaMPS/MTI for PSI and SBAS processing on Sentinel-1A/B data.

After reading all the Sentinel-related threads in the group I could not finally realize that how can I prepare the Sentinel-1data as input for StaMPS/MTI.

As far as I noticed from the manual only the SLC and IFGs format derived from the DORIS and ROI_PAC are readable by StaMPS/MTI, and if I am not mistaken, non of them has not been developed to handle such task so far.

Of course, I found some codes for reading and pre-processing the Sentinel-1 data (e.g. here and here) but I do not have any idea about that. I do not know whether the SLC and interferograms obtained from SNAP or SARscape could be used or not?   

It would be appreciated if someone could help or give me some hints.

Thanks for sharing your experience.

Mehdi

[David Bekaert]

Hi,

You will need to write a wrapper script that reads SNAP format.

Perhaps other users that have used SNAP or SARscape can let you know more details on how to do this?

If not mistaken, in the current version of stamps you will find such a script for gamma.

Cheers,

David

--
You received this message because you are subscribed to the Google Groups "MAINSAR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mainsar+unsubscribe@googlegroups.com.
To post to this group, send email to mai...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Katherine Smolyaninova]

Hi Mehdi,
To begin with you may go to the SNAP forum and read the topic Workflow between SNAP and StaMPS. http://forum.step.esa.int/t/workflow-between-snap-and-stamps. 
Good luck, 
Katherine

Среда, 15 марта 2017, 17:32 +04:00 от mehdidar...@gmail.com:

--
You received this message because you are subscribed to the Google Groups "MAINSAR" group.

To unsubscribe from this group and stop receiving emails from it, send an email to mainsar+u...@googlegroups.com.

To post to this group, send email to mai...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

All the best,
Katherine

[Katherine Smolyaninova]

Hi, Matteo,

 

Below  is  a  description of what I did to get  results for a case when a study area is within one IW swath. I also attach AndyHooper's modified scripts to prepare  data for StaMPS (mt_prep_gamma and ps_load_initial_gamma) which takes into account  geocoding  of pixels from SNAP. Andy kindly allowed to share these scripts. I think other people can use my experience to save their time.

 1/  S-1 images:   Split and Apply_Orbit_File. Split is done individually for each image because number of bursts covering the study area vary from image to image and you have to adjust this procedure manually. No batch processing can be used at this step. Apply orbit can be done using batch processing.

 2/Create stack using  Backgeocoding (   Radar-Coregistration-S-1Tops_Coregistration-S-1Backgeocoding. Master should be the first in the ProductSetReader). 12-13 images are OK for each Stack. If more images in the Stack - it is better to split it into several Stacks each having the same master. Master can be chosen by Radar-Interferometric-Insar Stack Overview.

 3/Deburst the result of 2/

 4/ Make subset of the result of 3/ (geographical coordinates) (Raster-Subset). Result file:  subset_master_Stack_Deb.

 5/  Apply Interferogram formation to the result of 4/ -.  Result file:   subset_master_Stack_Deb_ifg­.

 6/  Apply TopoPhaseRemoval  to the result of 5/ . Result file:  subset_master_Stack_Deb_ifg_dinsar.

 7/  Add elevation band to the result of 6/ :  subset_master_Stack_Deb_ifg­_dinsar .

 8/  Prepare data for StaMPS:  Radar-Interferometric-Psi/sbas-StaMPS export

 I exported to StaMPS  the two files:  subset_master_Stack_Deb_ifg_dinsar  and  subset_master_Stack_Deb . After  the Stamps Export procedure I got four  directories /diff0, /geo,/dem,/rslc  which I placed into the directory INSAR_master_data. In case of several Stacks export each Stack and then copy results to the directories directories /diff0, /geo,/dem,/rslc.

 9/ To have PS geocoded add lon and lat bands to the result of 4/ ( subset_master_Stack_Deb). (Right click -band maths-edit expressions-constants-lon (then lat) (uncheck "virtual" to save expression). Export result to ENVI or GAMMA format. Save result, for example, to /geo directory.  Rename files with lon and lat to master_data.lon and master_data.lat. These are   binary files with latitudes and longitudes for all pixels of the master crop.  (One file contains values of lat for all pixels and the other - of lon) All the other files (results of export of the other bands)  can  be deleted. The files  master_data.lon and master_data.lat should be placed into /geo directory.

 10/ Place the new scripts mt_prep_gamma to the /bin directory of StaMPS and ps_load_initial_gamma.m to the StaMPS matlab directory.

 11/ run mt_prep_gamma from the INSAR_masterdata directory.

 12/ run matlab and  launch StaMPS.

Hope this helps
Katherine

 

 Good luck!

[Mehdi]

Dear Katherine and David,

Thank you so much for your guidance and useful advice.

Best,

Mehdi

[Antobel]

Hi Mehdi, 

I try to produce STAMPS output by SNAP, if it is possible I want some informations

1: Is it necessary to create a subset of images? or 

Can I process backgeocoding_stack to obtain ifg stack?

2: Is  ifg stack  created directly using 'Interferogram formation'  to backgeocoding_stack?

Best regards

Antonella

[rana virk]

Hello Katherine,

Please send me complete flow to install StaMPS. I am less familiar with linux.

I have performed DInSAR using SNAP and now wish to perform PSInSAR.

Regards

ASVirk