Phase Retrieval with Random Illumination
|
|
- Justina Price
- 5 years ago
- Views:
Transcription
1 Phase Retrieval with Random Illumination Albert Fannjiang, UC Davis Collaborator: Wenjing Liao NCTU, July 212 1
2 Phase retrieval Problem Reconstruct the object f from the Fourier magnitude Φf. Why do we care? X-ray crystallography, single-molecule imaging, astronomy etc Nobel Prize in Chemistry for Hauptman and Karle: partial solution of phasing problem. 2
3 Phasing problem formulation Discrete finite objects Let n =(n 1,n 2 ) Z 2 and z =(z 1,z 2 ) C 2. multi-index : z n = z n 1 1 zn 2 2 Let the object be represented by f(n), n N =(N, N). Binary objects: white = 1, black =. 3
4 f L = Lena F L (w) = F L (w) e iθ L(w) f B = Barbara F B (w) = F B (w) e iθ B(w) F 1 (w) = F B (w) e iθ L(w) F 2 (w) = F L (w) e iθ B(w) f 1 = Φ F 1 f 2 = Φ F 2 4
5 F 1 (w) = F B (w) e iθ L(w) F 2 (w) = F L (w) e iθ B(w) f 1 = Φ F 1 f 2 = Φ F 2 5
6 Phase retrieval Discrete finite objects Let n =(n 1,n 2 ) Z 2 and z =(z 1,z 2 ) C 2. multi-index : z n = z n 1 1 zn 2 2 Let the object be represented by f(n), n N =(N, N). Fourier transform describes wave propagation F (e i2πw 1,e i2πw 2)= n f(n)e i2πn w Analytic continuation = z-transform F (z) = n f(n)z n. Discrete phase retrieval problem: Determine f(n) from Fourier magnitude data F (w), w =(e i2πw 1,e i2πw 2) [, 1] 2 6
7 Fourier magnitude data: F (w) 2 = N n= N m f(m + n)f (m)e i2πn w = N n= N C f (n)e i2πn w where C f (n) = m f(m + n)f (m) is the autocorrelation function of f. Fourier magnitude data contain complete information about autocorrelation function. Sampling Theorem: supp(c f ) [ N, N] 2 = [, 1] 2 is reduced to the Nyquist grid M = (k 1,k 2 ):k j =, 1 2N +1, 2 2N +1,, 2N 2N +1 7
8 Oversampling ratio σ = Fourier magnitude data number unknown image pixel number Standard ratio: Compressed sensing: 8
9 Trivial ambiguities Autocorrelation: Invariant under: (i) global phase, C f (n) = m+n N f(m + n)f (m) f(n) e iθ f(n), for some θ [, 2π], (ii) spatial translation f(n) f(n m), n m = n + m mod(n), some m Z 2 (iii) conjugate inversion (twin image) f(n) f (N n). 9
10 Nontrivial ambiguity THEOREM (Hayes 82, Pitts-Greenleaf 3) Let the z-transform F (z) of a finite complex-valued sequence {f(n)} be given by F (z) =αz m p k=1 F k (z), m N 2, α C where F k,k =1,...,p are nontrivial irreducible polynomials. Let G(z) be the z-transform of another finite sequence g(n). Suppose F (w) = G(w), w [, 1] 2. Then G(z) must have the form G(z) = α e iθ z p F k (z) Fk (1/z ), p N 2, θ R k I k I c where I is a subset of {1, 2,...,p}. Nontrivial ambiguity: Partial conjugate inversion on factors. 1
11 Random illumination Random Illumination P.F. Almoro et al. / Optics and Lasers in En cs and Lasers in Engineering 49 (211) Coded aperture imaging expanding the laser beam and impinging on a plate of translucent perpex, which acts as an opal diffuser, with unnoticeable grain and nearly Lambertian scattering of the light. In the fluorescence experiments the sample is coated with a thin layer (~5 µm) of solution of fluorescein diacetate (FDA) that reemits incoherent light in the green wavelengths of the optical spectrum. The process requires a high resolution image of the speckle that acts as the encodingdecoding mask. We take these reference images prior to each experiment by focusing at a transparent region in the sample plane using a lens with high NA (.4). Figure 2 displays the Diffuser generated pattern: Garcia-Zalevsky-Fixler 5 reference and its speckle autocorrelation. The size of the autocorrelation peak is the expected d: (a) plane wave orimage uniform illumination; (b) diffused illumination and (c mental setups for the after multiple-plane phase retrieval method: (a) plane wave resolution the superresolution process when a low NA lens is used. (a) (b) Fig. 2. (a) Encoding speckle pattern. (b) Autocorrelation of the encoding pattern. Once the reference speckle pattern is captured, 11 the sample is set in place and the lens is replaced by a lens with a low NA in the horizontal direction. Then the sample position is laterally scanned and the image set is captured. Note that instead of displacing the projected
12 X-ray holograhy with URA (Marchesini et al. 8). Coherent X-ray imaging will be a key technique for developing nanoscience and nanotechnology. One shot imaging: Bright and ultrashort X-ray pulse vaporizes the sample right after the pulse passes the sample. Random illumination amounts to replacing the original object f(n) by f(n) = f(n)λ(n) (illuminated object) 12
13 Random illumination f(n) =f(n)λ(n) (illuminated object) λ(n), representing the illumination field, is a known sequence of samples of random variables. Let λ(n) be continuous random variables with respect to the Lebesque measure on S 1 (the unit circle), R or C. Case of S 1 can be facilitated by a random phase modulator with λ(n) =e iφ(n) where φ(n) are continuous random variables on [, 2π]. Case of R: random amplitude modulator. Case of C: both phase and amplitude modulations. 13
14 Irreducibility THEOREM. Suppose that the support of the object {f(n)} has rank 2. Then the the z-transform of the illuminated object f(n)λ(n) is irreducible with probability one. False for rank 1 objects: fundamental thm of algebra σ = Fourier magnitude data number unknown image pixel number = 14
15 Absolute uniqueness Positivity THEOREM If f(n) is real and nonnegative for every n then, with probability one, f is determined absolutely uniquely by the Fourier magnitude measurement on the lattice L. Sector constraint THEOREM Suppose the phases of the object belong to [a, b] [, 2π]. Then the solution to the Fourier phasing problem has auniquesolutionwithprobabilityexponentiallyclosetounity (depending on the sparsity and the phase range b a.) 15
16 Complex objects w/o constraint THEOREM. Suppose that {λ 1 (n)} are i.i.d. continuous random variables with respect to the Lebesgue measure on S 1, R or C and in addition either one of the following conditions is true. (i) {λ 2 (n)} are i.i.d. continuous random variables with respect to the Lebesgue measure on S 1, R or C and {λ 2 (n)} are independent of {λ 1 (n)}. (ii) {λ 2 (n)} are deterministic. Then with probability one f(n) is uniquely determined, up to a constant phase factor, by the Fourier magnitude measurements with two illuminations λ 1 and λ 2. 16
17 Alternated projections Gerchberg-Saxton; Error Reduction (Fienup) f k+1 /f k Λ g k Φ G k P o T f k Λ 1 g k Φ 1 G k Oversampling method (Miao et al.) σ = image pixel number + zero-padding pixel number image pixel number. 17
18 Zero-padding cameraman 18
19 Multiple illuminations P 1 = Λ 1 1 Φ 1 T 1 ΦΛ 1 P 2 = Λ 1 2 Φ 1 T 2 ΦΛ 2. f k+1 = P o P 2 P 1 f k. 19
20 Error Reduction (Gerchberg-Saxton) Bregman 65: convex constraints = convergence to afeasible solution. Fourier magnitude data are a non-convex constraint! Nonconvexity or nonuniqueness? 2
21 Convergence Theorem 4. Let f C(N ) be an array with f() = and rank 2. Let λ(n) be i.i.d. continuous random variables on S 1. Let the Fourier magnitude be sampled on L. Let h be a fixed point of P o Pf θ such that Pf θ h satisfies the zero-padding condition. (a) If f is real-valued, h = ±f with probability one, (b) If f satisfies the sector condition of Theorem 2, then h = e iν f, for some ν, and satisfies the same sector constraint with probability at least 1 N (β α) S/2 (2π) S/2. 21
22 Hybrid-Input-Output (HIO) f k+1 /f k Λ g k Φ G k P o f k Λ 1 g k Φ 1 T G k (f k+1 (n)) = (f k(n)) (f k+1 (n)) = (f k (n)) β (f k(n)), Real-valued objects 22
23 Error metrics Relative error e( ˆf) = f ˆf/f min f ν [,2π) eiν ˆf/f Relative residual r( ˆf) = Y ΦΛP o{ ˆf} Y 23
24 (a) (b) 269 x x 2 24
25 at each at each (a) (b) (c) (f) (g) (h) at each at each at each.5 at each (i) 5 6 (d).5 (e) 4 (j) (k) (l) (e)-(h) Low resolution 4 x 4 block illumination with OR=2 Figure 1. (a)recovered cameraman by 5 ER s with single uniform illumina(i)-(l) High resolution illumination with OR=1 tion when ρ = 4. f 5 f4999 / f4999.2%, ε f (f5 ) 5.47%. (b)normalized error ε (fk ) at each. (c)recovered cameraman by 1 HIO +5 ER with one 25 uniform illumination when ρ = 4. f15 f149 / f149.3%, ε f (f15 ).93%. (d) ε (fk ) at each. (e)recovered cameraman by 2 ER iter-
26 .5 at each.5 at each (a) (b) (c) (d).5 at each 1 at each (e) (f) (g) (h).5 at each.5 at each (i) (j) (k) (l) (e) - (h) Low resolution 4 x 4 block illumination with OR=2 Figure 2. (a)recovered phantom by 5 ER s with single uniform illumination when (i) ρ - = (l) 4. High f 5 resolution f 4999 / f 4999 illumination.6%, ε f (f 5 with ) OR=1 14.7%. (b)normalized error ε(f k ) at each. (c)recovered phantom by 1 HIO +5 ER with single uniform illumination when ρ = 4. f f 149 / f 149.7%, ε f (f 15 ) 3.95%. (d) ε(f k ) at each. (e)recovered phantom by 5 ER s with single low resolution(block size: 4 4) random phase illumination when ρ = 2.
27 .5 relative residual at each.5 relative residual at each relative residual relative residual (a) (b) (c) (d).5 relative residual at each.5 relative residual at each relative residual relative residual (e) (f) (g) (h).5 relative residual at each.5 relative residual at each relative residual relative residual (i) (j) (k) (l) High resolution illumination with 5% Gaussian, Poisson and illuminator errors 27
28 .5 at each 1 at each (a) (b) (c) (d).5 at each 1 at each (e) (f) (g) (h).5 at each 1 at each (i) (j) (k) (l) Low resolution illumination with 5% Gaussian, Poisson and illuminator errors 28
29 Compressed measurement average relative error in 5 trials relative error versus oversampling rate high resolution & noise free low resolution & noise free high resolution & 5% gaussian nosie low resolution & 5% gaussian nosie high resolution & 5% poisson nosie low resolution & 5% poisson noise high resolution & 5% illuminator noise low resolution & 5% illuminator nosie average relative error in 5 trials relative error versus oversampling rate high resolution & noise free low resolution & noise free high resolution & 5% gaussian nosie low resolution & 5% gaussian nosie high resolution & 5% poisson nosie low resolution & 5% poisson noise high resolution & 5% illuminator noise low resolution & 5% illuminator nosie average relative error in 5 trials relative error versus oversampling rate high resolution & noise free low resolution & noise free high resolution & 5% gaussian nosie low resolution & 5% gaussian nosie high resolution & 5% poisson nosie low resolution & 5% poisson noise high resolution & 5% illuminator noise low resolution & 5% illuminator nosie oversampling rate oversampling rate oversampling rate (a) (b) (c) Figure 3 Figure 4. (a)average (a) relative real-valued error by adaptive HIO + adaptive ER in 5 trials versus oversampling rate for phantom with one random phase illumination. (b)average relative error by adaptive HIO (b) + adaptive positive ER in real 5 trials & versus imaginary oversampling parts rate for phantom with random phases between (c) no andconstraint π/2 with one random phase illumination. (c)average relative error by 2 HIO +3 ER in 5 trials versus oversampling rate for phantom with random phases between and 2π with one uniform illumination and one random phase illumination. 29
30 Noise stability 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 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 high resolution & gaussian noise low resolution & gaussian noise high resolution & poisson noise low resolution & poisson noise high resolution & illuminator noise low resolution & illuminator noise noise percentage noise percentage noise percentage (a) (b) (c) Figure 6 Figure 7. (a)average relative error by adaptive HIO + adaptive ER in 5 trials versus noise level for phantom with (a) real-valued one randomobjects phase illumination when ρ = 4. (b)average relative error by adaptive HIO + adaptive ER in 5 trials versus noise level for phantom with random phases between(b) and positive π/2 withreal one random & imaginary phase illumination partswhen ρ = 4. (c)average relative error by (c) 2 no HIO constraint +3 ER in 5 trials versus noise level for phantom with random phases between and 2π with one uniform illumination and one random phase illumination when ρ = 3. 3
31 Conclusions Random illumination as enabling tool for phase retrieval. Absolute uniqueness Fast convergence OR =1 (real) or 2 (complex) Proof of convergence: HIO? References: A. Fannjiang Absolute uniqueness in phase retrieval with random illumination Inverse Problems 212 ( arxiv: ) A. Fannjiang and W. Liao Phase retrieval with random phase illumination arxiv:
Fourier phase retrieval with random phase illumination
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 Classical
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 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 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 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 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 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 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 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 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 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 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 informationLaser Speckle and Applications in Optics
Laser Speckle and Applications in Optics M. FRANCON Optics Laboratory Faculty of Sciences University of Paris Paris, France Translated by HENRI H. ARSENAULT Department of Physics Laval University Quebec,
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 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 informationDensity Field Measurement by Digital Laser Speckle Photography
Density Field Measurement by Digital Laser Speckle Photography by M. Kawahashi and H. Hirahara Saitama University Department of Mechanical Engineering Shimo-Okubo 255, Urawa, Saitama, 338-8570, Japan ABSTRACT
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 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 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 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 informationWhat is a Quasicrystal?
July 23, 2013 Rotational symmetry An object with rotational symmetry is an object that looks the same after a certain amount of rotation. Rotational symmetry An object with rotational symmetry is an object
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 informationGeneral theory of diffraction
General theory of diffraction X-rays scatter off the charge density (r), neutrons scatter off the spin density. Coherent scattering (diffraction) creates the Fourier transform of (r) from real to reciprocal
More informationApplication of optimization technique to noncrystalline x-ray diffraction microscopy: Guided hybrid input-output method
Application of optimization technique to noncrystalline x-ray diffraction microscopy: Guided hybrid input-output method Chien-Chun Chen, 1 Jianwei Miao, 2 C. W. Wang, 3 and T. K. Lee 1 1 Institute of Physics,
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 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 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 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 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 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 informationPhase Retrieval for the Hubble Space Telescope and other Applications Abstract: Introduction: Theory:
Phase Retrieval for the Hubble Space Telescope and other Applications Stephanie Barnes College of Optical Sciences, University of Arizona, Tucson, Arizona 85721 sab3@email.arizona.edu Abstract: James R.
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 informationOn the FPA infrared camera transfer function calculation
On the FPA infrared camera transfer function calculation (1) CERTES, Université Paris XII Val de Marne, Créteil, France (2) LTM, Université de Bourgogne, Le Creusot, France by S. Datcu 1, L. Ibos 1,Y.
More informationContents Preface iii 1 Origins and Manifestations of Speckle 2 Random Phasor Sums 3 First-Order Statistical Properties
Contents Preface iii 1 Origins and Manifestations of Speckle 1 1.1 General Background............................. 1 1.2 Intuitive Explanation of the Cause of Speckle................ 2 1.3 Some Mathematical
More informationPolarization Shearing Interferometer (PSI) Based Wavefront Sensor for Adaptive Optics Application. A.K.Saxena and J.P.Lancelot
Polarization Shearing Interferometer (PSI) Based Wavefront Sensor for Adaptive Optics Application A.K.Saxena and J.P.Lancelot Adaptive Optics A Closed loop Optical system to compensate atmospheric turbulence
More informationOptics for Engineers Chapter 11
Optics for Engineers Chapter 11 Charles A. DiMarzio Northeastern University Nov. 212 Fourier Optics Terminology Field Plane Fourier Plane C Field Amplitude, E(x, y) Ẽ(f x, f y ) Amplitude Point Spread
More information2.71. Final examination. 3 hours (9am 12 noon) Total pages: 7 (seven) PLEASE DO NOT TURN OVER UNTIL EXAM STARTS PLEASE RETURN THIS BOOKLET
2.71 Final examination 3 hours (9am 12 noon) Total pages: 7 (seven) PLEASE DO NOT TURN OVER UNTIL EXAM STARTS Name: PLEASE RETURN THIS BOOKLET WITH YOUR SOLUTION SHEET(S) MASSACHUSETTS INSTITUTE OF TECHNOLOGY
More informationComputational Nanoscience
Computational Nanoscience Applications for Molecules, Clusters, and Solids KALMÄN VARGA AND JOSEPH A. DRISCOLL Vanderbilt University, Tennessee Щ CAMBRIDGE HP UNIVERSITY PRESS Preface Part I One-dimensional
More informationFoundations of Image Science
Foundations of Image Science Harrison H. Barrett Kyle J. Myers 2004 by John Wiley & Sons,, Hoboken, 0-471-15300-1 1 VECTORS AND OPERATORS 1 1.1 LINEAR VECTOR SPACES 2 1.1.1 Vector addition and scalar multiplication
More informationOptics for Engineers Chapter 11
Optics for Engineers Chapter 11 Charles A. DiMarzio Northeastern University Apr. 214 Fourier Optics Terminology Apr. 214 c C. DiMarzio (Based on Optics for Engineers, CRC Press) slides11r1 1 Fourier Optics
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 informationDiffraction Imaging with Coherent X-rays
Diffraction Imaging with Coherent X-rays John Miao Stanford Synchrotron Radiation Laboratory Stanford Linear Accelerator Center The Phase Problem: A Coherence Effect Detector Coherent X-rays Atoms The
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 informationhttp://topex.ucsd.edu/gmtsar amplitude and phase coherence and pixel matching resolution: optical vs. microwave D s = 2H sinθ r = 2H λ L H = 800km. Optical : L = 1m λ = 0.5µm D s = 0.8m Microwave : L
More informationToday. MIT 2.71/2.710 Optics 11/10/04 wk10-b-1
Today Review of spatial filtering with coherent illumination Derivation of the lens law using wave optics Point-spread function of a system with incoherent illumination The Modulation Transfer Function
More informationLecture 19 Optical MEMS (1)
EEL6935 Advanced MEMS (Spring 5) Instructor: Dr. Huikai Xie Lecture 19 Optical MEMS (1) Agenda: Optics Review EEL6935 Advanced MEMS 5 H. Xie 3/8/5 1 Optics Review Nature of Light Reflection and Refraction
More informationWeak-Beam Dark-Field Technique
Basic Idea recall bright-field contrast of dislocations: specimen close to Bragg condition, s î 0 Weak-Beam Dark-Field Technique near the dislocation core, some planes curved to s = 0 ) strong Bragg reflection
More informationJRE Group of Institutions ASSIGNMENT # 1 Special Theory of Relativity
ASSIGNMENT # 1 Special Theory of Relativity 1. What was the objective of conducting the Michelson-Morley experiment? Describe the experiment. How is the negative result of the experiment interpreted? 2.
More informationX-Ray Diffraction as a key to the Structure of Materials Interpretation of scattering patterns in real and reciprocal space
X-Ray Diffraction as a key to the Structure of Materials Interpretation of scattering patterns in real and reciprocal space Tobias U. Schülli, X-ray nanoprobe group ESRF OUTLINE 1 Internal structure of
More informationMERIS US Workshop. Instrument Characterization Overview. Steven Delwart
MERIS US Workshop Instrument Characterization Overview Steven Delwart Presentation Overview On-Ground Characterisation 1. Diffuser characterisation 2. Polarization sensitivity 3. Optical Transmission 4.
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 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 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 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 informationLab 3-4 : Confocal Microscope Imaging of Single-Emitter Fluorescence and Hanbury-Brown and Twiss Set Up, Photon Antibunching
Lab 3-4 : Confocal Microscope Imaging of Single-Emitter Fluorescence and Hanbury-Brown and Twiss Set Up, Photon Antibunching Mongkol Moongweluwan 1 1 Department of Physics and Astronomy, University of
More informationPhase Contrast. Zernike Phase Contrast (PC) Nomarski Differential Interference Contrast (DIC) J Mertz Boston University
Phase Contrast Zernike Phase Contrast (PC) Nomarski Differential Interference Contrast (DIC) J Mertz Boston University Absorption Scattering / n wave-number k 2 n = index of refraction / n wave-number
More informationThere and back again A short trip to Fourier Space. Janet Vonck 23 April 2014
There and back again A short trip to Fourier Space Janet Vonck 23 April 2014 Where can I find a Fourier Transform? Fourier Transforms are ubiquitous in structural biology: X-ray diffraction Spectroscopy
More informationMichelson Interferometer
Michelson Interferometer Objective Determination of the wave length of the light of the helium-neon laser by means of Michelson interferometer subsectionprinciple and Task Light is made to produce interference
More informationOPSE FINAL EXAM Fall 2016 YOU MUST SHOW YOUR WORK. ANSWERS THAT ARE NOT JUSTIFIED WILL BE GIVEN ZERO CREDIT.
CLOSED BOOK. Equation Sheet is provided. YOU MUST SHOW YOUR WORK. ANSWERS THAT ARE NOT JUSTIFIED WILL BE GIVEN ZERO CREDIT. ALL NUMERICAL ANSERS MUST HAVE UNITS INDICATED. (Except dimensionless units like
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 informationPhys 531 Lecture 27 6 December 2005
Phys 531 Lecture 27 6 December 2005 Final Review Last time: introduction to quantum field theory Like QM, but field is quantum variable rather than x, p for particle Understand photons, noise, weird quantum
More informationThe Plasma Phase. Chapter 1. An experiment - measure and understand transport processes in a plasma. Chapter 2. An introduction to plasma physics
The Plasma Phase Chapter 1. An experiment - measure and understand transport processes in a plasma Three important vugraphs What we have just talked about The diagnostics Chapter 2. An introduction to
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 informationHigh resolution traction force microscopy on small focal adhesions improved accuracy through optimal marker distribution and optical flow tracking
Scientific Reports High resolution traction force microscopy on small focal adhesions improved accuracy through optimal marker distribution and optical flow tracking Claude N. Holenstein 1,2,*, Unai Silvan
More informationElectricity & Optics
Physics 24100 Electricity & Optics Lecture 26 Chapter 33 sec. 1-4 Fall 2017 Semester Professor Koltick Interference of Light Interference phenomena are a consequence of the wave-like nature of light Electric
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 informationLecture 20. Wind Lidar (2) Vector Wind Determination
Lecture 20. Wind Lidar (2) Vector Wind Determination Vector wind determination Ideal vector wind measurement VAD and DBS technique for vector wind Coherent versus incoherent Detection Doppler wind lidar
More informationConstructive vs. destructive interference; Coherent vs. incoherent interference
Constructive vs. destructive interference; Coherent vs. incoherent interference Waves that combine in phase add up to relatively high irradiance. = Constructive interference (coherent) Waves that combine
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More informationThe Photon Counting Histogram: Statistical Analysis of Single Molecule Populations
The Photon Counting Histogram: Statistical Analysis of Single Molecule Populations E. Gratton Laboratory for Fluorescence Dynamics University of California, Irvine Transition from FCS The Autocorrelation
More informationCOMPUTATIONAL IMAGING. Berthold K.P. Horn
COMPUTATIONAL IMAGING Berthold K.P. Horn What is Computational Imaging? Computation inherent in image formation What is Computational Imaging? Computation inherent in image formation (1) Computing is getting
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 informationFourier Series. Combination of Waves: Any PERIODIC function f(t) can be written: How to calculate the coefficients?
Fourier Series Combination of Waves: Any PERIODIC function f(t) can be written: How to calculate the coefficients? What is the Fourier Transform It is the representation of an arbitrary NON-periodic signal
More informationControl of dispersion effects for resonant ultrashort pulses M. A. Bouchene, J. C. Delagnes
Control of dispersion effects for resonant ultrashort pulses M. A. Bouchene, J. C. Delagnes Laboratoire «Collisions, Agrégats, Réactivité», Université Paul Sabatier, Toulouse, France Context: - Dispersion
More informationX-ray Intensity Fluctuation Spectroscopy. Mark Sutton McGill University
X-ray Intensity Fluctuation Spectroscopy Mark Sutton McGill University McGill University Collaborators J-F. Pelletier K. Laaziri K. Hassani A. Fluerasu E. Dufresne G. Brown M. Grant Yale/MIT S. Mochrie
More informationSpeed of Light in Glass
Experiment (1) Speed of Light in Glass Objective:- This experiment is used to determine the speed of propagation of light waves in glass. Apparatus:- Prism, spectrometer, Halogen lamp source. Theory:-
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 informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/2//e50054/dc Supplementary Materials for Two-photon quantum walk in a multimode fiber Hugo Defienne, Marco Barbieri, Ian A. Walmsley, Brian J. Smith, Sylvain Gigan
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 informationSUPPLEMENTARY INFORMATION
scos SL PLB2 DOPC PL P BS PLB1 CL SF PL ND PLS O O BD PT IF DBS DD BS L3 CL SF PL ND BD BD f 0 OI L2 UST L1 S OC Quality assurance arm CCD f 0 + f US AO 50 Hz f 0 + f US AO 50 Hz PBS HWP BD RL RDD CL SF
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 informationFOURIER TRANSFORM METHODS David Sandwell, January, 2013
1 FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1. Fourier Transforms Fourier analysis is a fundamental tool used in all areas of science and engineering. The fast fourier transform (FFT) algorithm
More informationHigh resolution tomographic diffraction microscopy
High resolution tomographic diffraction microscopy J. Girard,Y. Ruan 1, E. Mudry 1, F. Drsek G. Maire 1, P. Chaumet 1, H. Giovannini 1,K. Belkebir 1, A. Talneau 2, A. Sentenac 1 1 Institut Fresnel (Marseille)
More informationComplex data mapping on a binary ferroelectric liquid crystal electrically addressed spatial light modulator for target recognition
Indian Journal of Pure & Applied Physics Vol. 45, August 2007, pp. 647-653 Complex data mapping on a binary ferroelectric liquid crystal electrically addressed spatial light modulator for target recognition
More informationNear Optimal Signal Recovery from Random Projections
1 Near Optimal Signal Recovery from Random Projections Emmanuel Candès, California Institute of Technology Multiscale Geometric Analysis in High Dimensions: Workshop # 2 IPAM, UCLA, October 2004 Collaborators:
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 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 informationHigh-Resolution. Transmission. Electron Microscopy
Part 4 High-Resolution Transmission Electron Microscopy 186 Significance high-resolution transmission electron microscopy (HRTEM): resolve object details smaller than 1nm (10 9 m) image the interior of
More informationMultifrequency Interferometric Imaging with Intensity-Only Measurements
SIAM J. IMAGING SCIENCES Vol., No. 3, pp. 5 32 c 7 Society for Industrial and Applied Mathematics Multifrequency Interferometric Imaging with Intensity-Only Measurements Miguel Moscoso, Alexei Novikov,
More informationPH 222-3A Spring 2010
PH -3A Spring 010 Interference Lecture 6-7 Chapter 35 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition) 1 Chapter 35 Interference The concept of optical interference is critical to understanding
More informationOptical Production of the Husimi Function of Two Gaussian Functions
Applied Mathematics & Information Sciences (3) (008), 309-316 An International Journal c 008 Dixie W Publishing Corporation, U. S. A. Optical Production of the Husimi Function of Two Gaussian Functions
More informationFourier phase analysis in radio-interferometry
Fourier phase analysis in radio-interferometry François Levrier Ecole Normale Supérieure de Paris In collaboration with François Viallefond Observatoire de Paris Edith Falgarone Ecole Normale Supérieure
More informationName :. Roll No. :... Invigilator s Signature :.. CS/B. Tech (New)/SEM-1/PH-101/ PHYSICS-I
Name :. Roll No. :..... Invigilator s Signature :.. CS/B. Tech (New)/SEM-1/PH-101/2011-12 2011 PHYSICS-I Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks. Candidates
More informationHighly Efficient and Anomalous Charge Transfer in van der Waals Trilayer Semiconductors
Highly Efficient and Anomalous Charge Transfer in van der Waals Trilayer Semiconductors Frank Ceballos 1, Ming-Gang Ju 2 Samuel D. Lane 1, Xiao Cheng Zeng 2 & Hui Zhao 1 1 Department of Physics and Astronomy,
More informationDynamic Light Scattering Employing One- and Two-Dimensional Detectors and Different Time-Correlation Approaches
Dynamic Light Scattering Employing One- and Two-Dimensional Detectors and Different Time-Correlation Approaches Ralf Klemt Summer Student Project FS-CXS September 10, 2015 Ralf Klemt Dynamic Light Scattering
More informationProbing the orbital angular momentum of light with a multipoint interferometer
CHAPTER 2 Probing the orbital angular momentum of light with a multipoint interferometer We present an efficient method for probing the orbital angular momentum of optical vortices of arbitrary sizes.
More informationCorrelation functions in optics and quantum optics, 4
Correlation functions in optics and quantum optics, 4 State Key Laboratory of Precision Spectroscopy East China Normal University, Shanghai, China Luis A. Orozco www.jqi.umd.edu The slides of the course
More informationSupporting Information
Supporting Information Devlin et al. 10.1073/pnas.1611740113 Optical Characterization We deposit blanket TiO films via ALD onto silicon substrates to prepare samples for spectroscopic ellipsometry (SE)
More informationInterference- Michelson Interferometer. Interference lecture by Dr. T.Vishwam
Interference- Michelson Interferometer Interference lecture by Dr. T.Vishwam * Measurement of the coherence length of a spectral line * Measurement of thickness of thin transparent flakes * Measurement
More information