Das hatte ich ich auch vermutet - eine Formel zu dem Thema würde mir sehr weiterhelfen...
Neben Inklination i und RAAN (right ascension of ascending node) brauchst du noch die Position von Sonne gegenüber Erde (also den Vektor Erde --> Sonne). Wenn du es genau haben willst musst du auch berücksichtigen dass JPL Ephemeriden und TLE nicht im gleichen Frame sind. JPL DE405 ist ICRF und DE200 ist in J2000, beide sind aber nahezu identisch. TLE ist auch fast ICRF, nur die Epoche ist nicht J2000.0 sondern im TLE Datensatz vermerkt und meist ein aktuelles Datum. Dadurch stimmen Präzession und Nutation der Erde zwischen den Frames nicht überein. Wenn du diese Abweichung vernachlässigst bzw. alles im gleichen (Erd-Äquator)Frame vorhanden ist, geht die Berechnung von Beta so:
Positionsvektor Sonne pVektor von Erde zur Sonne aus Ephemeriden
Orbitnormalenvektor nx = sin (i) * sin (RAAN)
y = - sin (i) * cos (RAAN)
z = cos (i)
Betabeta = 90° - arccos(p*n/|p|/|n|)
Beta ist 0° wenn die Sonne genau in der Orbitebene liegt. Steht die Sonne genau senkrecht auf dem Orbit ist Beta 90° (solch ein Orbit liegt immer in der Sonne, also kein Erdschatten).
Ich hoffe ich habe keinen Fehler gemacht.
Cosmo