Pith. sign in

REVIEW 3 major objections 7 minor 26 references

RatanSunPy: A robust preprocessing pipeline for RATAN-600 solar radio observations data

T0 review · 3 major / 7 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read RatanSunPy is an open-source Python package that automates the full processing chain for RATAN-600 solar radio scans: quiet-Sun-template calibration, local microwave source detection, NOAA active-region association, and Gaussian…

desk verdict Useful Python pipeline for RATAN-600 data, but the printed calibration formula is inverted and the paper needs a benchmark before I'd call it robust. read the letter →

arxiv 2412.08230 v1 pith:WG6WMBFI submitted 2024-12-11 astro-ph.SR astro-ph.IMcs.DLphysics.space-ph

I. Knyazeva , I. Lysov , E. Kurochkin , A. Shendrik , D. Derkach , N. Makarenko This is my paper · ORCID
classification astro-ph.SRastro-ph.IMcs.DLphysics.space-ph
keywords RATAN-600solarradioastronomyquietSuncalibrationactiveregionsmicrowavespectroscopyPythonpackageGaussiananalysiscycleobservations
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

RatanSunPy is a Python package that converts raw one-dimensional multi-frequency scans from the RATAN-600 radio telescope into calibrated solar spectra and structured tables of local sources and active regions. The authors are trying to establish that the decades-long RATAN-600 archive can be processed end-to-end without specialized standalone software: the package loads data, cleans it, calibrates it against a quiet-Sun template, detects local microwave sources from averaged circular polarization, and identifies them with NOAA active regions. If the package works as claimed, a 3-18 GHz record of the solar chromosphere and corona spanning more than two decades becomes usable by the broader solar-physics and machine-learning communities. The central bet is that a single quiet-Sun template built from minimum-activity observations can serve as the absolute calibration reference for every scan.

What carries the argument

The load-bearing mechanism is the quiet-Sun template calibration: a set of one-dimensional disk profiles at each frequency, built from minimum-activity observations, that is used to scale the background constant of every raw scan so that the observation matches a model quiet Sun. Supporting machinery includes the conversion of brightness temperatures to flux densities through the Rayleigh-Jeans law, symlet-wavelet denoising of the circular polarization signal for peak detection, a coordinate rotation that maps NOAA Solar Region Summary positions to the RATAN-600 scan geometry, and least-squares fitting of Gaussian components to active-region profiles to extract physical parameters.

What would settle it

Compare RatanSunPy-calibrated total solar fluxes, at several frequencies, with independent measurements from a small-antenna solar radio instrument such as the Nobeyama polarimeters across a range of activity levels. If the offsets between the two grow with solar activity or change with frequency in a way that tracks the difference between the template and the day's true quiet-Sun profile, the template-calibration assumption is falsified.

Watch

Extended reading notes

Core claim

The paper presents RatanSunPy as a complete, reproducible preprocessing pipeline for RATAN-600 solar observations. Its calibration method fits each raw scan so that its constant background level matches a quiet-Sun template: semi-profiles of the solar disk at each frequency, constructed from observations during minimal solar activity and normalized so that the area under the profile equals the total solar flux predicted by a historical brightness-temperature model via the Rayleigh-Jeans law. After calibration, the package denoises the frequency-averaged circular polarization signal with symlet wavelets, finds peaks corresponding to local sources, transforms Solar Region Summary positions into heliocentric coordinates, and matches sources to active regions. For each matched active region, Gaussian analysis at every frequency yields amplitudes, widths, fluxes, and brightness temperatures. The claimed deliverable is therefore both calibrated scans and physical source catalogs, ready for comparison with other wavelengths and for time-series studies.

Load-bearing premise

Every observation is calibrated by scaling its constant background to a quiet-Sun template built from minimum-activity observations and a historical brightness-temperature model; if that template is wrong for a given frequency, epoch, or solar-cycle phase, all calibrated fluxes, amplitudes, and brightness temperatures inherit the bias.

Editorial extensions

