python - 古代或 future 观察者的进动、自行和纪元

Question

我有一个关于纪元的问题。我想展示超古代观察者或遥远 future 的某个人看到的夜空。我是否应该将年份 和 纪元设置为相同？

Brandon Rhodes Halley Comet 示例中概述的三个场景如下。

1. my_star.compute(epoch='-100000')
2. my_star.compute(-100000',epoch='-100000')
3. my_star.compute('-100000')

对我来说，根据他对赤道望远镜的描述和第二个例子对1066年当代观测者的描述，我认为我应该选择第二个选项。也就是说，将纪元和年份设置为古代或 future 的年份。

halley.compute(epoch='1066')
This is probably useless: it computes the current position of halley, but returns coordinates relative to the direction the earth’s axis was pointing in the year 1066. Unless you use a Conquest-era star atlas, this is not useful.

halley.compute('1066', epoch='1066')
This is slightly more promising: it computes the position of halley in 1066 and returns coordinates for the orientation of the earth in that year. This might help you visualize how the object was positioned above contemporary observers, who considered it an ill omen in the imminent conflict between King Harold of England and William the Bastard. But to plot this position against a background of stars, you would first have to recompute each star’s position in 1066 coordinates.

halley.compute('1066')
This is what you will probably use most often; you get the position of halley in the year 1066 but expressed in the 2000 coordinates that your star atlas probably uses.

为什么我对我得到的描述如此不确定是因为 Ian Ridpath 和 Rick Pogge 得到了不同的结果。这些相关结果的链接发布在 github 问题中。 Rick 说他也得到了 StarryNight 的结果的证实。我不介意我是否错了，我只是看到我认为最好的哈雷例子与少数人的结果相冲突。

我还不能在 stackoverflow 上发布很多链接(或图像!)，但我在我的 github 问题中对此进行了非常详尽的记录 https://github.com/brandon-rhodes/pyephem/issues/61 .我只是想非常清楚我在寻找什么，并向您展示可能的情况。

如果您想测试两个工作的 IPython 笔记本，请参阅我的 notebooks 中的第 1 部分:proper-motion.ipynb 和第 2 部分:big-dipper.ipynb github 上的文件夹。

big-dipper.ipynb 示例中有很多代码，因此我不会将其发布到跟踪器中，但我会很快就此写一篇博文。请查看上面的 github 问题链接以获取所有代码和图像!

我相信不同的时代和在一起的描述是我想要的关于一位古代观察者在他的视角中看到夜空，以及一位 future 派观察者在他的视角中看到夜空的描述自己视角下的夜空。有时这很难理解，所以我非常希望社区能就此主题发表意见。

这是代码的较小部分。 big-dipper.ipynb 笔记本显示了更详尽的场景。

%pylab inline
import ephem

UMa = {
'dubhe': 'αUMa', #HIP 54061
'merak': 'βUMa', #HIP 53910
'phecda':'γUMa', #HIP 58001
'megrez':'δUMa', #HIP 59774
'alioth':'εUMa', #HIP 62956
'mizar': 'ζUMa', #HIP 65378
'alcor': '80UMa',#HIP 65477
'alcaid':'ηUMa', #HIP 67301
}

def const(const,year=None,epoch='2000',title=None):
    s = []

    for star in const.keys():
        s.append(ephem.star(star.capitalize()))
    
    for i in range(len(s)):
        if(year!=None):
            s[i].compute(year,epoch=epoch)
        else:
            s[i].compute(epoch=epoch)
    
    tau = 2.0 * pi
    degree = tau / 360.0
    hour = tau / 24.0

    ra_list = [star.a_ra / hour for star in s]
    dec_list = [star.a_dec / degree for star in s]

    mag_array = np.array([star.mag for star in s])
    size_array = (5 - mag_array) ** 1.5 * 4

    scatter(ra_list,dec_list,size_array)
    if(title!=None):
        pyplot.title(title)
    gca().xaxis.grid(True)
    gca().yaxis.grid(True)
    gca().invert_xaxis()
    return s

Answer 1

如果我想模仿地球上的古代观察者或 future 派观察者，我们应该同时改变纪元和时间。

已在 the issue tracker 中回答.

谢谢你的帮助:)