ELLIPSECT: A surface brightness analysis tool for GALFIT 3

El artículo presenta EllipSect, una herramienta en Python que complementa GALFIT 3 mediante la generación de gráficos de alta calidad, la exportación de datos y el cálculo de mediciones no paramétricas para evaluar modelos de brillo superficial.

Lo esencial

¡Claro que sí! Imagina que las galaxias son como gigantescas tartas de cumpleaños hechas de estrellas, gas y polvo. Los astrónomos quieren saber exactamente de qué están hechas estas tartas: ¿cuánto hay de "bizcocho" (el núcleo brillante)? ¿Cuánto de "crema" (el disco exterior)? ¿Y qué tan grande es la tarta en total?

Para hacer esto, usan un programa muy famoso y potente llamado GALFIT 3. Piensa en GALFIT como un chef matemático muy experto que toma una foto de la tarta y trata de adivinar la receta exacta, dibujando un modelo teórico sobre la foto para ver si coincide.

El problema es que GALFIT es como un chef que te da una lista de ingredientes y medidas en un cuaderno de notas muy técnico, pero no te muestra el pastel terminado para que tú lo veas y digas: "¡Ah, sí, eso se parece a la foto!". Además, GALFIT a veces se pierde cuando hay muchas tartas pegadas unas con otras (como en un grupo de galaxias) y no sabe bien dónde termina una y empieza la otra.

Aquí es donde entra en escena el héroe de esta historia: EllipSect.

¿Qué es EllipSect?

Imagina que EllipSect es el socio creativo y el inspector de calidad del chef GALFIT. Es una herramienta nueva (escrita en un lenguaje de programación llamado Python) que toma los resultados del chef y hace tres cosas mágicas:

1. El "Espejo Mágico" (Visualización):
   GALFIT te da los números, pero EllipSect te muestra las imágenes. Toma la foto original de la galaxia, el modelo que dibujó GALFIT y la "diferencia" entre ambas (lo que se le escapó al modelo).

   • La analogía: Es como si el chef te diera la receta, y tú, con EllipSect, pudieras ver una foto del pastel real, una foto del pastel ideal que el chef imaginó, y una foto de los trozos que sobró o faltó para que fueran idénticos. Esto ayuda a ver si el modelo es bueno o si hay que ajustarlo.

2. El "Cuchillo de Medición" (Mediciones no paramétricas):
   GALFIT mide partes específicas (como el radio del núcleo), pero a veces queremos saber cosas más generales que la receta no nos da directamente. EllipSect actúa como un cuchillo de cocina muy preciso que corta la galaxia en rebanadas (sectores) para medir cosas como:

   • El radio efectivo total: Si la galaxia tiene varias partes (un núcleo, un disco, una barra), EllipSect suma todo para decirte: "Oye, la mitad de toda la luz de esta galaxia está dentro de este círculo".
   • El "Radio de Petrosian" y "Radio de Kron": Son como reglas especiales para medir hasta dónde llega realmente la galaxia antes de que se desvanezca en la oscuridad del espacio.
   • La "Bumpiness" (Irregularidad): Imagina que pasas la mano por la superficie de la galaxia. ¿Es suave como la crema de una tarta perfecta, o tiene baches y protuberancias? EllipSect mide esos baches para decirte si la galaxia es tranquila o si ha tenido una pelea (colisiones con otras galaxias).

3. El "Detective de Fondo" (Cálculo del cielo):
   Para medir bien una galaxia, necesitas saber qué tan brillante es el fondo (el "cielo" negro detrás de la tarta). A veces, el chef GALFIT se equivoca al medir este fondo. EllipSect tiene dos trucos de detective:

   • Método de Gradiente: Mira anillos concéntricos alrededor de la galaxia (como los anillos de un árbol) para ver cómo cambia la luz hasta encontrar el punto donde la galaxia termina y el espacio vacío comienza.
   • Método de Cajas Aleatorias: Lanza "cajas" imaginarias en lugares aleatorios alrededor de la galaxia para tomar una muestra de la oscuridad del fondo y calcularla con precisión.

¿Por qué es importante esto?

