1,056 research outputs found
Source detection using a 3D sparse representation: application to the Fermi gamma-ray space telescope
The multiscale variance stabilization Transform (MSVST) has recently been
proposed for Poisson data denoising. This procedure, which is nonparametric, is
based on thresholding wavelet coefficients. We present in this paper an
extension of the MSVST to 3D data (in fact 2D-1D data) when the third dimension
is not a spatial dimension, but the wavelength, the energy, or the time. We
show that the MSVST can be used for detecting and characterizing astrophysical
sources of high-energy gamma rays, using realistic simulated observations with
the Large Area Telescope (LAT). The LAT was launched in June 2008 on the Fermi
Gamma-ray Space Telescope mission. The MSVST algorithm is very fast relative to
traditional likelihood model fitting, and permits efficient detection across
the time dimension and immediate estimation of spectral properties.
Astrophysical sources of gamma rays, especially active galaxies, are typically
quite variable, and our current work may lead to a reliable method to quickly
characterize the flaring properties of newly-detected sources.Comment: Accepted. Full paper will figures available at
http://jstarck.free.fr/aa08_msvst.pd
Post-Reconstruction Deconvolution of PET Images by Total Generalized Variation Regularization
Improving the quality of positron emission tomography (PET) images, affected
by low resolution and high level of noise, is a challenging task in nuclear
medicine and radiotherapy. This work proposes a restoration method, achieved
after tomographic reconstruction of the images and targeting clinical
situations where raw data are often not accessible. Based on inverse problem
methods, our contribution introduces the recently developed total generalized
variation (TGV) norm to regularize PET image deconvolution. Moreover, we
stabilize this procedure with additional image constraints such as positivity
and photometry invariance. A criterion for updating and adjusting automatically
the regularization parameter in case of Poisson noise is also presented.
Experiments are conducted on both synthetic data and real patient images.Comment: First published in the Proceedings of the 23rd European Signal
Processing Conference (EUSIPCO-2015) in 2015, published by EURASI
Wavelet transforms and their applications to MHD and plasma turbulence: a review
Wavelet analysis and compression tools are reviewed and different
applications to study MHD and plasma turbulence are presented. We introduce the
continuous and the orthogonal wavelet transform and detail several statistical
diagnostics based on the wavelet coefficients. We then show how to extract
coherent structures out of fully developed turbulent flows using wavelet-based
denoising. Finally some multiscale numerical simulation schemes using wavelets
are described. Several examples for analyzing, compressing and computing one,
two and three dimensional turbulent MHD or plasma flows are presented.Comment: Journal of Plasma Physics, 201
A Non-Local Structure Tensor Based Approach for Multicomponent Image Recovery Problems
Non-Local Total Variation (NLTV) has emerged as a useful tool in variational
methods for image recovery problems. In this paper, we extend the NLTV-based
regularization to multicomponent images by taking advantage of the Structure
Tensor (ST) resulting from the gradient of a multicomponent image. The proposed
approach allows us to penalize the non-local variations, jointly for the
different components, through various matrix norms with .
To facilitate the choice of the hyper-parameters, we adopt a constrained convex
optimization approach in which we minimize the data fidelity term subject to a
constraint involving the ST-NLTV regularization. The resulting convex
optimization problem is solved with a novel epigraphical projection method.
This formulation can be efficiently implemented thanks to the flexibility
offered by recent primal-dual proximal algorithms. Experiments are carried out
for multispectral and hyperspectral images. The results demonstrate the
interest of introducing a non-local structure tensor regularization and show
that the proposed approach leads to significant improvements in terms of
convergence speed over current state-of-the-art methods
Astronomical Data Analysis and Sparsity: from Wavelets to Compressed Sensing
Wavelets have been used extensively for several years now in astronomy for
many purposes, ranging from data filtering and deconvolution, to star and
galaxy detection or cosmic ray removal. More recent sparse representations such
ridgelets or curvelets have also been proposed for the detection of anisotropic
features such cosmic strings in the cosmic microwave background.
We review in this paper a range of methods based on sparsity that have been
proposed for astronomical data analysis. We also discuss what is the impact of
Compressed Sensing, the new sampling theory, in astronomy for collecting the
data, transferring them to the earth or reconstructing an image from incomplete
measurements.Comment: Submitted. Full paper will figures available at
http://jstarck.free.fr/IEEE09_SparseAstro.pd
- …