This IDL library is designed to be used on astronomical images. Its main aim is to stack data to allow a statistical detection of faint signal, using a prior. For instance, you can stack 160um data using the positions of galaxies detected at 24um or 3.6um, or use WMAP sources to stack Planck data. It can estimate error bars using bootstrap, and it can perform photometry (aperture photometry, or PSF fitting, or other that you can plug). The IAS Stacking Library works with gnomonic projections (RA---TAN), and also with HEALPIX projection.
https://ui.adsabs.harvard.edu/abs/2010A%26A...516A..43B and Bavouzet, 2008, PhD Thesis Univ. Paris Sud 11 https://tel.archives-ouvertes.fr/tel-00363975/fr/