Прецессия, собственное движение и эпохи для древних или будущих наблюдателей

Question

У меня вопрос по эпохам. Я хочу показать ночное небо, как его увидит супер-древний наблюдатель или кто-то в далеком будущем. Должен ли я установить год и эпоху одинаковыми?

Три сценария, изложенные в примере Кометы Брэндона Родса Галлея, заключаются в следующем.

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

Для меня, основываясь на его описаниях экваториальных телескопов и описании второго наблюдателя современного наблюдателя в 1066 году, я считаю, что мне следует выбрать второй вариант. То есть установить как эпоху, так и год в древние или будущие годы.

halley.compute (эпоха ='+1066')
Это, вероятно, бесполезно: он вычисляет текущее положение галеры, но возвращает координаты относительно направления, на которое указала земная ось в 1066 году. Если вы не используете звездный атлас эпохи Завоевания, это бесполезно.

halley.compute ('1066', эпоха ='1066')
Это немного более многообещающе: он вычисляет положение галеры в 1066 году и возвращает координаты ориентации Земли в этом году. Это может помочь вам представить, как объект был расположен над современными наблюдателями, которые считали его дурным предзнаменованием в неизбежном конфликте между королем Англии Гарольдом и Уильямом-ублюдком. Но чтобы построить это положение на фоне звезд, сначала нужно пересчитать положение каждой звезды в 1066 координатах.

halley.compute ('1066')
Это то, что вы, вероятно, будете использовать чаще всего; Вы получаете положение Галлея в 1066 году, но выражается в координатах 2000 года, которые, вероятно, использует ваш звездный атлас.

Почему я так не уверен в изображениях, которые я получил, потому что Йен Ридпат и Рик Погге получили разные результаты. Ссылки на эти релевантные результаты размещены в выпуске github. Рик говорит, что он также подтверждается результатами StarryNight. Я не возражаю, если я ошибаюсь, я просто вижу пример Галлея, который, как мне кажется, отборный, сталкивается с результатами нескольких человек.

Я пока не могу опубликовать много ссылок на stackru (или изображения!), Но я очень тщательно документировал эту проблему в github на https://github.com/brandon-rhodes/pyephem/issues/61. Я просто хотел быть очень ясным в том, что я ищу, и показать вам возможные сценарии.

Если вы хотите протестировать две рабочие тетради IPython, см. Часть 1: proper-motion.ipynb и часть2: big-dipper.ipynb в моей папке блокнотов на github.

Там много кода в big-dipper.ipynb Например, я не буду публиковать его в трекере, но скоро сделаю запись в блоге. Пожалуйста, смотрите ссылку на выпуск github выше для всего кода и изображений!

Я полагаю, что описание " Эпоха перемен" и " Год вместе" - это то, что я хочу, чтобы древний наблюдатель видел ночное небо в своей перспективе, а футуристический наблюдатель видел ночное небо в своей собственной перспективе. Иногда это довольно сложно понять, поэтому мне бы очень хотелось, чтобы сообщество внесло свой вклад по этому вопросу.

Вот меньшая часть кода. 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

Мы должны менять эпоху и время вместе, если я хочу подражать древнему или футуристическому наблюдателю на земле.

Ответил в номере трекера.

Спасибо вам за помощь:)