If this is right

  • A researcher can go from a raw RATAN-600 FITS file to a calibrated full-disk spectrum in one Python environment, with no local directory setup or standalone GUI.
  • Local microwave sources, including gyroresonance sources visible in circular polarization, are detected automatically and matched to NOAA active regions for the same day.
  • For each active region, the package produces frequency-resolved amplitudes, widths, fluxes, and brightness temperatures, enabling multi-day spectral evolution studies.
  • Because the archive spans solar cycles 23, 24, and part of 25, the pipeline makes cycle-length comparisons of microwave active-region properties practical.
  • The Python interface allows direct integration with machine-learning libraries, which the authors argue can support automated flare-forecasting systems.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable consequence of the template-calibration design is that all absolute fluxes inherit the quiet-Sun model's assumptions; comparing RatanSunPy outputs against independent total-flux measurements (for example, Nobeyama polarimeters) would reveal any activity-dependent bias.
  • Because source detection runs on circular polarization, active regions whose microwave emission is purely thermal and unpolarized will be underrepresented in automatically produced catalogs; combining with intensity-based detection would change the recoverable population.
  • If the template is stable over solar cycles, the package enables a self-consistent homogenized microwave database of active-region brightness temperatures across cycles 23-25; if not, cycle comparisons will need per-cycle calibration templates.
  • The package's design suggests a natural extension: mapping the one-dimensional detected sources back onto full-disk magnetograms could produce automated, daily microwave-magnetogram association statistics, which the paper does not itself evaluate.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. The paper presents RatanSunPy, an open-source Python package for accessing, calibrating, and analyzing one-dimensional multi-frequency solar observations from the RATAN-600 radio telescope. The package provides two data clients (RATANClient and SRSClient), a pipeline that calibrates raw scans using quiet-Sun templates, automatic detection of local sources from the averaged circular-polarization signal, matching of detected sources to NOAA active regions, and Gaussian fitting to estimate source parameters. The manuscript describes the package architecture, shows a usage example, and reports that the code is tested with pytest and doctest.

Significance. If the calibration and source-analysis steps are quantitatively validated, RatanSunPy would be a valuable community resource: it lowers the barrier to working with a multi-decade archive of microwave solar observations, is open source, integrates with the wider Python and SunPy ecosystem, and provides automated table outputs that could support active-region tracking and space-weather studies. The paper also deserves credit for shipping tests, documentation, and usage examples. However, the central claim of a 'robust preprocessing pipeline' is not yet supported by the evidence presented: there is no quantitative validation of the calibrated absolute flux scale, of source-detection completeness or false-positive rates, or of the active-region matching accuracy, and Section 8 explicitly acknowledges instability in source localization and Gaussian analysis.

major comments (3)
  1. [4.1, Eq. (2)] Equation (2) is dimensionally incorrect. The Rayleigh-Jeans law gives B_nu = 2 k_B T_b nu^2 / c^2, so the spectral flux density is proportional to nu^2/c^2. Equation (2) instead prints c^2/nu^2, which inverts the factor and would change the result by (c/nu)^4, roughly 8e-7 at 10 GHz. Because this conversion normalizes the quiet-Sun template and therefore sets the absolute flux scale for all calibrated scans and source parameters, this is a load-bearing step. Please correct the formula and, more importantly, verify the implemented code with a unit test and an independent calibration check such as comparison with Nobeyama or RSTN total-flux measurements.
  2. [6 and 8] The paper does not provide quantitative validation of the claimed robust preprocessing. Section 6 reports only that tests exist (pytest and doctest), which checks for coding bugs but not scientific accuracy. Section 8 admits that source localization and Gaussian analysis produce unstable results in certain cases. To support the central claim, please add quantitative metrics: calibrated flux densities compared with independent instruments, detection recall and precision against a labeled sample of active regions, and a characterization of when Gaussian analysis fails and by how much.
  3. [4.1] The calibration method assumes that the quiet-Sun template, based on solar-minimum observations and historical brightness temperatures from Shendrik et al. (2020) and Borovik (1997), provides the correct background level for every observation. Since the calibration procedure scales the observed background to match this template, any inaccuracy of the template at a given frequency, epoch, or solar-cycle phase would propagate directly into all source amplitudes and brightness temperatures. The manuscript gives no sensitivity analysis or independent benchmark for this assumption; please either validate the template against independent absolute calibration or quantify the resulting systematic uncertainty.
