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

arXiv comment18 Pages, 9 figures, published in MNRAS

References in corpus (37)

• 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
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Overview of the DESI Legacy Imaging Surveys
• Detailed Decomposition of Galaxy Images. II. Beyond Axisymmetric Models
• The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
• The DEEP Groth Strip Survey II. Hubble Space Telescope Structural Parameters of Galaxies in the Groth Strip
• The Absolute Magnitude of the Sun in Several Filters
• unWISE: unblurred coadds of the WISE imaging
• Imfit: A Fast, Flexible New Program for Astronomical Image Fitting
• Ultra Low Surface Brightness Imaging with the Dragonfly Telephoto Array
• MegaMorph - multi-wavelength measurement of galaxy structure: complete Sérsic profile information from modern surveys
• Galaxy Zoo DECaLS: Detailed Visual Morphology Measurements from Volunteers and Deep Learning for 314,000 Galaxies
• 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
• ProFound: Source Extraction and Application to Modern Survey Data
• The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s
• GALAPAGOS: From Pixels to Parameters
• Improvements to the Levenberg-Marquardt algorithm for nonlinear least-squares minimization
• ProFit: Bayesian Profile Fitting of Galaxy Images
• The Sloan Digital Sky Survey extended point spread functions
• Summarizing CPU and GPU Design Trends with Product Data
• PyMorph: Automated Galaxy Structural Parameter Estimation using Python
• Beyond Ellipse(s): Accurately Modelling the Isophotal Structure of Galaxies with $Isofit$ and $Cmodel$
• 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
• PSF--NET: A Non-parametric Point Spread Function Model for Ground Based Optical Telescopes
• DiskFit: a code to fit simple non-axisymmetric galaxy models either to photometric images or to kinematic maps
• Testing PSF Interpolation In Weak Lensing With Real Data
• Euclid: Nonparametric point spread function field recovery through interpolation on a graph Laplacian
• Posterior samples of source galaxies in strong gravitational lenses with score-based priors
• Pix2Prof: fast extraction of sequential information from galaxy imagery via a deep natural language 'captioning' model
• The Euclid Data Processing Challenges
• The Geometry of Hamiltonian Monte Carlo
• 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 (10)

• 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
• 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
• The Shape of FIREbox Galaxies and a Potential Tension with Low-mass Disks
• 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
• Separating the Optical to Near-Infrared Light of AGN and Their Host Galaxies