Как посчитать координаты Луны?

[Eddy_Em]

Казалось бы - совершенно разжеванная тема, однако, беру я libnova и пытаюсь с ее помощью определить, где же находится Луна.
Она выдает мне экваториальные координаты, совершенно отличающиеся от координат, которые я вижу в стеллариуме и разнообразных веб-сервисах с эфемеридами (ошибка больше, чем полградуса!).
Ковыряться в недрах стеллариума никакого желания нет: уж очень там много исходников.

Собственно вопрос: какую же взять библиотеку (в SOFA частично вычисляются координаты Луны, но функции, чтобы получить ее геоцентрические экваториальные координаты, нет), чтобы с точностью хотя бы в 5..10'' определить расстояние от избранного объекта на небе до Луны?

Язык - C, система - Gentoo Linux.

[mo]

Язык - C, система - Gentoo Linux.

К сожалению, не подскажу под твою платформу. На винде под ASCOM через https://ascom-standards.org/Help/Developer/html/M_ASCOM_Astrometry_AstroUtils_AstroUtils_EventTimes.htm , вот пример https://ascom-standards.org/Help/Developer/html/1c04e43b-7893-477c-be60-7eb8a14f2b85.htm

Понятно, что это не именно то, что ты ищешь. Но, быть может натолкнёт на нужную библиотеку.

P.S. Плиз, не надо про калькуляторы, игровые приставки и пр. Нет, так нет.

[GraY25]

(ошибка больше, чем полградуса!).

А реальные координаты наблюдателя по времени - точно учитываются?
Геоцентрических координат недостаточно, мы же крутимся по поверхности Земли а не просто смотрим на Луну из центра.

[xd]

Ну так горизонтальный параллакс Луны этот самый градус и даёт!

[Eddy_Em]

А реальные координаты наблюдателя по времени - точно учитываются?

Да, координаты не хуже секунды, время - доли миллисекунд.

Ну так горизонтальный параллакс Луны этот самый градус и даёт!

Ну и как тогда считать и кому верить?
Ведь в стеллариуме-то как-то рассчитывают это все с хорошей точностью...

[xd]

Информация из StarCalc

Код: [Выделить]

»Луна
Экват. координаты (на текущую эпоху):
Геоцентр.: Alp: 5ч 19м 38.1с   Del: +22° 58' 43"
Топоцентр.: Alp: 5ч 21м 39.5с   Del: +22° 19' 2"
r=395839 км;   Угловой размер: 30.19'
Возраст: 20.2 дней.;  Фаза: 0.71
Восход: 21ч 1м
Кульминация: 5ч 32м
Заход: 14ч 50м
08.10.2020 0:30:57
58.00N, 27.55E
GMT+3

[Eddy_Em]

Вот, что сейчас пишет libnova:

Код: [Выделить]

./lunar
JD 2459130.970771
lunar x 17833.446078  y 396102.605968  z 3351.711822
lunar long 87.422154  lat 0.484319
lunar RA 5:48.721563  Dec 23:53.741074
lunar az:  143.952048, lunar z: -13.684787
lunar distance km 396523.279634
lunar disk 0.655078
lunar phase 71.931376
lunar bright limb 88.724931
Стеллариум пишет: 5:47:02, +23:06:44

Даже для наших будущих роботов с поляной в полтора градуса это - критично! Эдак можно и Луну цепануть…

[xd]

Судя по том, что libnova не учитывает положение на Земле, она выдают геоцентрику.
Сделайте 3° от Луны и не ходите на те площадки...

[Eddy_Em]

Сделайте 3° от Луны и не ходите на те площадки...

Хотелось бы, все-таки, иметь более точные данные...
Народ же считает каким-то образом покрытия Луной всяких звезд и планет? Неужто они для этого сторонние сервисы используют?
Да и в том же стеллариуме нормально все считается.

Задача - совершенно избитая. И явно должно быть ее вменяемое решение, а не примерный алгоритм 1967 года…

[gasha]

http://thaec.rascanu.de/ пример программы на C++. Использует библиотеки calceph и iau sofa. Можно использоваать разные эфемериды: epm, inpop, jpl. Обеспечивает максимальную точность.

Код: [Выделить]

./thaecC
Searching for IERS files...
deltat.preds found (0.008808 days old)
UP: 0 of 0  DOWN: 0 of 0
UP: 0 of 0  DOWN: 0 of 0
tai-utc.dat found (0.008576 days old)
UP: 0 of 0  DOWN: 0 of 0
UP: 0 of 0  DOWN: 0 of 0
finals2000A.data found (4.560266 days old)
finals2000A.daily found (0.610162 days old)
eopc04_IAU2000.62-now found (4.706539 days old)
TTBIPM.act found (0.007905 days old)

Searching for Ephemeris files...
inpop13c_TDB_m100_p100_tt.dat found (0.008322 days old)

local:      13.10.2020 10:59:33
UT1:        13.10.2020 07:59:33
TT(=TDT):   13.10.2020 08:00:45.034539  (TT-UT1:  72.034539s)
TT(BIPMxx): 13.10.2020 08:00:45.034567  (TT-UT1:  72.034567s)
TCG:        13.10.2020 08:00:45.997480
TDB:        13.10.2020 08:00:45.032892  (TDB-TT:  -1.674827ms)
TCB:        13.10.2020 08:01:06.455733
JD (TT):    2459135.833855
MJD (UTC):    59135.333021
ERA:        09h 27m 51.940927s
CIP @ GCRS: X=409.367088" Y=0.715809" dX=0.000000" dY=0.000000"
CIO @ GCRS: s=-3.086146mas
CIP @ ITRS: x=0.000000" y=0.000000"
TIO @ ITRS: s'=-0.009768mas

Geocenter BCRS: x=+0.9304498 y=+0.3213587 z=+0.1393991
Observer  BCRS: x=+0.9304260 y=+0.3213721 z=+0.1394317

Sun: (8.297237 light-minutes)
BCRS: RA: 09.070365 dec: +17.341363, x=-0.0060884 y=+0.0058682 z=+0.0026405
CIRS: RA: 13.241198 dec: -07.986488
obs:  RA: 13.239649 dec: -07.946887, Az: 129.10525 Alt: +18.622129

Moon: (1.212628 light-seconds)
BCRS: RA: 01.276646 dec: +08.111339, x=+0.9282903 y=+0.3223512 z=+0.1400525
CIRS: RA: 10.358721 dec: +14.698996
obs:  RA: 10.358651 dec: +14.710114, Az: 172.12908 Alt: +54.386306


В исходнике можно пути поменять на ftp://cddis.gsfc.nasa.gov/pub/products/iers/, чтобы скачивались свежие данные. А то по приведённым там ссылкам глухо всё.

[Eddy_Em]

@gasha , спасибо!
Буду ковырять, как руки дойдут до этого.