Implementación de algoritmos de deconvolución en el software de reducción de datos astrofísicos Gnuastro

Abstract

La deconvolución es una operación matemática que permite recuperar imágenes que han sido distorsionadas por la respuesta óptica del instrumento de observación en presencia de ruido. En este trabajo, se han implementado diferentes algoritmos de deconvolución para el tratamiento de imágenes astronómicas. El objetivo final es compartir dichos algoritmos de manera libre a través de su incorporación en la librería open source Gnuastro. Concretamente, se han implementado los algoritmos de inversión directa, Wiener-Tikhonov, Richardson-Lucy y AWMLE. Además, se han desarrollado librerías de bajo nivel que permiten trabajar con matrices de números complejos y realizar transformaciones de Fourier y wavelet, las cuales también pueden ser aprovechadas en el futuro por cualquier desarrollador o investigador. Para validar los algoritmos se han utilizado tanto imágenes académicas con ruido acotado como imágenes astrofísicas más realistas y complejas que simulan diferentes escenarios astrofísicos. Por último, para asegurar la calidad del código y su correcta integración en Gnuastro, el trabajo ha sido supervisado por un administrador principal del proyecto de manera continua.