December 30th, 2019

Костерок

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

Решил я навсегда избавиться от устаревшей 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.
promo eddy_em august 17, 12:33 3
Buy for 10 tokens
Юра намедни напечатал корпус для хронометра. Для первого блина получилось неплохо: И еще немного фотографий: