BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michele Piana

Michele Piana appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRmath.NAastro-ph.IMTissues and OrgansArtificial IntelligenceApplicationsComputational Engineering, Finance, and ScienceNumerical Analysisphysics.med-phastro-ph.EPastro-ph.HEComputer VisionMachine Learningmath.DSMolecular NetworksNeural and Evolutionary Computingphysics.ao-phphysics.opticsphysics.plasm-phQuantitative Methods

Catalog footprint

What is connected

25works
20topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2024arXiv

AI-FLARES: Artificial Intelligence for the Analysis of Solar Flares Data

AI-FLARES (Artificial Intelligence for the Analysis of Solar Flares Data) is a research project funded by the Agenzia Spaziale Italiana and by the Istituto Nazionale di Astrofisica within the framework of the ``Attività di Studio per la Comunità Scientifica Nazionale Sole, Sistema Solare ed Esopianeti'' program. The topic addressed by this project was the development and use of computational methods for the analysis of remote sensing space data associated to solar flare emission. This paper overviews the main results obtained by the project, with specific focus on solar flare forecasting, reconstruction of morphologies of the flaring sources, and interpretation of acceleration mechanisms triggered by solar flares.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First hard X-ray imaging results by Solar Orbiter STIX

Context. The Spectrometer/Telescope for Imaging X-rays (STIX) is one of 6 remote sensing instruments on-board Solar Orbiter. It provides hard X-ray imaging spectroscopy of solar flares by sampling the Fourier transform of the incoming flux. Aims. To show that the visibility amplitude and phase calibration of 24 out of 30 STIX sub-collimators is well advanced and that a set of imaging methods is able to provide the first hard X-ray images of the flaring Sun from Solar Orbiter. Methods. We applied four visibility-based image reconstruction methods and a count-based one to calibrated STIX observations. The resulting reconstructions are compared to those provided by an optimization algorithm used for fitting the amplitudes of STIX visibilities. Results. When applied to six flares with GOES class between C4 and M4 which occurred in May 2021, the five imaging methods produce results morphologically consistent with the ones provided by the Atmospheric Imaging Assembly on-board the Solar Dynamic Observatory (SDO/AIA) in UV wavelengths. The $χ^2$ values and the parameters of the reconstructed sources are comparable between methods, thus confirming their robustness. Conclusions. This paper shows that the current calibration of the main part of STIX sub-collimators has reached a satisfactory level for scientific data exploitation, and that the imaging algorithms already available in the STIX data analysis software provide reliable and robust reconstructions of the morphology of solar flares.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Operational solar flare forecasting via video-based deep learning

Operational flare forecasting aims at providing predictions that can be used to make decisions, typically at a daily scale, about the space weather impacts of flare occurrence. This study shows that video-based deep learning can be used for operational purposes when the training and validation sets used for the network optimization are generated while accounting for the periodicity of the solar cycle. Specifically, the paper describes an algorithm that can be applied to build up sets of active regions that are balanced according to the flare class rates associated to a specific cycle phase. These sets are used to train and validate a Long-term Recurrent Convolutional Network made of a combination of a convolutional neural network and a Long-Short Memory network. The reliability of this approach is assessed in the case of two prediction windows containing the solar storm of March 2015 and September 2017, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bad and good errors: value-weighted skill scores in deep ensemble learning

In this paper we propose a novel approach to realize forecast verification. Specifically, we introduce a strategy for assessing the severity of forecast errors based on the evidence that, on the one hand, a false alarm just anticipating an occurring event is better than one in the middle of consecutive non-occurring events, and that, on the other hand, a miss of an isolated event has a worse impact than a miss of a single event, which is part of several consecutive occurrences. Relying on this idea, we introduce a novel definition of confusion matrix and skill scores giving greater importance to the value of the prediction rather than to its quality. Then, we introduce a deep ensemble learning procedure for binary classification, in which the probabilistic outcomes of a neural network are clustered via optimization of these value-weighted skill scores. We finally show the performances of this approach in the case of three applications concerned with pollution, space weather and stock prize forecasting.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Implementation paradigm for supervised flare forecasting studies: a deep learning application with video data

