REVIEW 3 major objections 4 minor 37 references
Structure and dynamics of erupting solar prominences using the Rolling Hough Transform: Toward a feature-oriented classification
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The Rolling Hough Transform sorts solar prominences by eruption phase in plane-of-sky images.
desk verdict A promising proof-of-concept for RHT-based prominence classification, but the polar-coordinate transform has a sign error that inverts the central radial/tangential labels as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the spatial Rolling Hough Transform (RHT), adapted for solar applications by Schad (2017) from the algorithm of Clark et al. (2014), which measures the orientation of curvilinear features by counting linearly aligned pixels in a circular kernel of width wlen = 51 around each pixel, applying an adaptive threshold of f = 0.25, and outputting a mean axial angle θ̄xy and a mean resultant length R̄, with only pixels satisfying R̄ ≥ 0.75 retained. The key modification for prominences is a polar conversion, θ_polar = α − θ̄xy, where α is the angle between the line from the pixel to the solar center and the eastward axis, which turns the classification into one about radial versus tangential alignment relative to the limb. The classification itself is a weighted sum of the angle histogram: bins at ±90° count fully as radial, bins at 0° count fully as tangential, and intermediate bins count in proportion, giving a ratio R/T whose crossing of 1 defines the label.
What would settle it
Recompute the classification for both events across a grid of cadences (1–10 minutes), kernel widths (21–101), adaptive thresholds f (0.10–0.50), and R̄ cutoffs (0.5–0.9); if any phase label changes between radial and tangential under small parameter changes, the claim that the classification tracks the dynamical phase is not robust.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the spatial RHT, applied over time to 304 Å images, returns per-pixel mean orientation angles for fine-scale prominence structures, and that these angles can be converted into a polar frame so that ±90° marks radial structure and 0° marks tangential structure. Weighted histograms of these angles yield a radial-to-tangential ratio R/T, with R/T ≥ 1 classifying a phase as radial and R/T < 1 as tangential, plus a strength |R/T − 1|. In the quiet-Sun eruption, the slow-rise phase is 64.93% tangential and the fast-rise phase is 58.26% radial; in the polar-swirl activation, the rise phase is 56.96% radial and the swirl phase is 50.59% tangential, with the two full-event classifications being tangential (54.95%) and radial (53.18%) respectively. The paper interprets this as evidence that the structural orientation carries dynamical information distinct from, but consistent with, motion measurements.
Load-bearing premise
The load-bearing premise is that the radial-versus-tangential percentages produced by the RHT are stable under the analyst's choices of time cadence, kernel width, and thresholds; if those choices shift, the phase labels could flip.
Editorial extensions
If this is right
- If the central claim is correct, the RHT gives a quantitative, automated way to detect the structural transition from tangential to radial alignment that accompanies eruption onset.
- The radial/tangential ratio can serve as a plane-of-sky structural complement to motion-based distinctions between eruptive and active prominences.
- The same pipeline can be extended to the temporal RHT to estimate plasma velocities, adding a kinematic term to the classification.
- Applied to a larger sample of events, the method could produce a feature-oriented catalogue of prominence dynamics that helps constrain eruption models.
Reading between the lines
- The paper does not quantify how sensitive the radial/tangential percentages are to the cadence, kernel width, and thresholds; a sensitivity scan on these two events would show whether the phase labels are stable or artefacts of the chosen parameters.
- Because RHT detects edges and boundaries, the classification may disproportionately reflect drainage paths and fine threads rather than bulk mass; relating the percentages to quantities like mass flux or twist would test this interpretation.
- The method could be transferred to on-disk filament observations, although projection effects and foreground emission would likely require different filtering and a re-definition of the polar reference frame.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a proof-of-concept application of the spatial Rolling Hough Transform (RHT) to SDO/AIA 304 Å movies of two limb prominences: a quiet-Sun eruption and a polar-crown swirl activation. The authors compute per-pixel orientation maps, convert the RHT line orientations into a polar frame relative to the solar limb, and introduce a radial/tangential classification based on weighted angle histograms (Eqs. 4-5). They report that the quiet-Sun event is predominantly tangential in the slow-rise phase and becomes more radial during/after eruption, while the polar-swirl event starts with a strong radial contribution and develops tangential structures during activation. The paper also provides timeslice kinematic measurements, a parameter table, and a cadence-sensitivity test for the detected-pixel histogram.
Significance. If the technical issues are resolved, this would be a useful proof-of-concept for a quantitative, feature-oriented descriptor of prominence morphology, complementing existing LOS and POS methods. Strengths of the manuscript include the validation of the RHT code against Schad's TRACE case, the deliberate choice of two contrasting dynamical events, and the explicit discussion of parameter choices and of the effect of cadence on detection. However, the central radial/tangential classification is not reproducible as written because of an angle-convention inconsistency between the RHT angle definition and the polar conversion, and the phase-level percentages are presented without uncertainty or sensitivity analysis. The paper's contribution is therefore potentially sound but requires corrective revision before the classification claims can be accepted.
major comments (3)
- [§3, §4.1, Eq. (3), Section 7] The polar-angle conversion in Eq. (3) is inconsistent with the RHT angle convention defined in Section 3. Section 3 states that the RHT output θ_xy is the line orientation, with 0° horizontal and ±90° vertical. Section 4.1 then defines θ_polar = α − θ_xy and states that θ_polar = 0° is tangential and θ_polar = ±90° is radial. For a pixel on the east limb (α = 0°), a horizontal feature is radial, yet Eq. (3) gives θ_polar = 0°, which the text defines as tangential; a vertical feature is tangential, yet Eq. (3) gives θ_polar = ±90°, which the text defines as radial. The correct limb-aligned conversion from line orientation φ is θ_polar = (φ − α − 90°) mod 180 (or an equivalent form), not α − φ. As written, a reader following the text would invert every radial/tangential label, and the values in Tables 3 and 4 and the narrative in Sections 5 and 7 would be reversed relative to the stated convention. Section 7 contains a related error: it says 'Radial (parallel to the solar limb) and tangential (perpendicular to the solar limb)', which is the reverse of the standard geometry and also contradicts Section 4.1's own wording. Please correct the equation and the wording, state explicitly which angle convention the code actually outputs (line orientation or normal angle), and verify the corrected conversion against a known feature, for example a radial drainage column on the east limb.
- [§6, Tables 3-4, Figure 12] The central classification percentages are reported to two decimal places without uncertainties, and the only sensitivity test provided shows that the detected-pixel histogram changes substantially with cadence (Figure 12). Several phase classifications sit very close to the 50/50 boundary: the swirl phase has R = 49.41% / T = 50.59%, the drainage phase R = 50.72% / T = 49.28%, and the plasma-ejection and collapse phases differ by only 3-4 percentage points. No analysis is given for how the R/T ratio or the resulting labels change under the choices of cadence, wlen = 51, f = 0.25, Rbar ≥ 0.75, smr_xy, med_xy, nlev, the limb-mask threshold, or the manually selected phase boundaries. Because the phase labels are the paper's central claim, please provide a sensitivity analysis for the classification percentages and the phase-boundary choices, or at minimum report the margins by which each classification exceeds the threshold and state how those margins are affected by reasonable parameter variation.
- [§3, Eq. (1)] Equation (1), written as ¯θxy = 0.5 arctan(y/x), is not a valid expression for the directional-statistics mean of an axial distribution; x and y are image coordinates, not the RHT histogram components. This appears to be a typographical compression of the standard circular mean, e.g., ¯θxy = 0.5 atan2(Σ w_i sin 2θ_i, Σ w_i cos 2θ_i), but as written the equation is not reproducible and cannot be the operation applied to the RHT angle bins. Please replace it with the explicit formula used in the code, including the handling of the 180° axial ambiguity, so that the central measurement is unambiguously defined.
minor comments (4)
- [§4.1] The description of the modulo/offset step after Eq. (3) is ambiguous: please specify the exact mapping to the interval [-90°, 90°] (for example, θ_polar = ((θ_polar + 90°) mod 180°) − 90°), since this affects the visualizations and the histograms.
- [§2.2, §2.3] The phase names are used inconsistently: Section 2.2 refers to a 'pre-eruptive phase' while the later analysis and Table 3 use 'slow-rise'; the polar-swirl event is labeled 'Rise' in Figure 2 and 'slow-rise' in the text. Please harmonize the terminology.
- [§6, Figure 12] The cadence test in Figure 12 shows only detected-pixel histograms, not the resulting angle distributions; showing how the RHT angle histograms and classification percentages change with cadence would make the test directly relevant to the paper's main claim.
- [General] There are several typographical and formatting issues, including 'T able' in the table captions, the inconsistent spelling 'an prominence activation', and duplicated words such as 'the the' in the text; a careful proofread is needed.
Circularity Check
No significant circularity: the RHT measurement is an independent image-processing output, the radial/tangential classifier is a threshold applied after the fact, and the only self-citation is ancillary.
full rationale
The derivation chain is self-contained. The spatial RHT angles are an independent image measurement produced by the published RHT code (Clark et al. 2014; Schad 2017), and the paper validates that code by reproducing Schad's TRACE result before applying it to the two prominence events. The polar-coordinate conversion in Equation (3) is a fixed coordinate convention, not a parameter fitted to the subsequent classification, and the radial/tangential classifier in Equations (4)-(5) is a deterministic threshold applied to the measured angle histograms. Phase boundaries are selected from light-curve intensity gradients before the RHT percentages are computed, so the phase-dependent classifications are not constructed so as to reproduce themselves. The only author self-citation, Régnier et al. (2011), appears in an ancillary kinematic comparison and is not load-bearing for the RHT-based structural claim. A possible sign inconsistency in Equation (3) would be a reproducibility or correctness concern rather than a circularity: it does not make the output equivalent to its input. The acknowledged dependence on user-chosen parameters (cadence, kernel width, thresholds) and the lack of a sensitivity analysis for the classification percentages are soundness limitations, not circular steps.
Assumptions & free parameters
free parameters (9)
- wlen (kernel diameter) =
51
- f (adaptive threshold) =
0.25
- Rbar threshold =
0.75
- smr_xy (spatial filter) =
30
- med_xy (median filter) =
3
- nlev (noise level) =
10
- Timeseries cadence =
1 min (quiet-Sun), 5 min (polar full), 10 min (polar phases)
- Phase boundary times =
16:30/17:45/18:00 UT (quiet-Sun); 12:15/20:30/23:45 UT (polar)
- Limb mask threshold =
not specified numerically
assumptions (6)
- standard math The Rolling Hough Transform and directional statistics formulas from Schad (2017) and Mardia & Jupp (1999) correctly compute mean orientation and coherence.
- domain assumption 304 Å He ii emission traces cool prominence plasma, so RHT orientation in this channel reflects prominence structure.
- domain assumption Plane-of-sky projection preserves enough structural orientation for classification despite projection effects.
- ad hoc to paper Linear weighting of 5-degree angle bins and the R/T >= 1 threshold provide a meaningful radial/tangential classification.
- domain assumption Phase boundaries derived from light-curve intensity gradients correspond to distinct dynamical phases.
- ad hoc to paper The limb-edge masking threshold removes the detected solar limb without losing significant prominence structure.
Cite this review
Pith. "Pith review of Structure and dynamics of erupting solar prominences using the Rolling Hough Transform: Toward a feature-oriented classification." pith.science (2026). https://pith.science/paper/AMVV3WMJ
@misc{pith2026250524515,
author = {Pith},
title = {Pith review of: Structure and dynamics of erupting solar prominences using the Rolling Hough Transform: Toward a feature-oriented classification},
year = {2026},
howpublished = {\url{https://pith.science/paper/AMVV3WMJ}},
note = {Machine review of arXiv:2505.24515}
}
read the original abstract
The classification of solar prominences has proven to be challenging due to their diverse morphologies and dynamical behaviour. Complexity is heightened when considering eruptive prominences, where the dynamics demand methods capable of capturing detailed structural information. While there exists a range of line-of-sight (LOS) and plane-of-sky (POS) techniques which have advanced our understanding of prominence motions, they are subject to limitations, emphasising the need for effective methods of extracting structural information from prominence dynamics. We present a proof-ofconcept for the spatial Rolling Hough Transform (RHT) algorithm, which identifies finescale structural orientation in the POS, applied to prominence structure and dynamics. We demonstrate the RHT approach using two contrasting prominence dynamics events using SDO/AIA 304 \r{A} observations: (1) a quiet-Sun eruption, (2) activation (swirl) of a polar-crown prominence. By analysing the light curves and movies from each event, we divide the events into distinct dynamical phases: from slow rise to drainage. The spatial RHT method enables us to extract structural information and localised dynamics for both events and the different evolution phases. We develop a classification to label the prominences as either radially or tangentially oriented structures. The quiet-Sun eruption has a predominately tangential structure in the slow-rise phase, but displays greater radial features during/after the eruption. The polar-swirl activation initially shows a strong radial contribution, which diminishes as more tangential structures appear during/after the activation. Our results demonstrate the successful application of the spatial RHT to prominences, leading to the classification of individual prominences and an insight into their dynamics.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Aschwanden, M. J. 2010a, SoPh, 262, 399, doi: 10.1007/s11207-010-9531-6 —. 2010b, SoPh, 262, 235, doi: 10.1007/s11207-009-9474-y Asensio Ramos, A., de la Cruz Rodr ´ ıguez, J., Mart ´ ınez Gonz´ alez, M. J., & Socas-Navarro, H. 2017a, A&A, 599, A133, doi: 10.1051/0004-6361/201629755 Asensio Ramos, A., Requerey, I. S., & Vitas, N. 2017b, A&A, 604, A11, doi...
-
[2]
Berger, T. E., Liu, W., & Low, B. C. 2012, ApJL, 758, L37, doi: 10.1088/2041-8205/758/2/L37
-
[3]
Berger, T. E., Shine, R. A., Slater, G. L., et al. 2008, ApJL, 676, L89, doi: 10.1086/587171
doi:10.1086/587171 2008
-
[4]
2008, ApJL, 689, L73, doi: 10.1086/595785
Sakurai, T. 2008, ApJL, 689, L73, doi: 10.1086/595785
doi:10.1086/595785 2008
-
[5]
2008, ApJ, 689, 593, doi: 10.1086/592761
Chae, J., & Sakurai, T. 2008, ApJ, 689, 593, doi: 10.1086/592761
-
[6]
E., Tripathi, D., Isobe, H., & Asai, A
Chifor, C., Mason, H. E., Tripathi, D., Isobe, H., & Asai, A. 2006, A&A, 458, 965, doi: 10.1051/0004-6361:20065687
-
[7]
Clark, S. E., Peek, J. E. G., & Putman, M. E. 2014, ApJ, 789, 82, doi: 10.1088/0004-637X/789/1/82 S ¸ahin, S., & Antolin, P. 2022, ApJL, 931, L27, doi: 10.3847/2041-8213/ac6fe9 S ¸ahin, S., Antolin, P., Froment, C., & Schad, T. A. 2023, ApJ, 950, 171, doi: 10.3847/1538-4357/acd44b
-
[8]
1976, SoPh, 49, 283, doi: 10.1007/BF00162453
Engvold, O. 1976, SoPh, 49, 283, doi: 10.1007/BF00162453
Show all 37 references
-
[9]
R., Alexander, D., & Liu, R
Gilbert, H. R., Alexander, D., & Liu, R. 2007, SoPh, 245, 287, doi: 10.1007/s11207-007-9045-z
2007 doi
-
[10]
Hundhausen, A. J. 2000, ApJ, 537, 503, doi: 10.1086/309030
2000 doi
-
[12]
2003, ApJ, 586, 562, doi: 10.1086/367614
Gopalswamy, N., Shimojo, M., Lu, W., et al. 2003, ApJ, 586, 562, doi: 10.1086/367614
2003 doi
-
[13]
Kucera, T. A. 2015, in Astrophysics and Space Science Library, Vol. 415, Solar Prominences, ed. J.-C. Vial & O. Engvold (Springer), 79, doi: 10.1007/978-3-319-10416-4 4
2015 doi
-
[14]
C., et al
Labrosse, N., Heinzel, P., Vial, J. C., et al. 2010, SSRv, 151, 243, doi: 10.1007/s11214-010-9630-6
2010 doi
-
[15]
R., Title, A
Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, SoPh, 275, 17, doi: 10.1007/s11207-011-9776-8
2012 doi
-
[16]
E., & Low, B
Liu, W., Berger, T. E., & Low, B. C. 2012, ApJL, 745, L21, doi: 10.1088/2041-8205/745/2/L21
2012 doi
-
[17]
2010, SSRv, 151, 333, doi: 10.1007/s11214-010-9628-0
Schmieder, B., & Aulanier, G. 2010, SSRv, 151, 333, doi: 10.1007/s11214-010-9628-0
2010 doi
-
[18]
V., & Jupp, P
Mardia, K. V., & Jupp, P. E. 1999, Directional Statistics, 1st edn., Wiley Series in Probability and Statistics (New York: Wiley), doi: 10.1002/9780470316979
1999 doi
-
[19]
I., Su, Y
McCauley, P. I., Su, Y. N., Schanche, N., et al. 2015, SoPh, 290, 1703, doi: 10.1007/s11207-015-0699-7
2015 doi
-
[20]
2014, SoPh, 289, 2945, doi: 10.1007/s11207-014-0523-9
Morgan, H., & Druckm¨ uller, M. 2014, SoPh, 289, 2945, doi: 10.1007/s11207-014-0523-9
2014 doi
-
[21]
R., & Woo, R
Morgan, H., Habbal, S. R., & Woo, R. 2006, SoPh, 236, 263, doi: 10.1007/s11207-006-0113-6
2006 doi
-
[22]
J., & Simon, G
November, L. J., & Simon, G. W. 1988, ApJ, 333, 427, doi: 10.1086/166758 O’Dwyer, B., Del Zanna, G., Mason, H. E.,
1988 doi
-
[23]
A., & Tripathi, D
Weber, M. A., & Tripathi, D. 2010, A&A, 521, A21, doi: 10.1051/0004-6361/201014872
2010 doi
-
[24]
2014, Living Reviews in Solar Physics, 11, 1, doi: 10.12942/lrsp-2014-1
Parenti, S. 2014, Living Reviews in Solar Physics, 11, 1, doi: 10.12942/lrsp-2014-1
2014 doi
-
[25]
D., Thompson, B
Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, SoPh, 275, 3, doi: 10.1007/s11207-011-9841-3 R´ egnier, S., Walsh, R. W., & Alexander, C. E. 2011, A&A, 533, L1, doi: 10.1051/0004-6361/201117381 22
2012 doi
-
[26]
Robbrecht, E., Berghmans, D., & Van der Linden, R. A. M. 2009, ApJ, 691, 1222, doi: 10.1088/0004-637X/691/2/1222
2009 doi
-
[27]
2017, SoPh, 292, 132, doi: 10.1007/s11207-017-1153-9
Schad, T. 2017, SoPh, 292, 132, doi: 10.1007/s11207-017-1153-9
2017 doi
-
[28]
2010, A&A, 514, A68, doi: 10.1051/0004-6361/200913477
Schmieder, B., Chandra, R., Berlicki, A., & Mein, P. 2010, A&A, 514, A68, doi: 10.1051/0004-6361/200913477
2010 doi
-
[29]
2017, A&A, 606, A30, doi: 10.1051/0004-6361/201730839
Schmieder, B., Zapi´ or, M., L´ opez Ariste, A., et al. 2017, A&A, 606, A30, doi: 10.1051/0004-6361/201730839
2017 doi
-
[30]
J., Elmore, C., Kliem, B., T¨ or¨ ok, T., & Title, A
Schrijver, C. J., Elmore, C., Kliem, B., T¨ or¨ ok, T., & Title, A. M. 2008, ApJ, 674, 586, doi: 10.1086/524294
2008 doi
- [31]
-
[32]
1995, The nature of solar prominences, Vol
Tandberg-Hanssen, E. 1995, The nature of solar prominences, Vol. 199, doi: 10.1007/978-94-017-3396-0
1995 doi
-
[33]
T., & Brekke, P
Thompson, W. T., & Brekke, P. 2000, SoPh, 195, 45, doi: 10.1023/A:1005203001242
2000 doi
-
[34]
2018, SoPh, 293, 57, doi: 10.1007/s11207-018-1276-7
Vincent, A. 2018, SoPh, 293, 57, doi: 10.1007/s11207-018-1276-7
2018 doi
-
[35]
Srivastava, A., & Joshi, N. C. 2013, ApJ, 778, 142, doi: 10.1088/0004-637X/778/2/142
2013 doi
-
[36]
T., Abbett, W
Welsch, B. T., Abbett, W. P., De Rosa, M. L., et al. 2007, ApJ, 670, 1434, doi: 10.1086/522422
2007 doi
-
[37]
Yashiro, S., Gopalswamy, N., Akiyama, S., & M¨ akel¨ a, P. A. 2020, Journal of Atmospheric and Solar-Terrestrial Physics, 205, 105324, doi: 10.1016/j.jastp.2020.105324
2020
-
[38]
Zhang, P., Buchlin, ´E., & Vial, J. C. 2019, A&A, 624, A72, doi: 10.1051/0004-6361/201834259
2019 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.