Поиск

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


Итак, пока сам не сделаешь анализ, ничего не прояснится…
Набросал файлик, сравнивающий sla, SOFA и NOVA.

Запускаем с дефолтными координатами:

ICRS (catalog), epoch J2000.0: 19 50 47.6000000 +08 52 12.300000
Date: (d/m/y +frac) 21/1/2020 +0.33831696
MJD by slalib: 58869.33831696
UTC by sofa: 58869, 0.33831696
date by sofa (TAI/TT): 58869/58869 & 0.338745/0.339118 ICRS -> CIRS: 19 50 42.5788574 +08 55 13.974206 CIRS -> ICRSc: 19 50 47.6000000 +08 52 12.300000 geocentric apparent: 19 51 43.2687797 +08 55 13.974206 geocentric apparent (sla): 19 51 43.2652993 +08 55 13.974165
apparent -> astrometric (sla): 19 50 47.6000000 +08 52 12.300000 ICRS -> observed: 19 50 41.8665206 +08 55 48.049071
A(bta)/Z: +155 39 35.29(+24 20 24.71)/+36 57 35.60 observed -> ICRS: 19 50 47.6000000 +08 52 12.299999 ICRS -> observed (sla): 19 51 42.5550336 +08 55 47.915150
A(bta)/Z: +155 39 35.65(+24 20 24.35)/+36 57 35.58 geocentric apparent (NOVA): 19 51 43.2678257 +08 55 13.981070

Apparent place (AP) все три вычислили с точностью до единиц mas. Однако, AP — не самое интересное. Ведь нужно навести телескоп на данное место. И тут начинается «веселуха». В NOVA нет однозначной процедуры, позволяющей вычислить observed place (OP), поэтому сравниваем лишь sla и SOFA. И видим, что sla на 1 минуту привирает по RA (или это я неправильно что-то ввел? но аргументы одни и те же использовались). Однако, на перевод в альт-азимутальную систему это почему-то особо не повлияло, и погрешность получилась на уровне долей секунд.
ОК, теперь возьмем приполярную зону: склонение 80, восхождение 11 и 23:

ICRS (catalog), epoch J2000.0: 11 00 00.0000000 +80 00 00.000000
Date: (d/m/y +frac) 21/1/2020 +0.34296686
MJD by slalib: 58869.34296686
UTC by sofa: 58869, 0.34296686
date by sofa (TAI/TT): 58869/58869 & 0.343395/0.343768 ICRS -> CIRS: 11 00 42.7476701 +79 53 21.761216 CIRS -> ICRSc: 11 00 00.0000000 +80 00 00.000000 geocentric apparent: 11 01 43.4376637 +79 53 21.761216 geocentric apparent (sla): 11 01 43.4342702 +79 53 21.760878
apparent -> astrometric (sla): 11 00 00.0000000 +80 00 00.000000 ICRS -> observed: 11 01 01.0254362 +79 53 57.912808
A(bta)/Z: -11 10 33. (-168 49 27. )/+51 54 15.42 observed -> ICRS: 10 59 59.9999962 +79 59 59.999992 ICRS -> observed (sla): 11 02 01.7014038 +79 53 58.201038
A(bta)/Z: -11 10 32.64(-168 49 27.36)/+51 54 15.39 geocentric apparent (NOVA): 11 01 43.4481354 +79 53 21.708150

./slalib_sofa_nova 23:0:0 80:0:0
ICRS (catalog), epoch J2000.0: 23 00 00.0000000 +80 00 00.000000
Date: (d/m/y +frac) 21/1/2020 +0.34438367
MJD by slalib: 58869.34438367
UTC by sofa: 58869, 0.34438367
date by sofa (TAI/TT): 58869/58869 & 0.344812/0.345184 ICRS -> CIRS: 22 59 16.2944361 +80 06 36.434407 CIRS -> ICRSc: 23 00 00.0000000 +80 00 00.000000 geocentric apparent: 23 00 16.9844514 +80 06 36.434407 geocentric apparent (sla): 23 00 16.9808922 +80 06 36.434817
apparent -> astrometric (sla): 23 00 00.0000000 +80 00 00.000000 ICRS -> observed: 22 59 00.9355098 +80 06 21.002559
A(bta)/Z: +12 46 30. 7(+167 13 29.93)/+41 59 17.70 observed -> ICRS: 23 00 00.0000011 +80 00 00.000001 ICRS -> observed (sla): 23 00 01.6365350 +80 06 20.725939
A(bta)/Z: +12 46 30.43(+167 13 29.57)/+41 59 17.68 geocentric apparent (NOVA): 23 00 16.9771567 +80 06 36.491739

Опять расхождение в 1 минуту по RA! Хотя AZ вычислены вполне одинаково.
И на экваторе примерно в кульминации расхождение по RA в 1 минуту (т.е. 15′!!!):

ICRS (catalog), epoch J2000.0: 19 00 00.0000000 +00 00 00.000000
Date: (d/m/y +frac) 21/1/2020 +0.34468266
MJD by slalib: 58869.34468266
UTC by sofa: 58869, 0.34468266
date by sofa (TAI/TT): 58869/58869 & 0.345111/0.345483 ICRS -> CIRS: 18 59 58.6672840 +00 01 40.392000 CIRS -> ICRSc: 19 00 00.0000000 -00 00 00.000000 geocentric apparent: 19 00 59.3573038 +00 01 40.392000 geocentric apparent (sla): 19 00 59.3538217 +00 01 40.391816
apparent -> astrometric (sla): 19 00 00.0000000 +00 00 00.000000 ICRS -> observed: 18 59 58.7074781 +00 02 25.529542
A(bta)/Z: -179 20 51.94(-00 39 08. 6)/+43 36 53.83 observed -> ICRS: 19 00 00.0000000 -00 00 00.000004 ICRS -> observed (sla): 19 00 59.3940403 +00 02 25.465919
A(bta)/Z: -179 20 51.58(-00 39 08.42)/+43 36 53.81 geocentric apparent (NOVA): 19 00 59.3574990 +00 01 40.387584

При том, что AZ-координаты опять более-менее одинаковы!
Сверим еще со Stellarium (у которого что-то свое используется, т.к. он ни от одной астрометрической библиотеки не зависит — не знаю, правда, зачем так авторы делали). Берем звезду 19:03:32.3 +1:49:6.2:

 ICRS -> observed: 19 03 30.2325669 +01 51 34.687907
A(bta)/Z: -178 19 25. 8(-01 40 34.92)/+41 48 20.60 ICRS -> observed (sla): 19 04 30.9197685 +01 51 34.625932
A(bta)/Z: -178 19 24.72(-01 40 35.28)/+41 48 20.58

Стеллариум показал A=-01:39:0, h=48:11:06. Т.е. с точки зрения вычисления OP в AZ-системе все прекрасно! Но откуда же берется 1 минута разницы во времени???

Вывод: для работы с AZ-монтировками однозначно можно хоть SOFA, хоть sla использовать (лучше SOFA, т.к. sla уже не поддерживается). Но вот с экваториальными все не так однозначно. Надо проводить следственный эксперимент. eddy_em.livejournal.com

Добавить комментарий