Top 10 Best Galaxies Software of 2026

DS9 is a desktop viewer focused on FITS display, WCS transformations, and region tools that can drive measurements like catalog cross-checks and pixel statistics. It supports browsing through image planes and spectral axes so that the same region can be inspected across channels. The integration story centers on FITS interoperability and scriptable operations that can be used in repeatable pipelines.

A key tradeoff is that DS9 is a visualization and analysis GUI rather than an end-to-end galaxy modeling system. It fits best when investigators need fast, WCS-correct visual QA of produced products like moment maps, lightcone slices, or extracted spectra before downstream statistical modeling.

CASA centers workflows around measurement-set data structures and the CASA task system, which makes calibration-to-imaging sequencing explicit. It includes interactive and scripted execution paths, which helps teams reproduce reductions by versioning task parameters in scripts. The imaging toolkit covers multi-frequency synthesis and deconvolution workflows used for radio continuum and spectral cubes.

A key tradeoff is that CASA is specialized for radio data models, so teams using optical-style tabular pipelines or non-radio formats often need preprocessing to convert into measurement sets. CASA fits best when datasets are already in radio interferometric form and the goal is calibrated imaging, spectral-line extraction, or derived cube products.

pynbody focuses on post-processing for galaxy formation simulations by letting users load snapshots, select particle subsets, and compute derived fields like centers of mass, velocities, and thermodynamic or kinematic quantities depending on the dataset. It supports workflow automation through Python scripts and reusable functions that operate on snapshot objects. The toolkit also includes routines that help build higher-level diagnostics such as radial profiles and phase-space style selections without forcing data export to external analysis stacks.

A key tradeoff is that pynbody analysis code assumes Python-based array semantics, so scaling to very large snapshots often depends on careful field selection and chunking patterns rather than automatic distributed execution. It fits best when iterative analysis is the bottleneck, such as refining subhalo centering or generating consistent mock observables across multiple runs inside a single Python environment.

CIGALE is a galaxies software solution focused on fitting galaxy SEDs using configurable star-formation and dust prescriptions. The core capability is producing synthetic spectral energy distributions from stellar population inputs and applying attenuation and emission components to match observed photometry and spectra.

CIGALE also supports workflow automation via parameter-grid runs, which helps standardize large model sweeps across many targets. Its configuration-driven approach makes it suitable for reproducible research pipelines that need consistent assumptions.

Astropy provides a Python framework for working with astronomical data, including FITS file handling and common coordinate transforms. It supplies a documented core API plus extensibility hooks so projects can build domain-specific modules and integrate them into analysis pipelines.

Core components include time and coordinate utilities, units-aware calculations, and tables for heterogeneous catalog data. Astropy also supports interoperability with the broader astronomy Python ecosystem through standardized data structures and I/O conventions.

Source Extractor suits astronomers who need repeatable catalog generation from large image sets, especially in command-line research workflows. Its distinct focus is detecting sources and measuring positions, fluxes, shapes, and morphology directly from astronomical images rather than modeling galaxy populations.

Configurable background subtraction, deblending, aperture photometry, and star-galaxy classification support survey preprocessing. Source Extractor writes catalogs in formats including FITS file format and supports dual-image mode for consistent detection and measurement across bands.

yt from yt-project.org focuses on programmatic analysis of large astrophysical simulation outputs using yt’s Python-native workflow. It converts simulation datasets into analysis objects that support derived fields, unit-aware calculations, and common visualization outputs like slices and projections.

The tool emphasizes reproducible scripts for tasks such as halo-oriented measurements, mock observations, and catalog-style exports rather than point-and-click analysis. yt also supports automation through its Python API, which enables batch processing of many snapshots with consistent field definitions and thresholds.

Photutils is a Python library for astronomical image analysis that focuses on measurable tasks like source detection, aperture photometry, and background estimation. It provides cohesive tools built around common data structures such as Astropy images and tables, which supports consistent preprocessing and downstream measurements.

The documented API covers segmentation, deblending, and fitting utilities that convert pixel data into catalogs of positions and fluxes. Photutils is distinct in its tight scope on measurement workflows rather than end-to-end survey pipelines.

Lenstronomy performs gravitational lensing and lens-galaxy forward modeling from parametric mass profiles to predicted images and observables. The library is built around a configurable model graph that connects lens mass, light profiles, PSF convolution, and likelihood evaluation in a single modeling loop.

It supports automated fitting via optimization and sampling workflows exposed through a Python API. Documentation centers on reproducible examples that map directly to common lensing use cases like quasar lenses and lens-galaxy light decomposition.

BAGPIPES is a documented Python framework for fitting galaxy spectral energy distributions with configurable stellar population models. It provides end-to-end SED fitting that covers photometry and spectra inputs with selectable priors and parameter constraints.

BAGPIPES runs locally from notebooks or scripts and exposes fit outputs suitable for downstream cataloging and analysis. It is designed for reproducible model fitting where the configuration files define the model, bands, and optimization workflow.