Solar flare forecasting can be realized by means of the analysis of magnetic data through artificial intelligence techniques. The aim is to predict whether a magnetic active region (AR) will originate solar flares above a certain class within a certain amount of time. A crucial issue is concerned with the way the adopted machine learning method is implemented, since forecasting results strongly depend on the criterion with which training, validation, and test sets are populated. In this paper we propose a general paradigm to generate these sets in such a way that they are independent from each other and internally well-balanced in terms of AR flaring effectiveness. This set generation process provides a ground for comparison for the performance assessment of machine learning algorithms. Finally, we use this implementation paradigm in the case of a deep neural network, which takes as input videos of magnetograms recorded by the Helioseismic and Magnetic Imager on-board the Solar Dynamics Observatory (SDO/HMI). To our knowledge, this is the first time that the solar flare forecasting problem is addressed by means of a deep neural network for video classification, which does not require any a priori extraction of features from the HMI magnetograms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gain and Loss of Function mutations in biological chemical reaction networks: a mathematical model with application to colorectal cancer cells

This paper studies a system of Ordinary Differential Equations modeling a chemical reaction network and derives from it a simulation tool mimicking Loss of Function and Gain of Function mutations found in cancer cells. More specifically, from a theoretical perspective, our approach focuses on the determination of moiety conservation laws for the system and their relation with the corresponding stoichiometric surfaces. Then we show that Loss of Function mutations can be implemented in the model via modification of the initial conditions in the system, while Gain of Function mutations can be implemented by eliminating specific reactions. Finally, the model is utilized to examine in detail the G1-S phase of a colorectal cancer cell.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mathematical modelling of nuclear medicine data

Positron Emission Tomography using 2-[18F]-2deoxy-D-glucose as radiotracer (FDG-PET) is currently one of the most frequently applied functional imaging methods in clinical applications. The interpretation of FDG-PET data requires sophisticated mathematical approaches able to exploit the dynamical information contained in this kind of data. Most of these approaches are formulated within the framework of compartmental analysis, which connects the experimental nuclear data with unknown tracer coefficients measuring the effectiveness of the tracer metabolism by means of Cauchy systems of ordinary differential equations. This paper provides a coincise overview of linear compartmental methods, focusing on the analytical solution of the forward compartmental problem and on the specific issues concerning the corresponding compartmental inverse problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MEM_GE: a new maximum entropy method for image reconstruction from solar X-ray visibilities

Maximum Entropy is an image reconstruction method conceived to image a sparsely occupied field of view and therefore particularly appropriate to achieve super-resolution effects. Although widely used in image deconvolution, this method has been formulated in radio astronomy for the analysis of observations in the spatial frequency domain, and an Interactive Data Language (IDL) code has been implemented for image reconstruction from solar X-ray Fourier data. However, this code relies on a non-convex formulation of the constrained optimization problem addressed by the Maximum Entropy approach and this sometimes results in unreliable reconstructions characterized by unphysical shrinking effects. This paper introduces a new approach to Maximum Entropy based on the constrained minimization of a convex functional. In the case of observations recorded by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), the resulting code provides the same super-resolution effects of the previous algorithm, while working properly also when that code produces unphysical reconstructions. Results are also provided of testing the algorithm with synthetic data simulating observations of the Spectrometer/Telescope for Imaging X-rays (STIX) in Solar Orbiter. The new code is available in the {\em{HESSI}} folder of the Solar SoftWare (SSW)tree.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Radiomics and artificial intelligence analysis of CT data for the identification of prognostic features in multiple myeloma

Multiple Myeloma (MM) is a blood cancer implying bone marrow involvement, renal damages and osteolytic lesions. The skeleton involvement of MM is at the core of the present paper, exploiting radiomics and artificial intelligence to identify image-based biomarkers for MM. Preliminary results show that MM is associated to an extension of the intrabone volume for the whole body and that machine learning can identify CT image features mostly correlating with the disease evolution. This computational approach allows an automatic stratification of MM patients relying of these biomarkers and the formulation of a prognostic procedure for determining the disease follow-up.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the two-step estimation of the cross--power spectrum for dynamical inverse problems

We consider the problem of reconstructing the cross--power spectrum of an unobservable multivariate stochatic process from indirect measurements of a second multivariate stochastic process, related to the first one through a linear operator. In the two--step approach, one would first compute a regularized reconstruction of the unobservable signal, and then compute an estimate of its cross--power spectrum from the regularized solution. We investigate whether the optimal regularization parameter for reconstruction of the signal also gives the best estimate of the cross--power spectrum. We show that the answer depends on the regularization method, and specifically we prove that, under a white Gaussian assumption: (i) when regularizing with truncated SVD the optimal parameter is the same; (ii) when regularizing with the Tikhonov method, the optimal parameter for the cross--power spectrum is lower than half the optimal parameter for the signal. We also provide evidence that a one--step approach would likely have better mathematical properties of the two--step approach. Our results apply particularly to the brain connectivity estimation from magneto/electro-encephalographic recordings and provide a formal interpretation of recent empirical results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Expectation Maximization and the retrieval of the atmospheric extinction coefficients by inversion of Raman lidar data

