La última noche astronómica de agosto

 Ayer terminaba mi periodo astronómico del verano, con pena he retirado el equipo astronómico de la azotea, pero la realidad manda.

En la última noche estuve siguiendo al asteroide 12069 que aún andaba cerca de la Luna (la noche anterior no pude dispararle porque tenía la Luna encima).

Tras mucho tratamiento de las imágenes , he optado por presentar las curva tras hacer stacking de dos imágenes (las tomas eran de 3 minutos cada una). He sustraído la mediana a la imagen para minimizar el ruido en la imagen.

La curva obtenida fue la siguiente:

Y lo que pasó a lo largo de la noche queda resumido en la animación siguiente:

Aún con todo este tratamiento, algo de ruido apareció dada la magnitud del asteroide y la cercanía de la Luna. Las condiciones no eran las mejores, pero ¿Quién tiene un cielo perfecto y oscuro?.

Luchando contra los elementos

 La pasada noche quise seguir un asteroide que ronda la magnitud 16 y para ello dejé el equipo disparando a tres minutos. El caso es que la luna estaba cerca pero no tanto (a mi parecer). Durante la noche la luna se fue acercando hasta que saturó todas las tomas. A media madrugada cuando le eché un vistazo a la tablet vi que las tomas estaban saturadas, así que cambié el plan a dos minutos por exposición para dejar de saturar las tomas.

Aún así las tomas eran penosas, pero no quería dejar tiradas en la papelera las tomas de esa noche, así que las que no pude salvar de ninguna manera sí fueron a la papelera. 

Las tomas originales, con tres minutos presentaban este histograma a inicio de la noche.

Pasé a sustraerles la mediana con un radio de 12 píxeles.

El fondo de cielo ya ha bajado a niveles muy bajos, me quedaba corregir los píxeles calientes para que no interfirieran en ninguna medición.

Y finalmente para aumentar la relación Señal-Ruido pasé las imágenes por el apilado (tres imágenes en esta ocasión) móvil de imágenes.

Ese puntito era el asteroide, "sacado" a la luz tras todo este tratamiento. La curva que pude medir sobre las imágenes finales ofreció el siguiente comportamiento.

Ya sé que lo mejor es tener unas imágenes perfectas de partida, pero no siempre se puede conseguir eso y por ello hay que pelearse en ocasiones para sacar la información que permanece oculta tras mucho ruido que se ha acoplado a las tomas.

Recordando gran campo

 Hacía tiempo que no juntaba objetivos fotográficos con CCD y en esta ocasión los he unido a la Star Adventurer. Para comprobar que hasta dos minutos con 135mm son factibles sin desplazamiento apreciable de estrellas en las tomas.

Con un buscador visual para localizar el campo buscado (aunque posiblemente se actualice a unos anillos coordenados en 3D, ya tenía unos hechos de cartón que daban buen resultado).

El campo obtenido era el siguiente:

Con ASI1600 a bin 2 se obtienen casi 8º por casi 6º de campo.

Y con dos minutos se puede llegar a magnitud 10ª en un cielo no muy contaminado (pequeña localidad).

Es cuestión de buscar un campo de observación acorde con estas magnitudes. El problema es el de siempre. ¿Cómo intentar manejar dos equipos simultáneamente?

Seguimos analizando V1209 Her

 En un intento de explicar que pasa en las estrellas pulsantes he creado esta animación sobre datos de V1209 Her puestos en fase.

Que ha salido de la noche del 1 de junio de 2026, donde se pudo observar la siguiente curva.

 Y produciendo el siguiente diagrama de fase.

Mi versión del eclipse de 12 de agosto de 2026 (II)

 Para terminar con el eclipse dejo una animación de las capturas que pude hacer desde la azotea con la canon 6D. Tras discutir con la IA muchas veces, pude conseguir un script que alineara perfectamente las imágenes (me toco hacerlo a mano, porque las propuestas de la IA fallaban y fallaban). 

De las tomas centradas pude medir el porcentaje eclipsado por geometría y obtuve las siguientes gráficas. El programa introduce un error hasta que comenzó el eclipse (no sé el por qué).

Y desde la hoja de cálculo (ya desechados los puntos iniciales erróneos).

Junto a los fallos geométricos habrá que añadir el minuto de adelanto que tiene la cámara respecto a la hora oficial.

Por optimización de Fourier se puede calcular el momento del máximo desde mi localidad:

Que acaba dando el máximo del eclipse a las 20h 41m 9s con un 93.06% de eclipse.

Las previsiones eran que el máximo ocurriría sobre las 20h38.8m con un máximo cercano al 93%. Este es el margen de error cometido con los cálculos geométricos sobre las imágenes capturadas. Quitando el minuto de diferencia que tenía la cámara, nos vamos a 20h40m09s , que está muy cercano a lo previsto (un minuto de diferencia aproximadamente).

Según la previsión...

Mi versión del eclipse de 12 de agosto de 2026

 Mi observación del eclipse consta de unas mediciones realizadas con los sensores DH22 y BMP280 conectados a una placa ESP32 y con el programa correspondiente que mostraba por pantalla led los datos cada 20 segundos. Como me fue imposible automatizar la recogida de datos, opté por apuntarlos manualmente cada 5 minutos.

Las diferentes gráficas muestran las variaciones medidas en el intervalo de las 19h15m a las 21h00 (hora oficial).

La primera gráfica, la temperatura que midió el sensor DH22.

La segunda gráfica es la temperatura medida desde el sensor BMP280.

El máximo del eclipse estaba previsto en mi localización sobre las 20h36 (hora oficial).

La siguiente gráfica pertenece a la presión atmosférica medida por el sensor BMP280.