Antes, si querías ver si el modelo de GALFIT era bueno, tenías que usar programas viejos y complicados que a veces se rompían o no funcionaban bien con galaxias que están muy juntas.

Con EllipSect, los astrónomos pueden:

• Ver rápidamente si su modelo es correcto (¡y generar gráficos hermosos para sus artículos!).
• Medir galaxias complejas que tienen varias partes (núcleo, disco, barra) sin confundirse.
• Automatizar el proceso: Imagina tener que analizar 100 tartas. EllipSect puede hacerlo todo por ti en segundos, dándote una lista de medidas precisas para cada una.

En resumen

EllipSect es el traductor y el artista que toma los datos fríos y matemáticos de GALFIT y los convierte en imágenes claras, gráficos bonitos y medidas útiles. Es la herramienta que permite a los astrónomos no solo "calcular" las galaxias, sino realmente "ver" y entender su estructura, como si finalmente pudieran tocar y medir la tarta estelar que están estudiando.

Es una herramienta amigable, potente y esencial para quien quiera entender la forma y la luz de las galaxias en nuestro universo.

Resumen técnico

Aquí tienes un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Resumen Técnico: EllipSect, una Herramienta de Análisis de Brillo Superficial para GALFIT 3

1. El Problema
El código GALFIT 3 es el estándar de la industria para el ajuste bidimensional (2D) de la distribución de brillo superficial (SB) de fuentes astronómicas, permitiendo modelar componentes como bulos, discos y barras simultáneamente. Sin embargo, GALFIT 3 presenta limitaciones significativas en el análisis posterior a la adaptación (post-fit):

• Falta de herramientas gráficas integradas: No incluye una herramienta nativa para comparar visualmente la imagen de entrada con el modelo resultante y los residuos de manera inmediata.
• Dependencia de flujos de trabajo obsoletos o complejos: Los usuarios tradicionalmente han dependido de tareas de 1D como ellipse de IRAF (que dejó de ser soportado en 2013) o requieren convertir formatos de salida de GALFIT a formatos legibles por otros códigos, un proceso propenso a errores y consumo de tiempo.
• Limitaciones en mediciones no paramétricas: GALFIT 3 no calcula automáticamente parámetros críticos derivados de modelos multicomponente, como el radio efectivo total ($r_e^t$) de una galaxia compuesta, ni métricas de calidad de ajuste no paramétricas (como el parámetro de "Bumpiness" o el radio de Petrosian).
• Dificultad en campos densos: En cúmulos o grupos compactos, el enmascaramiento manual es insuficiente; el ajuste simultáneo es necesario, pero el análisis de los componentes individuales y sus interacciones es difícil de visualizar y cuantificar con las herramientas actuales.

2. Metodología
EllipSect es una herramienta de análisis e implementación en Python diseñada para operar directamente sobre los archivos de salida de GALFIT 3 (específicamente el archivo galfit.nn y el cubo de datos FITS). Su metodología se basa en los siguientes pilares:

• Procesamiento de Datos: Lee los archivos de salida de GALFIT 3, incluyendo la imagen original, el modelo ajustado y la imagen de residuos. Utiliza la misma máscara que GALFIT para eliminar píxeles no deseados y resta el componente de cielo.
• Algoritmo de Seccionamiento (Sectors): En lugar de ajustar elipses isofóticas individuales (como hace IRAF), EllipSect divide la galaxia y su modelo en sectores angulares equidistantes alrededor del centro, utilizando la rutina sectors_photometry de la librería mgefit. Calcula el brillo superficial promedio en múltiples radios dentro de estos sectores.
• Parámetros Fijos vs. Variables: A diferencia de los métodos 1D que permiten que la elipticidad y el ángulo de posición varíen con el radio, EllipSect mantiene estos parámetros fijos (tomados del último componente del modelo GALFIT, típicamente exponencial). Esto es una ventaja en campos densos donde las galaxias están mezcladas en la línea de visión, evitando la necesidad de un "doble ajuste" (2D y 1D).
• Cálculo de Fotometría y Parámetros No Paramétricos:
  • Radio Efectivo Total ($r_e^t$): Calcula numéricamente el radio que contiene el 50% de la luz total sumando los modelos de todos los componentes (bulo, disco, barra, etc.), superando la aproximación de usar solo el $r_e$ de un componente individual.
  • Correcciones: Descarga automáticamente datos del NED para aplicar correcciones de extinción galáctica y atenuación de brillo superficial (SB dimming).
  • Cálculo del Cielo: Ofrece dos métodos propios para estimar el fondo de cielo (no incluidos en GALFIT): el método de gradiente (anillos concéntricos) y el método de cajas aleatorias, permitiendo al usuario refinar el modelo si la estimación inicial es incorrecta.
  • Métricas de Calidad: Calcula el Criterio de Información de Akaike (AIC) y el Criterio de Información Bayesiano (BIC) para comparar modelos, además de parámetros como "Bumpiness" (irregularidad estructural) y "Tidal parameter" (perturbaciones de marea).

