Obtenga una tabla de ángulo de fase lunar por día

JPL HORIZONTES puede hacer eso. Con ajustes como estos:

Tipo de efemérides: OBSERVADOR
Cuerpo objetivo: Luna [Luna] [301]
Ubicación del observador: Geocéntrico [500]
Intervalo de tiempo: Inicio = 2019-09-12 , Fin = 2019-09-19 , Paso = 1 d
Configuración de la tabla: CANTIDADES = 10 ,23,24

devuelve resultados como estos:

 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

donde Illu%esta la fraccion iluminada$k$, y S-T-Oes el ángulo de fase del observador del objetivo solar$i$. Estas cantidades están relacionadas por ( Meeus 41.1)

Una ubicación del observador en la superficie de la Tierra en lugar del geocentro da resultados ligeramente diferentes debido a la paralaje.

Las imágenes renderizadas de Tanner parecen estar indexadas por la diferencia entre la longitud de la eclíptica lunar y la solar . El ángulo de elongación del sol-observador-objetivo S-O-Tes una aproximación cercana en la mayoría de los casos. Cuando la /rcolumna sea /Tpara el seguimiento, use ese valor directamente. Donde /res /Lpara guía, reste S-O-Tde 360°.

Un solo ciclo de imágenes de la Luna no puede capturar los efectos de libración como lo hacen las imágenes SVS de la NASA (ver la respuesta de SpaceBread ), por lo que las imágenes de Tanner se ven un poco diferentes de la Luna real.

El Estudio de visualización científica de la NASA tiene una página que muestra la fase lunar y la libración en el transcurso de cada año. En la misma página puede encontrar un enlace a la tabla de datos que incluye, entre otras cosas, la fase de la luna para cada hora del año en UT. Puede descargar la tabla como un archivo TXT y JSON .
En caso de que solo desee una lista diaria y esté utilizando Linux, podría hacer algo como:

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

Finalmente solucioné mi problema al seguir usando las imágenes de Tanner (nombradas del 0 al 360) y al obtener las fases de la biblioteca SunCalc (que devuelve un número entre 0 y 1). De esta forma, podía asociar fácilmente imágenes de una fase a una fecha determinada.

Método 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";
}