pyephem - 计算次太阳点

Question

我刚刚开始使用 PyEphem。我的直接任务是，给定一个日期和时间，用 latitude-longitude 值计算地球上的次太阳点。我将深入研究 PyEphem 来解决这个问题，但如果有人已经这样做了，我将不胜感激示例代码。

Answer 1

我一直在寻找与 OP 相同的答案。许多帖子“提到”PyEphem 是如何走的路，但没有提供实际的例子。

这是我计算太阳下点的工作示例。将所有内容映射到 -180 到 + 180 度之间的经度。

greenwich = ephem.Observer()
greenwich.lat = "0"
greenwich.lon = "0"
greenwich.date = datetime.utcnow()
sun = ephem.Sun(greenwich)
sun.compute(greenwich.date)
sun_lon = math.degrees(sun.ra - greenwich.sidereal_time() )
if sun_lon < -180.0 :
  sun_lon = 360.0 + sun_lon 
elif sun_lon > 180.0 :
  sun_lon = sun_lon - 360.0
sun_lat = math.degrees(sun.dec)
print "Subsolar Point Sun Lon:",sun_lon, "Lat:",sun_lat

我不是 PyEphem 方面的专家，可能有更好的方法 - 但到目前为止，我的测试对我的目的有用。

附注是的.. 上面的格林威治实际上并没有设置为实际的纬度/经度...它实际上只需要 0.0 的经度来获得我们需要的适当的恒星时间。