Емельянов Эдуард Владимирович (eddy_em) wrote,
Емельянов Эдуард Владимирович
eddy_em

Categories:

Сравнение астрометрических библиотек

Решил я навсегда избавиться от устаревшей slalib (например, в apogee_control и прочих вещах, где нужен перевод координат между эпохами Jnow и J2000).
Существует несколько современных библиотек для этих целей: NOVAS, NOVA, SOFA и ERFA (правда, я так понял, что это — та же самая SOFA, но под другой лицензией).
Погуглив "nova sofa comparation" наткнулся на статью, но скачать ее нельзя — плати денежку... ОК, помимо статьи наткнулся на исходники. Пример есть и у автора в Readme на гитхабе, но я собрал и тоже запустил:
libNOVA UTC JNow  2458847.786754336208
NOVAS   TT  JNow  2458847.787555077113
                                                          d J2000             d NOVAS          d SOFA Atco13
                  ra deg           ,dec deg          ra''      ,dec''    ra''      ,dec''    ra''      ,dec''
J2000             345.986294056000,28.150389198100
libNOVA equ_prec2 346.228128882082,28.258449683874   870.605,  389.018    -0.009,    0.039   295.127, -160.093
NOVAS precession  346.228131342123,28.258438903511   870.614,  388.979                       295.136, -160.131
SOFA Atci13       346.222029679428,28.259039007764   848.648,  391.139   -21.966,    2.160   273.170, -157.971
SOFA Atco13       346.146149141993,28.302919840940   575.478,  549.110  -295.136,  160.131
SOFA Atco13 NO r  346.221947751140,28.259014432233   848.353,  391.051   -22.261,    2.072   272.875, -158.059
SOFA Atco13 NO eo 345.894173073050,28.302919840940  -331.636,  549.110 -1202.250,  160.131  -907.114,    0.000
ERFA Atco13       346.146149141993,28.302919840940   575.478,  549.110  -295.136,  160.131     0.000,    0.000

И кому верить?..

P.S. Почему-то поломался экспорт из dreamwidth в ЖЖ. Приходится руками копировать.
Что за чушь такая?
Failed to crosspost entry to eddy_em@LiveJournal: Client error: Your IP address is temporarily banned for exceeding the login failure rate.
Tags: c, астрофизика, всячина
Subscribe

  • Релюшки на CAN-шине

    Закончил с прошивкой для новой железяки. Как "наследница" USB-CAN переходника, она умеет все то же самое + несколько специфичных вещей (опрос…

  • Богатство нашего времени

    Вскрыл сегодня вот такой пакетик (лежит на работе уже довольно-таки давно, жрать не просит): Содержимое: Распаял модуль управления…

  • Продолжаю возиться с STM32F303

    Добавил работу с USART'ами: простейший вариант "почти эха" USART1 (чтение с двойной буферизацией в прерывании, блокирующая запись) и работу с тремя…

promo eddy_em august 17, 2019 12:33 3
Buy for 10 tokens
Юра намедни напечатал корпус для хронометра. Для первого блина получилось неплохо: И еще немного фотографий:
  • Post a new comment

    Error

    Anonymous comments are disabled in this journal

    default userpic
  • 0 comments