I'm making a program for predicting satellite passes. I'm trying to find out if the satellite is illuminated by the Sun and not in Earth's shadow. I need to know its phase angle: the angle between the observer on Earth, the satellite and the Sun.

Please explain it in simple terms if possible.

(Inserted by reviewer, extension to question, posted originally as answer)

I have TLE data for the satellite (contains right ascension). From that I got ECI position, azimuth, elevation, altitude. For observer I have latitude and longitude and for the sun: elevation and azimuth.


bumped to the homepage by Community 19 hours ago

This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.


Draw the triangle Sun-Earth-Satellite, we will first find the angle Sun-Earth-Sat. The angle between the Sun, the observer, and the satellite will be the angular separation between the Sun and the satellite on the observer's sphere. Here is how you could calculate that:
You have horizontal coordinates for both bodies so the easiest way to go from there would be to look at the spherical triangle Zenith-Sun-Sat, the angle at zenith will be the difference between azimuths, and the Zenith-Sun and Zenith-Sat lengths will be $90^{\circ}-h_{Sun}$ and $90^{\circ}-h_{Sat}$, respectively. Now using the cosine formula for spherical triangles, one may obtain the following formula:
Now use the sine theorem to find the angle E-Sun-Sat (the sine of this angle divided by the sine of the one we calculated will be equal to the ratio of distances from the Earth to the satellite, and from the satellite to the Sun, respectively), and to find the third angle, subtract the two from $180^{\circ}$.
Note: if you do not know the distance from the Sun to the satellite, I am certain you may use the distance from the Earth to the Sun as the error is probably negligible.

  • $\begingroup$ For most satellites, Earth-Sun-Sat is 1 arcminute or less and can be neglected. $\endgroup$ – Mike G 13 hours ago
  • $\begingroup$ You can estimate sinx = x, but I wouldn't say you can completely disregard it, that would change the picture altogether ... $\endgroup$ – Tosic 12 hours ago
  • $\begingroup$ Assuming Earth-Sun-Sat = 0 adds an arcminute of error to Sun-Sat-Earth. Whether you go to the additional trouble depends on the precision you require. $\endgroup$ – Mike G 12 hours ago

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.