Sparse Astronomical Data Analysis. Jean-Luc Starck. Collaborators: J. Bobin., F. Sureau. CEA Saclay
|
|
- Bathsheba Dennis
- 5 years ago
- Views:
Transcription
1 Sparse Astronomical Data Analysis Jean-Luc Starck Collaborators: J. Bobin., F. Sureau CEA Saclay
2 What is a good representation for data? A signal s (n samples) can be represented as sum of weighted elements of a given dictionary Dictionary (basis, frame) Atoms Ex: Haar wavelet coefficients Few large coefficients Many small coefficients Sorted index k Fast calculation of the coefficients Analyze the signal through the statistical properties of the coefficients Approximation theory uses the sparsity of the coefficients 2-2
3 Sparse Model: we consider a dictionary which has a fast transform/reconstruction operator: Local DCT Stationary textures Locally oscillatory Wavelet transform Piecewise smooth Isotropic structures Curvelet transform Piecewise smooth, edge
4 NGC2997
5 The STARLET Transform Isotropic Undecimated Wavelet Transform (a trous algorithm) " = B 3 # spline, 1 2 $(x 2 ) = 1 2 "( x 2 ) #"(x) h = [1,4,6,4,1]/16, g = % - h, h = g = % I(k,l) = c J,k,l + " J w j=1 j,k,l
6 Inverse Problem Regularization using Sparsity Y = AX + N Between all possible solutions, we want the one which has the sparsest representation in the dictionary Φ. It leads to the following optimization problem: 1 min α 1,,α T 2σ 2 Y AΦα2 + λ T α i p p, 0 p<2. i=1 X = Φα A sparse model can be interpreted in a Bayesian framework Assuming the coefficients α of the solution in the dictionary Φ follow a leptokurtic PDF with heavy tails such as the generalized Gaussian distribution form: T pdf α (α 1,...,α T ) exp λ α i p p 0 p<2. i=1
7 Denoising using a sparsity model Denoising using a sparsity prior on the solution: X is sparse in Φ, i.e.x = Φα where most of α are negligible. α arg min α 1 2 Y Φα 2 +t α p p, 0 p 1.
8 α arg min α p=0 1 2 Y Φα 2 + t2 2 α 0 ==> Solution via Iterative Hard Thresholding α (t+1) =HardThresh µt ( α (t) + µφ T (Y Φ α (t) )),µ=1/ Φ 2. 1st iteration solution: X = Φ HardThresh t (Φ T Y )= Φ,t (Y ) Exact for Φ orthonormal. Soft prior Noise Modeling
9 p=1 ==> Solution via iterative Soft Thresholding α (t+1) =SoftThresh µt ( α (t) + µφ T (Y Φ α (t) )),µ (0, 2/ Φ 2 ). 1st iteration solution: Exact for Φ orthonormal. X = Φ SoftThresh t (Φ T Y )= Φ,t (Y )
10 Sparsity and FERMI data S. Digel, Stanford I. Grenier, CEA J.-L. Starck, M.J. Fadili, S. Digel, B. Zhang and J. Chiang, "Source Detection Using a 3D Sparse Representation: Application to the Fermi Gamma-ray Space Telescope ", Astronomy and Astrophysics, 504, 2, pp , J. Schmitt, J.L. Starck, J.M. Casandjian, M.J. Fadili, I. Grenier, "Poisson Denoising on the Sphere: Application to the Fermi Gamma Ray Space Telescope", Astronomy and Astrophysics, 717, A26, 2010.
11 X E LAT Data Y underlying sources counts + src positions band: 16th band: 25th very low source intensity (~10-1 or even lower) multi-spectral images FERMI dedicated wavelets: 2D-1D Wavelets ψ(x, y, z) =ψ (xy) (x, y)ψ (z) (z) w j1,j 2 [k x,k y,k z ]=D ψ (xy) j 1 ψ (z) j 2 (x, y, z)
12 Results Multi-Spectral Image Restoration (MS-VST + [2D+1D] Wavelet) J.-L. Starck, M.J. Fadili, S. Digel, B. Zhang and J. Chiang, "Source Detection Using a 3D Sparse Representation: Application to the Fermi Gamma-ray Space Telescope ", Astronomy and Astrophysics, 504, 2, pp , A 35-source grid: 5 lines x 7 columns Thresh. Level = 4σ 16th band MS-VST 25th band MS-VST
13 Inverse Problems and Iterative Thresholding Minimizing Algorithm 1 T min α 1,,α T 2σ 2 Y AΦα2 + λ α i p p, 0 p<2. i=1 Iterative thresholding with a varying threshold was proposed in (Starck et al, 2004; Elad et al, 2005) for sparse signal decomposition in order to accelerate the convergence. The idea consists in using a different threshold at each iteration. For IST: For IHT: More Refs: Vonesch et al, 2007; Elad et al 2008; Wright et al., 2008; Nesterov, 2008 and Beck-Teboulle, 2009; Blumensath, 2008; Maleki et Donoho, 2009 ; etc.
14 DECONVOLUTION SIMULATION LUCY PIXON Wavelet
15 DECONVOLUTION - E. Pantin, J.-L. Starck, and F. Murtagh, "Deconvolution and Blind Deconvolution in Astronomy", in Blind image deconvolution: theory and applications, pp , J.-L. Starck, F. Murtagh, and M. Bertero, "The Starlet Transform in Astronomical Data Processing: Application to Source Detection and Image Deconvolution", Springer, Handbook of Mathematical Methods in Imaging, in press, 2011.
16 The Curse of Missing Data - Period detection in temporal series COROT: HD Bad pixels, cosmic rays, point sources in 2D images,...
17 Inp inting M. Elad, J.-L. Starck, D.L. Donoho, P. Querre, Simultaneous Cartoon and Texture Image Inpainting using Morphological Component Analysis (MCA)", ACHA, Vol. 19, pp , M.J. Fadili, J.-L. Starck and F. Murtagh, "Inpainting and Zooming using Sparse Representations", The Computer Journal, 52, 1, pp 64-79, Where M is the mask: M(i,j) = 0 ==> missing data M(i,j) = 1 ==> good data Iterative Hard Thresholding with a decreasing threshold. MCAlab available at:
18
19 Interpolation of Missing Data: Sparse Inpainting M. Elad, J.-L. Starck, D.L. Donoho, P. Querre, Simultaneous Cartoon and Texture Image Inpainting using Morphological Component Analysis (MCA)", ACHA, Vol. 19, pp , P. Abrial, Y. Moudden, J.L. Starck, M.J. Fadili, J. Delabrouille, and M. Nguyen, "CMB Data Analysis and Sparsity", Statistical Methodology, Vol 5, No 4, pp , Where M is the mask: M(i,j) = 0 ==> missing data M(i,j) = 1 ==> good data Sparse-Inpainting preserves the ISW and the Weak Lensing signal. L. Perotto, J. Bobin, S. Plaszczynski, J.-L. Starck, and A. Lavabre, "Reconstruction of the CMB lensing for Planck", A&A, F.-X. Dupe, A. Rassat, J.-L. Starck, M. J. Fadili, Measuring the Integrated Sachs-Wolfe Effect, in press, A&A, arxiv: Theoretical justification through the sampling theory of Compressed Sensing? Rauhut and Ward, Sparse Legendre expansion via l1 minimization, Constructive Approximation journal, 2010.
20 PLANCK and Sparsity Successor of WMAP (better resolution, better sensitivity, more channels) Launched on May 14, 2009 Two instruments LFI and HFI Nine Temperature maps at 30,44,70,100,143,217,353,545,857 GHz + Polarization Angular resolutions: 33, 24, 14, 10, 7.1, 5, 5, 5, 5 ==> DATA Released in 2013
21
22 The CMB exhibits Fluctuations WMAP The Cosmic Microwave Background (CMB) is a relic radiation (with a temperature equals to Kelvin) emitted 13 billion years ago when the Universe was about years old.
23 PLANCK SIMULATED POLARIZED DATA: Magnitude Orientation 2-23
24 CMB simulated map The importance of Source Separation Extra foregrounds are superimposed with the CMB!!! Point sources, galactic foregrounds,... etc IWCS - June, 1st 2011
25 30 GHz 44 GHz 70 GHz 100 GHz 143 GHz 857 GHz 545 GHz 353 GHz 217 GHz
26 Is Sparse models adapted to PLANCK Data? Components are sparse in a wavelet dictionary Spatial Domain Wavelet Domain Dust Dust CMB CMB
27 Isotropic Undecimated Wavelet on the Sphere Wavelets, Ridgelets and Curvelets on the Sphere, Astronomy & Astrophysics, 446, , Undecimated Wavelet Transform j=1 j=2 j=3 j=4 2-27
28 Curvelet functions on the sphere
29 Dictionaries Spherical Harmonics Wavelets Curvelets 29
30 Polarized Dictionary: E/B Polarized Wavelet J.-L. Starck, Y. Moudden and J. Bobin, "Polarized Wavelets and Curvelets on the Sphere", Astronomy and Astrophysics, 497, 3, pp , j=4 2-30
31 Curvelet and E/B Mode Decomposition 2-31
32 Polarized Data Denoising Where Hard thresholding corresponds to the following non linear operation:
33 Polarized Data Denoising
34 MRS Version V2.0 available since June 2010 Wavelet, Ridgelet and Curvelet on the Sphere : Software available at: J.-L. Starck, P. Abrial, Y. Moudden and M. Nguyen, Wavelets, Ridgelets and Curvelets on the Sphere, Astronomy & Astrophysics, 446, , Wavelet transforms Continuous Wavelet Transform (Mexican Hat) Orthogonal Wavelets Undecimated isotropic wavelet transform (Spline, Meyer and Needlet filters). Pyramidal wavelet transform 2. Ridgelet and Curvelet Transforms 3. Denoising using Wavelets and Curvelets 4. Gaussianity tests: Skewness, Kurtosis, Moment of order 5 and 6, Max, Higher Criticism 5. Astrophysical Component Separation (ICA on the Sphere): JADE, Fast ICA, GMCA. 6. Sparse Inpainting. Polarized Spherical Wavelets and Curvelets: SparsePol/Version 1.0 Software available at: J.-L. Starck, Y. Moudden and J. Bobin, "Polarized Wavelets and Curvelets on the Sphere", Astronomy and Astrophysics, 497, 3, pp , 2009.
35 Semi-Blind Source Separation Standard approaches amount to model each observation as a linear mixture of elementary components (i.e. CMB, SZ, Synchrotron, Free-Free, Dust...) : i; x i = j a ij s j + n j Which can be recast as: X = AS + N Blind source separation: The objective is to estimate both A and S simultaneously!! - CMB, SZ and Free-Free emission : their electromagnetic spectrum is well known (i.e. the related columns of A are known and fixed). - Synchrotron emission : rank-1 assumption / its electromagnetic spectrum is a power law with an unknown spectral index. We have nine channels and we search for nine sources: 4 sources are modeled and 5 are not modeled. 35
36 Sparse Component Separation: the GMCA Method ==> J. Bobin, J.-L. Starck, M.J. Fadili, and Y. Moudden, "Sparsity, Morphological Diversity and Blind Source Separation", IEEE Trans. on Image Processing, Vol 16, No 11, pp , J. Bobin, J.-L. Starck, M.J. Fadili, and Y. Moudden, "Blind Source Separation: The Sparsity Revolution", Advances in Imaging and Electron Physics, Vol 152, pp , Source: S =[s 1,s n ] We define a dictionary φ Data: X =[x 1,..., x n ]=AS + N GMCA searches a sparse solution S in the dictionary subject to the constraint that the norm X AS 2 is minimal. φ A and S are estimated alternately and iteratively in two steps : 1) Estimate S assuming A is fixed (iterative thresholding) : {S} = Argmin S j λ j s j W 1 + X AS 2 F,Σ 2) Estimate A assuming S is fixed (a simple least square problem) : {A} = Argmin A X AS 2 F,Σ 36
37 Planck - WG2 - Challenger 2 comparaisons - power spectrum GMCA has the lowest spectral residuals at low l. Plot shows C_l of (reconstructed CMB input CMB) evaluated at high galactic latitudes. Sam Leach (SISSA), June 19, 2008, WG2 meeting, Munich 37
38 Component Separation Planck provides data that do not share the same resolution : GHz 70GHz 44GHz 143GHz 100GHz Planck beam Spherical harmonic
39 More formally: Component Separation: more problems i; x i = b i j Globally: X = H (AS)+N where H a ij s j + n i is the multichannel convolution operator The mixture model no more holds! is singular! Spectral behavior varies spatially for some components (dust, synchroton): H 39
40 Component Separation => Use Wavelets to work at different resolutions: Planck beam H [1 9] 1 H [2 9] 2 H [3 9] 3 H [5 9] Spherical harmonic => Assume the mixing matrix varies smoothly Partitionning of the Wavelet Scales 40
41 LOCAL GMCA Undecimated Isotropic Wavelet Transform + Block Partitioning j=1 j=2 j=3 j=4 12*32 Blocks 12*16 Blocks 12*4 Blocks 12 Blocks Global
42 Channel 1 Channel 2 Channel 9 GMCA Scale
43 Inverse problem WVD y = Kf + n Wavelet-Vaguelette GMCA Decomposition f = j f = j Kf,Ψ j,k ψ j,k with K Ψ j,k = ψ j,k k (y, Ψ j,k )ψ j,k k Multi-channel WVD Y i = j K i(a j X) i + N i β j,k = Y i, Ψ (i) j,k = A j α j,k = X s, ψ j,k X s = j k à + j Y i, Ψ (i) j,k with K i Ψ(i) j,k = ψ j,k ψ j,k s GHz The sparse GMCA solution is obtained by minimizing: min α j,a j j 1 2σ 2 β j A j α j 2 s.t. α is sparse
44 CMB map estimation Input Input CMB Map
45 CMB map estimation lgmca Estimate CMB Map
46 CMB map estimation Residual at 60min Residual map for Global GMCA
47 CMB map estimation Residual at 60min Residual map for L-GMCA
48 CMB map estimation Residual per latitude Residual per band at 60min Standard deviation in mk Angle in degrees
49 Sparsity in Astrophysics Conclusions Sparsity is very efficient for Inverse problems (denoising, deconvolution, etc). Inpainting Component Separation. Perspectives PLANCK component separation. Euclid 49
Sparsity and Morphological Diversity in Source Separation. Jérôme Bobin IRFU/SEDI-Service d Astrophysique CEA Saclay - France
Sparsity and Morphological Diversity in Source Separation Jérôme Bobin IRFU/SEDI-Service d Astrophysique CEA Saclay - France Collaborators - Yassir Moudden - CEA Saclay, France - Jean-Luc Starck - CEA
More informationCompressed Sensing in Astronomy
Compressed Sensing in Astronomy J.-L. Starck CEA, IRFU, Service d'astrophysique, France jstarck@cea.fr http://jstarck.free.fr Collaborators: J. Bobin, CEA. Introduction: Compressed Sensing (CS) Sparse
More informationImage Processing by the Curvelet Transform
Image Processing by the Curvelet Transform Jean Luc Starck Dapnia/SEDI SAP, CEA Saclay, France. jstarck@cea.fr http://jstarck.free.fr Collaborators: D.L. Donoho, Department of Statistics, Stanford E. Candès,
More informationImage Processing in Astrophysics
AIM-CEA Saclay, France Image Processing in Astrophysics Sandrine Pires sandrine.pires@cea.fr NDPI 2011 Image Processing : Goals Image processing is used once the image acquisition is done by the telescope
More informationNew Multiscale Methods for 2D and 3D Astronomical Data Set
New Multiscale Methods for 2D and 3D Astronomical Data Set Jean Luc Starck Dapnia/SEDI SAP, CEA Saclay, France. jstarck@cea.fr http://jstarck.free.fr Collaborators: D.L. Donoho, O. Levi, Department of
More informationMorphological Diversity and Source Separation
Morphological Diversity and Source Separation J. Bobin, Y. Moudden, J.-L. Starck, and M. Elad Abstract This paper describes a new method for blind source separation, adapted to the case of sources having
More informationPoisson Denoising on the Sphere: Application to the Fermi Gamma Ray Space Telescope
Astronomy & Astrophysics manuscript no. AA3corrected2 March 4, 2010 c ESO Poisson Denoising on the Sphere: Application to the Fermi Gamma Ray Space Telescope J. Schmitt 1, J.L. Starck 1, J.M. Casandjian
More informationCMB Lensing Reconstruction on PLANCK simulated data
CMB Lensing Reconstruction on PLANCK simulated data L. Perotto, J. Bobin, S. Plaszczynski, J.-L. Starck, A. Lavabre Laurence Perotto - LAL Orsay ADA V, may 7-9 Outlines CMB lensing as a cosmological tools
More informationAuthor Name. Data Mining and Knowledge Discovery
Author Name Data Mining and Knowledge Discovery 2 Contents I Data Mining and Machine Learning for Astronomical Applications 1 1 Poisson Noise Removal in Spherical Multichannel Images : Application to Fermi
More informationAstronomical Data Analysis and Sparsity: from Wavelets to Compressed Sensing
ASTRONOMICAL DATA ANALYSIS 1 Astronomical Data Analysis and Sparsity: from Wavelets to Compressed Sensing Jean-Luc Starck and Jerome Bobin, Abstract Wavelets have been used extensively for several years
More informationarxiv: v1 [astro-ph.im] 13 Jun 2012
Astronomy & Astrophysics manuscript no. articlemc June 14, 2012 c ESO Multichannel Poisson denoising and deconvolution on the sphere : Application to the Fermi Gamma Ray Space Telescope J. Schmitt 1, J.L.
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationAstronomical Data Analysis and Sparsity: from Wavelets to Compressed Sensing
ASTRONOMICAL DATA ANALYSIS 1 Astronomical Data Analysis and Sparsity: from Wavelets to Compressed Sensing Jean-Luc Starck and Jerome Bobin, arxiv:0903.3383v1 [astro-ph.im] 19 Mar 2009 Abstract Wavelets
More informationSparse recovery for spherical harmonic expansions
Rachel Ward 1 1 Courant Institute, New York University Workshop Sparsity and Cosmology, Nice May 31, 2011 Cosmic Microwave Background Radiation (CMB) map Temperature is measured as T (θ, ϕ) = k k=0 l=
More informationDES Galaxy Clusters x Planck SZ Map. ASTR 448 Kuang Wei Nov 27
DES Galaxy Clusters x Planck SZ Map ASTR 448 Kuang Wei Nov 27 Origin of Thermal Sunyaev-Zel'dovich (tsz) Effect Inverse Compton Scattering Figure Courtesy to J. Carlstrom Observables of tsz Effect Decrease
More informationIN THE blind source separation (BSS) setting, the instantaneous
2662 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 16, NO. 11, NOVEMBER 2007 Sparsity and Morphological Diversity in Blind Source Separation Jérôme Bobin, Jean-Luc Starck, Jalal Fadili, and Yassir Moudden
More informationIN a series of recent papers [1] [4], the morphological component
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 16, NO. 11, NOVEMBER 2007 2675 Morphological Component Analysis: An Adaptive Thresholding Strategy Jérôme Bobin, Jean-Luc Starck, Jalal M. Fadili, Yassir Moudden,
More informationSNIa detection in the SNLS photometric analysis using Morphological Component Analysis
Prepared for submission to JCAP SNIa detection in the SNLS photometric analysis using Morphological Component Analysis arxiv:1501.02110v1 [astro-ph.im] 5 Jan 2015 A. Möller, a,b,1 V. Ruhlmann-Kleider,
More informationarxiv: v1 [astro-ph.co] 22 Feb 2012
Astronomy & Astrophysics manuscript no. paper1 c ESO 2018 October 29, 2018 True CMB Power Spectrum Estimation P. Paykari 1, J.-L. Starck 1, and M. J. Fadili 2 arxiv:1202.4908v1 [astro-ph.co] 22 Feb 2012
More informationMorphological Diversity and Sparsity for Multichannel Data Restoration
DOI 10.1007/s10851-008-0065-6 Morphological Diversity and Sparsity for Multichannel Data Restoration J. Bobin Y. Moudden J. Fadili J.-L. Starck Springer Science+Business Media, LLC 2008 Abstract Over the
More informationCMB Foreground and separation methods. Astro 448 Yuxi Zhao 12/5/2018
CMB Foreground and separation methods Astro 448 Yuxi Zhao 12/5/2018 Outline CMB foreground overview Thermal Dust AME (spin dust, spinning nanodiamonds? ) Synchrotron Free-free CO Zodiacal light, CIB Component
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationMultiscale Geometric Analysis: Thoughts and Applications (a summary)
Multiscale Geometric Analysis: Thoughts and Applications (a summary) Anestis Antoniadis, University Joseph Fourier Assimage 2005,Chamrousse, February 2005 Classical Multiscale Analysis Wavelets: Enormous
More information3D Sparse Representations
3D Sparse Representations Lanusse F. a Starck J.-L. a Woiselle A. c Fadili M.J. b a Laboratoire AIM, UMR CEA-CNRS-Paris 7, Irfu, Service d Astrophysique, CEA Saclay, F-91191 GIF-SUR-YVETTE Cedex, France.
More informationDictionary Learning for photo-z estimation
Dictionary Learning for photo-z estimation Joana Frontera-Pons, Florent Sureau, Jérôme Bobin 5th September 2017 - Workshop Dictionary Learning on Manifolds MOTIVATION! Goal : Measure the radial positions
More informationRecent developments on sparse representation
Recent developments on sparse representation Zeng Tieyong Department of Mathematics, Hong Kong Baptist University Email: zeng@hkbu.edu.hk Hong Kong Baptist University Dec. 8, 2008 First Previous Next Last
More informationInverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France
Inverse problems and sparse models (1/2) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Structure of the tutorial Session 1: Introduction to inverse problems & sparse
More informationPoisson noise removal in multivariate count data
Poisson noise removal in multivariate count data Jalal M. Fadili, Jean-Luc Starck, B. Zhang, Seth Digel To cite this version: Jalal M. Fadili, Jean-Luc Starck, B. Zhang, Seth Digel. Poisson noise removal
More informationSparse & Redundant Signal Representation, and its Role in Image Processing
Sparse & Redundant Signal Representation, and its Role in Michael Elad The CS Department The Technion Israel Institute of technology Haifa 3000, Israel Wave 006 Wavelet and Applications Ecole Polytechnique
More information34 Starlet Transform in Astronomical Data Processing
34 Starlet Transform in Astronomical Data Processing Jean-Luc Starck Fionn Murtagh Mario Bertero 34.1 Introduction...1491 34.1.1 Source Detection...1492 34.2 Standard Approaches to Source Detection...1493
More informationMultiresolution analysis & wavelets (quick tutorial)
Multiresolution analysis & wavelets (quick tutorial) Application : image modeling André Jalobeanu Multiresolution analysis Set of closed nested subspaces of j = scale, resolution = 2 -j (dyadic wavelets)
More informationTutorial: Sparse Signal Processing Part 1: Sparse Signal Representation. Pier Luigi Dragotti Imperial College London
Tutorial: Sparse Signal Processing Part 1: Sparse Signal Representation Pier Luigi Dragotti Imperial College London Outline Part 1: Sparse Signal Representation ~90min Part 2: Sparse Sampling ~90min 2
More informationSIGNAL SEPARATION USING RE-WEIGHTED AND ADAPTIVE MORPHOLOGICAL COMPONENT ANALYSIS
TR-IIS-4-002 SIGNAL SEPARATION USING RE-WEIGHTED AND ADAPTIVE MORPHOLOGICAL COMPONENT ANALYSIS GUAN-JU PENG AND WEN-LIANG HWANG Feb. 24, 204 Technical Report No. TR-IIS-4-002 http://www.iis.sinica.edu.tw/page/library/techreport/tr204/tr4.html
More informationCMB map restoration. Abstract
1 CMB map restoration J.Bobin, J.-L. Starck, F. Sureau and, J. Fadili Abstract arxiv:1111.3149v1 [astro-ph.co] 14 Nov 2011 Estimating the cosmological microwave background is of utmost importance for cosmology.
More informationTHE PLANCK MISSION The most accurate measurement of the oldest electromagnetic radiation in the Universe
THE PLANCK MISSION The most accurate measurement of the oldest electromagnetic radiation in the Universe Rodrigo Leonardi Planck Science Office ESTEC/ESA OVERVIEW Planck observational objective & science.
More informationSparsity in Underdetermined Systems
Sparsity in Underdetermined Systems Department of Statistics Stanford University August 19, 2005 Classical Linear Regression Problem X n y p n 1 > Given predictors and response, y Xβ ε = + ε N( 0, σ 2
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationThe microwave sky as seen by Planck
The microwave sky as seen by Planck Ingunn Kathrine Wehus Jet Propulsion Laboratory California Institute of Technology on behalf of the Planck Collaboration Bayesian component separation We adopt a parametric
More informationDeconvolution of confocal microscopy images using proximal iteration and sparse representations
Deconvolution of confocal microscopy images using proximal iteration and sparse representations François-Xavier Dupé, Jalal M. Fadili, Jean-Luc Starck To cite this version: François-Xavier Dupé, Jalal
More informationSparsity and Morphological Diversity in Blind Source Separation
Sparsity and Morphological Diversity in Blind Source Separation J. Bobin, J.-L. Starck, J. Fadili and Y. Moudden Abstract Over the last few years, the development of multi-channel sensors motivated interest
More informationApplication of deconvolution to images from the EGRET gamma-ray telescope
Application of deconvolution to images from the EGRET gamma-ray telescope Symeon Charalabides, Andy Shearer, Ray Butler (National University of Ireland, Galway, Ireland) ABSTRACT The EGRET gamma-ray telescope
More informationData analysis of massive data sets a Planck example
Data analysis of massive data sets a Planck example Radek Stompor (APC) LOFAR workshop, Meudon, 29/03/06 Outline 1. Planck mission; 2. Planck data set; 3. Planck data analysis plan and challenges; 4. Planck
More informationSignal Recovery, Uncertainty Relations, and Minkowski Dimension
Signal Recovery, Uncertainty Relations, and Minkowski Dimension Helmut Bőlcskei ETH Zurich December 2013 Joint work with C. Aubel, P. Kuppinger, G. Pope, E. Riegler, D. Stotz, and C. Studer Aim of this
More informationPlanck 2014 The Microwave Sky in Temperature and Polarisation Ferrara, 1 5 December The Planck mission
Planck 2014 The Microwave Sky in Temperature and Polarisation Ferrara, 1 5 December 2014 The Planck mission Marco Bersanelli Dipartimento di Fisica, Università degli Studi di Milano Planck-LFI Deputy PI
More informationInpainting for Compressed Images
Inpainting for Compressed Images Jian-Feng Cai a, Hui Ji,b, Fuchun Shang c, Zuowei Shen b a Department of Mathematics, University of California, Los Angeles, CA 90095 b Department of Mathematics, National
More informationPlanck and Virtual Observatories: Far Infra-red / Sub-mm Specificities
JF Sygnet (Planck HFI DPCM) Planck and Virtual Observatories: Far Infra-red / Sub-mm Specificities Planck Mission, Planck Data and its Specificities (4) Data Treatment and its Specificities (4) Planck
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationWAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II
WAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II Philipp Grohs 1 and Axel Obermeier 2 October 22, 2014 1 ETH Zürich 2 ETH Zürich, supported by SNF grant 146356 OUTLINE 3. Curvelets, shearlets and parabolic
More informationCorrelations between the Cosmic Microwave Background and Infrared Galaxies
Correlations between the Cosmic Microwave Background and Infrared Galaxies Brett Scheiner & Jake McCoy Based on work by Goto, Szapudi and Granett (2012) http://cdsads.u-strasbg.fr/abs/2012mnras.422l..77g
More informationSignal Restoration with Overcomplete Wavelet Transforms: Comparison of Analysis and Synthesis Priors
Signal Restoration with Overcomplete Wavelet Transforms: Comparison of Analysis and Synthesis Priors Ivan W. Selesnick a and Mário A. T. Figueiredo b a Polytechnic Institute of New York University, Brooklyn,
More informationPolarized Foregrounds and the CMB Al Kogut/GSFC
Polarized Foregrounds and the CMB Al Kogut/GSFC The Problem In A Nutshell Galactic foregrounds are 2 orders of magnitude brighter than primordial B-mode polarization from 30 to 200 GHz Polarized vs Unpolarized
More informationCosmic Microwave Background
Cosmic Microwave Background Following recombination, photons that were coupled to the matter have had very little subsequent interaction with matter. Now observed as the cosmic microwave background. Arguably
More informationGeneralized Orthogonal Matching Pursuit- A Review and Some
Generalized Orthogonal Matching Pursuit- A Review and Some New Results Department of Electronics and Electrical Communication Engineering Indian Institute of Technology, Kharagpur, INDIA Table of Contents
More informationHunting for Dark Matter in Anisotropies of Gamma-ray Sky: Theory and First Observational Results from Fermi-LAT
Hunting for Dark Matter in Anisotropies of Gamma-ray Sky: Theory and First Observational Results from Fermi-LAT Eiichiro Komatsu (Texas Cosmology Center, Univ. of Texas at Austin) MPA Seminar, September
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationConstraints on primordial abundances and neutron life-time from CMB
Constraints on primordial abundances and neutron life-time from CMB PhD Astronomy, Astrophysics and Space Science University of Sapienza and Tor Vergata Advisor: Alessandro Melchiorri Introduction Comparison
More informationNeutrinos in the era of precision Cosmology
Neutrinos in the era of precision Cosmology Marta Spinelli Rencontres du Vietnam Quy Nhon - 21 July 2017 The vanilla model: -CDM (Late times) cosmological probes Supernovae Ia standard candles fundamental
More informationhal , version 3-23 Feb 2008
Fast Poisson Noise Removal by Biorthogonal Haar Domain Hypothesis Testing B. Zhang a, a Quantitative Image Analysis Group URA CNRS 2582 of Institut Pasteur, 75724 Paris, France M. J. Fadili b b Image Processing
More informationDenoising and Compression Using Wavelets
Denoising and Compression Using Wavelets December 15,2016 Juan Pablo Madrigal Cianci Trevor Giannini Agenda 1 Introduction Mathematical Theory Theory MATLAB s Basic Commands De-Noising: Signals De-Noising:
More informationMMSE Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm
Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm Bodduluri Asha, B. Leela kumari Abstract: It is well known that in a real world signals do not exist without noise, which may be negligible
More informationIMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND
IMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND JIAN-FENG CAI, BIN DONG, STANLEY OSHER, AND ZUOWEI SHEN Abstract. The variational techniques (e.g., the total variation based method []) are
More informationFrequency-Domain Design and Implementation of Overcomplete Rational-Dilation Wavelet Transforms
Frequency-Domain Design and Implementation of Overcomplete Rational-Dilation Wavelet Transforms Ivan Selesnick and Ilker Bayram Polytechnic Institute of New York University Brooklyn, New York 1 Rational-Dilation
More informationMCALab: Reproducible Research in Signal and Image Decomposition and Inpainting
1 MCALab: Reproducible Research in Signal and Image Decomposition and Inpainting M.J. Fadili, J.-L. Starck, M. Elad and D. Donoho Abstract Morphological Component Analysis (MCA) of signals and images is
More informationEdge preserved denoising and singularity extraction from angles gathers
Edge preserved denoising and singularity extraction from angles gathers Felix Herrmann, EOS-UBC Martijn de Hoop, CSM Joint work Geophysical inversion theory using fractional spline wavelets: ffl Jonathan
More informationarxiv: v1 [astro-ph.co] 3 Feb 2016
Astronomische Nachrichten, 1 October 2018 Excess B-modes extracted from the Planck polarization maps H. U. Nørgaard - Nielsen National Space Institute (DTU Space), Technical University of Denmark, Elektrovej,
More informationA New Poisson Noisy Image Denoising Method Based on the Anscombe Transformation
A New Poisson Noisy Image Denoising Method Based on the Anscombe Transformation Jin Quan 1, William G. Wee 1, Chia Y. Han 2, and Xuefu Zhou 1 1 School of Electronic and Computing Systems, University of
More informationDeep Learning: Approximation of Functions by Composition
Deep Learning: Approximation of Functions by Composition Zuowei Shen Department of Mathematics National University of Singapore Outline 1 A brief introduction of approximation theory 2 Deep learning: approximation
More informationMultiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing
Multiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing ibeltal F. Alem, Daniel H. Chae, and Rodney A. Kennedy The Australian National
More informationINDEPENDENT COMPONENT ANALYSIS OF THE COSMIC MICROWAVE BACKGROUND
INDEPENDENT COMPONENT ANALYSIS OF THE COSMIC MICROWAVE BACKGROUND Jean-François Cardoso CNRS / ENST - TSI, Paris France Jacques Delabrouille, Guillaume Patanchon PCC - Collège de France, Paris, France
More information2.3. Clustering or vector quantization 57
Multivariate Statistics non-negative matrix factorisation and sparse dictionary learning The PCA decomposition is by construction optimal solution to argmin A R n q,h R q p X AH 2 2 under constraint :
More informationAstronomy. Astrophysics. Spherical 3D isotropic wavelets. F. Lanusse 1, A. Rassat 2,1, and J.-L. Starck Introduction
A&A 54, A9 (1 DOI: 1.151/4-6361/1118568 c ESO 1 Astronomy & Astrophysics Spherical 3D isotropic wavelets F. Lanusse 1, A. Rassat,1, and J.-L. Starck 1 1 Laboratoire AIM, UMR CEA-CNRS-Paris 7, Irfu, SEDI-SAP,
More informationAN INTRODUCTION TO COMPRESSIVE SENSING
AN INTRODUCTION TO COMPRESSIVE SENSING Rodrigo B. Platte School of Mathematical and Statistical Sciences APM/EEE598 Reverse Engineering of Complex Dynamical Networks OUTLINE 1 INTRODUCTION 2 INCOHERENCE
More informationDirectionlets. Anisotropic Multi-directional Representation of Images with Separable Filtering. Vladan Velisavljević Deutsche Telekom, Laboratories
Directionlets Anisotropic Multi-directional Representation of Images with Separable Filtering Vladan Velisavljević Deutsche Telekom, Laboratories Google Inc. Mountain View, CA October 2006 Collaborators
More informationSPARSE SIGNAL RESTORATION. 1. Introduction
SPARSE SIGNAL RESTORATION IVAN W. SELESNICK 1. Introduction These notes describe an approach for the restoration of degraded signals using sparsity. This approach, which has become quite popular, is useful
More informationThe cosmic microwave background radiation
The cosmic microwave background radiation László Dobos Dept. of Physics of Complex Systems dobos@complex.elte.hu É 5.60 May 18, 2018. Origin of the cosmic microwave radiation Photons in the plasma are
More informationSource detection using a 3D sparse representation: application to the Fermi gamma-ray space telescope
Astronomy & Astrophysics manuscript no. aa-0 c ESO 00 August, 00 Source detection using a D sparse representation: application to the Fermi gamma-ray space telescope J.-L. Starck,J.M.Fadili,S.Digel, B.
More informationIntroduction to Sparsity in Signal Processing
1 Introduction to Sparsity in Signal Processing Ivan Selesnick Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 212 2 Under-determined linear equations Consider a system
More informationInformation Field Theory. Torsten Enßlin MPI for Astrophysics Ludwig Maximilian University Munich
Information Field Theory Torsten Enßlin MPI for Astrophysics Ludwig Maximilian University Munich Information Field Theory signal field: degrees of freedom data set: finite additional information needed
More informationLINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING
LINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING JIAN-FENG CAI, STANLEY OSHER, AND ZUOWEI SHEN Abstract. Real images usually have sparse approximations under some tight frame systems derived
More informationCosmic Microwave Background
Cosmic Microwave Background Carlo Baccigalupi,, SISSA CMB lectures at TRR33, see the complete program at darkuniverse.uni-hd.de/view/main/winterschoollecture5 These lectures are available in pdf format
More informationAn Introduction to Sparse Representations and Compressive Sensing. Part I
An Introduction to Sparse Representations and Compressive Sensing Part I Paulo Gonçalves CPE Lyon - 4ETI - Cours Semi-Optionnel Méthodes Avancées pour le Traitement des Signaux 2014 Objectifs Part I The
More informationError Analysis for H 1 Based Wavelet Interpolations
Error Analysis for H 1 Based Wavelet Interpolations Tony F. Chan Hao-Min Zhou Tie Zhou Abstract We rigorously study the error bound for the H 1 wavelet interpolation problem, which aims to recover missing
More informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationImage Decomposition and Separation Using Sparse Representations: An Overview
INVITED PAPER Image Decomposition and Separation Using Sparse Representations: An Overview This overview paper points out that signal and image processing, as well as many other important areas of engineering,
More informationThe CMB sky observed at 43 and 95 GHz with QUIET
The CMB sky observed at 43 and 95 GHz with QUIET Hans Kristian Eriksen for the QUIET collaboration University of Oslo Munich, November 28th, 2012 QUIET (Q/U Imaging ExperimenT) QUIET is a ground-based
More informationThe Galactic diffuse gamma ray emission in the energy range 30 TeV 3 PeV
The Galactic diffuse gamma ray emission in the energy range 30 TeV 3 PeV Mount Rainier by Will Christiansen Silvia Vernetto & Paolo Lipari 35th ICRC 12-20 July 2017 - Busan - South Korea Gamma ray astronomy
More informationChallenges of foreground subtraction for (primordial) CMB B-modes
Challenges of foreground subtraction for (primordial) CMB B-modes Clive Dickinson Jodrell Bank Centre for Astrophysics University of Manchester Bennett et al. (2013) CMB2013, OIST, Okinawa, Japan, 10-14
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationPrimordial gravitational waves detected? Atsushi Taruya
21 May 2014 Lunch seminar @YITP Primordial gravitational waves detected? Atsushi Taruya Contents Searching for primordial gravitational waves from cosmic microwave background polarizations Gravitational-wave
More informationCompressive Sensing (CS)
Compressive Sensing (CS) Luminita Vese & Ming Yan lvese@math.ucla.edu yanm@math.ucla.edu Department of Mathematics University of California, Los Angeles The UCLA Advanced Neuroimaging Summer Program (2014)
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationSparse linear models and denoising
Lecture notes 4 February 22, 2016 Sparse linear models and denoising 1 Introduction 1.1 Definition and motivation Finding representations of signals that allow to process them more effectively is a central
More informationConstraints on inflationary models of the Universe based on CMB data
Astrophysics Group Cavendish Laboratory Constraints on inflationary models of the Universe based on CMB data Andrés Curto Martín Enrique Martínez González 8 th RES User Conference, Santander, 23/Sept/2014
More informationDark Matter searches with radio observations
Marco Taoso Dpt. of Physics and Astronomy UBC Vancouver Dark Matter searches with radio observations IDM 2012 Chicago, 23-27 July Search for DM with astrophysical observations Gamma- rays Microwave Radio
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationIs inflation really necessary in a closed Universe? Branislav Vlahovic, Maxim Eingorn. Please see also arxiv:
Is inflation really necessary in a closed Universe? Branislav Vlahovic, Maxim Eingorn North Carolina Central University NASA University Research Centers, Durham NC Please see also arxiv:1303.3203 Chicago
More informationConvex Optimization and l 1 -minimization
Convex Optimization and l 1 -minimization Sangwoon Yun Computational Sciences Korea Institute for Advanced Study December 11, 2009 2009 NIMS Thematic Winter School Outline I. Convex Optimization II. l
More informationPower spectrum exercise
Power spectrum exercise In this exercise, we will consider different power spectra and how they relate to observations. The intention is to give you some intuition so that when you look at a microwave
More information