Workflow between SNAP and StaMPS

My congratulation to you Mr. @falahfakhri, for the New Language and nationality Finnish.

1 Like

Thanks a lot,

1 Like

Hi Andy,
I’d like to ask you, in case of adding elevation by checking the box, is there any need to right click and add elevation band?

No, checking the box is enough.

1 Like

Hi ABraun, I have read your answer from some of the issue above that very helpful for me thank you so much,
I have a question for the Snap processing, after I Back-Geocoding create stack–> I go to subset but I can not see anything in the subset area tool ??
also I have that the same problem when I processing image in computing Interferogram I can not see the phase in the image (the phase show nothing)
Could you help me explain that thank you so much
I use five images with 1 master and 4 slave images


thank you so much for your time
Computing Interferogram image example

your first screenshot indicates that the file is too large for SNAP to produce a quick overview. But as shown in the second, the data is not correct as well.
Please go through all input files (_split_Orb) and visually check if all contain correct data.
As for the stacking, I recommend to only use VV because VH doesn’t bring you any benefit regarding phase information (I made an example here: Effect of polarization on interferometric products). Furthermore, it reduces the file size of your stack by 50% and for InSAR you can use only one polarization anyways.

Finally, you can create an RGB composite after back geocoding using intensities of different dates to see if all were included correctly.

Hi ABraun !! Thank you so much and really appreciate for your answer , I will try again about that.

Hi ABraun, It’s me again after I try to reduces the file size there is still that problem and also when I click to Intensity of Slave image that will not show anything could you explain that again
second you have mentioned visually check all contain data ?? How could I check that in my data thank you so much for your time


the left side I open Intensity of slave image and al there is viewing map of my raw satellite image:
If I try to subset image without quick overview there is error but I don’t know exactly what that imply ?

thank you so much “Tran Rick”

1 Like

Hi Tran,

I’d like to ask you two questions regarding to your processing,

The first one, Why you didn’t apply ESD?

and the second, why you don’t postpone datasubset after math Deburst step?

my suggestion was to open all intensity bands after you applied split and orbit files. If all show an image, you can proceed.
You can check if all bands of your back-geocoded stack contain information by making RGBs as described here: Fault marking by InSAR Data and difference between SLC_1SSV and SLC_1SDV?

If the images are not properly aligned, the interferogram will fail. Based on your information, it rather looks like some dates remain empty after back geocoding.

All the workfow needed by Trannick can be done automatically with the snap2stamps package.
See our latest publication using it on: https://www.mdpi.com/2072-4292/11/2/129
snap2stamps package available at: https://zenodo.org/record/1322353

All this is open and free just for helping scientist to re-invent the wheel, and go faster on their researchs.

I hope this helps!

2 Likes

Hi Falahfakhri, your first question (Enhance spectral Diversity ) provide a pixel or sub-pixel short to improve co-registration, but I’m not sure here"Just follow some step from this workflow between SNAP and Stamps in this forum
Second Question Datasubset, I think that will the same, because subset just help you choosing interested area for your study, but in my case I can not debursting for next step facing with error. I hope answer some point of your two questions
Thank you

Hi ABraun, firstly thank you so much for your answer, I did try about your instruction (my Stack image show just red color, and I try other image include some of images before stacking all show red that not look the same with image in your link provided) is data have problem or from my processing (I did try again four times are the same result)
what should I do for this problem ?


Thank you so much

Hi mdelgado, Thank you so much for your great document and I have read your paper (you mentioned the automatic processing chain process more than 160 sentinel 1A image ) that’s very big number is that we using this tool to process all image in 1 time ?? do you know where can I find a web can download image faster because in https://scihub.copernicus.eu/dhus/#/home take long time to download 1 image.
Second, in installation your package (in my computer I already install Python 3.6 configure with snap image below ) we do not need to install Linux system my computer run in window system? and also do we need

follow steps in this forum for Stamps method ??

The goal of my questions was to draw your attention to the processing sequences, for instance in your case to avoid any error, I’d suggest step by step process for each single product,
Split ----》apply orbit------》TOPSAR CORR-----》backgeocoding (S1 TOPS GEOCODING)------》Apply ESD------》DEBURST-----》EXPORT TO STAMPS

THE second item you need to be exported to Stamps, is
Start from ESD step to create INFGRM ADD LONG AND LAT BY CHECKING THE BOX------》DEBURST -----》EXPORT TO STAMP

In this case you would avoid any source of error, many steps could be done at once , but I’d prefer step by step.

Hi @Tranrick,

In the package user manual it is written as requirements python 2.7 (with pathlib package) and SNAP 6.0.0 or superior. Reading the manual you will see that you only need:
1)download all the S1 data and put it inside PROJECT/slaves folder compressed (zip files)
2) to set up the master (splitting it over the AOI)
3) set the parameters inside the file project.conf and ,
4) running the steps

@ABraun uses it on Windows, so there is no need to install linux inside the windows machine, at least for the snap2stamps processing. For the StaMPS part you should have other things in consideration (@ABraun can point you some threads on that)

For downloading S1 images you have alaska vertex, cnes-pep and scihub. Good luck with them! Normally I jump from one to another finding the data or the fastest download speed, but this is normally changing.

I hope this helps!

Of course you set all colors in the RGB to three different dates. This tells you which of them is maybe missing.
Alternatively, you can simply open all intensities after creating the stack and check if they all contain data.

As for snap2stamps on Windows, there are not much things different. I can provide you my config file to compare later this week (currently not in my office).

Hi Andy,
For step by step preparing data in SNAP 2 STAMPS, which is better to implement create STACK after DEBURST of ESD backgeocoded image and same to infgm?, or to create STACK after implemented backgeocoding for all pairs and then apply ESD and DEBURST?

ALSO TO check the box LAT and long for each single infgm?

well, if you use the snap2stamps scripts, you don’t need to make any decisions because it does all the work for you in an order that works.
If you want to do the same steps manually, I advise to stick to this order:

Capture
source: manual

as for the lat, lon and dem bands, snap2stamps creates them for each interferogram but in the end, the stamps export needs them only one time.

1 Like

Do you mean the python script?