3. Contribuciones Clave

• Automatización y Usabilidad: Proporciona una interfaz de línea de comandos y un módulo de Python fácil de usar que elimina la necesidad de conversión manual de formatos, generando gráficos de calidad para publicaciones y productos de datos exportables.
• Análisis de Componentes Individuales: Permite visualizar y analizar por separado cada componente del modelo multicomponente de GALFIT, facilitando la identificación de componentes irrelevantes o redundantes.
• Nuevas Mediciones Fotométricas: Introduce el cálculo sistemático de radios no estándar (Radio de Kron, Radio de Petrosian, radio al 90% de la luz) y parámetros de morfología (Bumpiness, Tidal) que GALFIT 3 no proporciona nativamente.
• Integración con Flujos de Trabajo Masivos: Diseñado para ser compatible con códigos de automatización como DGCG (Driver for GALFIT on Cluster Galaxies), permitiendo el análisis de cientos de galaxias en cúmulos densos.

4. Resultados y Ejemplos
El artículo presenta la aplicación de EllipSect en tres casos de estudio con diferentes instrumentos y complejidades:

• Holm 15A (KPNO): Una galaxia cD masiva en el cúmulo Abell 85. El modelo utilizó siete componentes gaussianos. EllipSect generó perfiles de brillo superficial detallados y gráficos de residuos, demostrando su capacidad para manejar perfiles complejos.
• M61 (2MASS): Una galaxia espiral en el banda Ks. Se utilizó un modelo con múltiples componentes Sérsic y exponenciales. Los gráficos multi-ángulo mostraron desviaciones donde el modelo no reproducía perfectamente el brillo, ilustrando la utilidad de la herramienta para identificar fallos en el ajuste.
• NGC 4261 (HST): Una galaxia con un núcleo activo, ajustada con un solo componente Sérsic. Los resultados mostraron la capacidad de la herramienta para generar visualizaciones claras de galaxia, modelo y residuo con alta resolución.

En todos los casos, la herramienta generó exitosamente gráficos de perfil de brillo superficial (SB), imágenes de cubo (galaxia, modelo, residuo) y archivos de datos con parámetros fotométricos calculados.

5. Significado e Impacto
EllipSect representa una evolución crucial en el análisis de morfología galáctica al cerrar la brecha entre el ajuste paramétrico complejo (GALFIT 3) y la interpretación física de los resultados.

• Precisión en Entornos Densos: Su enfoque de mantener parámetros geométricos fijos basados en el modelo 2D es superior a los métodos 1D tradicionales para galaxias en cúmulos, donde el enmascaramiento es difícil.
• Validación de Modelos: Al proporcionar métricas como $r_e^t$ y parámetros de calidad (AIC/BIC) específicos para el modelo completo, permite a los astrónomos validar si un modelo multicomponente es físicamente significativo o si es un sobreajuste.
• Reproducibilidad y Eficiencia: Al estar escrito en Python y ser de código abierto, facilita la integración en pipelines de análisis modernos, reduciendo el tiempo de procesamiento para grandes sondeos y mejorando la reproducibilidad de los estudios de estructura galáctica.

En conclusión, EllipSect no es solo un visualizador, sino una herramienta de análisis cuantitativo que extiende las capacidades de GALFIT 3, permitiendo una caracterización más robusta y detallada de la estructura y evolución de las galaxias.