minor comments (7)
  1. [Abstract] The sentence 'These levels remain some difficult to detect in the ultraviolet and X-ray ranges' is ungrammatical and should be rewritten.
  2. [3] The URL in the first paragraph contains a stray closing bracket: 'http://www.spbf.sao.ru/prognoz/]' should be 'http://www.spbf.sao.ru/prognoz/'.
  3. [4.2] The automatic peak detection relies on a discrete wavelet transform and the find_peaks function, but the wavelet parameters and detection threshold are not specified; please state them or refer to the source code with a versioned release.
  4. [4.4, Eq. (4)] Equation (4) uses c for the Gaussian width, which conflicts with c used for the speed of light in Eq. (2); consider renaming the width parameter to sigma.
  5. [2.2] The caption of Figure 2 (and the text referencing it) says '2024/31/07', which should be '2024/07/31' for consistency with Figure 1.
  6. [Throughout] The package name is spelled inconsistently: 'RatanSunPy' in the text versus 'RATANSunPy' in the repository URL and documentation links; please standardize the spelling.
  7. [6] The quality-assurance section would be stronger with a test-coverage metric and a statement that the tests include a regression test for the calibration formula, since that is the most safety-critical part of the pipeline.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the pipeline's calibrations and detections are not predictions derived from their own inputs; self-cited quiet-Sun data is an acknowledged external calibration input.

full rationale

RatanSunPy is a software pipeline, not a derivation of new physical predictions. Its calibration step scales raw RATAN-600 scans so that the background matches a quiet-Sun template; this matching is the calibration procedure itself rather than a predicted quantity. The template normalization is taken from historical brightness-temperature data attributed to previous work (Shendrik et al. 2020; Borovik 1997), so the absolute flux scale is anchored to an external, earlier data source, not to the outputs of this paper. Local-source detection uses circular-polarization extrema and Gaussian fitting on the calibrated scans; these are standard data-analysis operations and do not reduce to parameters fitted earlier and then reported as predictions. The paper also explicitly states that users may implement their own calibration methods, acknowledging that the adopted quiet-Sun calibration is a choice rather than a consequence of the package's own derivation. The self-citations present are load-bearing only in the sense that they document the origin of the calibration data and the preprocessing approach; they do not establish a theorem or forbid alternatives. The suspicious unit conversion in Eq. (2) is a potential correctness defect in the Rayleigh-Jeans transposition, but it is not a circularity: an incorrect conversion factor does not make a result equivalent to its input by construction. No fitted parameter is renamed as a prediction, and no claimed result is defined in terms of the very output it purports to produce. Therefore no circular steps are identified.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No new free parameters are fitted in this paper; the calibration table and processing steps are taken from prior SAO work. The assumptions are domain assumptions about the quiet-Sun background, Gaussian source shapes, polarization-based detection, and coordinate alignment. No new physical entities are introduced.

assumptions (4)
  • domain assumption The quiet-Sun template built from minimum-activity observations is an accurate background level for every scan.
    Section 4.1 scales the constant background of each observation to match the quiet-Sun template; if the template is biased, all calibrated fluxes inherit the bias.
  • domain assumption Local source one-dimensional profiles are approximately Gaussian after convolution with the antenna beam.
    Section 4.4 and Equation 4 model sources as Gaussians and fit them by least squares; non-Gaussian structures such as halos would bias fitted parameters.
  • domain assumption Averaged circular polarization V is a sufficient signal for detecting local sources.
    Section 4.2 uses peaks in frequency-averaged Stokes V for detection and acknowledges that thermal sources with weak polarization are missed.
  • domain assumption NOAA active-region positions from 2400 UTC can be rotated to the RATAN-600 observation time using differential solar rotation and FITS header parameters.
    Section 4.3 assumes the position angle and solar rotation corrections accurately align one-dimensional scans with two-dimensional active-region maps.

