pyephem - 使用 PyEphem 中的观察点为 Equinox/Solstice 日出/设置方位角

Question

基本上我想做的是:

• 获取观察点(使用纬度/经度)
• 在给定起始日期的情况下计算下一个春分和至日的日期
• 找出每个人的日落方位角
• 找出每个人的日出方位角

*请注意，我在 ArcGIS 中工作，所以我的一些值是从外部表中提取的

这是我所拥有的:

sun = ephem.Sun()
final = ephem.Observer()
final.lon = row[1]
final.lat = row[2]
final.elevation = row[3]

equinoxDate = ephem.next_equinox('0001/01/01')
equinoxSetDate = final.next_setting(ephem.Sun(), start=equinoxDate, use_center=True)
final.date = equinoxSetDate
sun.compute(final)

print sun.az

我一直卡在“下一个设置”部分。我收到 NeverUpError...如果我将其切换为 previous_setting、next_rising、previous_rising...没关系。我总是得到 NeverUpError 或 AlwaysUpError。

如果有人可以帮助我找到春分日落(在任何日期)的方位角，那么我可以确定剩下的。

如果有什么不清楚的地方，请告诉我。

谢谢!

Answer 1

我想通了。我收到错误是因为:final.lon、final.lat 被 pyephem 视为字符串。所以我通过说 final.lon = str(row[1]) 来切换它，然后从那里开始。现在效果很好! Brandon，您的值(value)观走在正确的轨道上。

sun = ephem.Sun()
final = ephem.Observer()
final.lon = str(row[1])
final.lat = str(row[2])
final.elevation = row[3]

equinoxDate = ephem.next_equinox('0001/01/01')
equinoxSetDate = final.next_setting(ephem.Sun(), start=equinoxDate, use_center=True)
final.date = equinoxSetDate
sun.compute(final)

print sun.az