kite icon indicating copy to clipboard operation
kite copied to clipboard

look vectors problem

Open CUMT530849827 opened this issue 1 year ago • 6 comments

I processed gmtsar data with kite and found that the look vectors for lifting the rail were incorrect. I processed two earthquakes in a row and found that it was the same in both cases. Could you tell me how this happened? How should I solve this problem https://github.com/pyrocko/kite/issues/73#issuecomment-1648016837
Here are the look vectors for the two earthquakes I worked with image image The following is the direction of the lift vector given in the case list image

CUMT530849827 avatar Nov 16 '24 12:11 CUMT530849827

Yes, so please make sure that the los vector is correct. You need to rotate your heading angle by 180 degree.

hvasbath avatar Nov 17 '24 19:11 hvasbath

@hvasbath @miili I calculated the enu of the ascending and descending orbits and found that there was no error here. I searched for many reasons and did not find it, but I found https://github.com/hvasbath/beat/issues/93 that the ascending vector in this problem was the same, and I wondered if there was a problem with the drawing script image image image https://github.com/hvasbath/beat/issues/93 The picture in question is below image

CUMT530849827 avatar Nov 18 '24 12:11 CUMT530849827

First of all there is no point in showing plots of BEAT. Please open the kite scenes with spool and check the look vectors there. BEAT does nothing but reading from Kite, if there is sth wrong then happens in the import from GMTSAR. I dont understand what the first column of your screenshot is supposed to be?

hvasbath avatar Nov 18 '24 12:11 hvasbath

I used spool to open the track lifting data, but I don't know how to check their look vectors image

CUMT530849827 avatar Nov 18 '24 12:11 CUMT530849827

You can select phi here, but already from the formula in scene_io we can see there is sth wrong. image

image

For ascending track this arrow you marked must look into the exact other direction- so rotated by 180 deg.

So the error is in the equation you marked just here: image

Until that is fixed I suggest you load your scene through a python script and add pi.

from kite.sceene import Scene
from numpy import pi

scn = Scene.load("your_scene_name")
scn.phi += pi
scn.save('kite_scene_outname')

Then please check again in spool. And once that is correct- please reimport to BEAT.

hvasbath avatar Nov 18 '24 13:11 hvasbath

Thank you very much for helping me find the mistake so kindly and carefully

CUMT530849827 avatar Nov 18 '24 13:11 CUMT530849827