arxiv: 2604.07374 · v1 · submitted 2026-04-07 · 🌌 astro-ph.IM

Recognition: no theorem link

PowerSpectR: An R Package for Radial Power Spectrum Estimation

Rafael S. de Souza

Authors on Pith no claims yet

Pith reviewed 2026-05-10 17:56 UTC · model grok-4.3

classification 🌌 astro-ph.IM

keywords power spectrumradial Fourier analysisimage morphologyR packageazimuthal medianspatial scalesFourier transformedge effects

0 comments

The pith

PowerSpectR is an R package that estimates radial Fourier power spectra from images by using azimuthal medians to produce slopes less biased by localized features.

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

The paper presents PowerSpectR as a tool for computing and visualizing median-based radial Fourier power spectra from imaging data. It describes a workflow that applies Hann windowing to reduce edge effects, transforms the data in the Fourier domain, and bins the resulting power radially while taking the median across angles instead of the mean. The author states that this median step helps limit the impact of bright compact sources, masking artifacts, and similar localized features that can distort results from standard estimators. The resulting power spectrum slopes are positioned as compact summaries of morphological complexity across different images. The package is made available under an open license for use in the R environment.

Core claim

PowerSpectR provides a workflow for estimating these slopes, combining edge-effect mitigation through Hann windowing, Fourier-domain analysis, and radial binning with azimuthal median statistics. The use of median aggregation helps to reduce sensitivity to bright compact sources, masking artifacts, and other localized features that can bias standard estimators.

What carries the argument

Azimuthal median statistics applied to radially binned Fourier power values after Hann windowing, which produces radial profiles by aggregating across angles with the median rather than the mean.

If this is right

Power spectrum slopes become usable as low-dimensional summaries of morphological complexity in imaging data.
The workflow supports visualization of the spectra alongside the derived slopes.
The package runs in R and is distributed under the MIT license.
Radial binning with medians is presented as directly applicable to any 2D image dataset.

Where Pith is reading between the lines

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

The median-based approach could be tested for consistency on simulated images where the true underlying power-law slope is known in advance.
The same radial-median technique might apply to non-image data arranged on a 2D grid, such as spatial maps from other scientific domains.
Users could extend the package by adding options for alternative window functions or binning schemes to compare robustness across choices.

Load-bearing premise

Azimuthal median statistics meaningfully reduce sensitivity to localized features compared with standard mean-based estimators.

What would settle it

A side-by-side comparison of slope values obtained from the same images with and without injected compact bright sources, computed once with azimuthal medians and once with means, to measure which aggregation method shows smaller shifts in the fitted slopes.

Figures

Figures reproduced from arXiv: 2604.07374 by Rafael S. de Souza.

**Figure 1.** Figure 1: Comparison of Hubble Space Telescope image cutouts of M101 (spiral) and M60 (elliptical). Left: images. Right: median radial power spectra computed with PowerSpectR. Credit: M60—NASA, ESA, and the Hubble Heritage (STScI/AURA)-ESA/Hubble Collaboration; M101—NASA, ESA, K. Kuntz (JHU), F. Bresolin (University of Hawaii), J. Trauger (JPL), J. Mould (NOAO), Y.-H. Chu (University of Illinois, Urbana), and STScI … view at source ↗

read the original abstract

