Con Astronomical Algorithms de J. Meeus (1998) he obtenido la fracción iluminada relevante del disco y el ángulo de posición.
Pero son números (¡todo en griego para mí!). Lo sé, Meeus explica que son proporciones de las longitudes de las extremidades y áreas iluminadas.
Bien, entonces tengo una proporción. ¿Cómo implemento eso junto con el ángulo de posición para calcular la dirección de la extremidad iluminada y la cantidad iluminada?
Supongo que la relación 1.0 está totalmente iluminada (es decir, Luna Llena) y 0.0 no es nada (es decir, Luna Nueva).
y luego está el ángulo de posición, asumo que la dirección en grados comienza en 0 apuntando hacia arriba. ¿Es esa la dirección de la forma de media luna, es decir, los cuernos apuntan en sentido contrario al ángulo direccional?
Me gustaría poder trazar esto para poder mostrar visualmente la fracción iluminada de la luna y la dirección para cualquier día en un año determinado.
Como dije, tengo los números, estoy tratando de traducir eso en un resultado visualmente más atractivo y significativo.
Meeus dice que este es el alargamiento selenocéntrico de la Tierra desde el Sol (ángulo de fase) - i.
la fórmula para la fracción iluminada k es:
k = ( 1 + cos i ) / 2
tan i = (R sen l) / (D - (R cos l))
R = distancia Tierra-Sol D = distancia Tierra-Luna l = elongación geocéntrica
Meeus llama k a la relación de las longitudes BC:AC lo mismo que la relación de las áreas: NBSC:NASC
Donde
Su comprensión de la iluminación es correcta. k es la relación entre la longitud iluminada BC y el diámetro AC. La Luna Nueva tiene una iluminación de 0, la Luna Llena una iluminación de 1 y el cuarto de fase una iluminación de 0,5.
El ángulo de posición se mide en sentido contrario a las agujas del reloj desde el norte (norte celestial, a lo largo de una línea de ascensión recta) hasta el extremo brillante C. Debe calcular el ángulo paraláctico, q, que es el ángulo desde "arriba" hasta el norte celestial. Del capítulo 14 (en la versión 2009)
Calcular H (ángulo horario local) requiere conocer el tiempo sidéreo local θ y la RA de la luna. H = θ - RA. Calcular el LST no es demasiado difícil, pero obtener el RA de la luna es más complicado, como se describe en el capítulo 47 de Algoritmos astronómicos de Meeus. Podría obtenerse más rápidamente simplemente extrayendo el valor de algún lugar usando cURL. Una vez que tenga estos valores, calcular q es sencillo.
Acabo de terminar mi blog detallando 3 formas de obtener la fase lunar vista desde tu ubicación. Las codificaciones de trabajo reales están incluidas. 3 formas de obtener la fase lunar como se ve en su ubicación
mike g
danny f
mike g
danny f
mike g
greg molinero
danny f
greg molinero
greg molinero