Помогите разобраться: точность положения планеты

[drevana]

Всем доброго времени суток
Я хочу получить точное положение планеты с помощью библиотеки Swiss Ephemeris.
Я пишу на nodejs, но я думаю код будет всем понятен.


Для проверки я хочу найти положение Солнца на 19 марта 0:00 UTC


В итоге функция библиотеки мне выдает longitude равным 358.10070052916114
что соответствует 28° 6' 02"


В бумажных эфемеридах на это время 28° 7'06


Я прикинула, что "мое" солнце будет в этой точке через час с хвостом.
Вопрос: что я не так делаю?


Код:

julianDate = jday(2019, 3, 19, 0, 0, 0);
result = swisseph.swe_calc_ut(julianDate, planet, swisseph.SEFLG_GEOCENTRIC | swisseph.SEFLG_SPEED | swisseph.SEFLG_NONUT)

[Единорог]

А почему бумажные эфемериды? Есть астропроги. С ними сравнить можно. Разница 1 минута угловой дуги это менише часа интервал для движения солнца.

[Единорог]

https://sotis-online.ru/?chr=ct:1;dt...me:New%20chart
Правильно посчитано?

[drevana]

Сообщение от Единорог

А почему бумажные эфемериды? Есть астропроги. С ними сравнить можно. Разница 1 минута угловой дуги это менише часа интервал для движения солнца.

Под бумажными я имела в виду pdf файлы Swiss Ephemeris.
Вы хотите сказать, что расчетные данные в pdf у них менее точные чем дает их же программа?
Я думала наоборот, что в файле точное положение, а значение функции зависит от настройки параметров.


Если, например, убрать флаг swisseph.SEFLG_NONUT, то полученное значение будет 28° 5' 47". Это еще дальше от указанного значения в pdf файле.
Да, вы правы, градус это меньше часа. В данном случае 29 мин 35 сек. Получится 28° 7' 0.57"
const julianDate = jday(2019, 3, 19, 0, 29, 35);


Но это все равно по моему много. Быстрые планеты за это время далеко убегут...

[drevana]

Сообщение от Единорог

https://sotis-online.ru/?chr=ct:1;dt:20190319000000;sx:;cid:1591 397;name:New%20chart
Правильно посчитано?

Если я оставлю флаг вычисления позиции без учета поправок на нутацию, то расхождения с Вашим вариантом в секундах. Это по-моему допустимые расхождения.


Солнце 28° 6' 2.52"
Луна 28° 56' 54.53"
Меркурий 20° 31' 31.62"
Венера 20° 37' 2.35"


Но без флага будет не так. Что правильно: с флагом, без флага или то что в pdf я не знаю...


Получается нельзя программно получить данные аналогично pdf файлам?

[DoReMi]

В приложенном коде много чего не видно. Например, какие эфемериды используются.

SEFLG_JPLEPH /* use JPL ephemeris */
SEFLG_SWIEPH /* use Swiss Ephemeris */
SEFLG_MOSEPH /* use Moshier ephemeris */

По умолчанию, если файлы эфемерид не используются или неверно указан к ним путь (не вижу в вашем коде swe_set_ephe_path), то используются эфемериды Мошьера, они "прибрехлительные"

Что у вас за PDF я вообще не представляю :)

[Единорог]

Сообщение от drevana

Под бумажными я имела в виду pdf файлы Swiss Ephemeris.
Вы хотите сказать, что расчетные данные в pdf у них менее точные чем дает их же программа?
Я думала наоборот, что в файле точное положение, а значение функции зависит от настройки параметров.


Если, например, убрать флаг swisseph.SEFLG_NONUT, то полученное значение будет 28° 5' 47". Это еще дальше от указанного значения в pdf файле.
Да, вы правы, градус это меньше часа. В данном случае 29 мин 35 сек. Получится 28° 7' 0.57"
const julianDate = jday(2019, 3, 19, 0, 29, 35);


Но это все равно по моему много. Быстрые планеты за это время далеко убегут...

Сообщение от drevana

Если я оставлю флаг вычисления позиции без учета поправок на нутацию, то расхождения с Вашим вариантом в секундах. Это по-моему допустимые расхождения.


Солнце 28° 6' 2.52"
Луна 28° 56' 54.53"
Меркурий 20° 31' 31.62"
Венера 20° 37' 2.35"


Но без флага будет не так. Что правильно: с флагом, без флага или то что в pdf я не знаю...


Получается нельзя программно получить данные аналогично pdf файлам?

Вам выше рекомендации дал специалист. Я никогда не имел дело со швейцарскими эфемеридами.

[drevana]

Сообщение от DoReMi

В приложенном коде много чего не видно. Например, какие эфемериды используются.

SEFLG_JPLEPH /* use JPL ephemeris */
SEFLG_SWIEPH /* use Swiss Ephemeris */
SEFLG_MOSEPH /* use Moshier ephemeris */

По умолчанию, если файлы эфемерид не используются или неверно указан к ним путь (не вижу в вашем коде swe_set_ephe_path), то используются эфемериды Мошьера, они "прибрехлительные"

Что у вас за PDF я вообще не представляю :)

Вот здесь лежат pdf эфемириды по годам:
https://www.astro.com/swisseph/swepha_e.htm


Здесь для примера эфемериды 2023 года:
https://www.astro.com/swisseph/ae/2000/ae_2023.pdf


По поводу библиотек: я прочитала на этом форуме что нужно указать папку где будут лежать три файла
semo_18.se1, sepl_18.se1, seas_18.se1


const swisseph = require('swisseph');
// Инициализация Swiss Ephemeris
swisseph.swe_set_ephe_path('../ephe'); // Установите путь к файлам эфемерид Swiss Ephemeris


Но указания упомянутых Вами констант у меня в коде нет.

SEFLG_JPLEPH /* use JPL ephemeris */
SEFLG_SWIEPH /* use Swiss Ephemeris */
SEFLG_MOSEPH /* use Moshier ephemeris */


А где их указывать?

[drevana]

Сообщение от DoReMi

В приложенном коде много чего не видно. Например, какие эфемериды используются.

SEFLG_JPLEPH /* use JPL ephemeris */
SEFLG_SWIEPH /* use Swiss Ephemeris */
SEFLG_MOSEPH /* use Moshier ephemeris */

По умолчанию, если файлы эфемерид не используются или неверно указан к ним путь (не вижу в вашем коде swe_set_ephe_path), то используются эфемериды Мошьера, они "прибрехлительные"

Что у вас за PDF я вообще не представляю :)

Ваши вопросы показали мне ошибку!! Спасибо! )))
Я сравнивала сначала 23 год, потом перешла на проверку 19, но при этом пользовалась файлом pdf 2023 года.

Переклинило


В pdf 2019 года 28° 5'47


Это ровно то значение, которое выдает функция БЕЗ параметра SEFLG_NONUT


Еще раз всем спасибо

[DoReMi]

result = swisseph.swe_calc_ut(julianDate, planet, swisseph.SEFLG_GEOCENTRIC | swisseph.SEFLG_SPEED | swisseph.SEFLG_NONUT | swisseph.SEFLG_SWIEPH)

Здесь добавьте.