I present here PowerSpectR, an R package for computing and visualizing median-based radial Fourier power spectra from imaging data. Power spectra provide a representation of spatial structure by decomposing contributions across spatial scales, and the resulting slopes can serve as compact, low-dimensional summaries of morphological complexity across images. PowerSpectR provides a workflow for estimating these slopes, combining edge-effect mitigation through Hann windowing, Fourier-domain analysis, and radial binning with azimuthal median statistics. The use of median aggregation helps to reduce sensitivity to bright compact sources, masking artifacts, and other localized features that can bias standard estimators. PowerSpectR is released under the MIT license at \href{https://github.com/RafaelSdeSouza/PowerSpectR}{this repository}.

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.

Desk Editor's Note private letter to a colleague

PowerSpectR is a basic R package for radial power spectra with median binning, but the manuscript supplies no tests or comparisons to support the claim that the median reduces bias.

read the letter

PowerSpectR is an R package that computes radial Fourier power spectra from images by applying Hann windowing, taking the FFT, binning radially, and using azimuthal medians rather than means. The abstract presents this as a way to get slopes less affected by compact sources and artifacts, yet the text contains no validation whatsoever for that advantage. No synthetic tests, no side-by-side comparisons with mean-based spectra, no bias or variance numbers, and no worked examples on real data appear anywhere in the manuscript. The package itself is the only concrete output: it bundles standard operations into an R workflow with some visualization helpers and is released under MIT on GitHub. For R users who already do image analysis in astrophysics and want a ready-made function for power-spectrum slopes, this could be a minor convenience. The underlying steps (Hann window, radial binning, median aggregation) are all drawn from existing literature, so nothing here is conceptually new. The median choice is a plausible robustness tweak, but without evidence it remains just a choice. The paper is therefore a software announcement rather than a research result. It is aimed at a narrow slice of R-using astronomers who need this exact workflow and are willing to trust the untested claim. Most readers looking for methodological progress or reproducible improvements will find little of value. I would not bring it to a reading group, would not cite it, and would not send it for peer review; it belongs in a tools note or repository readme if it belongs anywhere.

Referee Report

1 major / 2 minor

Summary. The manuscript presents PowerSpectR, an R package for estimating radial power spectra from images using a workflow that incorporates Hann windowing, Fourier transformation, radial binning, and azimuthal median statistics. The key innovation claimed is that the median aggregation in the azimuthal direction reduces the impact of localized features such as bright compact sources and masking artifacts on the estimated power spectrum slopes.

Significance. Should the median-based approach prove effective in reducing bias as claimed, the package would provide a practical tool for researchers in astrophysics and image analysis to obtain more robust summaries of spatial structure in imaging data. The open release of the code under MIT license with a public GitHub repository supports reproducibility and community adoption.

major comments (1)

Abstract: The claim that 'the use of median aggregation helps to reduce sensitivity to bright compact sources, masking artifacts, and other localized features that can bias standard estimators' is stated without any supporting validation, synthetic tests, comparisons to mean-based radial spectra, bias/variance quantification, or worked examples. This assertion is load-bearing for the package's central advantage and remains unverified in the text.

minor comments (2)

The manuscript is brief and would benefit from a dedicated usage section or vignette with at least one concrete example on sample imaging data to illustrate the workflow and output.
Consider adding installation instructions, required dependencies, and basic code snippets for core functions to improve accessibility for users.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their report and the opportunity to respond. We address the single major comment below and describe the planned revisions.

read point-by-point responses

Referee: Abstract: The claim that 'the use of median aggregation helps to reduce sensitivity to bright compact sources, masking artifacts, and other localized features that can bias standard estimators' is stated without any supporting validation, synthetic tests, comparisons to mean-based radial spectra, bias/variance quantification, or worked examples. This assertion is load-bearing for the package's central advantage and remains unverified in the text.

Authors: We agree that the manuscript as submitted does not contain synthetic tests, direct comparisons to mean-based estimators, or quantitative bias/variance results to support the abstract claim. This is a substantive gap. In the revised manuscript we will add a dedicated validation section that includes: (i) controlled simulations of power-law images with injected compact sources and artificial masks, (ii) side-by-side mean versus median radial spectra, (iii) Monte-Carlo quantification of slope bias and variance, and (iv) worked examples on both synthetic and real data. These additions will be referenced from the abstract and will directly substantiate the claimed robustness of the median approach. revision: yes

Circularity Check

0 steps flagged

No circularity: software implementation of standard Fourier operations

full rationale

The manuscript describes an R package that applies well-known operations (Hann windowing, 2D FFT, radial binning, and azimuthal median aggregation) to imaging data. No derivation chain, equations, fitted parameters, or self-referential definitions appear in the provided text. The statement that median statistics reduce sensitivity to compact sources is presented as a design choice without any supporting equations or reductions to prior results within the paper itself. No self-citations, uniqueness theorems, or ansatzes are invoked as load-bearing steps. The work is therefore self-contained as a software tool rather than a theoretical derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The description introduces no free parameters, background axioms, or postulated entities; it relies on standard Fourier analysis and median statistics already available in the prior literature.

pith-pipeline@v0.9.0 · 5414 in / 1010 out tokens · 63251 ms · 2026-05-10T17:56:59.855643+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

12 extracted references · 11 canonical work pages

[1]

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

Akeson, R., Armus, L., Bachelet, E., et al. 2019, arXiv:1902.05569. https://arxiv.org/abs/1902.05569

work page Pith review arXiv 2019
[2]

2025, imager: Image Processing Library Based on ’CImg’, doi: 10.32614/CRAN.package.imager

Barthelme, S. 2025, imager: Image Processing Library Based on ’CImg’, doi: 10.32614/CRAN.package.imager

work page doi:10.32614/cran.package.imager 2025
[3]

B., & Tukey, J

Blackman, R. B., & Tukey, J. W. 1958, Bell System Technical Journal, 37, 185, doi: https://doi.org/10.1002/j.1538-7305.1958.tb03874.x

work page doi:10.1002/j.1538-7305.1958.tb03874.x 1958
[4]

2013, The Astrophysical Journal, 771, 123, doi: 10.1088/0004-637X/771/2/123 CSST Collaboration, Gong, Y., Miao, H., et al

Burkhart, B., Lazarian, A., Ossenkopf, V., & Stutzki, J. 2013, The Astrophysical Journal, 771, 123, doi: 10.1088/0004-637X/771/2/123 CSST Collaboration, Gong, Y., Miao, H., et al. 2026, Science China Physics, Mechanics, and Astronomy, 69, 239501, doi: 10.1007/s11433-025-2809-0 de Souza, R. S. 2026, PowerSpectR: Robust Median-Based Power Spectrum Analysis ...

work page doi:10.1088/0004-637x/771/2/123 2013
[5]

G., & Scalo, J

Elmegreen, B. G., & Scalo, J. 2004, Annual Review of Astronomy and Astrophysics, 42, 211, doi: 10.1146/annurev.astro.41.011802.094859 Euclid Collaboration, Blanchard, A., Camera, S., et al. 2020, A&A, 642, A191, doi: 10.1051/0004-6361/202038071 Ivezi´ c, v., et al. 2019, The Astrophysical Journal, 873, 111, doi: 10.3847/1538-4357/ab042c

work page doi:10.1146/annurev.astro.41.011802.094859 2004
[6]

W., & Rosolowsky, E

Koch, E. W., & Rosolowsky, E. W. 2016, MNRAS, 452, 3435, doi: 10.1093/mnras/stv1370

work page doi:10.1093/mnras/stv1370 2016
[7]

W., Rosolowsky, E

Koch, E. W., Rosolowsky, E. W., Boyden, R. D., et al. 2019, AJ, 158, 1, doi: 10.3847/1538-3881/ab1cc0

work page doi:10.3847/1538-3881/ab1cc0 2019
[8]

A., de Souza, R

Kuhn, M. A., de Souza, R. S., Krone-Martins, A., et al. 2021, ApJS, 254, 33, doi: 10.3847/1538-4365/abe465

work page doi:10.3847/1538-4365/abe465 2021
[9]

2000, The Astrophysical Journal, 537, 720, doi: 10.1086/309040

Lazarian, A., & Pogosyan, D. 2000, The Astrophysical Journal, 537, 720, doi: 10.1086/309040

work page doi:10.1086/309040 2000
[10]

J., others, & de Souza, R

Pessi, P. J., others, & de Souza, R. S. 2024, Astronomy & Astrophysics, 691, A181, doi: 10.1051/0004-6361/202450535 R Core Team. 2025, R: A Language and Environment for Statistical Computing, R Foundation for Statistical

work page doi:10.1051/0004-6361/202450535 2024
[11]

G., Krone-Martins, A., et al

Stern, D., Djorgovski, S. G., Krone-Martins, A., et al. 2021, The Astrophysical Journal, 921, 42, doi: 10.3847/1538-4357/ac0f04 van der Schaaf, A., & van Hateren, J. 1996, Vision Research, 36, 2759, doi: https://doi.org/10.1016/0042-6989(96)00002-8

work page doi:10.3847/1538-4357/ac0f04 2021
[12]

2016, ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag New York)

Wickham, H. 2016, ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag New York). https://ggplot2.tidyverse.org

2016