yumorishita / LiCSBAS2

GNU General Public License v3.0
26 stars 4 forks source link

How to convert LOS to vertical displacement #24

Closed ClementGold closed 1 year ago

ClementGold commented 1 year ago

Hey yumorishita,

I wanted to know if velocity map correspond to the LOS displacement or the vertical displacement ?

In the case of LOS i need to convert with the angle, did LICSBAS2 doing this procedure or i need to project GNSS 3D (EW, NS, UD) deformation data to a LOS direction by multiplying LOS unit vector components (E.geo, N.geo, U.geo) and summing them.

I wish you a very nice day, Clément

yumorishita commented 1 year ago

It is LOS displacement. GNSS 3D data should be projected to the LOS as you mentioned to compare with the LOS displacement.

ClementGold commented 1 year ago

Thanks a lot, i made the hypothesis in my area only vertical displacement occured.

Do you know where I can find good data / paper to help me to achieve this Computation.

Best, Clément

yumorishita commented 1 year ago

Do you mean the GNSS data? It highly depends on your AOI, but NGL data might be useful.

ClementGold commented 1 year ago

Yes i find a gps near this station but I also know i can have the vertical displacement with : DisLOS / cos(i)

I suppose LICSBAS put the angle of the sensors somewhere in the output file ? I think it's the best to determine with one track the displacement

yumorishita commented 1 year ago

GEOCml*/U.geo is the unit vector of LOS and corresponds to cos(i).

ClementGold commented 1 year ago

I have find the file U.geo which is located in every big folder but doesn't find the GEOml file.

The best is to doing this step in the console and output the result but If i export the tiff file and doing the operation in QGIS it's fine too. U.geo file seems to have a special operation to read it

blueyeti1994 commented 4 months ago

Hey yumorishita,

I wanted to know if velocity map correspond to the LOS displacement or the vertical displacement ?

In the case of LOS i need to convert with the angle, did LICSBAS2 doing this procedure or i need to project GNSS 3D (EW, NS, UD) deformation data to a LOS direction by multiplying LOS unit vector components (E.geo, N.geo, U.geo) and summing them.

I wish you a very nice day, Clément

do you know how to read .geo format. I used notepad and QGIS both but notepad is showing weird characters and QGIS is giving me error of not valid raster file.