Information about vel file, vel filt etc and exporting in tiff

[ClementGold]

Dear Yumorishita,

I had a interrogation about the data which i extract from my result file. I wanted to put the subsidence in a area in mm/year so I did the following command to get the tiff : LiCSBAS_flt2geotiff.py -i TS_GEOCml1GACOSmaskclip/results/vel -p TS_GEOCml1GACOSmaskclip/info/EQA.dem_par

Because the vel.png show me as well a displacement by year.

But in QGIS with the tiff file it seems I have the same map as this one :

As I understant it's the velocity from the reference date and not a displacement by year.

Compared to this picture which is the displacement every year (the one who interrested me)

What is the input i need in the LiCSBAS_flt2geotiff ?

I wish you a great day. Clément

[ClementGold]

It seems I made a mistake the vel is the displacement / year so it seems do have the same unit and data

[yumorishita]

So, what is your question now? All solved?

[ClementGold]

I just wondering is the vel I export in tiff was the velocity filtered (blue value in time series) and I just now using vel.filt because it's more robust with less noice.

So yes you can reply and confirm my wondering as well closing issue :)

one extra question do you have any link to the deramp removing and why it's important to study urban ground deformation I find one paper only explain orbital error is removing.

[yumorishita]

vel and vel.filt are the velocity before and after filtering (step1-6), respectively. I have no specific reference about the deramp because it is quite a simple and standard method. I don't think deramping is necessary for Sentinel-1 in general because its orbital error is quite small. What paper do you mention?

[ClementGold]

It's the paper Characterization of land subsidence in Tabriz basin (NW Iran) using InSAR and watershed analyse but you're right it's from ENVISAT data.

Do you know where in the code I can apply the velocity * cos(i) In my area the mean of the angle it's like 38.6° corresponding to 0.66 radian I just try to find a better way and more properly.

[yumorishita]

gdal_calc.py would be useful.