有没有办法计算卫星(ISS)的视觉幅度?

时间:2013-11-02 07:24:26

标签: pyephem

刚刚在我的Raspberry Pi上运行了PyEphem的基础知识。工作得很好。

获取对象的next_pass(在我的情况下是ISS)的方法非常有用......但我真的也希望能够考虑观察者的实际视觉量。这样我只能列出国际空间站的可见通行证,而不是所有通行证。

是否有一些方法可以通过考虑观察者,国际空间站和太阳的位置来计算视觉量?

感谢您的任何建议。

3 个答案:

答案 0 :(得分:0)

看起来PyEphem不包含任何卫星幅度模型。卫星的亮度不仅取决于其大小和材料的反射程度,而且取决于卫星在特定时刻的特定方向 - 因为同一卫星上的不同表面通常具有明显不同的亮度。

事实上,许多卫星因为明亮的热屏蔽或发光的太阳能电池板在特定观察者的角度下捕捉太阳而定期闪烁。但是,当然,几英里外的观察者可能看不到那些闪光,这取决于所涉及的确切角度。如果卫星倾斜几度不同,那么地球观测者可能看不到闪光。

当然,所有这一切只适用于卫星在其海拔高度的阳光下;如果卫星远离地球夜晚的阴影,那么它将是一个非常黑暗的物体。明亮和黑暗之间的截止不是绝对的,而是渐进的,因为来自环形日出/日落的光线在地球周围的永久环中开始使卫星从实际进入充足阳光的位置亮几度。

而且,遗憾的是,PyEphem并没有尝试进行3D建模,甚至是粗略的近似,这对于准确生成卫星幅度是必要的。

答案 1 :(得分:0)

编辑:ANSWER 找到我自己对这个问题的回答。一定要阅读帖子的其余部分以获取详细信息。快速总结是我计算太阳和iss之间的分离角度。使用两个已知距离(sun.earth_distance)和iss.range)然后我解决三角形以获得相位角。这会被输入下面的幅度等式......和VIOLA。这返回的值似乎接近于Heavens-above.com上的值。这就是我想要的。

编辑:参见this post for the actual Python code我用来计算相位角。

注意:我真的不确定如何处理我的编辑。我看到之前的注释,我的编辑将在审核后显示 - 但我希望仍然可以看到我在帖子中显示的编辑内容。现在什么都没有了。在任何情况下,我的编辑实际上都是对我自己的问题的回答。主持人:请不要删除这个答案。它包含了我自己的问题的适当答案,我已经进一步研究,我认为是相关的。谢谢!!

许多提供卫星预测的在线服务实际上也提供了幅度预测。 Heavens-Above.com是一个这样的网站,所以肯定有技术可以做到这一点。

在天堂 - 上方可用的国际空间站的预测通常是非常可靠的,并考虑它何时将传递到地球的阴影。 编辑:事实上PyEphem确实显示卫星是否黯然失色。 body。黯然失色

虽然众所周知ISS偶尔闪现比预期更亮但仍然非常罕见,而且我对预测它并不太感兴趣(我同意尝试这样做似乎毫无意义)。铱价在视觉幅度上当然是变化很大的 - 但是根据观察者/卫星/其高度抛光天线和太阳的角度的位置,它们是非常可预测的。

如果可能的话,我真的只是想看一些使用PyEphem库来计算它的方法。

再次感谢任何实际执行此操作的指示。

编辑:我在这里找到了关于这一点的讨论 - > 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我用于ISS的是-1.3(1000km的固有亮度) - 一些在线资料引用了这一点。例如http://satobs.org/seesat/Aug-2005/0114.html和此Quicksat intrinsic magnitudes file: qsmag.zip

我有一种感觉,PyEphem有能力在给定起始幅度(std.mag)数的情况下进行此计算。互联网上有“标准”幅度数字(不确定它们是如何计算的......但是它们都是可用的)。

所以..我不是100%肯定的是如何获得那个角度B.我会看下一个(注意:做到这一点 - 见帖子顶部)

答案 2 :(得分:0)

由于B是卫星到观察者的太阳,不是国际空间站的高度加上太阳的高度吗?

或body.alt + abs(sun.alt)