Reconstruction from Digital Holograms by Statistical Methods
|
|
- Jesse Todd
- 5 years ago
- Views:
Transcription
1 Reconstruction from Digital Holograms by Statistical Methods Saowapak Sotthivirat Jeffrey A. Fessler EECS Department The University of Michigan 2003 Asilomar Nov. 12, 2003 Acknowledgements: Brian Athey, Emmett Leith, Kurt Mills 0
2 Outline Background on holography Image model Conventional numerical reconstruction Statistical holographic reconstruction Simulation results Conclusions and future work 1
3 Conventional Holography Record both amplitude and phase of a wave field Hologram = recorded interference pattern between an object beam and a reference beam object object beam recording medium Recording reference beam primary image hologram conjugate image Reconstructing reference beam 2
4 Image Plane Holography object P 1 P 2 incoherent source. object beam x 1 x 2.. hologram x 3. reference beam x 1 x 2. P 1 P 1 Optical sectioning property as in confocal microscopy No xy scanning Larger source decreases coherence and improves spatial resolution 3
5 Challenges Finite source size limits spatial resolution Optical reconstruction inconvenient Use statistical image reconstruction Hologram Numerical reconstruction on a computer Statistical image restoration Hologram Statistical image reconstruction 4
6 Measurement Model Source Object f( r )= Σ χ ( r ) j x j j Optical system h Object beam u o Reference beam u ref Record samples of + u + 2 o u ref on a digital detector Object beam, superposition: u o ( r )= h( r, r) f ( r)d r Discretized object model u o ( r i )=[Ax] i, Hologram intensity at recording plane I( r) = u o ( r)+u ref ( r) 2 a ij = h( r i, r)χ j ( r)d r = u o ( r) 2 + u ref ( r) 2 + u o ( r)u ref( r)+u o( r)u ref ( r) Mean of ith measured sample E[Y i ]=I( r i )+b i background dark current = u o ( r i )+u ref ( r i ) 2 + b i = [Ax] i + u i 2 + b i, i = 1,...,N 5
7 Conventional Reconstruction Method Plane-wave assumption: u ref ( r)=u ref e ı2π r α Hologram intensity I( r)= u o ( r) 2 +U 2 ref +U ref u o ( r)e ı2π r α +U ref u o( r)e ı2π r α Fourier transform of hologram intensity I ( ν)=i o ( ν)+u 2 refδ( ν)+u ref U o ( ν α)+u ref U o( ν α) r 2 ν 2 r 2 Zero order spectrum r 1 IFFT ν 1 r 1 Recorded Hologram First order spectra α α Magnitude of Fourier Transform Reconstructed Holographic Image 6
8 Statistical Model Poisson model of the hologram measurement data: } Y i Poisson { [Ax] i + u i 2 + b i, i = 1,...,N [Ax] i = P j=1 a ij x j Known: imaging system A, data {Y i }, reference beam {u i }, dark current offset {b i } Cost function: Ψ(x)=L(x)+V (x) Negative likelihood function: L(x)= N i=1 h i ([Ax] i ) h i (z)=( z + u i 2 + b i ) y i log( z + u i 2 + b i ) Penalty: V (x)=β r i=1 ψ([cx] i ) Penalized-likelihood (aka MAP) reconstruction: x = argminψ(x) x No closed-form solution need an iterative algorithm 7
9 Optimization Transfer Illustration 8
10 Optimization Transfer Principle S-step: find a majorizing surrogate function φ (n) 1. φ (n)( x (n)) = Ψ M-step: minimize φ (n) instead of Ψ (x (n)) 2. φ (n) (x) Ψ(x), x C P x (n+1) = argminφ (n) (x) x Fact: algorithm monotonically decreases cost function (cf. Newton) Ψ(x (n+1) ) Ψ(x (n) ) φ(x;x n) Φ(x) x ^ x n+1 x n x 9
11 Log-likelihood Surrogate Functions h (l R ;0) h R (l R ;l n ) q R (l R ;l n ) z (n) i N L(x)= h i ([Ax] i ) (independence) i=1 h i (z) =( z + u i 2 + b i ) y i log( z + u i 2 + b i ) q i (z;z (n) i )=h i (z (n) i )+ q i (z (n) i )(z z (n) i )+ 1 2 c(n) i z z (n) i l R =[Ax (n) ] i. Curvature: c (n) i = c(y i,u i,b i,[ax (n) ] i ) (JOSAA, in review) l n 2 10
12 S-step: L(x)= Quadratic Surrogate Function N h i ([Ax] i ) i=1 N q i ([Ax] i ;[Ax (n) ] i ) = Q(x; x (n) ) i=1 Q(x; x (n) ) = N q i ([Ax] i ;[Ax (n) ] i ) i=1 = L(x (n) )+ L(x (n) )(x x (n) )+ 1 2 (x x(n) ) A DA(x x (n) ) { } D = diag M-step: x (n+1) = argminq(x; x (n) )+V(x) x C P Surrogate Q is quadratic, so minimize by PCG algorithm. Quadratic surrogates for regularizing penalty function also available. (cf. half quadratic methods) c (n) i 11
13 Separable Quadratic Surrogate Algorithm Convexity of q i allows separable surrogate (De Pierro, 1995): ( P [ q i ([Ax] i ;z (n) [aij i )=q i π (x j x (n) ]) j )] ij +[Ax (n) ] i ; z (n) i j=1 π ij ) π ij 0 and P j=1 π ij = 1 P j=1 ( [aij π ij q (x j x (n) j )] i +[Ax (n) ] i ; z (n) i π ij Separable quadratic surrogate function: Q (x; x (n) P )= Q j (x j ; x (n) ) j=1 Q j (x j ; x (n) )= N i=1 Parallelizable update: x (n+1) j for j = 1,...,P ( π ij q [aij (x j x (n) j )] i +[Ax (n) ] i ; z (n) i π ij = argmin x j C ) Q j (x j ; x (n) )+V j (x j ; x (n) ), 12
14 How Many Holograms? Number of measurement elements N is fixed by the recorder (e.g., CCD camera pixels) Number of reconstructed pixels P is chosen by algorithm designer. Natural choice is P = N, but this is under-determined! (Measured data values are real, reconstructed field is complex.) Need P N/2 to avoid an under-determined problem. In conventional digital holography, FFT zero-padding is used. Possible strategies: 1 recorded hologram; reconstruct half-size image: P = N/2 1 recorded hologram; reconstruct full-size image: P = N (Must rely on regularization.) 2 recorded holograms with different reference beams; reconstruct full-size holographic image P = N/2 13
15 Simulation Data: Complex Object Original Object Magnitude Phase x 10 4 x Holograms Data #1 Data #2 0 14
16 Reconstruction Results: Complex Object Filtering Method 0.2 One Data Set, P=N/2 0.2 One Data Set, P=N 0.2 Two Data Sets, P=N/2 (NMSE = 0.16) (NMSE = 0.03) (NMSE = 0.03) (NMSE = 0.02) 0.2 Nonquadratic, edge-preserving regularization ψ(t)= t/δ log(1 + t/δ ), with δ chosen empirically 200 iterations of separable quadratic surrogate (SQS) algorithm 15
17 Contours of Marginal Objective Functions 300 x x x x x I x I x R One Data Set x R Two Data Sets h i (z) for z C, and h i ([Ax (n) ] i ) vs n (Unregularized.) 16
18 Real-Valued Object (Constrained) Hologram Data #1 Filtering Method 1 Data Set, P=N (NMSE = 0.19) (NMSE = 0.05) Original image Hologram Data #2 1 Data Set, P=N/2 (NMSE = 0.04) 2 Data Sets, P=N/2 (NMSE = 0.03) 17
19 Summary Statistical method for reconstructing a complex-valued object field from real-valued hologram intensity data Poisson statistical model. Extendable to others, e.g., Poisson+Gaussian (Snyder et al.) Optimization transfer monotonically decreases the cost function Preliminary simulations suggest improved image quality compared with conventional digital holography reconstruction method. Does not assume the reference beam is a plane-wave! But requires reasonably accurate reference beam model... Can incorporate constraints such as real-valued object field Accommodates hologram reference beam diversity, reducing the problem of multiple minima 18
20 Future Work? The usual story Use a space-variant model Regularization parameter selection Faster converging algorithms Test with real data Specific to this topic Different types of digital holography (e.g., Fresnel, Fourier) Phase retrieval problems 19
Part 3. Algorithms. Method = Cost Function + Algorithm
Part 3. Algorithms Method = Cost Function + Algorithm Outline Ideal algorithm Classical general-purpose algorithms Considerations: nonnegativity parallelization convergence rate monotonicity Algorithms
More informationOptimization transfer approach to joint registration / reconstruction for motion-compensated image reconstruction
Optimization transfer approach to joint registration / reconstruction for motion-compensated image reconstruction Jeffrey A. Fessler EECS Dept. University of Michigan ISBI Apr. 5, 00 0 Introduction Image
More informationChapter 16 Holography
Chapter 16 Holography Virtually all recording devices for light respond to light intensity. Problem: How to record, and then later reconstruct both the amplitude and phase of an optical wave. [This question
More informationRelaxed linearized algorithms for faster X-ray CT image reconstruction
Relaxed linearized algorithms for faster X-ray CT image reconstruction Hung Nien and Jeffrey A. Fessler University of Michigan, Ann Arbor The 13th Fully 3D Meeting June 2, 2015 1/20 Statistical image reconstruction
More informationA Study of Numerical Algorithms for Regularized Poisson ML Image Reconstruction
A Study of Numerical Algorithms for Regularized Poisson ML Image Reconstruction Yao Xie Project Report for EE 391 Stanford University, Summer 2006-07 September 1, 2007 Abstract In this report we solved
More informationTransmission Electron Microscopy
L. Reimer H. Kohl Transmission Electron Microscopy Physics of Image Formation Fifth Edition el Springer Contents 1 Introduction... 1 1.1 Transmission Electron Microscopy... 1 1.1.1 Conventional Transmission
More informationStatistical Methods for Image Reconstruction
Statistical Methods for Image Reconstruction Jeffrey A. Fessler EECS Department The University of Michigan Johns Hopkins University: Short Course May 11, 2007 0.0 Image Reconstruction Methods (Simplified
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationRegularized Multiplicative Algorithms for Nonnegative Matrix Factorization
Regularized Multiplicative Algorithms for Nonnegative Matrix Factorization Christine De Mol (joint work with Loïc Lecharlier) Université Libre de Bruxelles Dept Math. and ECARES MAHI 2013 Workshop Methodological
More informationOptics.
Optics www.optics.rochester.edu/classes/opt100/opt100page.html Course outline Light is a Ray (Geometrical Optics) 1. Nature of light 2. Production and measurement of light 3. Geometrical optics 4. Matrix
More informationHolographic interferometry for the study of liquids
Holographic interferometry for the study of liquids J. Colombani Prof. J. Bert - Laboratoire PMCN Université Claude Bernard Lyon 1 - France 1 Etymology latin inter : between greek ολοs : all latin metrum
More informationLight Propagation in Free Space
Intro Light Propagation in Free Space Helmholtz Equation 1-D Propagation Plane waves Plane wave propagation Light Propagation in Free Space 3-D Propagation Spherical Waves Huygen s Principle Each point
More informationStatistical Methods for Image Reconstruction
Statistical Methods for Image Reconstruction Jeffrey A. Fessler EECS Department The University of Michigan NSS-MIC Oct. 19, 2004 0.0 Image Reconstruction Methods (Simplified View) Analytical (FBP) Iterative
More informationScan Time Optimization for Post-injection PET Scans
Presented at 998 IEEE Nuc. Sci. Symp. Med. Im. Conf. Scan Time Optimization for Post-injection PET Scans Hakan Erdoğan Jeffrey A. Fessler 445 EECS Bldg., University of Michigan, Ann Arbor, MI 4809-222
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 informationMetrology and Sensing
Metrology and Sensing Lecture 5: Interferometry I 08--6 Herbert Gross Winter term 08 www.iap.uni-jena.de Schedule Optical Metrology and Sensing 08 No Date Subject Detailed Content 6.0. Introduction Introduction,
More informationA Limited Memory, Quasi-Newton Preconditioner. for Nonnegatively Constrained Image. Reconstruction
A Limited Memory, Quasi-Newton Preconditioner for Nonnegatively Constrained Image Reconstruction Johnathan M. Bardsley Department of Mathematical Sciences, The University of Montana, Missoula, MT 59812-864
More informationConvex Optimization. Lecture 12 - Equality Constrained Optimization. Instructor: Yuanzhang Xiao. Fall University of Hawaii at Manoa
Convex Optimization Lecture 12 - Equality Constrained Optimization Instructor: Yuanzhang Xiao University of Hawaii at Manoa Fall 2017 1 / 19 Today s Lecture 1 Basic Concepts 2 for Equality Constrained
More informationChapter 7. Interference of Light
Chapter 7. Interference of Light Last Lecture Superposition of waves Laser This Lecture Two-Beam Interference Young s Double Slit Experiment Virtual Sources Newton s Rings Film Thickness Measurement by
More informationSparsity Regularization
Sparsity Regularization Bangti Jin Course Inverse Problems & Imaging 1 / 41 Outline 1 Motivation: sparsity? 2 Mathematical preliminaries 3 l 1 solvers 2 / 41 problem setup finite-dimensional formulation
More informationGraphical Models for Collaborative Filtering
Graphical Models for Collaborative Filtering Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012 Sequence modeling HMM, Kalman Filter, etc.: Similarity: the same graphical model topology,
More informationChapter 10. Interference of Light
Chapter 10. Interference of Light Last Lecture Wave equations Maxwell equations and EM waves Superposition of waves This Lecture Two-Beam Interference Young s Double Slit Experiment Virtual Sources Newton
More informationMetrology and Sensing
Metrology and Sensing Lecture 5: Interferometry I 06--09 Herbert Gross Winter term 06 www.iap.uni-jena.de Preliminary Schedule No Date Subject Detailed Content 8.0. Introduction Introduction, optical measurements,
More informationImage restoration: numerical optimisation
Image restoration: numerical optimisation Short and partial presentation Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS BINP / 6 Context
More informationImage restoration. An example in astronomy
Image restoration Convex approaches: penalties and constraints An example in astronomy Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS
More informationInverse problem and optimization
Inverse problem and optimization Laurent Condat, Nelly Pustelnik CNRS, Gipsa-lab CNRS, Laboratoire de Physique de l ENS de Lyon Decembre, 15th 2016 Inverse problem and optimization 2/36 Plan 1. Examples
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationPH880 Topics in Physics
PH880 Topics in Physics Modern Optical Imaging (Fall 2010) Monday Fourier Optics Overview of week 3 Transmission function, Diffraction 4f telescopic system PSF, OTF Wednesday Conjugate Plane Bih Bright
More informationPRINCIPLES OF PHYSICAL OPTICS
PRINCIPLES OF PHYSICAL OPTICS C. A. Bennett University of North Carolina At Asheville WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION CONTENTS Preface 1 The Physics of Waves 1 1.1 Introduction
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 informationAN EFFICIENT COMPUTATIONAL METHOD FOR TOTAL VARIATION-PENALIZED POISSON LIKELIHOOD ESTIMATION. Johnathan M. Bardsley
Volume X, No. 0X, 200X, X XX Web site: http://www.aimsciences.org AN EFFICIENT COMPUTATIONAL METHOD FOR TOTAL VARIATION-PENALIZED POISSON LIKELIHOOD ESTIMATION Johnathan M. Bardsley Department of Mathematical
More informationSPATIAL RESOLUTION PROPERTIES OF PENALIZED WEIGHTED LEAST-SQUARES TOMOGRAPHIC IMAGE RECONSTRUCTION WITH MODEL MISMATCH
SPATIAL RESOLUTION PROPERTIES OF PENALIZED WEIGHTED LEAST-SQUARES TOMOGRAPHIC IMAGE RECONSTRUCTION WITH MODEL MISMATCH Jeffrey A. Fessler COMMUNICATIONS & SIGNAL PROCESSING LABORATORY Department of Electrical
More informationOptimized first-order minimization methods
Optimized first-order minimization methods Donghwan Kim & Jeffrey A. Fessler EECS Dept., BME Dept., Dept. of Radiology University of Michigan web.eecs.umich.edu/~fessler UM AIM Seminar 2014-10-03 1 Disclosure
More informationSparsity Regularization for Image Reconstruction with Poisson Data
Sparsity Regularization for Image Reconstruction with Poisson Data Daniel J. Lingenfelter a, Jeffrey A. Fessler a,andzhonghe b a Electrical Engineering and Computer Science, University of Michigan, Ann
More informationA Majorize-Minimize subspace approach for l 2 -l 0 regularization with applications to image processing
A Majorize-Minimize subspace approach for l 2 -l 0 regularization with applications to image processing Emilie Chouzenoux emilie.chouzenoux@univ-mlv.fr Université Paris-Est Lab. d Informatique Gaspard
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
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 informationPhotometric Stereo: Three recent contributions. Dipartimento di Matematica, La Sapienza
Photometric Stereo: Three recent contributions Dipartimento di Matematica, La Sapienza Jean-Denis DUROU IRIT, Toulouse Jean-Denis DUROU (IRIT, Toulouse) 17 December 2013 1 / 32 Outline 1 Shape-from-X techniques
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More informationOptimization by General-Purpose Methods
c J. Fessler. January 25, 2015 11.1 Chapter 11 Optimization by General-Purpose Methods ch,opt Contents 11.1 Introduction (s,opt,intro)...................................... 11.3 11.1.1 Iterative optimization
More informationAnalytical Approach to Regularization Design for Isotropic Spatial Resolution
Analytical Approach to Regularization Design for Isotropic Spatial Resolution Jeffrey A Fessler, Senior Member, IEEE Abstract In emission tomography, conventional quadratic regularization methods lead
More informationOptical Systems Program of Studies Version 1.0 April 2012
Optical Systems Program of Studies Version 1.0 April 2012 Standard1 Essential Understand Optical experimental methodology, data analysis, interpretation, and presentation strategies Essential Understandings:
More informationCOMPRESSIVE SENSING IN COHERENT IMAGING AND ANTENNA SYNTHESIS
COMPRESSIVE SENSING IN COHERENT IMAGING AND ANTENNA SYNTHESIS A Thesis Presented by Hongyi Zhang to The Department of Electrical and Computer Engineering In partial fulfillment of the requirements for
More informationNature of Light Part 2
Nature of Light Part 2 Fresnel Coefficients From Helmholts equation see imaging conditions for Single lens 4F system Diffraction ranges Rayleigh Range Diffraction limited resolution Interference Newton
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 informationSTAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song
STAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song Presenter: Jiwei Zhao Department of Statistics University of Wisconsin Madison April
More informationGravitational Lensing of the CMB
Gravitational Lensing of the CMB SNAP Planck 1 Ω DE 1 w a.5-2 -1.5 w -1 -.5 Wayne Hu Leiden, August 26-1 Outline Gravitational Lensing of Temperature and Polarization Fields Cosmological Observables from
More informationAccelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems)
Accelerated Dual Gradient-Based Methods for Total Variation Image Denoising/Deblurring Problems (and other Inverse Problems) Donghwan Kim and Jeffrey A. Fessler EECS Department, University of Michigan
More informationMulti-stage phase retrieval algorithm based upon the gyrator transform
Multi-stage phase retrieval algorithm based upon the gyrator transform José A. Rodrigo, Hamootal Duadi, Tatiana Alieva 3 and Zeev Zalevsky Instituto de Óptica (CSIC), Imaging and Vision Department. Serrano,
More informationEffect of additive noise on phase measurement in digital holographic microscopy
Topical Editor: Dr. Tristan Colomb, 3D Res. 0, 01006 (011) 10.1007/01(011)6 3DR EXPRESS w Effect of additive noise on phase measurement in digital holographic microscopy Nitesh Pandey Bryan Hennelly Received:
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 informationOptimization with nonnegativity constraints
Optimization with nonnegativity constraints Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 30, 2007 Seminar: Inverse problems 1 Introduction Yves van Gennip February 21 2 Regularization strategies
More informationA memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration
A memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration E. Chouzenoux, A. Jezierska, J.-C. Pesquet and H. Talbot Université Paris-Est Lab. d Informatique Gaspard
More informationRegularizing inverse problems using sparsity-based signal models
Regularizing inverse problems using sparsity-based signal models Jeffrey A. Fessler William L. Root Professor of EECS EECS Dept., BME Dept., Dept. of Radiology University of Michigan http://web.eecs.umich.edu/
More informationarxiv: v4 [math.oc] 1 Apr 2018
Generalizing the optimized gradient method for smooth convex minimization Donghwan Kim Jeffrey A. Fessler arxiv:607.06764v4 [math.oc] Apr 08 Date of current version: April 3, 08 Abstract This paper generalizes
More informationLinear & nonlinear classifiers
Linear & nonlinear classifiers Machine Learning Hamid Beigy Sharif University of Technology Fall 1396 Hamid Beigy (Sharif University of Technology) Linear & nonlinear classifiers Fall 1396 1 / 44 Table
More informationProbabilistic Latent Semantic Analysis
Probabilistic Latent Semantic Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationMetrology and Sensing
Metrology and Sensing Lecture 5: Interferometry I 017-11-16 Herbert Gross Winter term 017 www.iap.uni-jena.de Preliminary Schedule No Date Subject Detailed Content 1 19.10. Introduction Introduction, optical
More informationOptical Sciences Center, Rm 704 University of Arizona Tucson, AZ Office Hours: Call for appointment or see after class
Term: Spring 2000 Course #: OPTI 505 Course Title: Diffraction and Interferometry Instructor: James C. Wyant Optical Sciences Center, Rm 704 University of Arizona Tucson, AZ 85721 Phone: 520-621-2448 E-Mail:
More informationReconstruction from Fourier Samples (Gridding and alternatives)
Chapter 6 Reconstruction from Fourier Samples (Gridding and alternatives) ch,four The customer of sampling and reconstruction technology faces a large gap between the practice of (non-uniform) sampling
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 informationAdvanced holographic nondestructive testing system for residual stress analysis
Advanced holographic nondestructive testing system for residual stress analysis Anatoli Kniazkov, Yuri Onischenko, George Dovgalenko, Gregory Salamo Laser Laboratory, Physics Department, University of
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationCOMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, August 23, 2004, 12:14 PM)
COMPUTER GENERATED HOLOGRAMS Optical Sciences 67 W.J. Dallas (Monday, August 3, 4, 1:14 PM) PART IV: CHAPTER FOUR OPTICAL TESTING Part IV: Chapter Four Page 1 of 1 Introduction In optical testing an element
More informationTwo-Material Decomposition From a Single CT Scan Using Statistical Image Reconstruction
/ 5 Two-Material Decomposition From a Single CT Scan Using Statistical Image Reconstruction Yong Long and Jeffrey A. Fessler EECS Department James M. Balter Radiation Oncology Department The University
More information2.71/2.710 Optics. MIT 2.71/ /5/01 wk1-b-1
2.71/2.710 Optics Units: 3-0-9, Prerequisites: 8.02, 18.03 2.71: meets the Department Elective requirement 2.710: H-Level, meets the MS requirement in Design 9/5/01 wk1-b-1 Class objectives Cover the fundamental
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 informationOutline. Scientific Computing: An Introductory Survey. Optimization. Optimization Problems. Examples: Optimization Problems
Outline Scientific Computing: An Introductory Survey Chapter 6 Optimization 1 Prof. Michael. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationKarush-Kuhn-Tucker Conditions. Lecturer: Ryan Tibshirani Convex Optimization /36-725
Karush-Kuhn-Tucker Conditions Lecturer: Ryan Tibshirani Convex Optimization 10-725/36-725 1 Given a minimization problem Last time: duality min x subject to f(x) h i (x) 0, i = 1,... m l j (x) = 0, j =
More informationCoherent Superposition of Mueller-Jones States: Synthesizing Mueller Matrices with an Ellipsometer
Coherent Superposition of Mueller-Jones States: Synthesizing Mueller Matrices with an Ellipsometer Ertan Kuntman, Mehmet A. Kuntman Oriol Arteaga SPECTROSCOPIC ELLIPSOMETRY FOCUS TOPIC Young s double slit
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 informationUniform Quadratic Penalties Cause Nonuniform Spatial Resolution
Uniform Quadratic Penalties Cause Nonuniform Spatial Resolution Jeffrey A. Fessler and W. Leslie Rogers 3480 Kresge 111, Box 0552, University of Michigan, Ann Arbor, MI 48109-0552 ABSTRACT This paper examines
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 informationNumerical Linear Algebra Primer. Ryan Tibshirani Convex Optimization
Numerical Linear Algebra Primer Ryan Tibshirani Convex Optimization 10-725 Consider Last time: proximal Newton method min x g(x) + h(x) where g, h convex, g twice differentiable, and h simple. Proximal
More informationNumerical Optimization: Basic Concepts and Algorithms
May 27th 2015 Numerical Optimization: Basic Concepts and Algorithms R. Duvigneau R. Duvigneau - Numerical Optimization: Basic Concepts and Algorithms 1 Outline Some basic concepts in optimization Some
More informationSupplementary Note 1 Description of the sample and thin lamella preparation Supplementary Figure 1 FeRh lamella prepared by FIB and used for in situ
Supplementary Note 1 Description of the sample and thin lamella preparation A 5nm FeRh layer was epitaxially grown on a go (1) substrate by DC sputtering using a co-deposition process from two pure Fe
More informationA quantative comparison of two restoration methods as applied to confocal microscopy
A quantative comparison of two restoration methods as applied to confocal microscopy Geert M.P. van Kempen 1, Hans T.M. van der Voort, Lucas J. van Vliet 1 1 Pattern Recognition Group, Delft University
More informationApplication of imaging ellipsometry: graphene - pinpointing and ellipsometric characterization ULRICH WURSTBAUER CHRISTIAN RÖLING PETER H.
Application of imaging ellipsometry: graphene - pinpointing and ellipsometric characterization ULRICH WURSTBAUER CHRISTIAN RÖLING PETER H. THIESEN Introduction to imaging ellipsometry Ellipsometry in general
More information, b = 0. (2) 1 2 The eigenvectors of A corresponding to the eigenvalues λ 1 = 1, λ 2 = 3 are
Quadratic forms We consider the quadratic function f : R 2 R defined by f(x) = 2 xt Ax b T x with x = (x, x 2 ) T, () where A R 2 2 is symmetric and b R 2. We will see that, depending on the eigenvalues
More informationVariance Images for Penalized-Likelihood Image Reconstruction
IEEE TRANSACTIONS ON MEDICAL IMAGING, VERSION September 29, 2 1 Variance Images for Penalized-Likelihood Image Reconstruction Jeffrey A. Fessler 424 EECS Bldg., University of Michigan, Ann Arbor, MI 4819-2122
More informationOptics. n n. sin c. sin
Optics Geometrical optics (model) Light-ray: extremely thin parallel light beam Using this model, the explanation of several optical phenomena can be given as the solution of simple geometric problems.
More information10. Unconstrained minimization
Convex Optimization Boyd & Vandenberghe 10. Unconstrained minimization terminology and assumptions gradient descent method steepest descent method Newton s method self-concordant functions implementation
More informationStatic and Dynamic Optimization (42111)
Static and Dynamic Optimization (421) Niels Kjølstad Poulsen Build. 0b, room 01 Section for Dynamical Systems Dept. of Applied Mathematics and Computer Science The Technical University of Denmark Email:
More informationMotion Estimation (I)
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationCoordinate descent. Geoff Gordon & Ryan Tibshirani Optimization /
Coordinate descent Geoff Gordon & Ryan Tibshirani Optimization 10-725 / 36-725 1 Adding to the toolbox, with stats and ML in mind We ve seen several general and useful minimization tools First-order methods
More informationSparsifying Transform Learning for Compressed Sensing MRI
Sparsifying Transform Learning for Compressed Sensing MRI Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and Coordinated Science Laborarory University of Illinois
More informationSuperiorized Inversion of the Radon Transform
Superiorized Inversion of the Radon Transform Gabor T. Herman Graduate Center, City University of New York March 28, 2017 The Radon Transform in 2D For a function f of two real variables, a real number
More informationSparse Regularization via Convex Analysis
Sparse Regularization via Convex Analysis Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering New York University Brooklyn, New York, USA 29 / 66 Convex or non-convex: Which
More informationOptical Metrology and Sensing
Optical Metrology and Sensing Lecture 1: Introduction 2016-10-18 Herbert Gross Winter term 2016 www.iap.uni-jena.de 2 Preliminary Schedule No Date Subject Detailed Content 1 18.10. Introduction Introduction,
More informationFish embryo multimodal imaging by laser Doppler digital holography
Fish embryo multimodal imaging by laser Doppler digital holography Nicolas Verrier, Daniel Alexandre, Pascal Picart, Michel Gross To cite this version: Nicolas Verrier, Daniel Alexandre, Pascal Picart,
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationA Digital Holographic Approach for Co-Phasing of Segmented Telescopes: Proof of Concept Using Numerical Simulations
PUBLICATIONS OF THE ASTRONOMICAL SOCIETY OF THE PACIFIC, 126:280 286, 2014 March 2014. The Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A. A Digital Holographic Approach for
More informationNonlinear Programming
Nonlinear Programming Kees Roos e-mail: C.Roos@ewi.tudelft.nl URL: http://www.isa.ewi.tudelft.nl/ roos LNMB Course De Uithof, Utrecht February 6 - May 8, A.D. 2006 Optimization Group 1 Outline for week
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 informationProvable Alternating Minimization Methods for Non-convex Optimization
Provable Alternating Minimization Methods for Non-convex Optimization Prateek Jain Microsoft Research, India Joint work with Praneeth Netrapalli, Sujay Sanghavi, Alekh Agarwal, Animashree Anandkumar, Rashish
More informationChristophe De Beule. Graphical interpretation of the Schrödinger equation
Christophe De Beule Graphical interpretation of the Schrödinger equation Outline Schrödinger equation Relation between the kinetic energy and wave function curvature. Classically allowed and forbidden
More informationRegularization Parameter Selection Methods for Ill-Posed Poisson Maximum Likelihood Estimation
Regularization Parameter Selection Methods for Ill-Posed Poisson Maximum Likelihood Estimation Johnathan M. Bardsley and John Goldes Department of Mathematical Sciences University of Montana Missoula,
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationNoise, Image Reconstruction with Noise!
Noise, Image Reconstruction with Noise! EE367/CS448I: Computational Imaging and Display! stanford.edu/class/ee367! Lecture 10! Gordon Wetzstein! Stanford University! What s a Pixel?! photon to electron
More information