I tried to make DEM by using TanDEM-X products and then plan to compare with one DEM like ‘ASTER 1sec GDEM’.
I read some steps in STEP forum and even follow this utube toturial
I did like below:
Making interferogram from TDM-CoSSC-Experimental
Phase to height
Geocoding (map projection: UTM zone, parameter projection:4north)
I do not know which kind of DEM, I should use in geocoding section?
Actually I used one product from ‘ASTER 1sec GDEM’ and then compared my result with it over land but results were not the same (figure1).
I would suggest that you try “phase to elevation” instead and use the ASTER DEM as reference so both have the same base.
If both are in one product you can also create a scatter plot and see if the shift between both rasters is constant or just affecting values within a certain range.
I used this ’phase to elevaton’ option in SNAP and then I used ‘External DEM’(figure1) but the problem is for ASTER data, there is some error values over sea (figure2 and 3) that I removed them by MATLAB (figure 4) and then used it in SNAP but when I made a DEM from TanDEM-X, I also chose ‘DEM’ option there (figure5) but I do not know why DEM values over sea are (-1,…). (Figure 6)
It is weird, it should be zero and I do not know this minus values are coming from where?
Yes I can remove them but my problem is not only with this DEM. It is easy to remove it in DEM but my problem is with DEM that I made by TanDEM-X images by using this DEM. When I made two DEMs for two different days and subtract from each other then I got minus value over sea.(fig1)
I discussed with some of my colleagues about this DEM produces with SNAP and they guess that this minus values over water that actually are sea ice is here with high coherence could be a error in SNAP. So, I asked SNAP developers care about this issue and tell us their idea about it.
I do not know who are SNAP developers but I add some people that I think they are. If you know more people, please tag their names here. Thank you
If you download ASTER DEM in that area I mentioned (ASTGTM2_N71W157), then you can see, it has some error values over sea(figure1) but we do not have any DEM over sea and it should be zero or maybe NAN. That’s why I used MATLAB and convert all ASTER DEM values over sea to zero (figure2) or NAN (figure3) and I used them in ‘phase to elevation’ and ‘geocoding’ but I do not know why I got minus values over sea. In my idea, it should be zero, NOT minus.
These minus value are coming from where?
I did not use sea level here? how should I use it in calculation?
I only convert ‘unwrapped phae’ to ‘phase-to-elevation’ and when I checked my result with ASTER DEM over land, they are same but they are not same over sea.