Fourier phase retrieval with random phase illumination
|
|
- Charity Wilcox
- 5 years ago
- Views:
Transcription
1 Fourier phase retrieval with random phase illumination Wenjing Liao School of Mathematics Georgia Institute of Technology Albert Fannjiang Department of Mathematics, UC Davis. IMA August 5, 27
2 Classical phase retrieval Fourier Transform: Let {f (n)} be a discrete array with n = (n,..., n d ) Z d. We assume n N and d 2. F (ω) = N f (n)e 2πiω n n= ω = (ω,..., ω d ) [, ] d Phase retrieval: recover f from the Fourier intensity measurements Y = Φf, where Φ represents Fourier transform. Chanllenge: non-uniqueness non-convexity 2 / 26
3 Non-uniqueness Sampling: L = F (ω) 2 = N n= N f (m + n)f (m) m }{{} auto-correlation function of f { (ω,..., ω d ) : ω k =, 2N k +,..., Oversampling ratio = 2 d e 2πin ω } 2N k 2N k +, k =,..., d Trivial ambiguity: the auto-correlation function is invariant under: spatial translation: f ( ) f ( + t) conjugate reflection: f ( ) f (N ) constant global phase change: f ( ) e iθ f ( ) 2 Nontrivial ambiguity: exists for few images when d 2 [Bruck & Sodin 979, Bates 982, Hayes 982] 3 / 26
4 Non-convexity Standard phasing algorithms: [Fienup 982, Marchesini 27] Find f such that f n (n)e 2πiω n = n f (n)e 2πiω n, ω L 2 f satisfies the same constraint as f, such as the support or positivity constraint (a) cameraman (b) phantom Figure : Images with a loose support 4 / 26
5 relative residual at each iteration iteration relative residual at each iteration iteration relative residual at each iteration iteration relative residual at each iteration iteration Numerical stagnation [Fienup and Wackerman 986] Residual: r(f k ) = Φf k Y Algorithms: Error Reduction (ER) [Gerchberg & Saxton 972, Fienup 982] and Hybrid Input Output (HIO) [Fienup 982] relative residual relative residual (a) ER (b) Residual (c) HIO + ER (d) Residual relative residual relative residual (e) ER (f) Residual (g) HIO + ER (h) Residual 5 / 26
6 Outline Our approach: random phase illumination 2 Recovery with exact mask information Uniqueness Numerical results 3 Recovery with phase-uncertain mask Uniqueness Phasing algorithms Numerical results 6 / 26
7 Random phase illumination Structured illuminations for phase retrieval: Ptychography [Thibault et al. 28, Dierolf et al. 2, Maiden and Rodenburg 29, Marchesini et al ] Masking [Marchesini et al. 28, Candès et al. 2,... ] Inverse Problems 28 (22) 758 Random phase mask: Inverse Problems 28 (22) 758 A Fannjiang A Fannjiang (a) (a) (b) (b) Figure. Illumination Figureof. aillumination partially transparent of a partiallyobject transparent (the blue object oval) (thewith blue aoval) deterministic with a deterministic (a) or (a) or random field λ created randomby field a diffuser λ created(b) byfollowed a diffuser by (b) followed intensity by an measurement intensity measurement of the diffraction of the diffraction pattern. In the case pattern. of wavefront In the casereconstruction, of wavefront reconstruction, therandom modulator theillumination randomis modulator placed at isthe placed exitatpupil the exit pupil instead of the entrance insteadpupil of theas entrance in (b). pupil as in (b). deterministic illumination λ(n) = λ(n) = e iφ(n) : φ(n) random, known or roughly known problem of stagnation problem of may stagnation be duemay to the be possibility due to the possibility of the iterative of the process iterative process to approach to approach the object and the its object twin or and shifted its twin image, or shifted the support image, the not support tight enough not tight or enough the boundary or the boundary not not sharp enough [3, sharp 4, enough 8]. [3, Besides 4, 8]. the Besides uniqueness the uniqueness issue, phase issue, retrieval phase retrieval is also inherently is also inherently nonconvex and many nonconvex and many researchers have researchers have believed the lack of convexity in the Fourier magnitude believed the lack of convexity in the Fourier magnitude constraint to be a main, constraint to be a main, if not the if not the dominant, source of numerical problems with the standard dominant, source of numerical problems with the standard phasing algorithms [4, 2, 3]. While there have been dazzling advances in applications of phasing algorithms [4, 2, 3]. While there f (n)λ(n)e 2πin ω have been dazzling advances in applications of phase retrieval in the past decades [22], we still do not know just how much of the error and phase retrieval in the past decades [22], we still do not know just how much of the error and stagnation problems is attributable to to the lack of uniqueness or convexity. stagnation problems is attributable to to the lack of uniqueness or convexity. We propose here ton refocus on the issue of uniqueness as uniqueness is undoubtedly We propose here to refocus on the issue of uniqueness as uniqueness is undoubtedly the first foundational issue of any inverse problem, including phase retrieval. Specifically the first foundational issue of any inverse problem, including phase retrieval. Specifically we will first establish uniqueness in the absolute sense with random illumination under we will first establish general, physically uniqueness reasonable in the absolute object constraints sense with (figure random ) and illumination secondlydemonstrate under that7 / 26 Measurement: Y (ω) =
8 How does the problem change? f physical process Φλf Classical phase retrieval Φf λ(n) = phase retrieval f Phase retrieval with random phase mask } Φλf phase retrieval f known or roughly known λ(n) 8 / 26
9 Outline Our approach: random phase illumination 2 Recovery with exact mask information Uniqueness Numerical results 3 Recovery with phase-uncertain mask Uniqueness Phasing algorithms Numerical results 9 / 26
10 Uniqueness with exact mask information Real-valued images: Theorem [Fannjiang 22] Suppose f is real-valued and its support has rank 2. Then with probability one, f is determined by one random phase illumination up to a ± sign. Complex-valued images: Theorem [Fannjiang 22] Suppose the support of f has rank 2. Then with probability one, f is determined by two random phase illuminations up to a global phase. The rank of a support set is the dimension of its convex hull. / 26
11 relative residual at each iteration iteration relative residual at each iteration iteration relative residual at each iteration iteration relative residual at each iteration iteration Numerical results relative residual relative residual (a) ER (b) Residual (c) HIO + ER (d) Residual relative residual relative residual (e) ER (f) Residual (g) HIO + ER (h) Residual We observe global convergence as the residual decreases to. / 26
12 Stability to noise Error versus noise-to-signal ratio average relative error in 5 trials relative error versus noise level high resolution & gaussian noise low resolution & gaussian noise high resolution & poisson noise low resolution & poisson noise high resolution & illuminator noise low resolution & illuminator noise average relative error in 5 trials relative error versus noise level one uniform + one high resolution & gaussian noise one uniform + one low resolution & gaussian noise one uniform + one high resolution & poisson noise one uniform + one low resolution & poisson noise one uniform + one high resolution & illuminator noise one uniform + one low resolution & illuminator noise noise percentage noise percentage Left: non-negative phantom with one random phase mask and oversampling rate = 2 Right: complex-valued image whose intensity is the phantom with two random phase masks and oversampling rate = 3 2 / 26
13 Outline Our approach: random phase illumination 2 Recovery with exact mask information Uniqueness Numerical results 3 Recovery with phase-uncertain mask Uniqueness Phasing algorithms Numerical results 3 / 26
14 Phase-uncertain mask Exact mask: λ(n) = e iφ(n) Mask estimate: λ (n) = e iφ (n) φ(n) φ (n) δπ, φ (n) + δπ φ (n) ± δπ, n, φ (n) φ(n) δπ Fourier magnitude data: Y = Φλf Fourier phase retrieval with a roughly known mask: } Φλf phase retrieval f λ 4 / 26
15 Uniqueness for real-valued images with one phase-uncertain mask Theorem [Fannjiang and L. 23] Suppose f is a real-valued image with sparsity S and its support has rank 2. The exact mask phases {φ(n)} are i.i.d. samples in [, 2π) and δ < /2. If we can find real-valued f and λ(n) = exp (i φ(n)) such that then Φ λ f = Φλf ; λ satisfies mask constraint: φ(n) φ (n) ± δπ, f (n) = ±f (n) n λ(n) = e iθ λ(n) on n where f (n) with probability no less than (2δ) S/2 d j= (N j + ). 5 / 26
16 Uniqueness for complex-valued images with two phase-uncertain masks Theorem [Fannjiang and L. 23] Suppose f is a complex-valued image with sparsity S and its support has rank 2. The first mask λ () = λ has i.i.d. random phases and δ <. The second mask λ (2) is exactly known. If we can find f and λ(n) = exp (i φ(n)) such that then Φ λ f = Φλf and Φλ (2) f = Φλ (2) f ; λ satisfies mask constraint: φ(n) φ (n) ± δπ, f (n) = e iα f (n) n λ(n) = e iα2 λ(n) on n where f (n), with probability no less than δ S/2 d j= (N j + ). 6 / 26
17 Algorithm: alternating updates between image and mask Example: recovery of real-valued images with one phase-uncertain mask Objective: Find f and λ such that Φ λ f = Φλf ; f is real-valued; λ satisfies mask constraint. Alternating updates: (f k, λ k ) standard phasing method of (f k+, λ k ) (f k+, λ k+ ). algorithms alternating projection Image update: (f k, λ k ) (f k+, λ k ) by standard phasing algorithms [Fienup 982, Marchesini 27] Mask update: (f k+, λ k ) (f k+, λ k+ ) by alternating projection 7 / 26
18 standard phasing Image update (f k, λ k ) (f k+, λ k ) algorithms Two operators: Object-constraint fitting operator P o Fourier intensity data fitting operator P f,k = λ k Φ T Φλ k where { Y (ω) exp (i G(ω)) if G(ω) > T G(ω) = Y (ω) if G(ω) =. where Y = Φλf. Two phasing algorithms: Error Reduction (ER) [Gerchberg & Saxton 972, Fienup 982]: f k+ = P o P f,k f k Douglas-Rachford (DR) [Douglas & Rachford 956]: f k+ = I + (2P o I )(2P f,k I ) f k. 2 8 / 26
19 method of Mask udpate (f k+, λ k ) (f k+, λ k+ ) alternating projection Two operators: Fourier intensity data fitting operator Q f,k { ] [Φ λ k = Q f,kλ k (n) = T ΦΛ k f k+ (n)/f k+ (n) if f k+ (n). λ k (n) else. Mask-constraint fitting operator Q m φ (n) δπ { Q m λ exp(i λ k (n) = k (n)) orthogonal projection λ k (n) φ (n) ± δπ else.. Mask update by alternating projection: λ k+ = Q m Q f,k λ k 9 / 26
20 Alternating updates between image and mask (f k, λ k ) standard phasing method of (f k+, λ k ) (f k+, λ k+ ) algorithms alternating projection Alternating Error Reduction (AER): (f k+, λ k+ ) = (P o P f,k f k, Q m Q f,k λ k ) (AER) Theorem Let r(f k, λ k ) = ΦΛ k f k Y. ER iteration satisfies r(f k+, λ k+ ) r(f k, λ k ). Alternating Douglas Rachford and Error Reduction (DRER): ( ) (f k+, λ k+ ) = 2 [I + (2P o I )(2P f,k I )]f k, Q m Q f,k λ k (DRER) 2 / 26
21 Test images and masks Test images: cameraman mandrill 2 2 phantom High-Resolution Mask (HRM) and Low-Resolution Mask (LRM): HRM LRM 2 / 26
22 Non-negative images and one phase-uncertain mask Illumination: One low-resolution phase-uncertain mask with uncertainty δ =.3 Algorithm: DRER + AER recovered images mask error residual.9 Relative residual of each iterate iteration.9 Relative residual of each iterate iteration.9 Relative residual of each iterate iteration 22 / 26
23 How large can the mask uncertainty δ be? Image error versus δ with % gaussian or poisson noise average relative error in 5 trials relative error versus mask uncertainty HRM, % gaussian nosie LRM, % gaussian nosie HRM, % poisson noise LRM, % poisson noise average relative error in 5 trials relative error versus mask uncertainty HRM, % gaussian nosie LRM, % gaussian nosie HRM, % poisson noise LRM, % poisson noise uncertainty level (i) non-negative cameraman with one phase-uncertain mask uncertainty level (j) complex-valued image whose intensity is cameraman with one uniform illumination and one phaseuncertain mask 23 / 26
24 Summary Phase retrieval with random phase illumination uniqueness global convergence and stability of the standard phasing algorithms 2 Phase-uncertain mask uniqueness of the image and the mask algorithm: alternating update between the image and the mask good reconstruction of the image and the mask 3 Convergence analysis for non-convex problems Wirtinger flow [Candès, Li and Soltanolkotabi 25] Geometric analysis of phase retrieval [Sun, Qu and Wright 26] 24 / 26
25 Continuing: sensor calibration in imaging N sensors s point sources location: ω j [, ) amplitudes: x j C Measurements: {y n : n =,..., N } where Unknowns: y n = g n s j= x je 2πinω j + e n Calibration parameters {g n } N n= Source locations {ω j } s j= and source amplitudes {x j} s j= Connection to phase retrieval: equivalent to phase retrieval if g n = Y. C. Eldar, W. Liao and S. Tang, Sensor calibration for off-the-grid spectral estimation, arxiv: / 26
26 Thank you for your attention! Wenjing Liao Georgia Institute of Technology 26 / 26
Phase Retrieval with Random Illumination
Phase Retrieval with Random Illumination Albert Fannjiang, UC Davis Collaborator: Wenjing Liao NCTU, July 212 1 Phase retrieval Problem Reconstruct the object f from the Fourier magnitude Φf. Why do we
More informationFOURIER PHASING WITH PHASE-UNCERTAIN MASK
FOURIER PHASING WITH PHASE-UNCERTAIN MASK ALBERT FANNJIANG AND WENJING LIAO Abstract. Uniqueness, up to a global phase, is proved for phasing with a random phaseuncertain mask (PUM). Phasing algorithms
More informationPhase retrieval with random phase illumination
A. Fannjiang and W. Liao Vol. 29, No. 9 / September 212 / J. Opt. Soc. Am. A 1847 Phase retrieval with random phase illumination Albert Fannjiang* and Wenjing Liao Department of Mathematics, University
More informationABSOLUTE UNIQUENESS OF PHASE RETRIEVAL WITH RANDOM ILLUMINATION
ABSOLUTE UNIQUENESS OF PHASE RETRIEVAL WITH RANDOM ILLUMINATION ALBERT FANNJIANG Abstract. Random illumination is proposed to enforce absolute uniqueness and resolve all types of ambiguity, trivial or
More informationNonlinear Analysis with Frames. Part III: Algorithms
Nonlinear Analysis with Frames. Part III: Algorithms Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD July 28-30, 2015 Modern Harmonic Analysis and Applications
More informationarxiv: v1 [cs.it] 26 Oct 2015
Phase Retrieval: An Overview of Recent Developments Kishore Jaganathan Yonina C. Eldar Babak Hassibi Department of Electrical Engineering, Caltech Department of Electrical Engineering, Technion, Israel
More informationarxiv: FOURIER PHASE RETRIEVAL WITH A SINGLE MASK BY THE DOUGLAS-RACHFORD ALGORITHM
arxiv:59.888 FOURIER PHASE RETRIEVAL WITH A SINGLE MASK BY THE DOUGLAS-RACHFORD ALGORITHM PENGWEN CHEN AND ALBERT FANNJIANG Abstract. Douglas-Rachford (DR) algorithm is analyzed for Fourier phase retrieval
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 informationPhase Retrieval by Linear Algebra
Phase Retrieval by Linear Algebra Pengwen Chen Albert Fannjiang Gi-Ren Liu December 12, 2016 arxiv:1607.07484 Abstract The null vector method, based on a simple linear algebraic concept, is proposed as
More informationMismatch and Resolution in CS
Mismatch and Resolution in CS Albert Fannjiang, UC Davis Coworker: Wenjing Liao, UC Davis SPIE Optics + Photonics, San Diego, August 2-25, 20 Mismatch: gridding error Band exclusion Local optimization
More informationFOURIER PHASE RETRIEVAL WITH A SINGLE MASK BY DOUGLAS-RACHFORD ALGORITHMS
FOURIER PHASE RETRIEVAL WITH A SINGLE MASK BY DOUGLAS-RACHFORD ALGORITHMS PENGWEN CHEN AND ALBERT FANNJIANG Abstract. The Fourier-domain Douglas-Rachford (FDR) algorithm is analyzed for phase retrieval
More informationPHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION. A Thesis MELTEM APAYDIN
PHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION A Thesis by MELTEM APAYDIN Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the
More informationGlobally Convergent Levenberg-Marquardt Method For Phase Retrieval
Globally Convergent Levenberg-Marquardt Method For Phase Retrieval Zaiwen Wen Beijing International Center For Mathematical Research Peking University Thanks: Chao Ma, Xin Liu 1/38 Outline 1 Introduction
More informationFrequency Subspace Amplitude Flow for Phase Retrieval
Research Article Journal of the Optical Society of America A 1 Frequency Subspace Amplitude Flow for Phase Retrieval ZHUN WEI 1, WEN CHEN 2,3, AND XUDONG CHEN 1,* 1 Department of Electrical and Computer
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 informationarxiv: v1 [physics.optics] 5 Mar 2012
Designing and using prior knowledge for phase retrieval Eliyahu Osherovich, Michael Zibulevsky, and Irad Yavneh arxiv:1203.0879v1 [physics.optics] 5 Mar 2012 Computer Science Department, Technion Israel
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 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 informationPhase recovery with PhaseCut and the wavelet transform case
Phase recovery with PhaseCut and the wavelet transform case Irène Waldspurger Joint work with Alexandre d Aspremont and Stéphane Mallat Introduction 2 / 35 Goal : Solve the non-linear inverse problem Reconstruct
More informationCoherent imaging without phases
Coherent imaging without phases Miguel Moscoso Joint work with Alexei Novikov Chrysoula Tsogka and George Papanicolaou Waves and Imaging in Random Media, September 2017 Outline 1 The phase retrieval problem
More informationPhase retrieval with random Gaussian sensing vectors by alternating projections
Phase retrieval with random Gaussian sensing vectors by alternating projections arxiv:609.03088v [math.st] 0 Sep 06 Irène Waldspurger Abstract We consider a phase retrieval problem, where we want to reconstruct
More informationJournal of Fourier Analysis and Applications (doi: /s )
Journal of Fourier Analysis and Applications (doi:10.1007/s00041-017-9536-8) PHASE RETRIEVAL WITH ONE OR TWO DIFFRACTION PATTERNS BY ALTERNATING PROJECTIONS WITH THE NULL INITIALIZATION PENGWEN CHEN, ALBERT
More informationSolving Corrupted Quadratic Equations, Provably
Solving Corrupted Quadratic Equations, Provably Yuejie Chi London Workshop on Sparse Signal Processing September 206 Acknowledgement Joint work with Yuanxin Li (OSU), Huishuai Zhuang (Syracuse) and Yingbin
More informationThe Iterative and Regularized Least Squares Algorithm for Phase Retrieval
The Iterative and Regularized Least Squares Algorithm for Phase Retrieval Radu Balan Naveed Haghani Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD April 17, 2016
More informationA Geometric Analysis of Phase Retrieval
A Geometric Analysis of Phase Retrieval Ju Sun, Qing Qu, John Wright js4038, qq2105, jw2966}@columbiaedu Dept of Electrical Engineering, Columbia University, New York NY 10027, USA Abstract Given measurements
More informationGreedy Signal Recovery and Uniform Uncertainty Principles
Greedy Signal Recovery and Uniform Uncertainty Principles SPIE - IE 2008 Deanna Needell Joint work with Roman Vershynin UC Davis, January 2008 Greedy Signal Recovery and Uniform Uncertainty Principles
More informationSelf-Calibration and Biconvex Compressive Sensing
Self-Calibration and Biconvex Compressive Sensing Shuyang Ling Department of Mathematics, UC Davis July 12, 2017 Shuyang Ling (UC Davis) SIAM Annual Meeting, 2017, Pittsburgh July 12, 2017 1 / 22 Acknowledgements
More informationPhase Retrieval from Local Measurements in Two Dimensions
Phase Retrieval from Local Measurements in Two Dimensions Mark Iwen a, Brian Preskitt b, Rayan Saab b, and Aditya Viswanathan c a Department of Mathematics, and Department of Computational Mathematics,
More informationINEXACT ALTERNATING OPTIMIZATION FOR PHASE RETRIEVAL WITH OUTLIERS
INEXACT ALTERNATING OPTIMIZATION FOR PHASE RETRIEVAL WITH OUTLIERS Cheng Qian Xiao Fu Nicholas D. Sidiropoulos Lei Huang Dept. of Electronics and Information Engineering, Harbin Institute of Technology,
More informationReshaped Wirtinger Flow for Solving Quadratic System of Equations
Reshaped Wirtinger Flow for Solving Quadratic System of Equations Huishuai Zhang Department of EECS Syracuse University Syracuse, NY 344 hzhan3@syr.edu Yingbin Liang Department of EECS Syracuse University
More informationMethoden moderner Röntgenphysik I. Coherence based techniques II. Christian Gutt DESY, Hamburg
Methoden moderner Röntgenphysik I Coherence based techniques II Christian Gutt DESY Hamburg christian.gutt@desy.de 8. January 009 Outline 18.1. 008 Introduction to Coherence 8.01. 009 Structure determination
More informationNonlinear optimization algorithm for retrieving the full complex pupil function
Nonlinear optimization algorithm for retrieving the full complex pupil function Gregory R. Brady and James R. Fienup The Institute of Optics, University of Rochester, Rochester, New York 14627-0186 gbrady@optics.rochester.edu
More informationOn The 2D Phase Retrieval Problem
1 On The 2D Phase Retrieval Problem Dani Kogan, Yonina C. Eldar, Fellow, IEEE, and Dan Oron arxiv:1605.08487v1 [cs.it] 27 May 2016 Abstract The recovery of a signal from the magnitude of its Fourier transform,
More informationHybrid projection reflection method for phase retrieval
Bauschke et al. Vol. 20, No. 6/June 2003/J. Opt. Soc. Am. A 1025 Hybrid projection reflection method for phase retrieval Heinz H. Bauschke Department of Mathematics and Statistics, University of Guelph,
More informationMitigating the effect of noise in the hybrid input output method of phase retrieval
Mitigating the effect of noise in the hybrid input output method of phase retrieval Russell Trahan III* and David Hyland Department of Aerospace Engineering, Texas A&M University, HRBB 70, Ross Street,
More informationInformation and Resolution
Information and Resolution (Invited Paper) Albert Fannjiang Department of Mathematics UC Davis, CA 95616-8633. fannjiang@math.ucdavis.edu. Abstract The issue of resolution with complex-field measurement
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 informationIncremental Reshaped Wirtinger Flow and Its Connection to Kaczmarz Method
Incremental Reshaped Wirtinger Flow and Its Connection to Kaczmarz Method Huishuai Zhang Department of EECS Syracuse University Syracuse, NY 3244 hzhan23@syr.edu Yingbin Liang Department of EECS Syracuse
More informationIntroduction to Sparsity. Xudong Cao, Jake Dreamtree & Jerry 04/05/2012
Introduction to Sparsity Xudong Cao, Jake Dreamtree & Jerry 04/05/2012 Outline Understanding Sparsity Total variation Compressed sensing(definition) Exact recovery with sparse prior(l 0 ) l 1 relaxation
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 informationPHASE RETRIEVAL FOR IMAGING PROBLEMS.
PHASE RETRIEVAL FOR IMAGING PROBLEMS. FAJWEL FOGEL, IRÈNE WALDSPURGER, AND ALEXANDRE D ASPREMONT ABSTRACT. We study convex relaxation algorithms for phase retrieval on imaging problems. We show that exploiting
More informationarxiv: v1 [math.na] 16 Apr 2019
A direct solver for the phase retrieval problem in ptychographic imaging arxiv:1904.07940v1 [math.na] 16 Apr 2019 Nada Sissouno a,b,, Florian Boßmann d, Frank Filbir b, Mark Iwen e, Maik Kahnt c, Rayan
More informationMitigating the Effect of Noise in Iterative Projection Phase Retrieval
Mitigating the Effect of Noise in Iterative Projection Phase Retrieval & David Hyland Imaging and Applied Optics: Optics and Photonics Congress Introduction Filtering Results Conclusion Outline Introduction
More informationSeparating Laguerre-Gaussian Radial Modes Via Projective Measurements RACHEL SAMPSON
Separating Laguerre-Gaussian Radial Modes Via Projective Measurements RACHEL SAMPSON Laguerre-Gaussian (LG) Beams Cylindrically symmetric and form a basis set for paraxial light beams Properties of beams
More informationAuthor s Accepted Manuscript
Author s Accepted Manuscript Combining ptychographical algorithms with the Hybrid Input-Output (HIO) algorithm A.P. Konijnenberg, W.M.J. Coene, S.F. Pereira, H.P. Urbach www.elsevier.com/locate/ultramic
More informationConvex Phase Retrieval without Lifting via PhaseMax
Convex Phase etrieval without Lifting via PhaseMax Tom Goldstein * 1 Christoph Studer * 2 Abstract Semidefinite relaxation methods transform a variety of non-convex optimization problems into convex problems,
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 informationarxiv: v3 [astro-ph.im] 9 Nov 2018
Research Article Journal of the Optical Society of America A Wavefront retrieval through random pupil-plane phase probes: Gerchberg-Saxton approach EUGENE PLUZHNIK,2,*, DAN SIRBU,2, RUSLAN BELIKOV, EDUARDO
More informationPhaseless Imaging Algorithms and Applications
Phaseless Imaging Algorithms and Applications James R. Fienup Robert E. Hopkins Professor of Optics University of Rochester Institute of Optics fienup@optics.rochester.edu Presented to: Institute for Mathematics
More informationPhase Retrieval using the Iterative Regularized Least-Squares Algorithm
Phase Retrieval using the Iterative Regularized Least-Squares Algorithm Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD August 16, 2017 Workshop on Phaseless
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 informationStatic telescope aberration measurement using lucky imaging techniques
DOI 10.1007/s10686-012-9291-4 ORIGINAL ARTICLE Static telescope aberration measurement using lucky imaging techniques Marcos López-Marrero Luis Fernando Rodríguez-Ramos José Gil Marichal-Hernández José
More informationFundamental Limits of Weak Recovery with Applications to Phase Retrieval
Proceedings of Machine Learning Research vol 75:1 6, 2018 31st Annual Conference on Learning Theory Fundamental Limits of Weak Recovery with Applications to Phase Retrieval Marco Mondelli Department of
More informationFundamental Limits of PhaseMax for Phase Retrieval: A Replica Analysis
Fundamental Limits of PhaseMax for Phase Retrieval: A Replica Analysis Oussama Dhifallah and Yue M. Lu John A. Paulson School of Engineering and Applied Sciences Harvard University, Cambridge, MA 0238,
More informationSolving PhaseLift by low-rank Riemannian optimization methods for complex semidefinite constraints
Solving PhaseLift by low-rank Riemannian optimization methods for complex semidefinite constraints Wen Huang Université Catholique de Louvain Kyle A. Gallivan Florida State University Xiangxiong Zhang
More informationarxiv: v2 [cs.it] 20 Sep 2011
Phase Retrieval via Matrix Completion Emmanuel J. Candès, Yonina C. Eldar, Thomas Strohmer and Vladislav Voroninski August 2011 Abstract arxiv:1109.0573v2 [cs.it] 20 Sep 2011 This paper develops a novel
More informationCombining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation
UIUC CSL Mar. 24 Combining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation Yuejie Chi Department of ECE and BMI Ohio State University Joint work with Yuxin Chen (Stanford).
More informationBM3D-prGAMP: Compressive Phase Retrieval Based on BM3D Denoising
BM3D-prGAMP: Compressive Phase Retrieval Based on BM3D Denoising Chris Metzler, Richard Baraniuk Rice University Arian Maleki Columbia University Phase Retrieval Applications: Crystallography Microscopy
More informationarxiv: v2 [stat.ml] 12 Jun 2015
Phase Retrieval using Alternating Minimization Praneeth Netrapalli Prateek Jain Sujay Sanghavi June 15, 015 arxiv:1306.0160v [stat.ml] 1 Jun 015 Abstract Phase retrieval problems involve solving linear
More informationPHASE retrieval seeks to recover a signal from the magnitudes
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 64, NO. 3, DECEMBER 1, 016 6105 Phase Retrieval from 1D Fourier Measurements: Convexity, Uniqueness, and Algorithms Kejun Huang, Student Member, IEEE, Yonina
More informationOn the recovery of measures without separation conditions
Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Applied and Computational Mathematics Seminar Georgia Institute of Technology October
More informationThe Analysis Cosparse Model for Signals and Images
The Analysis Cosparse Model for Signals and Images Raja Giryes Computer Science Department, Technion. The research leading to these results has received funding from the European Research Council under
More informationEric G. Johnson: Diffractive and Micro-Optical Systems. Justin Coleman. College of Optical Sciences. The University of Arizona
Eric G. Johnson: Diffractive and Micro-Optical Systems Justin Coleman College of Optical Sciences The University of Arizona jcoleman@email.arizona.edu Abstract The material discussed in this paper relates
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More informationThe Projected Power Method: An Efficient Algorithm for Joint Alignment from Pairwise Differences
The Projected Power Method: An Efficient Algorithm for Joint Alignment from Pairwise Differences Yuxin Chen Emmanuel Candès Department of Statistics, Stanford University, Sep. 2016 Nonconvex optimization
More informationFourier Phase Retrieval: Uniqueness and Algorithms
Fourier Phase Retrieval: Uniqueness and Algorithms Tamir Bendory 1, Robert Beinert 2 and Yonina C. Eldar 3 arxiv:1705.09590v3 [cs.it] 6 Nov 2017 1 The Program in Applied and Computational Mathematics,
More informationSparse Phase Retrieval Using Partial Nested Fourier Samplers
Sparse Phase Retrieval Using Partial Nested Fourier Samplers Heng Qiao Piya Pal University of Maryland qiaoheng@umd.edu November 14, 2015 Heng Qiao, Piya Pal (UMD) Phase Retrieval with PNFS November 14,
More informationThe uniform uncertainty principle and compressed sensing Harmonic analysis and related topics, Seville December 5, 2008
The uniform uncertainty principle and compressed sensing Harmonic analysis and related topics, Seville December 5, 2008 Emmanuel Candés (Caltech), Terence Tao (UCLA) 1 Uncertainty principles A basic principle
More informationRecovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm
Recovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm J. K. Pant, W.-S. Lu, and A. Antoniou University of Victoria August 25, 2011 Compressive Sensing 1 University
More informationSpectral Initialization for Nonconvex Estimation: High-Dimensional Limit and Phase Transitions
Spectral Initialization for Nonconvex Estimation: High-Dimensional Limit and hase Transitions Yue M. Lu aulson School of Engineering and Applied Sciences Harvard University, USA Email: yuelu@seas.harvard.edu
More informationPhaseLift: A novel framework for phase retrieval. Vladislav Voroninski
PhaseLift: A novel framework for phase retrieval. by Vladislav Voroninski A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Mathematics in the
More informationarxiv: v1 [physics.optics] 27 May 2011
Iterative Algorithms for Ptychographic Phase Retrieval arxiv:1105.5628v1 [physics.optics] 27 May 2011 Chao Yang Computational Research Division, Lawrence Bereley National Laboratory, Bereley, CA 94720.
More informationPhase Retrieval: Hubble and the James Webb Space Telescope
Phase Retrieval: Hubble and the James Webb Space Telescope James R. Fienup Robert E. Hopkins Professor of Optics University of Rochester Institute of Optics fienup@optics.rochester.edu March 2003 JRF 2/03-1
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 informationComposite nonlinear models at scale
Composite nonlinear models at scale Dmitriy Drusvyatskiy Mathematics, University of Washington Joint work with D. Davis (Cornell), M. Fazel (UW), A.S. Lewis (Cornell) C. Paquette (Lehigh), and S. Roy (UW)
More informationSolving Large-scale Systems of Random Quadratic Equations via Stochastic Truncated Amplitude Flow
Solving Large-scale Systems of Random Quadratic Equations via Stochastic Truncated Amplitude Flow Gang Wang,, Georgios B. Giannakis, and Jie Chen Dept. of ECE and Digital Tech. Center, Univ. of Minnesota,
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 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 informationOptimal Spectral Initialization for Signal Recovery with Applications to Phase Retrieval
Optimal Spectral Initialization for Signal Recovery with Applications to Phase Retrieval Wangyu Luo, Wael Alghamdi Yue M. Lu Abstract We present the optimal design of a spectral method widely used to initialize
More informationResearch Statement Qing Qu
Qing Qu Research Statement 1/4 Qing Qu (qq2105@columbia.edu) Today we are living in an era of information explosion. As the sensors and sensing modalities proliferate, our world is inundated by unprecedented
More informationQuantized Iterative Hard Thresholding:
Quantized Iterative Hard Thresholding: ridging 1-bit and High Resolution Quantized Compressed Sensing Laurent Jacques, Kévin Degraux, Christophe De Vleeschouwer Louvain University (UCL), Louvain-la-Neuve,
More informationAlgorithms for phase retrieval problems
Algorithms for phase retrieval problems Part I Irène Waldspurger CNRS et CEREMADE (Université Paris Dauphine) Équipe MOKAPLAN (INRIA) 11 décembre 2017 Journées de géométrie algorithmique Aussois Introduction
More informationINVERSE PROBLEMS IN X-RAY SCIENCE STEFANO MARCHESINI
INVERSE PROBLEMS IN X-RAY SCIENCE STEFANO MARCHESINI overview of photon science inverse problems linear (under-determined): tomography non-linear: diffraction methods Inverse problems and the data deluge
More informationThe Phase Problem: A Mathematical Tour from Norbert Wiener to Random Matrices and Convex Optimization
The Phase Problem: A Mathematical Tour from Norbert Wiener to Random Matrices and Convex Optimization Thomas Strohmer Department of Mathematics University of California, Davis Norbert Wiener Colloquium
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 informationMulti-frequency phase retrieval from noisy data
Multi-frequency phase retrieval from noisy data Vladimir Katkovnik and Karen Egiazarian Laboratory of Signal Processing, Tampere University of Technology, P.O. Box 553 FI-33101 Tampere, Finland e-mails:
More informationA Tutorial on Compressive Sensing. Simon Foucart Drexel University / University of Georgia
A Tutorial on Compressive Sensing Simon Foucart Drexel University / University of Georgia CIMPA13 New Trends in Applied Harmonic Analysis Mar del Plata, Argentina, 5-16 August 2013 This minicourse acts
More informationNonconvex Methods for Phase Retrieval
Nonconvex Methods for Phase Retrieval Vince Monardo Carnegie Mellon University March 28, 2018 Vince Monardo (CMU) 18.898 Midterm March 28, 2018 1 / 43 Paper Choice Main reference paper: Solving Random
More informationMedian-Truncated Gradient Descent: A Robust and Scalable Nonconvex Approach for Signal Estimation
Median-Truncated Gradient Descent: A Robust and Scalable Nonconvex Approach for Signal Estimation Yuejie Chi, Yuanxin Li, Huishuai Zhang, and Yingbin Liang Abstract Recent work has demonstrated the effectiveness
More informationThree Recent Examples of The Effectiveness of Convex Programming in the Information Sciences
Three Recent Examples of The Effectiveness of Convex Programming in the Information Sciences Emmanuel Candès International Symposium on Information Theory, Istanbul, July 2013 Three stories about a theme
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationIEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 11, NOVEMBER On the Performance of Sparse Recovery
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 11, NOVEMBER 2011 7255 On the Performance of Sparse Recovery Via `p-minimization (0 p 1) Meng Wang, Student Member, IEEE, Weiyu Xu, and Ao Tang, Senior
More informationPhase Retrieval Meets Statistical Learning Theory: A Flexible Convex Relaxation
hase Retrieval eets Statistical Learning Theory: A Flexible Convex Relaxation Sohail Bahmani Justin Romberg School of ECE Georgia Institute of Technology {sohail.bahmani,jrom}@ece.gatech.edu Abstract We
More informationPHASE RETRIEVAL FROM LOCAL MEASUREMENTS: IMPROVED ROBUSTNESS VIA EIGENVECTOR-BASED ANGULAR SYNCHRONIZATION
PHASE RETRIEVAL FROM LOCAL MEASUREMENTS: IMPROVED ROBUSTNESS VIA EIGENVECTOR-BASED ANGULAR SYNCHRONIZATION MARK A. IWEN, BRIAN PRESKITT, RAYAN SAAB, ADITYA VISWANATHAN Abstract. We improve a phase retrieval
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 informationPhase Retrieval via Matrix Completion
SIAM J. IMAGING SCIENCES Vol. 6, No. 1, pp. 199 225 c 2013 Society for Industrial and Applied Mathematics Phase Retrieval via Matrix Completion Emmanuel J. Candès, Yonina C. Eldar, Thomas Strohmer, and
More informationEstimating Unknown Sparsity in Compressed Sensing
Estimating Unknown Sparsity in Compressed Sensing Miles Lopes UC Berkeley Department of Statistics CSGF Program Review July 16, 2014 early version published at ICML 2013 Miles Lopes ( UC Berkeley ) estimating
More informationVectorial Phase Retrieval of 1-D Signals
1632 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 61, NO. 7, APRIL 1, 2013 Vectorial Phase Retrieval of 1-D Signals Oren Raz, Nirit Dudovich, Boaz Nadler, Member, IEEE Abstract Reconstruction of signals
More informationAn algebraic approach to phase retrieval
University of Michigan joint with Aldo Conca, Dan Edidin, and Milena Hering. I learned about frame theory from... AIM Workshop: Frame theory intersects geometry http://www.aimath.org/arcc/workshops/frametheorygeom.html
More information