Subsidence map in 3d view

thank you, @heinzollerketchup.
Generally, the longer the temporal baseline (=the time between two images) the higher is the decorrelation which leads to small coherence values. It also depends on the cover of the analyzed surface. Bare ground mostly shows slower decorrelation than one with high vegetation cover.

As you can see here (tested with S1 data), coherence decreases from 0.7 to nearly 0.1 within 6 12-day-cycles (12, 24, 36, 48, 60 and 72 days, blue to red)


Please also have a look at the ESA tutorials:



thanks again, your post just explains my second question. My study area (in Mongolia) and has a low coherence even in a 12 day time gap (mean 0.4). So according to your response this result is influenced by the surface (shrubs, woodland, grasland, bare soil) and not by the time interval.

Is there a possibility to increase the coherence in areas with covered earth surface?

The wavelength of Sentinel-1 is about 5 cm and strongly interacts with smaller vegetation cover such as shrubs. Ich these shrubs slightly change between two image acquisitions coherence rapidly decreases.

It cannot be increased by methods of data processing. Higher coherence can be achieved through precise orbit parameters, shorter time intervals and longer wavelengths (ALOS PALSAR with L-band, ca. 21 cm is less sensitive to smaller surface cover changes).

So all you can do with sentinel is applying the orbit files before your co-registration (also increase the number of GCPs) and hope that the coherence is of better quality.

1 Like

Ok, I will give it a try.
Thank you!

"Select a point where you assume no change and read the value. Let’s say its 45 cm

Subtract 45 from your imagein order to have zero values at those areas with no change. The remaining variation should be due to subsidence between two overpasses."

can you explain it briefly … you say “45 cm”… from where I get this value for my case… is it from phase or coherence …from coherence map> properties> pixel information tab>band… little confused

this value is hypothetical, it replaces any value you measure at your location where you don’t expect change. You measure it from the unwrapped phase.

ok I understood … but ABraun can you show me procedure of subtraction step wise … that you mentioned

"Subtract 45 from your imagein order to have zero values at those areas with no change. The remaining variation should be due to subsidence between two overpasses:

Assuming that your unwrapped phase is “Unwrapped Phase_ifg_srd_VV_05Dec2008_24Apr2009”
and your reference height is 15:

Right-click on your unwrapped phase open Band maths and enter
Unwrapped Phase_ifg_srd_VV_05Dec2008_24Apr2009 - 15

If the measured value is negative, add it instead.



in this picture pixel value show 2o8 … what does it mean`

if you expect no subsidence at this location it means that your subsidence map is 208 mm too high. If there was no change at this pixel there should be a zero. After unwrapping your product a product of relative changes within the raster, given as absolute values (mm in this case if you applied the formula correctly).

So make your raster 208 smaller and you should theoretically have the absolute height changes during the two overpasses.

However, there’s a strong ramp in your image wich distorts the result massively. I don’t think you did anything wrong so far but maybe there’s too much time between your images.

1 Like

is there any option in sentinel 1 toolbox by which i can get the information of incidence angle of satellite of any radar radar image.

There is the incidence angle tie point grid and also the incidence near and far in the abstracted metadata.

i am trying to unwrapped the phase on VMware … it seems a never ending process… it shows treesize 33876567 , pivots 1376897546, improvements 11567898… what does these term means… is there any way to know how much time it takes…

1 Like

Going back to the post of Sept. 8-th.

How to find the “incident_angle” in order to process the data?
Is it in the metadata?

hello ABraun
is there any option in sentinel 1 toolbox from where i can get the wavelength information of my terrasar image…

you said that the crop possess of eliminating the pixel that have low coherence will be done on unwrapped phase … but when i do this process on unwrapped phase … a error message is shown…

but when i did it on coherence image only the required pixel remain no error message shown… is it correct…

I can’t help you when I don’t know the error message.
As you work with the information of the unwrapped phase you need to mask out the locations with low coherence there, of course. It’s kind of map algebra with two rasters: If coherence is low, set unwrapped phase to nan/invalid.

Wavelength is 3.1 cm - I didn’t find it on the rush in the metadata, sorry.

and where i can get the information of incident angle in sentinel 1

incidence angle is dependent from the distance to the sensor.
So there is always a difference between the far and the near range.:
You find it under Abstracted_Metadata:

If you want the incident angle of a specific pixel, you can open it as a Tie-Point-Grid:

can u please check weather … the steps I did is right??

  1. this is my unwrapped phase

  1. then I did the band math …

  1. here is new image after doing the band math calculation

  1. then I am opened up the coherence image … and did the crop out process of eliminating pixel that have low coherence value…

then I am identifying the stable pixel, the highlighted one…

what will I do now…