Сравнение астрометрических библиотек - 2
Jan. 9th, 2020 08:28 am![[personal profile]](https://www.dreamwidth.org/img/silk/identity/user.png){kind=small}
eddy_emКак я уже писал, получается, что разные библиотеки для вычисления apparent place дают разные данные. С расхождением в минуты!
Один из авторов любезно поделился со мной статьей, где таких расхождений нет: там говорится о 5-6mas, а не минутах!!!
Берем все тот же comparelibNOVA-NOVAS-SOFA-ERFA. Вводим туда координаты БТА, сегодняшнюю поправку DUT1 и смотрим, что поменяется, если в первом случае атмосферное давление сделать равным 100кПа, а во втором — 78кПа (обычное давление на горе). Давление как на уровне моря:
Давление как на БТА:
Запустил однострочником, чтобы как можно меньше было расхождение по JD.
В принципе, судя по коду, сравнение делается некорректно. В libnova вызывается функция ln_get_equ_prec2, никакого отношения к apparent place не имеющая, да и julian date не вычисляется нормально. В NOVAS julian date вычисляется, исходя из поправок для TT/TAI, не учитывая DUT1, но все равно идет просто пересчет координат на новую эпоху, без apparent place.
В общем, сравнивать их надо корректно. Но, судя по тому, что NOVA не использует DUT1, ее результаты таки должны быть менее правильными. Хотя, понятное дело: даже для БТА ошибка вычисления координат в ±1'' — мелочь по сравнению с ошибками наведения и сопровождения.
Один из авторов любезно поделился со мной статьей, где таких расхождений нет: там говорится о 5-6mas, а не минутах!!!
Берем все тот же comparelibNOVA-NOVAS-SOFA-ERFA. Вводим туда координаты БТА, сегодняшнюю поправку DUT1 и смотрим, что поменяется, если в первом случае атмосферное давление сделать равным 100кПа, а во втором — 78кПа (обычное давление на горе). Давление как на уровне моря:
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.228458943796,28.258597189533 871.794, 389.549 -0.009, 0.039 324.298, -266.516
NOVAS precession 346.228461429350,28.258586387865 871.803, 389.510 324.306, -266.555
SOFA Atci13 346.221443558472,28.258637867687 846.538, 389.695 -25.264, 0.185 299.042, -266.370
SOFA Atco13 346.138376291446,28.332629488577 547.496, 656.065 -324.306, 266.555
SOFA Atco13 NO r 346.221423118462,28.258608573980 846.465, 389.590 -25.338, 0.080 298.969, -266.475
SOFA Atco13 NO eo 345.886031202046,28.332629488577 -360.946, 656.065 -1232.749, 266.555 -908.442, 0.000
ERFA Atco13 346.138376291446,28.332629488577 547.496, 656.065 -324.306, 266.555 0.000, 0.000
Давление как на БТА:
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.228458943945,28.258597189600 871.794, 389.549 -0.009, 0.039 258.673, -208.067
NOVAS precession 346.228461429500,28.258586387931 871.803, 389.510 258.682, -208.106
SOFA Atci13 346.221443558383,28.258637867562 846.538, 389.695 -25.264, 0.185 233.418, -207.921
SOFA Atco13 346.156605302540,28.316393568600 613.120, 597.616 -258.682, 208.106
SOFA Atco13 NO r 346.221423120371,28.258608573613 846.465, 389.590 -25.338, 0.080 233.344, -208.026
SOFA Atco13 NO eo 345.904260212808,28.316393568600 -295.322, 597.616 -1167.124, 208.106 -908.442, 0.000
ERFA Atco13 346.156605302540,28.316393568600 613.120, 597.616 -258.682, 208.106 0.000, 0.000
Запустил однострочником, чтобы как можно меньше было расхождение по JD.
В принципе, судя по коду, сравнение делается некорректно. В libnova вызывается функция ln_get_equ_prec2, никакого отношения к apparent place не имеющая, да и julian date не вычисляется нормально. В NOVAS julian date вычисляется, исходя из поправок для TT/TAI, не учитывая DUT1, но все равно идет просто пересчет координат на новую эпоху, без apparent place.
В общем, сравнивать их надо корректно. Но, судя по тому, что NOVA не использует DUT1, ее результаты таки должны быть менее правильными. Хотя, понятное дело: даже для БТА ошибка вычисления координат в ±1'' — мелочь по сравнению с ошибками наведения и сопровождения.
