Расчет положения восходящих и нисходящих узлов Луны из файлов эфемерид

Вот способ сделать это полностью онлайн, если вы не хотите загружать библиотеки CSPICE. Посетите https://naif.jpl.nasa.gov/naif/webgeocalc.html -> щелкните WebGeocalc и щелкните поиск позиции:

Заполните данные следующим образом:

Обратите внимание, что планетоцентрическая широта, равная 0, означает, что эклиптическая широта равна 0, поскольку мы используем систему отсчета ECLIPDATE (эклиптика даты).

Нажмите рассчитать:

Результаты — это все случаи пересечения Луной эклиптики между 1970 и 2040 годами. Сюда входят как восходящие, так и нисходящие узлы.

Конечно, вам нужны не времена, а фактические позиции, поэтому нажмите «Сохранить все интервалы»:

Вернитесь в меню расчета и нажмите «вектор состояния»:

Заполните его следующим образом:

Обратите внимание, что вам может потребоваться перетащить список сохраненных интервалов в поле «Список интервалов», иначе он может быть заполнен автоматически.

После заполнения нажмите «Рассчитать»:

И вот оно: список положений Луны, когда Луна находится в нисходящем и восходящем узлах.

Обратите внимание, что прямое восхождение узла движется несколько линейно:

но не совсем.

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

Ответ с использованием SPICE превосходен!

Самый короткий путь — использовать интерфейс JPL Horizons в режиме ORBITAL ELEMENTS, загрузить таблицу и использовать OMили$\Omega$, долгота восходящего узла в градусах.

  Symbol meaning:

    JDTDB    Julian Day Number, Barycentric Dynamical Time
      EC     Eccentricity, e
      QR     Periapsis distance, q (km)
      IN     Inclination w.r.t X-Y plane, i (degrees)
      OM     Longitude of Ascending Node, OMEGA, (degrees)
      W      Argument of Perifocus, w (degrees)
      Tp     Time of periapsis (Julian Day Number)
      N      Mean motion, n (degrees/sec)
      MA     Mean anomaly, M (degrees)
      TA     True anomaly, nu (degrees)
      A      Semi-major axis, a (km)
      AD     Apoapsis distance (km)
      PR     Sidereal orbit period (sec)

#!/usr/bin/env python
# -*- coding: utf-8 -*-

import numpy as np
import matplotlib.pyplot as plt

with open('datamoon.txt', 'r') as infile:
    linez = infile.readlines()

lines = [line.split(',') for line in linez]

Omega = np.array([float(line[5]) for line in lines])

JD = np.array([float(line[0]) for line in lines])
JD -= JD[0]

plt.plot(JD/365.2564 + 10, Omega)
plt.xlim(10, 30)
plt.xlabel('years since 2000')
plt.ylabel('geocentric lunar ascending node Ω (deg)')
plt.show()

Индийские родные альманахи указывают положение узлов. Восходящий узел — Рааху, нисходящий — Кету. https://srirangaminfo.com/vakya-panchangam-srirangam.php