We consider the problem of retrieving the aerosol extinction coefficient from Raman lidar measurements. This is an ill--posed inverse problem that needs regularization, and we propose to use the Expectation--Maximization (EM) algorithm to provide stable solutions. Indeed, EM is an iterative algorithm that imposes a positivity constraint on the solution, and provides regularization if iterations are stopped early enough. We describe the algorithm and propose a stopping criterion inspired by a statistical principle. We then discuss its properties concerning the spatial resolution. Finally, we validate the proposed approach by using both synthetic data and experimental measurements; we compare the reconstructions obtained by EM with those obtained by the Tikhonov method, by the Levenberg-Marquardt method, as well as those obtained by combining data smoothing and numerical derivation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Radon transform and the Hough transform: a unifying perspective

The Radon transform is a linear integral transform that mimics the data formation process in medical imaging modalities like X-ray Computerized Tomography and Positron Emission Tomography. The Hough transform is a pattern recognition technique, which is mainly used to detect straight lines in digital images and which has been recently extended to the automatic recognition of algebraic plane curves. Although defined in very different ways, in numerical applications both transforms ultimately take an image as an input and provide, as an output, a function defined on a parameter space. The parameters in this space describe a family of curves, which represent either the integration domains considered in the (generalized) Radon transform, or the curves to be detected by means of the Hough transform. In both cases, the 2D plot of the intensity values of the output function is the so-called (Radon or Hough) sinogram. While the Hough sinogram is produced by an algorithm whose implementation requires that the parameter space be discretized in cells, the Radon sinogram is mathematically defined on a continuous parameter space, which in turn may need to be discretized just for physical or numerical reasons. In this paper, by considering a more general and n-dimensional setting, we prove that, whether the input image is described as a set of points (possibly with different intensity values) or as a piecewise constant function, its (rescaled) Hough sinogram converges to the corresponding Radon sinogram as the discretization step in the parameter space tends to zero. We also show that this result may have a notable impact on the image reconstruction problem of inverting the Radon sinogram recorded by a medical imaging scanner, and that the description of the Hough transform problem within the framework of regularization theory for inverse problems is worth investigating.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

DESAT: an SSW tool for SDO/AIA image de-saturation

Saturation affects a significant rate of images recorded by the Atmospheric Imaging Assembly on the Solar Dynamics Observatory. This paper describes a computational method and a technological pipeline for the de-saturation of such images, based on several mathematical ingredients like Expectation Maximization, image correlation and interpolation. An analysis of the computational properties and demands of the pipeline, together with an assessment of its reliability are performed against a set of data recorded from the Feburary 25 2014 flaring event.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse diffraction for the Atmospheric Imaging Assembly in the Solar Dynamics Observatory

The Atmospheric Imaging Assembly in the Solar Dynamics Observatory provides full Sun images every 1 seconds in each of 7 Extreme Ultraviolet passbands. However, for a significant amount of these images, saturation affects their most intense core, preventing scientists from a full exploitation of their physical meaning. In this paper we describe a mathematical and automatic procedure for the recovery of information in the primary saturation region based on a correlation/inversion analysis of the diffraction pattern associated to the telescope observations. Further, we suggest an interpolation-based method for determining the image background that allows the recovery of information also in the region of secondary saturation (blooming).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Systematic Approach to the Reconstruction of Saturated SDO/AIA Images

EUV images of solar flares provided by the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory (SDO/AIA) are often affected by saturation effects in their core, physically most interesting region. We introduce an image reconstruction procedure that allows recovering information in the primary saturation domain using the secondary images produced by the diffraction fringes as input data. Such a procedure is based on standard image-processing tools like correlation, convolution, and back-projection. Its effectiveness is tested in the case of SDO/AIA observations of the July 8 2013 flaring event.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sliding window and compressive sensing for low-field dynamic magnetic resonance imaging

We describe an acquisition/processing procedure for image reconstruction in dynamic Magnetic Resonance Imaging (MRI). The approach requires sliding window to record a set of trajectories in the k-space, standard regularization to reconstruct an estimate of the object and compressed sensing to recover image residuals. We validated this approach in the case of specific simulated experiments and, in the case of real measurements, we showed that the procedure is reliable even in the case of data acquired by means of a low-field scanner.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The process of data formation for the Spectrometer/Telescope for Imaging X-rays (STIX) in Solar Orbiter

The Spectrometer/Telescope for Imaging X-rays (STIX) is a hard X-ray imaging spectroscopy device to be mounted in the Solar Orbiter cluster with the aim of providing images and spectra of solar flaring regions at different photon energies in the range from a few keV to around 150 keV. The imaging modality of this telescope is based on the Moire pattern concept and utilizes 30 sub-collimators, each one containing a pair of co-axial grids. This paper applies Fourier analysis to provide the first rigorous description of the data formation process in STIX. Specifically, we show that, under first harmonic approximation, the integrated counts measured by STIX sub-collimators can be interpreted as specific spatial Fourier components of the incoming photon flux, named visibilities. Fourier analysis also allows the quantitative assessment of the reliability of such interpretation. The description of STIX data in terms of visibilities has a notable impact on the image reconstruction process, since it fosters the application of Fourier-based imaging algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Compartmental analysis of nuclear imaging data for the quantification of FDG liver metabolism

This paper utilizes compartmental analysis and a statistical optimization technique in order to reduce a compartmental model describing the metabolism of labelled glucose in liver. Specifically, we first design a compartmental model for the gut providing as output the tracer concentration in the portal vein. This quantity is then used as one of the two input functions in a compartmental model for the liver. This model, in turn, provides as output the tracer coefficients quantitatively describing the effectiveness with which the labelled glucose is transported between the different compartments. For both models, the computation of the solutions for the inverse problems is performed by means of an Ant Colony Optimization algorithm. The validation of the whole process is realized by means of synthetic data simulated by solving the forward problem of the compartmental system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Expectation Maximization for Hard X-ray Count Modulation Profiles

This paper is concerned with the image reconstruction problem when the measured data are solar hard X-ray modulation profiles obtained from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI)} instrument. Our goal is to demonstrate that a statistical iterative method classically applied to the image deconvolution problem is very effective when utilized for the analysis of count modulation profiles in solar hard X-ray imaging based on Rotating Modulation Collimators. The algorithm described in this paper solves the maximum likelihood problem iteratively and encoding a positivity constraint into the iterative optimization scheme. The result is therefore a classical Expectation Maximization method this time applied not to an image deconvolution problem but to image reconstruction from count modulation profiles. The technical reason that makes our implementation particularly effective in this application is the use of a very reliable stopping rule which is able to regularize the solution providing, at the same time, a very satisfactory Cash-statistic (C-statistic). The method is applied to both reproduce synthetic flaring configurations and reconstruct images from experimental data corresponding to three real events. In this second case, the performance of Expectation Maximization, when compared to Pixon image reconstruction, shows a comparable accuracy and a notably reduced computational burden; when compared to CLEAN, shows a better fidelity with respect to the measurements with a comparable computational effectiveness. If optimally stopped, Expectation Maximization represents a very reliable method for image reconstruction in the RHESSI context when count modulation profiles are used as input data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Return currents and energy transport in the solar flaring atmosphere

According to a standard ohmic perspective, the injection of accelerated electrons into the flaring region violates local charge equilibrium and therefore, in response, return currents are driven by an electric field to equilibrate such charge violation. In this framework, the energy loss rate associated to these local currents has an ohmic nature and significantly shortens the acceleration electron path. In the present paper we adopt a different viewpoint and, specifically, we study the impact of the background drift velocity on the energy loss rate of accelerated electrons in solar flares. We first utilize the Rutherford cross-section to derive the formula of the energy loss rate when the collisional target has a finite temperature and the background instantaneously and coherently moves up to equilibrate the electron injection. We then use the continuity equation for electrons and imaging spectroscopy data provided by RHESSI to validate this model. Specifically, we show that this new formula for the energy loss rate provides a better fit of the experimental data with respect to the model based on the effects of standard ohmic return currents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Specific Acceleration Rate in Loop-structured Solar Flares -- Implications for Electron Acceleration Models

We analyze electron flux maps based on RHESSI hard X-ray imaging spectroscopy data for a number of extended coronal loop flare events. For each event, we determine the variation of the characteristic loop length $L$ with electron energy $E$, and we fit this observed behavior with models that incorporate an extended acceleration region and an exterior "propagation" region, and which may include collisional modification of the accelerated electron spectrum inside the acceleration region. The models are characterized by two parameters: the plasma density $n$ in, and the longitudinal extent $L_0$ of, the acceleration region. Determination of the best-fit values of these parameters permits inference of the volume that encompasses the acceleration region and of the total number of particles within it. It is then straightforward to compute values for the emission filling factor and for the {\it specific acceleration rate} (electrons s$^{-1}$ per ambient electron above a chosen reference energy). For the 24 events studied, the range of inferred filling factors is consistent with a value of unity. The inferred mean value of the specific acceleration rate above $E_0=20$ keV is $\sim10^{-2}$ s$^{-1}$, with a 1$σ$ spread of about a half-order-of-magnitude above and below this value. We compare these values with the predictions of several models, including acceleration by large-scale, weak (sub-Dreicer) fields, by strong (super-Dreicer) electric fields in a reconnecting current sheet, and by stochastic acceleration processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Determination of the Acceleration Region Size in a Loop-structured Solar Flare

In order to study the acceleration and propagation of bremsstrahlung-producing electrons in solar flares, we analyze the evolution of the flare loop size with respect to energy at a variety of times. A GOES M3.7 loop-structured flare starting around 23:55 on 2002 April 14 is studied in detail using \textit{Ramaty High Energy Solar Spectroscopic Imager} (\textit{RHESSI}) observations. We construct photon and mean-electron-flux maps in 2-keV energy bins by processing observationally-deduced photon and electron visibilities, respectively, through several image-processing methods: a visibility-based forward-fit (FWD) algorithm, a maximum entropy (MEM) procedure and the uv-smooth (UVS) approach. We estimate the sizes of elongated flares (i.e., the length and width of flaring loops) by calculating the second normalized moments of the intensity in any given map. Employing a collisional model with an extended acceleration region, we fit the loop lengths as a function of energy in both the photon and electron domains. The resulting fitting parameters allow us to estimate the extent of the acceleration region which is between $\sim 13 \rm{arcsec}$ and $\sim 19 \rm{arcsec}$. Both forward-fit and uv-smooth algorithms provide substantially similar results with a systematically better fit in the electron domain.The consistency of the estimates from these methods provides strong support that the model can reliably determine geometric parameters of the acceleration region. The acceleration region is estimated to be a substantial fraction ($\sim 1/2$) of the loop extent, indicating that this dense flaring loop incorporates both acceleration and transport of electrons, with concurrent thick-target bremsstrahlung emission.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Properties of the Acceleration Regions in Several Loop-structured Solar Flares

Using {\em RHESSI} hard X-ray imaging spectroscopy observations, we analyze electron flux maps for a number of extended coronal loop flares. For each event, we fit a collisional model with an extended acceleration region to the observed variation of loop length with electron energy $E$, resulting in estimates of the plasma density in, and longitudinal extent of, the acceleration region. These quantities in turn allow inference of the number of particles within the acceleration region and hence the filling factor $f$ -- the ratio of the emitting volume to the volume that encompasses the emitting region(s). We obtain values of $f$ that lie mostly between 0.1 and 1.0; the (geometric) mean value is $f = 0.20 \times ÷3.9$, somewhat less than, but nevertheless consistent with, unity. Further, coupling information on the number of particles in the acceleration region with information on the total rate of acceleration of particles above a certain reference energy (obtained from spatially-integrated hard X-ray data) also allows inference of the specific acceleration rate (electron s$^{-1}$ per ambient electron above the chosen reference energy). We obtain a (geometric) mean value of the specific acceleration rate $η(20$ keV) $ = (6.0 \times / ÷3.4) \times 10^{-3}$ electrons s$^{-1}$ per ambient electron; this value has implications both for the global electrodynamics associated with replenishment of the acceleration region and for the nature of the particle acceleration process.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Regularization of constrained maximum likelihood iterative algorithms by means of statistical stopping rule

In this paper we propose a new statistical stopping rule for constrained maximum likelihood iterative algorithms applied to ill-posed inverse problems. To this aim we extend the definition of Tikhonov regularization in a statistical framework and prove that the application of the proposed stopping rule to the Iterative Space Reconstruction Algorithm (ISRA) in the Gaussian case and Expectation Maximization (EM) in the Poisson case leads to well defined regularization methods according to the given definition. We also prove that, if an inverse problem is genuinely ill-posed in the sense of Tikhonov, the same definition is not satisfied when ISRA and EM are optimized by classical stopping rule like Morozov's discrepancy principle, Pearson's test and Poisson discrepancy principle. The stopping rule is illustrated in the case of image reconstruction from data recorded by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). First, by using a simulated image consisting of structures analogous to those of a real solar flare we validate the fidelity and accuracy with which the proposed stopping rule recovers the input image. Second, the robustness of the method is compared with the other classical stopping rules and its advantages are shown in the case of real data recorded by RHESSI during two different flaring events.

0 citations0 reviewsNo code linkedOpen access