how reviews work

0 comments
Cite this review

Pith. "Pith review of RatanSunPy: A robust preprocessing pipeline for RATAN-600 solar radio observations data." pith.science (2026). https://pith.science/paper/WG6WMBFI

@misc{pith2026241208230,
  author       = {Pith},
  title        = {Pith review of: RatanSunPy: A robust preprocessing pipeline for RATAN-600 solar radio observations data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WG6WMBFI}},
  note         = {Machine review of arXiv:2412.08230}
}
read the original abstract

The advancement of observational technologies and software for processing and visualizing spectro-polarimetric microwave data obtained with the RATAN-600 radio telescope opens new opportunities for studying the physical characteristics of solar plasma at the levels of the chromosphere and corona. These levels remain some difficult to detect in the ultraviolet and X-ray ranges. The development of such methods allows for more precise investigation of the fine structure and dynamics of the solar atmosphere, thereby deepening our understanding of the processes occurring in these layers. The obtained data also can be utilized for diagnosing solar plasma and forecasting solar activity. However, using RATAN-600 data requires extensive data processing and familiarity with the RATAN-600. This paper introduces RatanSunPy, an open-source Python package developed for accessing, visualizing, and analyzing multi-band radio observations of the Sun from the RATAN-600 solar complex. The package offers comprehensive data processing functionalities, including direct access to raw data, essential processing steps such as calibration and quiet Sun normalization, and tools for analyzing solar activity. This includes automatic detection of local sources, identifying them with NOAA (National Oceanic and Atmospheric Administration) active regions, and further determining parameters for local sources and active regions. By streamlining data processing workflows, RatanSunPy enables researchers to investigate the fine structure and dynamics of the solar atmosphere more efficiently, contributing to advancements in solar physics and space weather forecasting.

Figures

Figures reproduced from arXiv: 2412.08230 by the authors.

