Phase Unwrapping to Vertical Displacement

Hi, why I cannot see the tie-point-grids folder in my ‘‘unwrapped phase product’’ (the last step, after I import it from snaphu). I want to convert the unwrapped phase to vertical displacement by right click on the product> band math> band math nexpresion editor>
can anyone help me please? thanks

Tiepoints

the product you show in the screenshot contains displacement_vv and elevation. That means you already applied “phase to displacement” on this data. Use the product of the previous step, directly after the snaphu import. It contains tie-point grids (as copied from the reference product defined in the first tab of the snaphu import (1-Read-Phase").

Dear ABraun,
Again, Many thanks for your swift reply!
I have one thing else, I did not pin any GCP in my product through the processing of DinSAR, does this mean I missed one big step in the procedure? Do I need necessarily to identify a stable point of X,y coordinates in the image before prior to proceeding? if so, where shall I do that step?

Thank you so much Dear ABraun

no, this is not mandatory. If you manually calculate the displacement, you just need to define a reference height later to define where is zero displacement in your raster.

1 Like

Thanks… where I can define the the reference height in the raster, please? specifically on which product…
Thanks Again Mr ABraun:grinning::grinning::grinning::grinning:

in the band maths, please see here:

1 Like

Ok thanks I understand, but where I shall start to do this procedure? I do not know as I am new user of SNAP. Do you have any guidelines or steps to do so, if so, please forward it to me… With regards

At which point do you struggle?
InSAR is not the easiest technique to start with. Did you see the SNAP tutorials?
http://step.esa.int/main/doc/tutorials/

ASF also offers great ones:
https://www.asf.alaska.edu/asf-tutorials/sar-and-insar-guides/

ESA provides excellent materials on the theories behind it: http://www.esa.int/About_Us/ESA_Publications/InSAR_Principles_Guidelines_for_SAR_Interferometry_Processing_and_Interpretation_br_ESA_TM-19

I know InSAR techniques but I never used the SNAP! I mean the step of Selecting a good stable point in SNAP software? I hope you understand my point… I am so sorry If I bothered you with my questions ?

SNAP cannot do this for you, this is simply a step of logic. If you derive deformation with the equation named above (using the incidence angle), you get a raster showing the relative displacement between both dates. It has no absolute meaning unless you subtract that number from the whole raster.
I tried to explain it here: Subsidence map in 3d view

1 Like

So, my question is, how I do that subtraction you mentioned earlier from the whole raster unwrapped phase map, please? can you show me that by taking a screenshot in SNAP software? I do really appreciate your help with me
Many thanks Dear ABraun :grinning:

First you derive the displacement from the unwrapped phase, either with the module “phase to displacement” or manually entering the equation on the band maths. If the latter is the case, you need to reference it to zero displacement manually.

Please simply read the first couple of posts in this topic, everything is explained there already very well:

A tutorial is given here: Proccessing SAR images
Please also read the notes two posts below, some parts of the tutorial are outdated and no longer required with the current version of SNAP.

1 Like

Many thanks for your assistance. I will follow the topics you have sent it to me… I do not know how to thank you Dear ABraun!

Dear ABraun, I have successfully reached the deformation map of my interest area that we have been discussing on it from last week.

I wonder how I can get all the data (latitude, longitude along with uplift and downlift/cm ) in a text file or excel file? is there any way to import such thing to GIS software in order to classify the data and do some analysis

Many thanks

I have found the same topic you are talking about it our nice forum, I will follow it

Thanks

you apply Terrain Correction on your data (you can include lat/lon bands) and then you can open the img files in the data folder in any GIS.
Alternatively, you can place pins at significant locations and export them (including coordinates and the LOS subsidence) with the Pin Manager, as demonstrated here: Export values from Sentinel 2

The final product should be vertical displacement not LOS! isn’t it?
I have applied the equation as you sated before, (Unw_Phase * wavelength in mm) / (-4 * PI * cos(rad(incident_angle))

I am I right? Thanks

according to my understanding, the equation is just approximating the vertical elevation change. In fact, it is still largely representing the change of surfaces along the line of sight. There is no simple way retrieving it from a single InSAR image pair (confirmed by mengdahl here). johngan nicely explained the difference here.
A comparison between the formula and the phase to displacement module is given here - they are nearly the same.
This was discussed this in several topics:

Two posts naming formulas:

Hi. Can I ask how vertical displacement interprets? Is the line of sight direction considered in this situation? This is my unwrapped vertical displacement (vv) result.
kesilmemi%C5%9F-zoomyap%C4%B1lm%C4%B1%C5%9F100%20- table100

Greetings Braun!

I want to ask you about the formula. Are there any references that certify the reliability of the results that they are produced from this formula?

Thank you in advance :slight_smile: