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

Categories:

Сравнение астрометрических библиотек - 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 уже не поддерживается). Но вот с экваториальными все не так однозначно. Надо проводить следственный эксперимент.
Tags: c, астрофизика
Subscribe

  • Новая железяка - продолжение

    Шел третий день отпуска… Когда собирал железяку, возник насущный вопрос: а сможет ли такой слабенький движочек гонять такую дурищу (прежде всего —…

  • Новая железяка для БТА/Ц1000

    Уже больше полугода занимаюсь разработкой, вот, наконец-то в мастерских взялись за меня и начали выдавать первые детали. Сегодня сделал тестовую…

  • Пара зарядок с али

    3 августа заказал пару новых зарядок (нонче такое дело, что зарядок много не бывает). Сегодня получил. Обозреваю. Первая зарядка. Аж три USB и типа…

promo eddy_em september 3, 12:13 8
Buy for 10 tokens
Уже больше полугода занимаюсь разработкой, вот, наконец-то в мастерских взялись за меня и начали выдавать первые детали. Сегодня сделал тестовую сборку (как обычно, местами пришлось "доработать напильником"): Пока прибор без названия (да и как-то не лезет в голову ничего, у меня нет…
  • Post a new comment

    Error

    Anonymous comments are disabled in this journal

    default userpic
  • 0 comments