Seeing Through the RFI Mask: DT-FISTA Inpainting of Pulsar Dynamic Spectra via Fourier Sparsity
Abstract
Radio frequency interference (RFI) frequently contaminates radio observations, and excising affected samples leaves masked regions in sampled images. Replacing masked samples with zero or the mean causes spectral leakage in Fourier-domain power representations. This contamination is particularly severe for modern pulsar scintillation studies, which rely on identifying and characterizing faint parabolic scintillation arcs in secondary spectra, the 2D power spectra of dynamic spectra. Scintillation arcs probe discrete scattering screens in the interstellar medium (ISM) and constrain their physical properties. However, leakage can obscure faint arcs whose power is often <10−3 of the dominant pulsar signal power. Wiener inpainting can mitigate this leakage, but a dense implementation requires 𝒪(N2) memory and 𝒪(JN2) time for N pixels and J iterations, which is computationally expensive for full-resolution spectra containing millions of pixels. We introduce the decreasing-threshold fast iterative shrinkage-thresholding algorithm (DT-FISTA). Using a sparse Fourier representation, DT-FISTA initially estimates masked intensities from the strongest Fourier components and progressively refines them with weaker components as the soft threshold decreases. Nesterov momentum extrapolation accelerates these updates. Each iteration uses 2D fast Fourier transforms (FFT) and pointwise operations, requiring 𝒪(Nlog N) time and 𝒪(N) memory. Tests on 300×300 pixel patches from Five-hundred-meter Aperture Spherical radio Telescope (FAST) observations yield lower reconstruction errors than the dense Wiener implementation, with runtime and memory requirements reduced by factors >103. Applied to full-resolution dynamic spectra, DT-FISTA suppresses mask-induced artifacts and reveals previously obscured faint arcs, enabling efficient ISM studies.
Related
- Fast Radio Burst Classification on Dispersed Dynamic Spectra Using a Wavelet-based Convolutional Network
- Improved Understanding of Lunar Surface Hydration Signatures through Comparison of Thermally Corrected Spectra from the Moon Mineralogy Mapper (M 3)
- Deciphering the VHE gamma-rays spectra of Extreme High-Frequency peaked BL Lac Objects through a Photohadronic Model
- Determining elemental composition in laboratory meteorite ablation spectra through radiative transfer modeling
- Editorial: Dynamic exospheres of terrestrial bodies through the solar system
BibTeX
@misc{huang2026seeing,
title = {{Seeing Through the RFI Mask: DT-FISTA Inpainting of Pulsar Dynamic Spectra via Fourier Sparsity}},
author = {Yuanshang Huang and Xun Shi},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.23140},
eprint = {2610.23140},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.23140},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Seeing Through the RFI Mask: DT-FISTA Inpainting of Pulsar Dynamic Spectra via Fourier Sparsity
AU - Yuanshang Huang
AU - Xun Shi
PY - 2026
DA - 2026/10/16
PB - arXiv
DO - 10.48550/arxiv.2610.23140
UR - https://arxiv.org/abs/2610.23140
AB - Radio frequency interference (RFI) frequently contaminates radio observations, and excising affected samples leaves masked regions in sampled images. Replacing masked samples with zero or the mean causes spectral leakage in Fourier-domain power representations. This contamination is particularly severe for modern pulsar scintillation studies, which rely on identifying and characterizing faint parabolic scintillation arcs in secondary spectra, the 2D power spectra of dynamic spectra. Scintillation arcs probe discrete scattering screens in the interstellar medium (ISM) and constrain their physical properties. However, leakage can obscure faint arcs whose power is often <10^−3 of the dominant pulsar signal power. Wiener inpainting can mitigate this leakage, but a dense implementation requires 𝒪(N^2) memory and 𝒪(JN^2) time for N pixels and J iterations, which is computationally expensive for full-resolution spectra containing millions of pixels. We introduce the decreasing-threshold fast iterative shrinkage-thresholding algorithm (DT-FISTA). Using a sparse Fourier representation, DT-FISTA initially estimates masked intensities from the strongest Fourier components and progressively refines them with weaker components as the soft threshold decreases. Nesterov momentum extrapolation accelerates these updates. Each iteration uses 2D fast Fourier transforms (FFT) and pointwise operations, requiring 𝒪(Nlog N) time and 𝒪(N) memory. Tests on 300×300 pixel patches from Five-hundred-meter Aperture Spherical radio Telescope (FAST) observations yield lower reconstruction errors than the dense Wiener implementation, with runtime and memory requirements reduced by factors >10^3. Applied to full-resolution dynamic spectra, DT-FISTA suppresses mask-induced artifacts and reveals previously obscured faint arcs, enabling efficient ISM studies.
ER -