3

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

获取对象（在我的情况下是国际空间站）的 next_pass 的方法非常有用......但我真的也希望能够考虑到观察者的实际视觉幅度。这样我就可以只列出国际空间站的可见通行证，而不是全部。

有没有一些方法可以通过考虑观察者、国际空间站和太阳的位置来计算视觉大小？

感谢您的任何建议。

4

3 回答 3

1

编辑：答案找到了我自己对这个问题的答案。请务必阅读帖子的其余部分以获取详细信息。快速总结是我计算了太阳和国际空间站之间的分离角。使用两个已知距离（sun.earth_distance）和iss.range），然后我求解三角形以获得相位角。这被泵入下面的幅度方程......和VIOLA。此返回的值似乎接近 Heavens-above.com 上的值。这就是我想要的。

编辑：有关我用来计算相位角的实际 Python 代码，请参阅这篇文章。

注意：我真的不确定我的编辑是如何处理的。我看到之前的注释，我的编辑将在审核后显示 - 但是我希望仍然会在帖子中看到我的编辑。现在什么都没有了。无论如何，我的编辑实际上更多是对我自己问题的回答。主持人：请不要删除这个答案。它包含对我自己的问题的适当答案，我已经进一步研究并且我认为是相关的。谢谢你！！

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

Heavens-Above 上对国际空间站的预测通常非常可靠，并且考虑到它何时会进入地球阴影。 编辑：事实上，PyEphem 确实显示了卫星是否黯然失色。 身体。黯然失色

此外，虽然众所周知国际空间站偶尔会比预测的更亮，但这仍然是一种非常罕见的情况，我对预测这一点不太感兴趣（我同意尝试这样做似乎毫无意义）。铱星的票价在视觉幅度上当然是非常可变的 - 但仍然可以根据观察者/卫星/高度抛光天线和太阳的角度的位置非常可预测。

如果可能的话，我真的只是在寻找一些指向使用 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 是 -1.3（1000 公里处的固有亮度） - 一些在线资源引用了这个。例如http://satobs.org/seesat/Aug-2005/0114.html和这个Quicksat 内在震级文件：qsmag.zip

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

所以..我不是 100% 确定的一点是如何获得那个角度 B。我接下来会看那个（注意：解决了这个问题 - 见帖子顶部）

于 2013-11-03T22:10:52.870 回答
1

PyEphem 看起来不包含任何卫星震级模型。卫星的亮度不仅取决于其大小和材料的反射程度，还取决于卫星在特定时刻的特定方向——因为同一颗卫星上的不同表面通常具有显着不同的亮度。

而且，事实上，由于明亮的热屏蔽或闪亮的太阳能电池板恰好以特定观察者的角度捕捉太阳，许多卫星在它们从头顶飞过时会周期性地闪烁。但是，当然，几英里外的观察者可能看不到这些闪光，这取决于所涉及的确切角度。如果卫星倾斜几度不同，那么地球上的观察者可能根本看不到闪光。

当然，所有这些都只适用于卫星在其高度处于阳光下的情况。如果卫星远离地球夜间的阴影，那么它将是一个相当黑暗的物体。明暗之间的分界不是绝对的，而是渐进的，因为来自环绕地球的永久环中的环形日出/日落的光开始使卫星从实际进入全日照的位置变亮几度。

而且，不幸的是，PyEphem 没有尝试进行那种 3D 建模，甚至没有尝试对其进行粗略的近似，而这对于准确产生卫星震级是必要的。

于 2013-11-02T14:21:29.837 回答
0

Since B is the sun to the satellite to the observer, wouldn't that be the altitude of the ISS plus the sun's altititude?

or body.alt + abs(sun.alt)

于 2014-01-17T15:00:09.800 回答