Thank you for the help @mdelgado But could you please tell me how to do this?
Thank you for the suggestion @ABraun. I am short of 5 FBD images. Can I use FBD and FBS ones in that case?
You can combine FBD and FBS as long as both have one common polarization.
well, as you already have the unwrapped phase, you could open it using any mathematical software such Matlab or python and calculate the plane that minimises (best approximate) using least squares or so. Then once you got that plane, remove it from the unwrapped phase. Note that then you could have still other contributions in the phase and not only deformation, but you will be one step closer.
If you plan to use PSI or SBAS then this probably is not needed as that kind of processing normally correct those effects.
I hope this helps
Is there any progress regarding the issue with the shifted elevation and topo bands?
I have a similar problem of shifted bands, in the case of processing two ALOS SLC images (26/12/2006 – 31/12/2008, Bperp=255). SRTM 3arc was used in this case.
What drove my attention to check for this potential shift, in the first place, is the fact that before topography removal the interferogram was like that:
And after topo-removal:
Something seemed to have gone wrong with the topography removal tool. I then read @mdelgado post…
So in cases like this what should we do? Just discard the pair?
Please note that for this particular case, the two acquisitions refer to different reception paths. Could this be the problem?
Thank you all in advance
Have you applied ALOS Deskewing before the topographic phase removal?
I wonder if it makes a difference if you do it before or afterwards.
thank you for your quick response.
No, deskewing was not applied since the images were provided by ESA (Alos-1 DInSAR Deformation Map). ALOS support team confirmed that the products do not need deskewing.
I see. You are right then, of course.
I am also interested in the issue, but it seems to be of technical nature.
Just to exclude this as an error source: Did you compare the results of topo-phase removal performed during the interferogram processing (using the checkbox) with calculating it as a second individual step after the interferogram was written (using the separate module)?
Unfortunately the problem remains whether topography removal is performed in one step (checked box) or two steps. The following figure is from one step process (only elevation band can be calculated).
In fact, as already discussed in ALOS SLC L1.1 interferometry processing parameters, even trying with higher orbit interpolation degrees, results are awkward, don’t you think? I would have thought that the 5th “orbit interpolation degree” could be safely selected in every processing (results should be more accurate than using for eg 3rd orbit interpolation degree). Now, it seems that a major problem apart from the orbit precision, is related with the topography removal tool… Is this a SNAP processing error or I should try something else?
What is the actual problem - low quality of ALOS-1 orbits or something in SNAP? I’m not sure if going for the highest possible orbit interpolation degree is physically the correct solution. Perhaps this issue requires baseline-refinement, which we plan to add into SNAP in the future.
mdelgado summarized it in the post above. The integration of elevation data for topographic phase removal reduces the quality of the interferogram.
It seems that orbit interpolation degree affects the problem encountered with the shifted bands and the topography removal performance.
Using the highest orbit Interpolation degree 5, in the two- step processing (1st step for subtracting flat-earth phase, 2nd step for topography removal), bands are not shifted and results are like this:
However, would you trust such interferometry results?!
Would you like me to upload the respective images? is it possible to upload here?
I would like to share with you that in every ALOS pair I processed, I noticed that in the case of using 3rd orbit interpolation degree to remove topography (either in one step -checked box- or in two steps), topo-band and elevation band appear shifted with relation to the differential interferogram produced. On the other hand, when using 5th orbit interpolation degree, bands don’t appear shifted. However, when comparing the two differential interferograms produced (one with the 3rd and the other with the 5th degree) should not they appear shifted between each other as well? Because they do not… and they are not yet geocoded. Could this all thing with the shifted bands be just a visualization error of SNAP?
In the following figures, the left images are snapshots of a 2-year differential interferogram after using 5th orbit interpolation degree and the right ones are after using the default 3rd degree, for topography removal.
Is there a way to check if topography is safely removed in either cases? Do these products-especially the 5th degree- look normal to you? Would it be better to wait for future baseline-refinement?