Hola Jordi
Pues es un problema bastante difícil. MorphologicalTransform no te va a ayudar demasiado. La única solución que se me ocurre es trabajar en dos fases. Ante todo tienes que hacer una buena máscara de estrellas y actuar con ella activa, o destrozarás la imagen.
Yo no he usado máscara de estrellas en este ejemplo por brevedad, pero por supuesto que deberás hacerlo.
Como puedes ver, lo que he hecho es extraer una PSF a partir de una de las estrellas problemáticas. Esto es muy fácil: define un preview sobre la estrella y arrastra el view selector (la pestaña de la izquierda con el nombre en vertical) fuera de la ventana, con lo cual generarás una nueva ventana con el contenido del preview. Después lo pasas a escala de grises y le ajustas los histogramas
sin recortar ni un sólo píxel. Es muy importante que el preview lo definas
centrado sobre la estrella con precisión.
Después seleccionas esta PSF para aplicar una deconvolución. Muy pocas iteraciones, porque de lo contrario freirás las estrellas. Fíjate que esto corregirá todas las estrellas que tengan aproximadamente ese tamaño, pero casi no tocará a las estrellas mayores o menores. O sea que tendrás que repetir la operación con otra(s) PSF(s) adaptadas a diferentes tamaños de estrellas, posiblemente con otras tantas máscaras de estrellas.
La deconvolución te dejará un perfil estelar totalmente plano, que deberás corregir posteriormente, por ejemplo con el script del que estáis hablando, o con wavelets.
Edito: se me olvidaba decir que no deberías aplicar la deconvolución a la crominancia, sino sólo a la luminancia. Si deconvolucionas la crominancia perderás toda la información de color, que después no podrás recuperar.
Bueno, esta es la idea. No es fácil y es muy laborioso. Tú tendrás que valorar si te compensa el trabajo ;-)