Y la siguiente gráfica es la altura barométrica medida por el sensor BMP280.

De las gráficas se desprende que la caída de presión se adelanta respecto al momento del máximo eclipse, mientras la temperatura se retrasa muy poco respecto dicho momento.

También medí la humedad con el sensor DH22, pero no parece seguir un patrón claro, salvo un pico cerca del momento del máximo eclipse.

Un asteroide y dos supernovas

 La pasada noche pude observar durante unas horas el asteroide Isolda. Tras la sesión de asteroide me fui a buscar un par de supernovas que están al alcance actualmente.

Del asteroide dejo la curva obtenida de la noche en cuestión.

Y las capturas que pude realizar de las supernovas.

Y recientemente descubierta (el 9 de agosto de 2026)

Algunos datos del 2024

 Hace dos veranos estuve observando una delta Scuti de Bootes. Hasta hace poco tiempo no he encontrado los momentos para analizar parte de las noches observadas en verano.

La estrella V377 Boo era el objeto de estudio ante la posibilidad de saltos bruscos en el periodo. El caso es que usando el algoritmo MPMA pude obtener los siguientes resultados:

Y para el periodo encontrado la gráfica O-C tiene este comportamiento:

Aún me quedan observaciones de 2024 y 2025 por analizar, espero encontrar los ratitos necesarios para terminar el estudio sobre V377 Boo.

En la búsqueda de mejoras

 Siguiendo la búsqueda de mejoras en la fotometría he estado haciendo pruebas con tratamiento de imagen y métodos de fotometría.

Usando los valores de APASS para las estrellas:

11.77 sería 11.921

Var1  sería 14,318

Var2 sería 15,890

Var4 sería 15,907

Para diversos tratamientos obtuve:

Imagen original tras apilar una serie de tomas.

Imagen tras limpiar con la mediana (para las estrellas de calcula la mediana del fondo de cielo).

Imagen tras limpiar con la mediana el fondo de cielo y multiplicar los píxeles de las estrellas por un factor >1.

Se puede ver que los dos primeros tratamientos son muy parecidos, al multiplicar por un factor los píxeles de las estrellas se distorsiona la fotometría.

Posteriormente probé un algoritmo de fotometría PSF con las mismas imágenes, obteniendo los siguientes resultados:

Con la imagen original.

COMP 1 972.276 826.000 137690.82 -12.8473 11.9210 11.9210

TARGET 1 867.740 941.134 4497.50 -9.1324 15.6358 N/A

TARGET 2 1000.366 1002.779 3627.25 -8.8989 15.8693 N/A

Con la imagen de cielo limpiado por mediana.

COMP 1 972.272 825.719 138287.91 -12.8520 11.9210 11.9210

TARGET 1 1074.344 881.000 15208.87 -10.4552 14.3177 N/A

TARGET 2 867.752 941.141 4413.32 -9.1119 15.6610 N/A

TARGET 3 1000.372 1002.809 3515.48 -8.8650 15.9080 N/A

Tras lo cual diría que para estrellas débiles lo mejor puede ser usar fotometría PSF con el fondo lo más limpio posible (usando sustracción de la mediana).

Por otro lado estuve probando en otro campo de estrellas la linealidad de ambos métodos (sustracción de mediana solamente en el cielo y sustracción de mediana a toda la imagen). Pude obtener las siguientes gráficas que indican que la sustracción de mediana solamente al cielo es más lineal que la sustracción de toda la imagen.

Para medir estrellas débiles (supernovas, etc...) puede justificarse que se limpie el cielo usando la mediana y el uso de fotometría basada en PSF en vez de círculos fotométricos.

Adiós estrellas

 Buscando como mejorar fotometría se me ocurrió que podría eliminar los píxeles perteneciente a estrellas en el campo de fondo del cielo. De modo que calculase la mediana de la zona fondo de cielo sin contar con los píxeles de estrellas que en zonas muy densas es imposible de evitar.

Tras unas pruebas se pudo conseguir un algoritmo que realizase ese cálculo y después la fotometría.

De este modo nos aseguramos que realmente estamos midiendo fondo de cielo cuando lo necesitamos para realizar la fotometría. Para hacer las pruebas usé el cúmulo abierto NGC6811, en el Cisne.

Allí tenemos un buen número de estrellas del tipo Delta Scuti y parecía una buena oportunidad para calibrar la efectividad del algoritmo. Las variables señaladas se suponen que son de dicho tipo.

Por su versatilidad y seguridad en las medidas he querido comparar los resultados del algoritmo con los obtenidos por el programa Fotodif, usando los mismos radios fotométricos en ambos programas. Los resultados han sido los siguientes (muestro la noche a bin4, porque la noche a bin1 salieron unas curvas más ruidosas).

Como se puede observar el algoritmo parece fiable en cuanto la obtención de fotometría con estrellas en el fondo del cielo. En nuestro caso las estrellas 1, 3, 4, 11, 12, 14, 15, 16, y 17 tenían estrellas en la zona de cálculo del fondo del cielo.

Mi segunda opción fue pensar que podría modificar el programa Phasthon para no tener que escoger círculo de fondo del cielo anterior o posterior al asteroide. Así que implementé el algoritmo en Phasthon. Y obtuve el siguiente resultado para la noche del 5 de agosto donde el recorrido del asteroide fue el siguiente: 

La gráfica con los puntos rojos corresponde a eliminar estrellas de la zona que se usa de fondo de cielo alrededor del asteroide en cada posición.

Así que podemos dar por bueno el nuevo algoritmo para realizar fotometría de estrellas. Solamente nos queda desarrollar nuevas ideas fotométricas.