计算次太阳点

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

我去寻找与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 的经度即可获得我们需要的适当的恒星时间。