Normalized to: Barbey, N.
[1]
oai:arXiv.org:1103.5904 [pdf] - 1053049
TomograPy: A Fast, Instrument-Independent, Solar Tomography Software
Submitted: 2011-03-30
Solar tomography has progressed rapidly in recent years thanks to the
development of robust algorithms and the availability of more powerful
computers. It can today provide crucial insights in solving issues related to
the line-of-sight integration present in the data of solar imagers and
coronagraphs. However, there remain challenges such as the increase of the
available volume of data, the handling of the temporal evolution of the
observed structures, and the heterogeneity of the data in multi-spacecraft
studies.
We present a generic software package that can perform fast tomographic
inversions that scales linearly with the number of measurements, linearly with
the length of the reconstruction cube (and not the number of voxels) and
linearly with the number of cores and can use data from different sources and
with a variety of physical models: TomograPy
(http://nbarbey.github.com/TomograPy/), an open-source software freely
available on the Python Package Index. For performance, TomograPy uses a
parallelized-projection algorithm. It relies on the World Coordinate System
standard to manage various data sources. A variety of inversion algorithms are
provided to perform the tomographic-map estimation. A test suite is provided
along with the code to ensure software quality. Since it makes use of the
Siddon algorithm it is restricted to rectangular parallelepiped voxels but the
spherical geometry of the corona can be handled through proper use of priors.
We describe the main features of the code and show three practical examples
of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1
data. Static and smoothly varying temporal evolution models are presented.
[2]
oai:arXiv.org:1103.4481 [pdf] - 416792
Morphology, dynamics and plasma parameters of plumes and inter-plume
regions in solar coronal holes
Wilhelm, K.;
Abbo, L.;
Auchere, F.;
Barbey, N.;
Feng, L.;
Gabriel, A. H.;
Giordano, S.;
Imada, S.;
Llebaria, A.;
Matthaeus, W. H.;
Poletto, G.;
Raouafi, N. -E.;
Suess, S. T.;
Teriaca, L.;
Wang, Y. -M.
Submitted: 2011-03-23
Coronal plumes, which extend from solar coronal holes (CH) into the high
corona and - possibly - into the solar wind (SW), can now continuously be
studied with modern telescopes and spectrometers on spacecraft, in addition to
investigations from the ground, in particular, during total eclipses. Despite
the large amount of data available on these prominent features and related
phenomena, many questions remained unanswered as to their generation and
relative contributions to the high-speed streams emanating from CHs. An
understanding of the processes of plume formation and evolution requires a
better knowledge of the physical conditions at the base of CHs, in plumes and
in the surrounding inter-plume regions (IPR). More specifically, information is
needed on the magnetic field configuration, the electron densities and
temperatures, effective ion temperatures, non-thermal motions, plume
cross-sections relative to the size of a CH, the plasma bulk speeds, as well as
any plume signatures in the SW. In spring 2007, the authors proposed a study on
"Structure and dynamics of coronal plumes and inter-plume regions in solar
coronal holes" to the International Space Science Institute (ISSI) in Bern to
clarify some of these aspects by considering relevant observations and the
extensive literature. This review summarizes the results and conclusions of the
study. Stereoscopic observations allowed us to include three-dimensional
reconstructions of plumes. Multi-instrument investigations carried out during
several campaigns led to progress in some areas, such as plasma densities,
temperatures, plume structure and the relation to other solar phenomena, but
not all questions could be answered concerning the details of plume generation
process(es) and interaction with the SW.
[3]
oai:arXiv.org:1012.0497 [pdf] - 1042336
Feasibility and performances of compressed-sensing and sparse map-making
with Herschel/PACS data
Submitted: 2010-12-02
The Herschel Space Observatory of ESA was launched in May 2009 and is in
operation since. From its distant orbit around L2 it needs to transmit a huge
quantity of information through a very limited bandwidth. This is especially
true for the PACS imaging camera which needs to compress its data far more than
what can be achieved with lossless compression. This is currently solved by
including lossy averaging and rounding steps on board. Recently, a new theory
called compressed-sensing emerged from the statistics community. This theory
makes use of the sparsity of natural (or astrophysical) images to optimize the
acquisition scheme of the data needed to estimate those images. Thus, it can
lead to high compression factors.
A previous article by Bobin et al. (2008) showed how the new theory could be
applied to simulated Herschel/PACS data to solve the compression requirement of
the instrument. In this article, we show that compressed-sensing theory can
indeed be successfully applied to actual Herschel/PACS data and give
significant improvements over the standard pipeline. In order to fully use the
redundancy present in the data, we perform full sky map estimation and
decompression at the same time, which cannot be done in most other compression
methods. We also demonstrate that the various artifacts affecting the data
(pink noise, glitches, whose behavior is a priori not well compatible with
compressed-sensing) can be handled as well in this new framework. Finally, we
make a comparison between the methods from the compressed-sensing scheme and
data acquired with the standard compression scheme. We discuss improvements
that can be made on ground for the creation of sky maps from the data.
[4]
oai:arXiv.org:0802.0113 [pdf] - 9663
A Time-Evolving 3D Method Dedicated to the Reconstruction of Solar
plumes and Results Using Extreme Ultra-Violet Data
Submitted: 2008-02-01
An important issue in the tomographic reconstruction of the solar poles is
the relatively rapid evolution of the polar plumes. We demonstrate that it is
possible to take into account this temporal evolution in the reconstruction.
The difficulty of this problem comes from the fact that we want a 4D
reconstruction (three spatial dimensions plus time) while we only have 3D data
(2D images plus time). To overcome this difficulty, we introduce a model that
describes polar plumes as stationary objects whose intensity varies
homogeneously with time. This assumption can be physically justified if one
accepts the stability of the magnetic structure. This model leads to a bilinear
inverse problem. We describe how to extend linear inversion methods to these
kinds of problems. Studies of simulations show the reliability of our method.
Results for SOHO/EIT data show that we are able to estimate the temporal
evolution of polar plumes in order to improve the reconstruction of the solar
poles from only one point of view. We expect further improvements from
STEREO/EUVI data when the two probes will be separated by about 60 degrees.