pith. sign in

arxiv: 2606.20360 · v1 · pith:4JKKJO5Anew · submitted 2026-06-18 · 🌌 astro-ph.IM

Lightstack: A Python Package for Creating Photometric Data Cubes

Pith reviewed 2026-06-26 15:31 UTC · model grok-4.3

classification 🌌 astro-ph.IM
keywords lightstackphotometric data cubesmulti-band photometryPython packagePSF matchingimage stackingastronomy softwaredata preparation
0
0 comments X

The pith

Lightstack supplies a Python package whose three-step workflow turns separate multi-filter images into aligned photometric data cubes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents lightstack as a tool that takes standalone telescope images across different filters and converts them into data cubes ready for multi-band photometry. It does so through a fixed sequence of cropping a region of interest, stacking the cropped frames, and matching their point-spread functions. A sympathetic reader would care because upcoming large surveys will generate enormous volumes of multi-wavelength imaging that require consistent, reproducible preparation before scientific analysis can begin. The package is released under an MIT license with example notebooks, lowering the barrier for researchers who need to work with data from facilities such as Hubble, JWST, Roman, or Rubin.

Core claim

We present lightstack, a Python package for combining standalone images into photometric data cubes. The workflow consists of three main steps: cropping a region of interest from a mosaic across all available filters; stacking the images to construct the data cube; and performing PSF matching on the cube. This package is intended for preparing data for studies involving multi-band photometry.

What carries the argument

The lightstack package and its three-step workflow of cropping a region of interest, stacking the images, and performing PSF matching on the resulting cube.

If this is right

  • Researchers gain a reproducible way to prepare multi-filter mosaics from HST, JWST, Roman, or Rubin data for photometry.
  • Multi-band studies tracing physical processes across wavelengths can start from consistently processed cubes.
  • The open MIT-licensed code allows direct inspection and modification of each step in the pipeline.
  • Jupyter tutorials supplied with the package demonstrate end-to-end use on real survey data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The package may be extended to accept time-domain or spectroscopic inputs without changing its core cropping-stacking-PSF sequence.
  • Direct comparison benchmarks against other astronomy image-combination tools would clarify where lightstack sits relative to existing alternatives.
  • Integration hooks for downstream machine-learning photometry pipelines could be added once the cubes are produced.

Load-bearing premise

The three-step workflow of cropping, stacking, and PSF matching is sufficient and appropriate for preparing data for the intended multi-band photometry studies.

What would settle it

A side-by-side test in which lightstack cubes yield photometric measurements that differ systematically from those produced by an independent, established pipeline on the same input images would show the workflow is not adequate.

Figures

Figures reproduced from arXiv: 2606.20360 by Alberto Krone-Martins, Ana L. Chies-Santos, Andressa Wille, Emille E. O. Ishida, Rafael S. de Souza, Thallis Pessi.

