I have used external DEM 10x10m in the SNAP S1 Toolbox for InSAR processing, but it always failed.
The data that I use (Sentinel-1A images and external DEM 10x10m) already same coordinate system.
Is anyone know why I always failed to process it?
Is it because the data that I use have different spatial resolution?
steps that I followed:
open product (master and slaves), split, deburst, subset area, coregistration, interferogram formation and topographic phase removal. In topographic phase removal, I choose external DEM. Then, the process was error.
Dear colleagues, I face the same problem with external DEM. Working with external pieces from SRTM and ASTER the results is fine, but with the external one I have the result is black images in all bands. Before start I followed the recommendations to use WGS84 - no luck. Any comment/experience is welcome.
Thank you for reply. For the case I report here I’d like to use DEM having 10x10m spatial resolution. The file was reprojected in QGIS and works fine there and with other software too.
One possibility for this issue could be the small area covered by this DEM. Could you advise please.
Dear ABraun, sorry for the delayed reply. In QGIS I used Raster->Projections->Warp dialog and setting appropriate parameters -tr switch. After processing I got 30x30m file and worked fine with SNAP.
Actually the problem was the size of DEM tiles (4x4 km). I started with only two tiles (32 sq.km) and it didn’t work, but getting more (12 - 192 sq.km) it worked fine even without reprojection.
My guess is that the reprojection of the 2m DEM data to WGS84 has not been correct and now SNAP treats it somehow wrong.
As the elevation data is a asc file with no metadata you have to make sure that you select assign the correct coordinate reference system (UTM 35 N) before you reproject.
Edit: It took 135 minutes (16 GB RAM) but it worked. The upper image is with topographic phase removal based on the 2m DEM, the second is without.
no matter what QGIS thinks about the projection of the data (it assumes WGS84 after loading, but it is not), you have to select the correct UTM zone in the reprojection module (here it is 35N) as the source CRS and WGS84 as the target reference system.
I did exactly same before and now , but still the slave is empty, I’ll send you the DEM, and both S1 identifier, as private message, hope you have time, many thanks in advance. I also tried up 10 m *asc,. but still similar empty slave. But using GETASS30 autodownload, there is no problem.
just for clarification: After which step is the slave empty?
First I thought you wanted to use the 2m DEM for the topographic phase removal (because this was the problem in the topic you linked).
I don’t think it makes sense to use a 2m DEM for the back-geocoding of S1 SLC data. Especially if the covered area is much smaller than the S1 data. Use GETASSE30 and it works.
Both tiles of the DEM covers the only first burst, which covers the AOI, The slave is empty after TOPSAR corr. GETASSE30 works, but doesn’t implement the demand, I explained to you, the goal only to test the DEM for the further processes.