So, basically forming the interferogram after subsetting works perfectly!
Now, I have these results (Image attached). Now, what can I do next? How can I convert phase to displacement? Any further processing should be done here?
I tried to use “ps_plot(‘v-do’)” command but I am not getting all the images as I was expecting. What can I do to get all images? (Image attached of v-do)
Good evening, I finally managed to finish the PSI analysis in which I worked with 56 images, take SNAP’s recommendation as a master image. When visualizing the resustados they do not generate confidence, the time intervals between one image and another is 1 month, it is my time series graphs I have very pronounced changes.
Any recommendation?
It depends on your research question how you want to to analyze the results.
Extract from the manual
'w' for wrapped phase
'w-d' for wrapped phase minus smoothed dem error
'w-o' for wrapped phase minus orbital ramps ('w-dm', 'w-do', 'w-dmo')
'p' for spatially filtered wrapped phase
'u' for unwrapped phase
'u-d' for unwrapped phase minus dem error
'u-m' for unwrapped phase minus and master AOE
'u-o' for unwrapped phase minus orbital ramps
'u-a' for unwrapped phase minus topo-correlated atmosphere
('u-dm', 'u-do', 'u-da', 'u-dmo', 'u-dma', 'u-dms', 'u-dmao', 'u-dmos')
'usb' for unwrapped phase of small baseline ifgs
('usb-d', 'usb-o', 'usb-a' also 'usb-do','usb-da', 'usb-dao')
'rsb' residual between unwrapped phase of sb ifgs and inverted
'd' for spatially correlated DEM error (rad/m)
'm' for AOE phase due to master
'o' for orbital ramps
's' for atmosphere and orbit error (AOE) phase due to slave
'v' mean LOS velocity (MLV) in mm/yr ('v-d', 'v-o', 'v-a', 'v-do,
'v-da', 'v-dao')
'vs' standard deviation of MLV (mm/yr) ('vs-d', 'vs-o', also 'vs-do')
'vdrop' MLV calculated from all but current ifg (mm/yr)
'vdrop-d'
'vdrop-o' (also 'vdrop-do')
good job! I think intervals of one months are not bad because they allow you to cover large periods with a reasonable amount of images. And despite of single spikes in your graph, you can still detect trends with the third drop-down menu (exemplified here)
You can also exclude large outliers (if the same date is faulty for all points), by removing it from the StaMPS processing by the drop tags described in the manual.
When I plotted ps_plot(‘u’), I got all the unwrapped interferograms (17 for my data) in one image.
I need all the interferograms in one image just like ps_plot(‘u’), but instead of phase angle, I need all values in mm or m. I am trying to get these images just for display purpose. Analysis will be based on time series.
thank you ,I’v solve the problem successfully even the result is no good.however i have trouble with configuring snappy of SNAP - StaMPS Workflow - 1.1 code(@thho etc.recommanded):when I do step 2.3.1 “Configure snappy SNAP - Python interface” after “1/sudo apt-get update,2/sudo apt-get install python-jpy,3/cd /home/yzx/snap/bin, 4/./snappy-conf /home/yzx/anaconda2/bin/python2,5/cp,6/spyder”,terminal seems frozen ?
Ask for help!
I have finished all the steps of SNAP and StaMPS,and get the plot.But i want use these results to other applications for further analysis.But i don’t know how to convert the format of the result.Please give me some hints,thanks!
Is ESD can increase the precision of PSI about StaMPS to work after pre-processing in SNAP(before deburst after stack)?I read through the forum finding whether or not implement ESD step is ambiguous ,someone is is ingored where other masters recommanded
Hi @ABraun, instead I have a doubt about TOPSAR Deburst in the PS processing. Is it better apply it before or after Interferogram formation? I read different points of view about this. Thanks in advance.
it would depend of the conditions, as @ABraun says, if you work with only 1 burst, you surely do not needed.
From once side, and following your above sentence,… I would ask who is this someone? If masters recommend it, I would trust masters, as they normally do science since decades, and this someone could be introducing error in his and your measurements if his assumptions are wrong.
From other side, if this someone, tried, tested, and published recognised works doing it in that way… then give a try. Again… it depends if you then do any type of check or do post-processing to verify that the assumption taken is valid.
From my understanding when you use several burst and you want a good coregistration for InSAR using Sentinel-1 TOPSAR, you need at 1/1000 pixel accuracy, that can only be reached using ESD over areas with medium/high coherence level. This I have listened from JPL, DLR, ESA and other so-called ‘masters’. So I follow their advice.
However, and depending of what is your goal, you can always try both options and evaluate their implication in the final results, so you can draw your own conclusions. And then share it with the community.