Figure 1
Figure 1. Figure 1: Workflow of the code. The example illustrates the steps for constructing a PSF matched photometric data cube, in this case using JWST images. In the first step, a region of interest is cropped: a spiral galaxy observed in the F090W filter. This cropping must be repeated for all available filters. In the second step, all cutouts are stacked in order of wavelengths, resulting in a data cube with dimensions (… view at source ↗
read the original abstract

Multi-band photometry traces diverse physical processes across a wide range of wavelengths. In recent decades, this field has been driven by the rapid growth of multi-imaging datasets, from high-resolution observation from Hubble Space Telescope and James Webb Space Telescope to the forthcoming large-scale surveys enabled by the Roman Space Telescope and Rubin Observatory, for example. In this work, we present lightstack, a Python package for combining standalone images into photometric data cubes. The workflow consists of three main steps: cropping a region of interest from a mosaic across all available filters; stacking the images to construct the data cube; and performing PSF matching on the cube. This package is intended for preparing data for studies involving multi-band photometry. The code is released under an MIT license and is available on GitHub together with a Jupyter tutorial notebook. The version used for this publication (v0.2.1) is archived on Zenodo.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript presents lightstack, a Python package for constructing photometric data cubes from standalone multi-band images. The described workflow consists of three steps: cropping a region of interest across filters from mosaics, stacking the cropped images into a data cube, and applying PSF matching to the resulting cube. The package is released under an MIT license, with source code and a Jupyter tutorial notebook on GitHub and version v0.2.1 archived on Zenodo; it is positioned as a tool to prepare data for multi-band photometry studies in the context of growing imaging datasets from HST, JWST, Roman, and Rubin.

Significance. If the package implements the stated workflow correctly, it supplies a convenient, open-source utility for data preparation in multi-wavelength photometry, a task that is increasingly relevant given the volume of imaging data from current and upcoming facilities.

minor comments (2)
  1. [Abstract] Abstract: the closing phrase 'for example' after the list of telescopes is redundant and disrupts sentence flow; rephrasing the clause would improve readability.
  2. [Abstract] Abstract: no references are provided to comparable existing tools or packages for multi-band cube construction, which would help situate the contribution.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive review of the manuscript and their recommendation to accept.

Circularity Check

0 steps flagged

Tool-release description with no derivations or load-bearing claims

full rationale

The manuscript is a software package announcement. It states the purpose of lightstack, lists the three workflow steps (cropping, stacking, PSF matching), and notes the MIT license and availability on GitHub/Zenodo. No equations, fitted parameters, predictions, uniqueness theorems, or self-citations appear. The central claim is simply that the package implements the described workflow; this is not derived from prior results within the paper and does not reduce to any input by construction. The text is self-contained as a tool description and contains no circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

This is a software-package description paper rather than a derivation of new scientific results; therefore the ledger contains no free parameters, axioms, or invented entities.

pith-pipeline@v0.9.1-grok · 5708 in / 1113 out tokens · 20374 ms · 2026-06-26T15:31:14.910756+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

17 extracted references · 15 canonical work pages · 10 internal anchors

  1. [1]

    The Propagation of Uncertainties in Stellar Population Synthesis Modeling. I. The Relevance of Uncertain Aspects of Stellar Evolution and the Initial Mass Function to the Derived Physical Properties of Galaxies. , keywords =. doi:10.1088/0004-637X/699/1/486 , archivePrefix =. 0809.4261 , primaryClass =

  2. [2]

    J-PAS: The Javalambre-Physics of the Accelerated Universe Astrophysical Survey

    J-PAS: The Javalambre-Physics of the Accelerated Universe Astrophysical Survey. arXiv e-prints , keywords =. doi:10.48550/arXiv.1403.5237 , archivePrefix =. 1403.5237 , primaryClass =

  3. [3]

    doi:10.1111/j.1365-2966.2005.09839.x , eprint =

    Evolutionary population synthesis: models, analysis of the ingredients and application to high-z galaxies. , keywords =. doi:10.1111/j.1365-2966.2005.09270.x , archivePrefix =. astro-ph/0410207 , primaryClass =

  4. [4]

    UV-extended E-MILES stellar population models: young components in massive early-type galaxies

    UV-extended E-MILES stellar population models: young components in massive early-type galaxies. , keywords =. doi:10.1093/mnras/stw2231 , archivePrefix =. 1612.01187 , primaryClass =

  5. [5]

    , keywords =

    The miniJPAS survey: A preview of the Universe in 56 colors. , keywords =. doi:10.1051/0004-6361/202038841 , archivePrefix =. 2007.01910 , primaryClass =

  6. [6]

    The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s

    The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s. arXiv e-prints , keywords =. doi:10.48550/arXiv.1902.05569 , archivePrefix =. 1902.05569 , primaryClass =

  7. [7]

    LSST: from Science Drivers to Reference Design and Anticipated Data Products

    LSST: From Science Drivers to Reference Design and Anticipated Data Products. , keywords =. doi:10.3847/1538-4357/ab042c , archivePrefix =. 0805.2366 , primaryClass =

  8. [8]

    , keywords =

    JADES Initial Data Release for the Hubble Ultra Deep Field: Revealing the Faint Infrared Sky with Deep JWST NIRCam Imaging. , keywords =. doi:10.3847/1538-4365/acf44d , archivePrefix =. 2306.02466 , primaryClass =

  9. [9]

    SAGUI: SED-based Segmentation of Multi-band Galaxy Images -- Application to JADES in GOODS-South

    SAGUI: SED-based Segmentation of Multi-band Galaxy Images -- Application to JADES in GOODS-South. arXiv e-prints , keywords =. doi:10.48550/arXiv.2604.18812 , archivePrefix =. 2604.18812 , primaryClass =

  10. [10]

    Overview of the JWST Advanced Deep Extragalactic Survey (JADES)

    Overview of the JWST Advanced Deep Extragalactic Survey (JADES). , keywords =. doi:10.3847/1538-4365/ae3163 , archivePrefix =. 2306.02465 , primaryClass =

  11. [11]

    , keywords =

    The JADES Origins Field: A New JWST Deep Field in the JADES Second NIRCam Data Release. , keywords =. doi:10.3847/1538-4365/ae1137 , archivePrefix =. 2310.12340 , primaryClass =

  12. [12]

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

    The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package. , keywords =. doi:10.3847/1538-4357/ac7c74 , archivePrefix =. 2206.14220 , primaryClass =

  13. [13]

    The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package

    The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package. , keywords =. doi:10.3847/1538-3881/aabc4f , archivePrefix =. 1801.02634 , primaryClass =

  14. [14]

    Astropy: A Community Python Package for Astronomy

    Astropy: A community Python package for astronomy. , keywords =. 2013. doi:10.1051/0004-6361/201322068 , archivePrefix =. 1307.6212 , primaryClass =

  15. [15]

    reproject: Python-based astronomical image reprojection

  16. [16]

    , year = 2007, month = may, journal =

    Matplotlib: A 2D Graphics Environment. Computing in Science and Engineering , keywords =. doi:10.1109/MCSE.2007.55 , adsurl =

  17. [17]

    2026 , doi =

    Wille, Andressa and Pessi, Thallis , title =. 2026 , doi =