dbekaert / StaMPS

Stanford Method for Persistent Scatterers
http://homepages.see.leeds.ac.uk/~earahoo/stamps/
GNU General Public License v3.0
226 stars 119 forks source link

GEOLOCATION CORRECTION - STAMPS #47

Open yurioliveira-ctec opened 1 year ago

yurioliveira-ctec commented 1 year ago

Dear,

Any ideas on how to update the lon/lat coordinates using the residual topographic information from the estimated PS during the StaMPS-steps? I found a good approach here: https://www.mdpi.com/2072-4292/11/11/1258 I am even using that same article to validate the methodology of my work. The first image below is the result obtained in StaMPS without the lon/lat coordinate correction approach. The second is with the approach, how can I improve?

Image 1 Captura de tela de 2023-07-05 09-22-05

Image 2 Captura de tela de 2023-07-05 09-22-25

The abortion adopted in the article uses the equations below: Captura de tela de 2023-07-05 09-34-49

The parameters I used were: theta - angle of view stored in la2.mat as la phi - heading stored in params.mat deltaH - DEM error that is provided by in rad/m, being stored in scla2.mat as K_ps_uw, to convert to meters I used the K2q.m function. However, due to the non-existence of parameters for Sentinel-1, I adopted the following values:

Captura de tela de 2023-07-05 09-45-44

For the range_pixel parameter present in the K2q.m function, I used the values stored in ps2.mat as ij(:,3).

It is worth mentioning that for the InSAR processing, the DEM SRTM 1Sec HGT was not used, as it presented problems in the download, as shown in the image below. Thus, the Copernicus 30m Global DEM was used.

Captura de tela de 2023-07-05 09-57-10

Thanks in advance if you can help me.

GengPeiLin commented 10 months ago

Hi @yurioliveira-ctec , would you mind to let me know how you get S1-ASC paramters for K2q.m? I would like to do that for decending and other platforms.

cheers,

GP

yurioliveira-ctec commented 10 months ago

Dear @GengPeiLin, how are you doing? I obtained these parameters from a file within the "rslc" folder, which is created during the StaMPS-SNAP process, and from the master image metadata. Slant_range was the most complex, and I still have doubts about it, but I adopted the "range_spacing" parameter from the metadata. The others I obtained from the ".rslc.par" file, and for the height "h", I used the difference between "sar_to_earth_center" and "earth_radius_below_sensor". In summary, it was from these sources that I extracted the parameters for "K2q.m".

Unfortunately, it's been a while since I've been able to study in-depth to improve these parameters. However, if you have any developments or improvements regarding these parameters, could you please share them?

Best regards, Yuri de Oliveira Lima

GengPeiLin commented 10 months ago

Thank you Yuri @yurioliveira-ctec, I am using ISCE2, but your experience is very helpful, I will try to find the equivalent.

regards,

GP