Получить таблицу фазового угла Луны за сутки

JPL HORIZONS может это сделать. С такими настройками:

Тип эфемерид: НАБЛЮДАТЕЛЬ
Целевое тело: Луна [Луна] [301]
Местоположение наблюдателя: геоцентрическое [500]
Промежуток времени: начало = 12.09.2019 , окончание = 19.09.2019 , шаг = 1
Параметры таблицы d: КОЛИЧЕСТВА = 10 ,23,24

он возвращает такие результаты:

 Date__(UT)__HR:MN       Illu%    S-O-T /r    S-T-O
***************************************************
 2019-Sep-12 00:00      95.688 155.9731 /T  23.9640
 2019-Sep-13 00:00      98.620 166.4750 /T  13.4888
 2019-Sep-14 00:00      99.814 175.0444 /T   4.9422
 2019-Sep-15 00:00      99.231 169.9061 /L  10.0670
 2019-Sep-16 00:00      96.874 159.5743 /L  20.3723
 2019-Sep-17 00:00      92.791 148.7668 /L  31.1541
 2019-Sep-18 00:00      87.074 137.7510 /L  42.1469
 2019-Sep-19 00:00      79.858 126.5419 /L  53.3370

где Illu%освещенная доля$k$, а S-T-O- фазовый угол солнце-цель-наблюдатель$i$. Эти величины связаны соотношением ( Meeus 41.1 )

Расположение наблюдателя на поверхности Земли вместо геоцентра дает немного другие результаты из-за параллакса.

Визуализированные изображения Таннера, по-видимому, индексируются по разнице между лунной и солнечной эклиптической долготой . Угол элонгации солнце-наблюдатель-цель S-O-Tв большинстве случаев является близким приближением. Если /rстолбец предназначен /Tдля трейлинга, используйте это значение напрямую. Где /rдля /Lинтерлиньяжа вычесть S-O-Tиз 360°.

Один цикл изображений Луны не может зафиксировать эффекты либрации , как это делают изображения НАСА SVS (см. ответ SpaceBread ), поэтому изображения Таннера выглядят немного иначе, чем реальная Луна.

В студии научной визуализации НАСА есть страница , показывающая фазу Луны и либрацию в течение каждого года. На этой же странице вы можете найти ссылку на таблицу данных, включающую, среди прочего, фазы луны для каждого часа года по всемирному времени. Вы можете скачать таблицу в виде файла TXT и JSON .
Если вам нужен только ежедневный список и вы используете Linux, вы можете сделать что-то вроде:

cat mooninfo_2019.txt | grep "00:00 UT"  > mooninfo_2019_daily.txt

Я, наконец, решил свою проблему, по-прежнему используя изображения Таннера (с именами от 0 до 360) и получая фазы из библиотеки SunCalc (которая возвращает число от 0 до 1). Таким образом, я мог легко связать изображения фазы с заданной датой.

Метод Javascript:

getMoonImageName (year, month, day) {
  var ill = SunCalc.getMoonIllumination(new Date(year, month-1, day))
  return Math.round(ill['phase'] * 360).toString().padStart(3, "0") + ".png";
}