IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM
|
|
- Regina Waters
- 6 years ago
- Views:
Transcription
1 IMAGE RECONSTRUCTION FROM UNDERSAMPLED FOURIER DATA USING THE POLYNOMIAL ANNIHILATION TRANSFORM Anne Gelb, Rodrigo B. Platte, Rosie Renaut Development and Analysis of Non-Classical Numerical Approximation Tools AFOSR ( ) Developing Fast, Accurate, and Robust Numerical Algorithms for Extracting Actionable Information from Acquired Sensing Data AFOSR ( ) AFOSR Computational Mathematics Program Review Meeting (August 215)
2 PROJECT: DEVELOPMENT AND ANALYSIS OF NON-CLASSICAL NUMERICAL APPROXIMATION TOOLS THEME: Sampling (direct and indirect) and function (or feature) reconstruction. Consider {(T N [f (x)]) j : j = 1,..., N} is given. Design P M, such that P M (T N [f ])(x) g(f (x)). RELATED QUESTIONS: Restraints on the collected data, {(T N [f ](x)) j : j = 1,..., N}? Convergence rates and stability? What regularization techniques can be applied to ensure numerical accuracy and stability?
3 APPLICATIONS RADAR IMAGING Data from synthetic aperture radar can be seen as incomplete Fourier data.
4 APPLICATIONS RADAR IMAGING Data from synthetic aperture radar can be seen as incomplete Fourier data. Medical imaging: MRI and tomography
5 APPLICATIONS RADAR IMAGING Data from synthetic aperture radar can be seen as incomplete Fourier data. Medical imaging: MRI and tomography Solution of PDEs on irregular domains
6 Function Reconstruction from Undersampled Fourier Data Using the Polynomial Annihilation Transform In collaboration with Richard Archibald (Oak Ridge)
7 THE FOURIER RECONSTRUCTION PROBLEM Let Given 1 ˆf (ω) = f (x)e iπω x dx. 2 R d ˆf (ωj ), j =,..., N, We want to recover f (x j ), j =,... M. DIFFICULITES: 1 If f is not smooth (Gibbs phenomenon). 2 ˆf is not sampled on a structured grid (no FFT). 3 The data may be undersampled (N < M). Remark: NUFFT addresses 2 to some extent, but not 1 and 3. Regularization is needed to address lack of uniqueness.
8 1D EXAMPLES EXAMPLE f a (x) = { 1 x if 1 x < 1 x otherwise ; f b(x) = cos πx 2 if 1 x < 1 2 cos 3πx 2 if 1 2 x < 1 2 cos 7πx 2 if 1 2 x 1 Test Function 1 Uniform Grid Test Function Uniform Grid data points are shown
9 FOURIER DATA, 1D EXAMPLE 1 Target function.5 f(x) x.3 Fourier data.15 R(ˆf(ω)) I(ˆf(ω)) ω ω Wave numbers ω j were drawn from a normal distribution with standard deviation 2N/6, where 2N + 1 is the number of recovered function values.
10 FOURIER RECONSTRUCTION FROM PARTIAL SUM, 1D EXAMPLE 1 Test Function Uniform Grid 1 Fourier Reconstruction Test Function 1 Filtered Fourier Reconstruction Test Function f a (x) S N f a (x) S σ N f a(x) Test Function 1 Fourier Reconstruction 1 Filtered Fourier Reconstruction Uniform Grid Test Function Test Function f b (x) S N f b (x) S σ N f b(x) S N f (x) = N ˆf (ωj )e iπωjx and SNf σ (x) = j= N σ jˆf (ωj )e iπωjx j=
11 l 1 REGULARIZATION [Candès, Tao, Romberg, Donoho, Lustig, Goldstein, Osher,...] min f J(f) such that Ff ˆf 2 =, ˆf consists of samples of the Fourier transform of the target function/image, f. F Is the Fourier transform operator (matrix). J is an appropriate l 1 regularization term, (e.g. Lf 1 ). Typically for measured data the related (TV) denoising problem, min f J(f) such that Ff ˆf 2 < σ,
12 THE POLYNOMIAL ANNIHILATION TRANSFORM [Archibald, Gelb, Yoon, 25] Let [f ](x) = f (x + ) f (x ) be the jump function in 1D. Edge map: f [f ]. POLYNOMIAL ANNIHILATION TRANSFORM L m f (x) = 1 c j (x)f (x j ), q m (x) x j S x S Ij = {x j m,, x j+ m }, S I 2 2 j = {x j m+1,, x 2 j+ m 1 } 2 for m even and odd respectively (centered case). c j (x)p l (x j ) = p (m) l (x), j = 1,..., m + 1, x j S x q m (x) = x j S + x c j (x), where S + x is the set of points x j S x such that x j x.
13 THE POLYNOMIAL ANNIHILATION TRANSFORM [Archibald, Gelb, Yoon, 25] Let [f ](x) = f (x + ) f (x ) be the jump function in 1D. Edge map: f [f ]. POLYNOMIAL ANNIHILATION TRANSFORM For example, not assuming periodicity, for m = 4, L m =
14 FOURIER RECONSTRUCTION USING THE PA TRANSFORM IN 1D min f L m f 1 such that Ff ˆf 2 =, 1 PA m=1 Test Function 1 PA m=2 Test Function 1 PA m=3 Test Function PA (m = 1) for f a PA (m = 2) for f a PA (m = 3) for f a PA m=1 1 PA m=2 1 PA m=3 Test Function Test Function Test Function PA (m = 1) for f b PA (m = 2) for f b PA (m = 3) for f b
15 PHASE TRANSITION DIAGRAMS Recovery of piecewise constant functions m = 1 (TV) m = 2 m = 3 The colormap shows the fraction of successful recoveries in 2 trials. A trial is deemed successful when the relative l 2 error is below 1 2.
16 PHASE TRANSITION DIAGRAMS Recovery of piecewise linear functions m = 1 (TV) m = 2 m = 3 The colormap shows the fraction of successful recoveries in 2 trials. A trial is deemed successful when the relative l 2 error is below 1 2.
17 PHASE TRANSITION DIAGRAMS Recovery of piecewise quadratic polynomials m = 1 (TV) m = 2 m = 3 The colormap shows the fraction of successful recoveries in 2 trials. A trial is deemed successful when the relative l 2 error is below 1 2.
18 THE SPLIT BREGMAN ALGORITHM [Goldstein and Osher, 29] min f L m f 1 such that Ff ˆf 2 < σ, ALGORITHM While Ff ˆf 2 > σ 1 f k+1 = arg min Ff ˆf 2 + λ d k L m f b k 2 2 f 2 d k+1 = arg min d 1 + λ d L m f k+1 b k 2 2 d 3 b k+1 = b k + L m f k+1 d k+1 Remarks: (1) is a linear least-squares problem (2) is solved using the shrink operator (3) is a feedback term (the residual is added back into the regularization)
19 2D EXAMPLES Masks used for sampling in Fourier space Gaussian Sampling Tomographic Sampling
20 PA (m = 2) PA (m = 2) 2D EXAMPLE - NO NOISE, GAUSSIAN SAMPLE (5%) TV TGVSH
21 2D EXAMPLE: CROSS SECTION ERROR COMPARISON Fourier TV 6 TGVSHCS Spa m=2 Spa m=
22 Noisy Fourier data (1 db SNR) and using tomographic sampling (5%) TV TGVSH PA (m = 2) PA (m = 3)
23 Noisy Fourier data (5 db SNR) and using Gaussian sampling (5%) TV TGVSH PA (m = 2) PA (m = 3)
24 SAR image: Noisy Fourier data (5 db SNR) undersampled by 5% Original SAR Image TV TGVSH Rel. l error =.14 Rel. l2 error =.13 2 PA m = 2, Rel. l2 error =.96 PA m = 3, Rel. l2 error =.9
25 REMARKS The Polynomial Annihilation transform improves accuracy away from discontinuities. The Split Bregman Algorithm can be adapted to use the PA transform. The method is especially effective when data are undersampled and noisy. The computational cost is comparable to TV regularized recovery.
26 A frame theoretic approach to the nonuniform fast Fourier transform (Gelb and Song, 214) Frame operator Frame approximation S(f )(x) = j ˆf (ωj )e iωjx f (x) = j ˆf (ωj )S 1 e iωjx Change of Frames S 1 e iωjx l k b l,j e ijx New frame approximation A frm (f ) = ˆf (ωj )b l,j e ijx Convolutional gridding A cg (f )(x) = j n l m j n l ω j q α jˆf (ωj )ĝ(l ω j) e ilx g(x)
27 A frame theoretic approach to the nonuniform fast Fourier transform (Gelb and Song, 214)
28 Fourier reconstruction of univariate piecewise-smooth functions from non-uniform spectral data with exponential convergence rates (Platte, Gutierrez and Gelb, 215)
29 y 3 A windowed Fourier method for approximation of non-periodic functions on equispaced nodes (Platte 215) u uw 1 w x
30 y 3 A windowed Fourier method for approximation of non-periodic functions on equispaced nodes (Platte 215) u uw 1 w x
31 STARTING GRANT ( ) Developing Fast, Accurate, and Robust Numerical Algorithms for Extracting Actionable Information from Acquired Sensing Data AFOSR Supported by: Dr. Arje Nachman Sensing, Surveillance and Navigation Program and Dr. Jean-Luc Cambier Computational Mathematics Program
32 OBJECTIVES Construct fast, accurate, and robust numerical algorithms that extract actionable information from acquired sensing data in collaboration with Air Force laboratories. Study synthetic aperture radar (SAR) as the prototype sensing system, with the focus on the detection, classification, tracking of targets, and the extraction of important image information. Thoroughly analyze current localized feature detection and image reconstruction methods so that it is possible to develop new algorithms that are effective under non-idyllic circumstances (flight path variations, clutter, jamming). Expand current localized feature detection and imaging techniques, when possible, to a framework amenable to other relevant representations (such as wavelets and short time Fourier transforms). Develop high-order algorithms for solving partial differential equations (PDEs) where non-standard algorithms may be more suitable for computations.
33 WINDOWED FOURIER METHODS FOR PDES Combining the benefits of Fourier spectral methods with domain decomposition (partition of unit) schemes WF Fourier WFadapt RBF ε =1 RBF ε =2 RBF ε =1 y.5.5 error (log scale) x degrees of freedom
34 WINDOWED FOURIER METHODS FOR PDES Combining the benefits of Fourier spectral methods with domain decomposition (partition of unit) schemes.
35 PUBLICATIONS D. Denker, R. Archibald and A. Gelb, (215) An Adaptive Fourier Filter for Relaxing Time Stepping Constraints for Explicit Solvers, Lecture Notes in Computational Science and Engineering (ICOSAHOM Proceedings), Springer, to appear. R. Archibald, A. Gelb, and R.B. Platte (215) Image Reconstruction from Undersampled Fourier Data Using the Polynomial Annihilation Transform, J.Sci. Comput., in revision. R.B. Platte, A. Gutierrez, and A. Gelb, (215) Edge informed Fourier reconstruction from non-uniform spectral data with exponential convergence rates, Appl. Comput. Harmon. Anal., DOI:1.116/j.acha G. Wasserman, R. Archibald and A. Gelb,(214) Image Reconstruction from Fourier Data Using Sparsity of Edges, J. Sci. Comput., DOI 1.17/s A. Gelb and G. Song, (214) A Frame Theoretic Approach to the Non-Uniform Fast Fourier Transform, SIAM J.Numer. Anal, 52: A. Martinez, A. Gelb and A. Gutierrez, (214) Edge Detection from Non-Uniform Fourier Data Using the Convolutional Gridding Algorithm, J. Sci. Comput, 61: G. Song and A. Gelb, (213) Approximating the Inverse Frame Operator from Localized Frames, Appl. Comput. Harmon. Anal., 35:1, Q. Huang, R. Eubank and R. A. Renaut, (215), Functional Partial Canonical Correlation, Bernoulli, 21:2, S. Vatankhah, V. E Ardestani and R. A Renaut, (214) Automatic estimation of the regularization parameter in 2-D focusing gravity inversion: an application to the Safo manganese mine in northwest of Iran, J. Geophysics and Engineering, 11:4,
36 PUBLICATIONS, CONT. S. Vatankhah, R. A. Renaut and V. E. Ardestani, (214), Regularization Parameter Estimation for Underdetermined problems by the χ 2 principle with application to 2D focusing gravity inversion, Inverse Problems,3: J. Hansen, J. Hogue, G. Sander, R. A. Renaut, and S. C. Popat (214), Non-negatively constrained least squares and parameter choice by the residual periodogram for the inversion of electrochemical impedance spectroscopy, submitted. R.A. Renaut, R. Baker, M. Horst, C. Johnson and D. Nasir (213) Stability and error analysis of the polarization estimation inverse problem for microbial fuel cells, Inverse Problems 29: R.B. Platte (215) C compactly supported and positive definite radial kernels, SIAM J. Sci. Comput. (to appear). R.B. Platte (215) A windowed Fourier method for approximation of non-periodic functions on equispaced nodes, Lecture Notes in Computational Science and Engineering (ICOSAHOM Proceedings), Springer (to appear). A.A. Mitrano, R.B. Platte (215) A numerical study of divergence-free kernel approximations, Appl. Numer. Math., 96, R.B. Platte, A.J. Gutierrez, A. Gelb (214) Fourier reconstruction of univariate piecewise-smooth functions from non-uniform spectral data with exponential convergence rates, Appl. Comput. Harmon. Anal., DOI 1.116/j.acha J.M. Martel, R.B. Platte (215) Stability of radial basis function methods for convection problems on the circle and sphere, J. Sci. Comput. (submitted). B. Adcock, R.B. Platte (215) A mapped polynomial method for high-accuracy approximations on arbitrary grids, SIAM J. Numer. Anal. (submitted). Thank you!
Numerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationTowards a Mathematical Theory of Super-resolution
Towards a Mathematical Theory of Super-resolution Carlos Fernandez-Granda www.stanford.edu/~cfgranda/ Information Theory Forum, Information Systems Laboratory, Stanford 10/18/2013 Acknowledgements This
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
More informationOn the coherence barrier and analogue problems in compressed sensing
On the coherence barrier and analogue problems in compressed sensing Clarice Poon University of Cambridge June 1, 2017 Joint work with: Ben Adcock Anders Hansen Bogdan Roman (Simon Fraser) (Cambridge)
More informationDirect Methods for Reconstruction of Functions and their Edges from Non-Uniform Fourier Data
Direct Methods for Reconstruction of Functions and their Edges from Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu ICERM Research Cluster Computational Challenges in Sparse and Redundant
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
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 informationGoing off the grid. Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison
Going off the grid Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison Joint work with Badri Bhaskar Parikshit Shah Gonnguo Tang We live in a continuous world... But we work
More informationIntroduction to Compressed Sensing
Introduction to Compressed Sensing Alejandro Parada, Gonzalo Arce University of Delaware August 25, 2016 Motivation: Classical Sampling 1 Motivation: Classical Sampling Issues Some applications Radar Spectral
More informationBeyond incoherence and beyond sparsity: compressed sensing in the real world
Beyond incoherence and beyond sparsity: compressed sensing in the real world Clarice Poon 1st November 2013 University of Cambridge, UK Applied Functional and Harmonic Analysis Group Head of Group Anders
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationGREEDY SIGNAL RECOVERY REVIEW
GREEDY SIGNAL RECOVERY REVIEW DEANNA NEEDELL, JOEL A. TROPP, ROMAN VERSHYNIN Abstract. The two major approaches to sparse recovery are L 1-minimization and greedy methods. Recently, Needell and Vershynin
More informationFourier Reconstruction from Non-Uniform Spectral Data
School of Electrical, Computer and Energy Engineering, Arizona State University aditya.v@asu.edu With Profs. Anne Gelb, Doug Cochran and Rosemary Renaut Research supported in part by National Science Foundation
More informationCompressed sensing. Or: the equation Ax = b, revisited. Terence Tao. Mahler Lecture Series. University of California, Los Angeles
Or: the equation Ax = b, revisited University of California, Los Angeles Mahler Lecture Series Acquiring signals Many types of real-world signals (e.g. sound, images, video) can be viewed as an n-dimensional
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationFourier reconstruction of univariate piecewise-smooth functions from non-uniform spectral data with exponential convergence rates
Fourier reconstruction of univariate piecewise-smooth functions from non-uniform spectral data with exponential convergence rates Rodrigo B. Platte a,, Alexander J. Gutierrez b, Anne Gelb a a School of
More informationExponential decay of reconstruction error from binary measurements of sparse signals
Exponential decay of reconstruction error from binary measurements of sparse signals Deanna Needell Joint work with R. Baraniuk, S. Foucart, Y. Plan, and M. Wootters Outline Introduction Mathematical Formulation
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 informationRecovery of Compactly Supported Functions from Spectrogram Measurements via Lifting
Recovery of Compactly Supported Functions from Spectrogram Measurements via Lifting Mark Iwen markiwen@math.msu.edu 2017 Friday, July 7 th, 2017 Joint work with... Sami Merhi (Michigan State University)
More informationRobust Uncertainty Principles: Exact Signal Reconstruction from Highly Incomplete Frequency Information
1 Robust Uncertainty Principles: Exact Signal Reconstruction from Highly Incomplete Frequency Information Emmanuel Candès, California Institute of Technology International Conference on Computational Harmonic
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 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 informationA Wavelet-based Method for Overcoming the Gibbs Phenomenon
AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 A Wavelet-based Method for Overcoming the Gibbs Phenomenon NATANIEL GREENE Department of Mathematics Computer
More informationSupport Detection in Super-resolution
Support Detection in Super-resolution Carlos Fernandez-Granda (Stanford University) 7/2/2013 SampTA 2013 Support Detection in Super-resolution C. Fernandez-Granda 1 / 25 Acknowledgements This work was
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationSignal Recovery From Incomplete and Inaccurate Measurements via Regularized Orthogonal Matching Pursuit
Signal Recovery From Incomplete and Inaccurate Measurements via Regularized Orthogonal Matching Pursuit Deanna Needell and Roman Vershynin Abstract We demonstrate a simple greedy algorithm that can reliably
More informationRui ZHANG Song LI. Department of Mathematics, Zhejiang University, Hangzhou , P. R. China
Acta Mathematica Sinica, English Series May, 015, Vol. 31, No. 5, pp. 755 766 Published online: April 15, 015 DOI: 10.1007/s10114-015-434-4 Http://www.ActaMath.com Acta Mathematica Sinica, English Series
More informationChebfun and equispaced data
Chebfun and equispaced data Georges Klein University of Oxford SIAM Annual Meeting, San Diego, July 12, 2013 Outline 1 1D Interpolation, Chebfun and equispaced data 2 3 (joint work with R. Platte) G. Klein
More informationINVERSION OF ELECTROCHEMICAL IMPEDANCE SPECTROSCOPIC MEASUREMENTS
BIOFUEL CELL POLARIZATION ESTIMATION INVERSION OF ELECTROCHEMICAL IMPEDANCE SPECTROSCOPIC MEASUREMENTS IMPORTANCE OF MODEL FORMULATION Rosemary Renaut, Mathematics, ASU Jakob Hansen, Jarom Hogue and Grant
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationEnhanced Compressive Sensing and More
Enhanced Compressive Sensing and More Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Nonlinear Approximation Techniques Using L1 Texas A & M University
More informationCompressive sampling meets seismic imaging
Compressive sampling meets seismic imaging Felix J. Herrmann fherrmann@eos.ubc.ca http://slim.eos.ubc.ca joint work with Tim Lin and Yogi Erlangga Seismic Laboratory for Imaging & Modeling Department of
More informationRecovery of high order accuracy in radial basis function approximation for discontinuous problems
Recovery of high order accuracy in radial basis function approximation for discontinuous problems Chris L. Bresten, Sigal Gottlieb 1, Daniel Higgs, Jae-Hun Jung* 2 Abstract Radial basis function(rbf) methods
More informationExact Signal Recovery from Sparsely Corrupted Measurements through the Pursuit of Justice
Exact Signal Recovery from Sparsely Corrupted Measurements through the Pursuit of Justice Jason N. Laska, Mark A. Davenport, Richard G. Baraniuk Department of Electrical and Computer Engineering Rice University
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationA Non-sparse Tutorial on Sparse FFTs
A Non-sparse Tutorial on Sparse FFTs Mark Iwen Michigan State University April 8, 214 M.A. Iwen (MSU) Fast Sparse FFTs April 8, 214 1 / 34 Problem Setup Recover f : [, 2π] C consisting of k trigonometric
More informationCompressed Sensing and Neural Networks
and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationRandom projections. 1 Introduction. 2 Dimensionality reduction. Lecture notes 5 February 29, 2016
Lecture notes 5 February 9, 016 1 Introduction Random projections Random projections are a useful tool in the analysis and processing of high-dimensional data. We will analyze two applications that use
More informationCOMPUTATIONAL ISSUES RELATING TO INVERSION OF PRACTICAL DATA: WHERE IS THE UNCERTAINTY? CAN WE SOLVE Ax = b?
COMPUTATIONAL ISSUES RELATING TO INVERSION OF PRACTICAL DATA: WHERE IS THE UNCERTAINTY? CAN WE SOLVE Ax = b? Rosemary Renaut http://math.asu.edu/ rosie BRIDGING THE GAP? OCT 2, 2012 Discussion Yuen: Solve
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 information1-Bit Matrix Completion
1-Bit Matrix Completion Mark A. Davenport School of Electrical and Computer Engineering Georgia Institute of Technology Yaniv Plan Mary Wootters Ewout van den Berg Matrix Completion d When is it possible
More informationReconstructing conductivities with boundary corrected D-bar method
Reconstructing conductivities with boundary corrected D-bar method Janne Tamminen June 24, 2011 Short introduction to EIT The Boundary correction procedure The D-bar method Simulation of measurement data,
More informationImproved Radon Based Imaging using the Shearlet Transform
Improved Radon Based Imaging using the Shearlet Transform Glenn R. Easley a, Flavia Colonna b, Demetrio Labate c a System Planning Corporation, Arlington, Virginia b George Mason University, Fairfax, Virginia
More informationTRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS
TRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München, München, Arcistraße
More informationCS 229r: Algorithms for Big Data Fall Lecture 19 Nov 5
CS 229r: Algorithms for Big Data Fall 215 Prof. Jelani Nelson Lecture 19 Nov 5 Scribe: Abdul Wasay 1 Overview In the last lecture, we started discussing the problem of compressed sensing where we are given
More informationAN INVERTIBLE DISCRETE AUDITORY TRANSFORM
COMM. MATH. SCI. Vol. 3, No. 1, pp. 47 56 c 25 International Press AN INVERTIBLE DISCRETE AUDITORY TRANSFORM JACK XIN AND YINGYONG QI Abstract. A discrete auditory transform (DAT) from sound signal to
More informationNoisy Signal Recovery via Iterative Reweighted L1-Minimization
Noisy Signal Recovery via Iterative Reweighted L1-Minimization Deanna Needell UC Davis / Stanford University Asilomar SSC, November 2009 Problem Background Setup 1 Suppose x is an unknown signal in R d.
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationSparse Approximation: from Image Restoration to High Dimensional Classification
Sparse Approximation: from Image Restoration to High Dimensional Classification Bin Dong Beijing International Center for Mathematical Research Beijing Institute of Big Data Research Peking University
More informationSolving Poisson s Equations Using Buffered Fourier Spectral Method
Solving Poisson s Equations Using Buffered Fourier Spectral Method Yinlin Dong Hassan Abd Salman Al-Dujaly Chaoqun Liu Technical Report 2015-12 http://www.uta.edu/math/preprint/ Solving Poisson s Equations
More informationUsing Frame Theoretic Convolutional Gridding for Robust Synthetic Aperture Sonar Imaging
Using Frame Theoretic Convolutional Gridding for Robust Synthetic Aperture Sonar Imaging John McKay, Anne Gelb, Vishal Monga, Raghu G. Raj Department of Electrical Engineering, Pennsylvania State University,
More informationA Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms
Universal Journal of Integral Equations 4 (2016), 54-62 www.papersciences.com A Singular Integral Transform for the Gibbs-Wilbraham Effect in Inverse Fourier Transforms N. H. S. Haidar CRAMS: Center for
More informationStatistical Issues in Searches: Photon Science Response. Rebecca Willett, Duke University
Statistical Issues in Searches: Photon Science Response Rebecca Willett, Duke University 1 / 34 Photon science seen this week Applications tomography ptychography nanochrystalography coherent diffraction
More informationNear Ideal Behavior of a Modified Elastic Net Algorithm in Compressed Sensing
Near Ideal Behavior of a Modified Elastic Net Algorithm in Compressed Sensing M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas M.Vidyasagar@utdallas.edu www.utdallas.edu/ m.vidyasagar
More informationCompressive Sensing and Beyond
Compressive Sensing and Beyond Sohail Bahmani Gerorgia Tech. Signal Processing Compressed Sensing Signal Models Classics: bandlimited The Sampling Theorem Any signal with bandwidth B can be recovered
More informationSIGNAL AND IMAGE RESTORATION: SOLVING
1 / 55 SIGNAL AND IMAGE RESTORATION: SOLVING ILL-POSED INVERSE PROBLEMS - ESTIMATING PARAMETERS Rosemary Renaut http://math.asu.edu/ rosie CORNELL MAY 10, 2013 2 / 55 Outline Background Parameter Estimation
More informationSeismic compression. François G. Meyer Department of Electrical Engineering University of Colorado at Boulder
Seismic compression François G. Meyer Department of Electrical Engineering University of Colorado at Boulder francois.meyer@colorado.edu http://ece-www.colorado.edu/ fmeyer IPAM, MGA 2004 Seismic Compression
More informationInterpolation via weighted l 1 minimization
Interpolation via weighted l 1 minimization Rachel Ward University of Texas at Austin December 12, 2014 Joint work with Holger Rauhut (Aachen University) Function interpolation Given a function f : D C
More informationA new method on deterministic construction of the measurement matrix in compressed sensing
A new method on deterministic construction of the measurement matrix in compressed sensing Qun Mo 1 arxiv:1503.01250v1 [cs.it] 4 Mar 2015 Abstract Construction on the measurement matrix A is a central
More informationFast and Robust Phase Retrieval
Fast and Robust Phase Retrieval Aditya Viswanathan aditya@math.msu.edu CCAM Lunch Seminar Purdue University April 18 2014 0 / 27 Joint work with Yang Wang Mark Iwen Research supported in part by National
More informationRapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization
Rapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization Shuyang Ling Department of Mathematics, UC Davis Oct.18th, 2016 Shuyang Ling (UC Davis) 16w5136, Oaxaca, Mexico Oct.18th, 2016
More informationInverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France.
Inverse problems and sparse models (6/6) Rémi Gribonval INRIA Rennes - Bretagne Atlantique, France remi.gribonval@inria.fr Overview of the course Introduction sparsity & data compression inverse problems
More informationThresholds for the Recovery of Sparse Solutions via L1 Minimization
Thresholds for the Recovery of Sparse Solutions via L Minimization David L. Donoho Department of Statistics Stanford University 39 Serra Mall, Sequoia Hall Stanford, CA 9435-465 Email: donoho@stanford.edu
More informationA Survey of Compressive Sensing and Applications
A Survey of Compressive Sensing and Applications Justin Romberg Georgia Tech, School of ECE ENS Winter School January 10, 2012 Lyon, France Signal processing trends DSP: sample first, ask questions later
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 informationModel-Based Compressive Sensing for Signal Ensembles. Marco F. Duarte Volkan Cevher Richard G. Baraniuk
Model-Based Compressive Sensing for Signal Ensembles Marco F. Duarte Volkan Cevher Richard G. Baraniuk Concise Signal Structure Sparse signal: only K out of N coordinates nonzero model: union of K-dimensional
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationMinimizing Isotropic Total Variation without Subiterations
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Minimizing Isotropic Total Variation without Subiterations Kamilov, U. S. TR206-09 August 206 Abstract Total variation (TV) is one of the most
More informationGauge optimization and duality
1 / 54 Gauge optimization and duality Junfeng Yang Department of Mathematics Nanjing University Joint with Shiqian Ma, CUHK September, 2015 2 / 54 Outline Introduction Duality Lagrange duality Fenchel
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 informationA discretized Newton flow for time varying linear inverse problems
A discretized Newton flow for time varying linear inverse problems Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München Arcisstrasse
More informationA Lower Bound Theorem. Lin Hu.
American J. of Mathematics and Sciences Vol. 3, No -1,(January 014) Copyright Mind Reader Publications ISSN No: 50-310 A Lower Bound Theorem Department of Applied Mathematics, Beijing University of Technology,
More informationVariational Methods in Signal and Image Processing
Variational Methods in Signal and Image Processing XU WANG Texas A&M University Dept. of Electrical & Computer Eng. College Station, Texas United States xu.wang@tamu.edu ERCHIN SERPEDIN Texas A&M University
More informationPhase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms
Phase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms Aditya Viswanathan Department of Mathematics and Statistics adityavv@umich.edu www-personal.umich.edu/
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 informationSimultaneous Sparsity
Simultaneous Sparsity Joel A. Tropp Anna C. Gilbert Martin J. Strauss {jtropp annacg martinjs}@umich.edu Department of Mathematics The University of Michigan 1 Simple Sparse Approximation Work in the d-dimensional,
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 informationLecture 2: Reconstruction and decomposition of vector fields on the sphere with applications
2013 Dolomites Research Week on Approximation : Reconstruction and decomposition of vector fields on the sphere with applications Grady B. Wright Boise State University What's the problem with vector fields
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationIntroduction to the Mathematics of Medical Imaging
Introduction to the Mathematics of Medical Imaging Second Edition Charles L. Epstein University of Pennsylvania Philadelphia, Pennsylvania EiaJTL Society for Industrial and Applied Mathematics Philadelphia
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
More informationSolving DC Programs that Promote Group 1-Sparsity
Solving DC Programs that Promote Group 1-Sparsity Ernie Esser Contains joint work with Xiaoqun Zhang, Yifei Lou and Jack Xin SIAM Conference on Imaging Science Hong Kong Baptist University May 14 2014
More informationVariational methods for restoration of phase or orientation data
Variational methods for restoration of phase or orientation data Martin Storath joint works with Laurent Demaret, Michael Unser, Andreas Weinmann Image Analysis and Learning Group Universität Heidelberg
More informationNoncommutative Uncertainty Principle
Noncommutative Uncertainty Principle Zhengwei Liu (joint with Chunlan Jiang and Jinsong Wu) Vanderbilt University The 12th East Coast Operator Algebras Symposium, Oct 12, 2014 Z. Liu (Vanderbilt) Noncommutative
More information1-Bit Matrix Completion
1-Bit Matrix Completion Mark A. Davenport School of Electrical and Computer Engineering Georgia Institute of Technology Yaniv Plan Mary Wootters Ewout van den Berg Matrix Completion d When is it possible
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Error Estimates for Filtered Back Projection Matthias Beckmann and Armin Iske Nr. 2015-03 January 2015 Error Estimates for Filtered Back Projection Matthias
More informationSparse and Fast Fourier Transforms in Biomedical Imaging
Sparse and Fast Fourier Transforms in Biomedical Imaging Stefan Kunis 1 1 joint work with Ines Melzer, TU Chemnitz and Helmholtz Center Munich, supported by DFG and Helmholtz association Outline 1 Trigonometric
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 5: Correlation and Power-Spectrum CCNY
1 BME 50500: Image and Signal Processing in Biomedicine Lecture 5: Correlation and Power-Spectrum Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More informationSUMMARY. H (ω)u(ω,x s ;x) :=
Interpolating solutions of the Helmholtz equation with compressed sensing Tim TY Lin*, Evgeniy Lebed, Yogi A Erlangga, and Felix J Herrmann, University of British Columbia, EOS SUMMARY We present an algorithm
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
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 information