简而言之，我计算了太阳和国际空间站之间的分离角。使用两个已知距离（sun.earth_distance）和 iss.range），然后我求解三角形以获得相位角。这被输入到下面的幅度方程......和 VIOLA。返回的值似乎与 Heavens-above.com 上的值接近。这就是我想要的。
See 这篇文章提供了实际的Python代码 https://stackoverflow.com/questions/19759501/calculating-the-phase-angle-between-the-sun-iss-and-an-observer-on-the-earth我用来计算相位角。
许多提供卫星预测的在线服务实际上也提供震级预测。 Heavens-Above.com 就是这样一个网站，因此肯定有一些技术可以做到这一点。
《天堂之上》对国际空间站的预测通常相当可靠，并且考虑到了它何时会进入地球的阴影。编辑：事实上 PyEphem 确实显示卫星是否被食。 body.eclipsed
此外，虽然众所周知国际空间站偶尔会比预测的更亮，但这仍然是一种非常罕见的情况，而且我对预测这种情况不太感兴趣（而且我同意尝试这样做似乎毫无意义）。当然，铱星的费用在视觉幅度上变化很大，但仍然可以根据观察者/卫星的位置/高度抛光的天线与太阳的角度来预测。
如果可能的话，我实际上只是在寻找一些使用 PyEphem 库计算此值的方法的指针。
我在这里找到了关于这一点的讨论->http://www.satobs.org/seesat/Apr-2001/0313.html http://www.satobs.org/seesat/Apr-2001/0313.html
>Mag = Std. Mag - 15 + 5*LOG(Range) -
> 2.5*LOG(SIN(B) + (pi-B)*COS(B))
>
>where Range is in km, and B is in radians and measures
>the angle from the sun to the satellite to the observer.
>At full phase, B is 0; at new phase, B is pi (i.e.
>satellite transiting the sun).
enter code here
我在国际空间站使用的 Standard.Mag 是 -1.3（1000 公里处的固有亮度）——一些在线资源引用了这一点。例如http://satobs.org/seesat/Aug-2005/0114.html http://satobs.org/seesat/Aug-2005/0114.html和这个Quicksat 固有震级文件：qsmag.zip http://prismnet.com/%7Emmccants/programs/index.html
我感觉 PyEphem 有能力在给定起始震级 (std.mag) 的情况下进行此计算。互联网上有“标准”震级数字（不确定它们是如何计算的......但它们仍然可用）。
所以..我不是100%确定的是如何获得那个角度B。接下来我会看看这个（注意：解决了这个问题 - 请参阅帖子顶部）
