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.

Adiós darks

 Siguiendo con la creación de algoritmos, he estado desarrollando un algoritmo que "mata" a los píxeles calientes de una imagen, de modo que si nos olvidamos de los darks o nos salen regulares, podemos acudir al algoritmo a limpiar de píxeles calientes nuestras imágenes.

Os dejo una visión aclaratoria de dos imágenes, la de la derecha es la imagen original con ruido generado por un tratamiento anterior de la toma, la imagen de la izquierda ha sido procesada por el algoritmo en cuestión (Pixelhot, lo llamo yo). Los niveles de brillo se han igualado en las dos imágenes para no introducir variaciones entre las imágenes mostradas.

Mi idea es usarlo para corregir tomas donde un píxel caliente se mete en los círculos de medición fotométrica y arruina unas mediciones sin querer. Se trabaja con dos parámetros que el usuario puede modificar a su gusto, de hecho se puede apreciar como en las estrellas brillantes introduce un error, que para nuestros propósitos es despreciable.

En fin, que seguimos creando scripts que nos permitan exprimir nuestras imágenes astronómicas la máximo. Seguiremos pensando nuevas utilidades para el astrónomo amateur.