AstroPhot: Fitting Everything Everywhere All at Once in Astronomical Images

Abstract

We present AstroPhot, a fast, powerful, and user-friendly Python based astronomical image photometry solver. AstroPhot incorporates automatic differentiation and GPU (or parallel CPU) acceleration, powered by the machine learning library PyTorch. Everything: AstroPhot can fit models for sky, stars, galaxies, PSFs, and more in a principled Chi^2 forward optimization, recovering Bayesian posterior information and covariance of all parameters. Everywhere: AstroPhot can optimize forward models on CPU or GPU; across images that are large, multi-band, multi-epoch, rotated, dithered, and more. All at once: The models are optimized together, thus handling overlapping objects and including the covariance between parameters (including PSF and galaxy parameters). A number of optimization algorithms are available including Levenberg-Marquardt, Gradient descent, and No-U-Turn MCMC sampling. With an object-oriented user interface, AstroPhot makes it easy to quickly extract detailed information from complex astronomical data for individual images or large survey programs. This paper outlines novel features of the AstroPhot code and compares it to other popular astronomical image modeling software. AstroPhot is open-source, fully Python based, and freely accessible here: https://github.com/Autostronomy/AstroPhot

18 Pages, 9 figures, published in MNRAS

References in corpus (16)

• Array Programming with NumPy
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• unWISE: unblurred coadds of the WISE imaging
• Imfit: A Fast, Flexible New Program for Astronomical Image Fitting
• The mass relations between supermassive black holes and their host galaxies at 1<z<2 with HST-WFC3
• lenstronomy II: A gravitational lensing software ecosystem
• The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s
• GALAPAGOS: From Pixels to Parameters
• ProFit: Bayesian Profile Fitting of Galaxy Images
• The Sloan Digital Sky Survey extended point spread functions
• The Automated Photometry Of Transients (AutoPhOT) pipeline
• AutoProf -- I. An automated non-parametric light profile pipeline for modern galaxy surveys
• ProFuse: Physical Multi-Band Structural Decomposition of Galaxies and the Mass-Size-Age Plane
• Posterior samples of source galaxies in strong gravitational lenses with score-based priors
• pysersic: A Python package for determining galaxy structural properties via Bayesian inference, accelerated with jax
• $\mathrm{morphofit}$: An automated galaxy structural parameters fitting package

Cited by in corpus (9)

• Euclid: Early Release Observations -- Dwarf galaxies in the Perseus galaxy cluster
• Euclid: Early Release Observations -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions
• Euclid: Early Release Observations -- The surface brightness and colour profiles of the far outskirts of galaxies in the Perseus cluster
• The Shape of FIREbox Galaxies and a Potential Tension with Low-mass Disks
• Disentangling transients and their host galaxies with Scarlet2: A framework to forward model multi-epoch imaging
• Investigating the Growth of Little Red Dot Descendants at z<4 with the JWST
• Tackling the Problem of Distributional Shifts: Correcting Misspecified, High-Dimensional Data-Driven Priors for Inverse Problems
• Zangetsu: A Candidate of Isolated, Quiescent, and Backsplash Ultra-Diffuse Galaxy in the COSMOS Field
• A review on fundamental bounds and estimators for photometry and astrometry of celestial point sources using array detectors, from first principles