With phase jump I mean strong changes within two pixels, e.g. caused by height changes above the wavelength scale (5 cm for Sentinel-1).
With phase ramp I mean the overall trend superimposing the fringes which is often introduced after unwrapping.
If you have worked with StaMPS before, in step 7, the spatially correlated DEM error is calculated and additionally the phase ramp is calculated when setting the parameter “scla_deramo, y”, how would the phase ramp that may have been introduced during the unwrapping affect the estimated mean LOS displacement velocity?
PS analysis is less prone to unwrapping errors because it constructs a network of stable scatterers. The mentioned ramp is a bigger problem for traditional raster based DInSAR where large parts of phase noise are unwrapped.
I had (stamps_tsexport.csv) and I want to calculate vertical displacement or subsidence from relation :
Based on the reference( Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang, Indonesia by SBAS DInSAR analyses using Envisat-ASAR, ALOS-PALSAR, and Sentinel-1A SAR data)
How can I do this using snsp or what ? please help and thank you in advance
This equation is not very precise, but an approximation of actual displacement.
I don’t know how do implement it in Matlab (would probably be the easiest way), but if you have already exported it to CSV, you can load it into a GIS (e.g. QGIS) together with a raster of the local incidence angle (exported from SNAP) and extract the angles to each PS point. In the attribute table of the points you could then apply the equation based on both variables.
hello ABraun … I want to calculate vertical displacement …Please help me understand this relationship. Is it intended that I always replace LOS displacement with a positive value even if it is negative?
Based on the reference( Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang, Indonesia by SBAS DInSAR analyses using Envisat-ASAR, ALOS-PALSAR, and Sentinel-1A SAR data)
Thank you dear in advance
hello ABraun … I want to calculate vertical displacement …Please help me understand this relationship. Is it intended that I always replace LOS displacement with a positive value even if it is negative?
Based on the reference( Monitoring of long-term land subsidence from 2003 to 2017 in coastal area of Semarang, Indonesia by SBAS DInSAR analyses using Envisat-ASAR, ALOS-PALSAR, and Sentinel-1A SAR data)
Thank you dear in advance
if coh_IW3_VV_09Dec2016_31May2016 >= 0.97 then coh_IW3_VV_09Dec2016_31May2016 else 0
and the result new coherence band had shown only the pixels whose value is a equal or above the condition parameter. but when I applied range doppler terrain correction to the new coherence band, few more pixels have appeared whose value is below 0.97 which doesn’t match the condition.
here’s the data statistics of the input coherence map before terrain correction and the output after the correction
A) What could be the problem?
B) Should i do terrain correction before coherence masking ?
Edit, I have tried changing the image and DEM resampling method from bilinear interpolation to the nearest neighbor in order to preserve the original pixel values from the source image grid instead of averaging the values of multiple pixels and creating a new value into the output corrected grid and yes the problem seems to have been partially solved as there are no more pixels with coherence values that are different the conditional case in the source band. However the minimum and maximum values don’t exactly match the source band
What makes you think that they change? The histogram indicates that the data ranges between -0.11 and +0.18. The color legend in Google Earth might not cover the entire range of values in favour of better contrasts.
There’s a difference between the displacement map produced from Phase-to-displacement operator (bottom one) and the one extracted from bandmath (top one), Also difference in statistics
What is even weirder, is that when transfer the phase to displacement band that was created using the operator, into the unwrappaed product package, the statistics also change and hence the color pattern (using the same exported and imported color palette) also changes, even adding a the color palette itself changes the statistcs. Here take a look please
be careful, these numbers only affect how the data is displayed, e.g. which color is assigned to each value. you can manually modify it or move the slider to change it. But the pixel values stay the same. You can use the statistics tool to calculate the actual range of pixel values.
@ABraun Thank you for your constant assistance in all my inquiries and I hope Excuse me for my ignorance sometimes
I obtained the displacement from an ascending orbit and a descending orbit in order to calculate the horizontal and vertical displacements resulting from the earthquake using the equations
, but when I compare with the values of the measurements of the displacements of two GPS points available in the study area, I find a difference. What is the reason and what should I do to get the same values
GPS is point-based temporal resolution is very high and spatial resolution is very less, in the case of InSAR it is the reverse.
I think you are using only two SAR images which is very difficult to match with GPS. Better try with PS-InSAR (at least three years of SAR data is required) to match InSAR vertical with GPS vertical.