Figure 2
Figure 2. RATAN-600 solar data for 2024/31/07, same date as for [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Overlay of RATAN-600 Observation on SDO HMI Magnetogram [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Workflow of RatanSunPy package [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figures from the paper (7 more)
Figure 5
Figure 5. Figure 5: Entity structure in RatanSunPy expand the range of calibration approaches available in future updates to accommodate broader scientific applications. The idea behind the method is as follows. A series of solar observations with minimal solar activity during the minimum…
Figure 7
Figure 7. Figure 7: Obtained profiles of the quiet Sun der the curve matches the model of the quiet Sun. This model is based on the historical observations of the Sun in terms of brightness temperatures at different frequencies, as exemplified in the table presented in Shendrik et al. (20…
Figure 9
Figure 9. Figure 9: Example of raw scan from Ratan-600 4.2. Automatic detection of local sources A local source (LS) on the Sun corresponds to a local en￾hancement of the magnetic field and is visually characterized by a peak on the intensity graph. For automatic peak detection, the find …
Figure 10
Figure 10. Figure 10: Example of calibrated scan from Ratan-600 [PITH_FULL_IMAGE:figures/full_fig_p007_10.png]
Figure 11
Figure 11. Figure 11: Denoised signal 4.3. Identification of local sources with NOAA active regions To identify local sources of radiation (e.g., interspot sources, sources above the leading sunspot, sources at the base of loops, etc.), it is necessary to correlate the data with observatio…
Figure 13
Figure 13. Figure 13: Three day of evolution of AR 2673 5. Usage example 5.1. Usage Example ”Comprehensive usage examples can be found in the exam￾ple notebooks available on the project’s github repository 3 and within the official documentation 4 . The main functionality and operations wi…
Figure 12
Figure 12. Figure 12: Gaussian analysis for AR 2674 on September 3, 2017 [PITH_FULL_IMAGE:figures/full_fig_p008_12.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Dynamics of microwave sources associated with the neutral line and the magnetic-field parameters of sunspots as a factor in predicting large flares

    Abramov-Maximov, V ., Borovik, V ., Opeikina, L., Tlatov, A., 2015. Dynamics of microwave sources associated with the neutral line and the magnetic-field parameters of sunspots as a factor in predicting large flares. Solar Physics 290, 53–77

    work page 2015

  2. [2]

    The measure- ment of magnetic fields in the solar atmosphere above sunspots using gy- roresonance emission

    Akhmedov, S., Gelfreikh, G., Bogod, V ., Korzhavin, A., 1982. The measure- ment of magnetic fields in the solar atmosphere above sunspots using gy- roresonance emission. Solar Physics 79, 41–58

    work page 1982

  3. [3]

    Preface: Radio heliophysics: Science and forecasting

    Leibacher, J., 2015. Preface: Radio heliophysics: Science and forecasting. Solar Physics 290, 2393–2396

    work page 2015

  4. [4]

    Ratan-600 radio telescope in the 24th solar-activity cycle

    Bogod, V ., 2011. Ratan-600 radio telescope in the 24th solar-activity cycle. i. new opportunities and tasks. Astrophysical bulletin 66, 190–204

    work page 2011

  5. [5]

    Ratan-600 radio telescope in the 24th solar-activity cycle

    Bogod, V ., Alesin, A., Pervakov, A., 2011. Ratan-600 radio telescope in the 24th solar-activity cycle. ii. multi-octave spectral and polarization high- resolution solar research system. Astrophysical Bulletin 66, 205–214

    work page 2011

  6. [6]

    Study of the magneto- spheres of active regions on the sun by radio astronomy techniques

    Bogod, V ., Kal’tman, T., Peterova, N., Yasnov, L., 2017. Study of the magneto- spheres of active regions on the sun by radio astronomy techniques. Cosmic research 55, 1–11

    work page 2017

  7. [7]

    A method of forecasting solar activity based on radio astronomical observa- tions

    Bogod, V ., Svidskiy, P., Kurochkin, E., Shendrik, A., Everstov, N., 2018. A method of forecasting solar activity based on radio astronomical observa- tions. Astrophysical Bulletin 73, 478–486

    work page 2018

  8. [8]

    Study of the quiet Sun in the microwave range on the RATAN-600 and BPR radio telescopes in the XX-XXII solar activity cy- cles

    Borovik, V ., 1997. Study of the quiet Sun in the microwave range on the RATAN-600 and BPR radio telescopes in the XX-XXII solar activity cy- cles. Ph.D. thesis. SAO RAS

    work page 1997

Show all 26 references
  1. [9]

    Local maximum in the microwave spectrum of solar active regions as a factor in predicting powerful flares

    Borovik, V ., Grigor’eva, I.Y ., Korzhavin, A., 2012. Local maximum in the microwave spectrum of solar active regions as a factor in predicting powerful flares. Geomagnetism and Aeronomy 52, 1032–1043

    2012

  2. [10]

    Orthonormal bases of compactly supported wavelets

    Daubechies, I., 1988. Orthonormal bases of compactly supported wavelets. Communications on pure and applied mathematics 41, 909–996

    1988

  3. [11]

    Processing of one-dimensional data vector arrays in win- dows os, workscan version 2.3

    Garaimov, V ., 1997. Processing of one-dimensional data vector arrays in win- dows os, workscan version 2.3. Preprint . Grigor’eva, I.Y ., Livshits, M., 2014. Weak solar flares with a detectable flux of hard x rays: Specific features of microwave radiation in the corresponding...

    1997

  4. [12]

    Radioheliographs and radio reflectors in solar plasma studies

    Khaikin, V ., Storozhenko, A., Bogod, V ., 2019. Radioheliographs and radio reflectors in solar plasma studies. Astrophysical bulletin 74, 221–233

    2019

  5. [13]

    On the possibility of probing the flare productivity of an active region in the early stage of emergence

    Kutsenko, A.S., Abramenko, V .I., Kutsenko, O.K., 2021. On the possibility of probing the flare productivity of an active region in the early stage of emergence. Monthly Notices of the Royal Astronomical Society 501, 6076– 6082

    2021

  6. [14]

    The idl astronomy user’s library, in: Astronomical Data Analysis Software and Systems II, p

    Landsman, W., 1993. The idl astronomy user’s library, in: Astronomical Data Analysis Software and Systems II, p. 246

    1993

  7. [15]

    Stokes inversion techniques with neural networks: anal- ysis of uncertainty in parameter estimation

    Derkach, D., 2023. Stokes inversion techniques with neural networks: anal- ysis of uncertainty in parameter estimation. Solar Physics 298, 98

    2023

  8. [16]

    Two-dimensional mapping of the sun with the ratan-600

    Nindos, A., Alissandrakis, C., Gelfreikh, G., Borovik, V ., Korzhavin, A., Bo- god, V ., 1996. Two-dimensional mapping of the sun with the ratan-600. Solar Physics 165, 41–59

    1996

  9. [17]

    Three-dimensional structure of microwave sources from solar rotation stereoscopy versus magnetic extrapolations

    Yasnov, L., 2011. Three-dimensional structure of microwave sources from solar rotation stereoscopy versus magnetic extrapolations. The astrophysical journal 737, 82

    2011

  10. [18]

    Revisiting the estimation of solar radio emission flux densities

    Opeikina, L., Majorova, E., Korzhavin, A., 2015. Revisiting the estimation of solar radio emission flux densities. Astrophysical bulletin 70, 371–378

    2015

  11. [19]

    Ratan-600: the world’s biggest reflector at the’cross roads’

    Parijskij, Y .N., 1993. Ratan-600: the world’s biggest reflector at the’cross roads’. IEEE Antennas and Propagation Magazine 35, 7–12

    1993

  12. [20]

    Increased microwave radiation brightness as a sign of flare-producing active regions based on observations of noaa active region 12371

    Peterova, N., Topchilo, N., Kurochkin, E., 2021. Increased microwave radiation brightness as a sign of flare-producing active regions based on observations of noaa active region 12371. Geomagnetism and Aeronomy 61, S24–S35

    2021

  13. [21]

    Rfi-resistant decimeter band radiometer for the ratan-600 radio telescope

    Ripak, A., Bogod, V ., Grenkov, S., Lebedev, M., 2023. Rfi-resistant decimeter band radiometer for the ratan-600 radio telescope. Astrophysical Bulletin 78, 622–634

    2023

  14. [22]

    Spatial and temporal properties of the solar cycle from ratan-600 microwave observations

    Shendrik, A., Bogod, V ., Kuzanyan, K., Kaltman, T., Kurochkin, E., 2020. Spatial and temporal properties of the solar cycle from ratan-600 microwave observations. Geomagnetism and Aeronomy 60, 999–1006. Smol’kov, G., 2008. Radio features of initial factors of space weather. S...

    2020

  15. [23]

    Modeling of solar atmosphere parameters above sunspots using ratan-600 microwave observa- tions

    Stupishin, A., Yasnov, L., Kaltman, T., Bogod, V ., 2018. Modeling of solar atmosphere parameters above sunspots using ratan-600 microwave observa- tions. Solar physics 293, 13. 10

    2018

  16. [24]

    Ratan-600 radio telescope in the 24th solar activity cycle

    Tokhchukova, S.K., 2011. Ratan-600 radio telescope in the 24th solar activity cycle. iv. information system for ratan-600 solar observations. Astrophysical Bulletin 66, 379–388

    2011

  17. [25]

    southern sector with a flat reflector

    Tokhchukova, S.K., Korzhavin, A., Bogod, V ., Kurochkin, E., Shendrik, A., 2014. Computation of the horizontal size of the ratan-600 beam pattern for the “southern sector with a flat reflector” mode with allowance for the parameters of primary feeds. Astrophysical Bulletin 69, 356–367

    2014

  18. [26]

    Peculiarities of polarized radio emission of solar active regions

    Yasnov, L., Kal’Tman, T., Bogod, V ., 2011. Peculiarities of polarized radio emission of solar active regions. Astronomy reports 55, 82–90. 11

    2011

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.