Compte tenu des éléments suivants:
- Temps, t
- L'ensemble des données d'éphémérides IS-200, E, d'un satellite GPS correspondant au temps t
- La position ECEF du satellite GPS, P = (x, y, z), dérivée du temps et des éphémérides, (t, E).
- Supposons que la terre n'est que l'ellipsoïde WGS-84.
- Tous les points du WGS-84 ont l'angle de masque, m.
Trouvez les éléments suivants:
- l'anneau de couverture, R, sur WGS-84 du satellite GPS. c'est-à-dire la frontière qui distingue quels points WGS-84 sont en vue le satellite au point P = (x, y, z) et quels points WGS-84 ne sont pas en vue
Solutions acceptables:
- Une spline sur WGS-84 qui se rapproche de R.
- Un polygone sur WGS-84 qui se rapproche de R.
- Ou une formule (s) qui me donne R.
Ce que j'ai essayé jusqu'à présent:
- Soit e ^ 2 = 0,0066943799901264; excentricité au carré
Nous avons une position ECEF WGS-84 par latitude géodésique phi et longitude lambda:
r = 1 / (sqrt (1-e ^ 2 sin ^ 2 (phi))) * (cos (phi) * cos (lambda), cos (phi) * sin (lambda), (1-e ^ 2) * péché (phi))
Je convertis ensuite ECEF en cadre géographique est-nord en haut (ENU) avec phi et lambda en utilisant la matrice:
(-sin(lambda) cos(lambda) 0 )
C= (-cos(lambda)*sin(phi) -sin(lambda)*sin(phi) cos(phi))
( cos(lambda)*cos(phi) sin(lambda)*cos(phi) sin(phi))
- Soit G = C (P - r)
- Prenez la composante z de G. si la composante z de G est supérieure à sin (m), alors je sais que le point, r, est en vue. Mais cela ne suffit pas pour obtenir la solution que je recherche. Je pourrais juste trouver un tas de points qui sont en vue et prendre la coque convexe de ces points, mais ce n'est pas du tout efficace.