Deconvolution of the Fourier Spectrum

Total Page:16

File Type:pdf, Size:1020Kb

Deconvolution of the Fourier Spectrum Koninklijk Meteorologisch Instituut van Belgi¨e Institut Royal M´et´eorologique de Belgique Deconvolution of the Fourier spectrum F. De Meyer 2003 Wetenschappelijke en Publication scientifique technische publicatie et technique Nr 33 No 33 Uitgegeven door het Edit´epar KONINKLIJK METEOROLOGISCH l’INSTITUT ROYAL INSTITUUT VAN BELGIE METEOROLOGIQUE DE BELGIQUE Ringlaan 3, B -1180 Brussel Avenue Circulaire 3, B -1180 Bruxelles Verantwoordelijke uitgever: Dr. H. Malcorps Editeur responsable: Dr. H. Malcorps Koninklijk Meteorologisch Instituut van Belgi¨e Institut Royal M´et´eorologique de Belgique Deconvolution of the Fourier spectrum F. De Meyer 2003 Wetenschappelijke en Publication scientifique technische publicatie et technique Nr 33 No 33 Uitgegeven door het Edit´epar KONINKLIJK METEOROLOGISCH l’INSTITUT ROYAL INSTITUUT VAN BELGIE METEOROLOGIQUE DE BELGIQUE Ringlaan 3, B -1180 Brussel Avenue Circulaire 3, B -1180 Bruxelles Verantwoordelijke uitgever: Dr. H. Malcorps Editeur responsable: Dr. H. Malcorps 1 Abstract The problem of estimating the frequency spectrum of a real continuous function, which is measured only at a finite number of discrete times, is discussed in this tutorial review. Based on the deconvolu- tion method, the complex version of the one-dimensional CLEAN algorithm provides a simple way to remove the artefacts introduced by the sampling and the effects of missing data from the computed Fourier spectrum. The technique is very appropriate in the case of equally spaced data, as well as for data samples randomly distributed in time. An example of the application of CLEAN to a synthetic spectrum of a small number of harmonic components at discrete frequencies is shown. The case of a low signal-to-noise time sequence is also presented by illustrating how CLEAN recovers the spectrum of the declination component of the geomagnetic field, measured in the magnetic observatory of Dourbes. 1. Introduction Synthesis maps in radio astronomy are generated by Fourier transformation of interferome- ter visibility data to produce a map of the sky (Thompson et al., 1986). Because only a finite number of visibility points are gathered, the problems of sampling apply to interferometer maps. The CLEAN algorithm is widely used in two-dimensional image reconstruction and performs an approximate deconvolution of the ‘true map’ of the sky from the ‘dirty map’ procreated from the data, in essence, removing the false features introduced inherently by the finite sampling. By adapting the two-dimensional CLEAN procedure developed for use in aperture synthesis, Roberts et al. (1987) discussed a rather intuitive way of handling the difficulties associated with the spurious apparent responses that arise from the incompleteness of the discrete sampling of a continuous time function. A multitude of spectral estimation algorithms for discrete time series have been proposed (Kay and Marple, 1981). Comparisons among various competing techniques have been based on limited computer simulations, which can be misleading. Here we select the conventional Fourier spectral estimator because the effects of finite, discrete sampling on the performance of the classical periodogram are precisely known and can be almost completely rectified. Estimation of the power spectral density, or simply the spectrum, of discretely sampled deterministic and stochastic processes is usually based on procedures employing the Fast Fourier Transform or FFT (Brigham, 1974), which implicitly assumes that the data are periodic outside the observ- ing interval. Unfortunately the FFT technique is directly applicable only when the input data sequence is evenly spaced in time. Often observations cannot be controlled to the extent that observability constraints may cause an unequally spaced sample domain. One is then forced to extract spectral information from these irregularly sampled data. The FFT approach to spectrum analysis is computationally efficient and produces reason- able outcomes for a large class of signal processes (Briggs and Henson, 1995). In spite of these advantages, there are several inherent performance limitations of the FFT procedure. The most prominent shortcoming is that of frequency resolution, i.e., the ability to distinguish the spec- tral responses of distinct signals. The frequency resolution is roughly the reciprocal of the time interval over which sampled data are available. A second imperfection is due to the implicit windowing of the data that occurs when processing with the FFT. Windowing manifests itself as ‘leakage’ in the spectral domain, that is, energy in the main lobe of a spectral response ‘leaks’ into the sidelobes, obscuring and distorting other spectral responses that are present. In fact, weak signal spectral responses can be masked by higher sidelobes from stronger spectral com- ponents. These two performance limitations of the FFT approach are particularly troublesome when analyzing short data records. 2 In an attempt to alleviate the inherent shortcomings of the FFT approach, many alternative spectral estimation procedures have been proposed. For instance, Kay and Marple (1981) discuss the improvement that may result from non-traditional methods such as the autoregressive (AR) method or the maximum entropy method (MEM). Claims have been made concerning the de- gree of improvement obtained in the spectral resolution and the signal detectability when these numerical techniques are applied to numerical data. These performance advantages, though, strongly depend upon the signal-to-noise ratio (SNR) of the input data, as might be expected. In fact, for low enough SNR’s the ‘modern’ spectral estimators are often no better than those obtained with conventional FFT processing. Even in those cases where improved spectral fidelity is achieved by use of a different spectral estimation procedure, the computational requirements of that alternate method may be substantially higher than FFT computation. This may make some spectral estimators unattractive for real-time implementation. Schwarz (1978) has shown that CLEAN is a statistically correct method of least squares fit- ting sinusoidal functions to the observations which is the conventional framework of the discrete Fourier transform. Fourier inversion of a finite representation of a continuous function leads to a frequency spectrum distorted by (1) the limited frequency resolution due to the finite time span of the data sample, and (2) spurious apparent responses that are caused by the incompleteness of the sampling. There is an intuitive way of handling the difficulties associated with the latter effect. With a view to practical applications heuristic arguments are used in this discussion to make the discourse more easily accessible. The CLEAN algorithm performs, in essence, a non-linear deconvolution of the Fourier spectrum on the frequency axis, equivalent to a least squares interpolation in the time domain. The numerical method is particularly suitable for functions whose spectra are dominated by a few harmonic components at discrete frequencies. The CLEAN technique for time series spectral analysis (Roberts et al., 1987) works by subtracting from the noisy or ‘dirty’ spectrum, which is the convolution of the ‘true’ spectrum with the ‘dirty beam’ generated by the finite number of data, the response expected from a sinusoidal component with frequency corresponding to the maximum of the spectrum. In this way a residual spectrum is produced from which both the component and the spurious features due to sampling have been removed. This procedure of ‘CLEANing’ the spectrum is repeated on successive residual spectra until nothing is left but noise. The set of resulting ‘clean com- ponents’ is then used as a model for the time function. This model is convolved with a ‘clean beam’ which has the same resolution as the original ‘spectral window’ but no sidelobes. The convolution of the clean components model with the clean beam serves to weight down the interpolated information at time instants outside those sampled. Finally, to preserve the noise level of the spectrum and to incorporate any components which were not resolved by CLEAN, the final residual spectrum is added to the convolved clean components to produce the ‘clean spectrum’. The one-dimensional version of the CLEAN deconvolution technique is especially useful for spectral analysis of unequally spaced observations and provides a simple way to dispose of the artefacts of missing data, a situation which is frequently encountered in practical time series. As a first example, the application of CLEAN will be illustrated by means of the frequency analysis of a synthetic spectrum, that is generated by a small number of harmonic components at discrete frequencies. The method is also used in the search for periodic variability in the time series of the hourly values of the magnetic declination, measured in the observatory of Dourbes, Belgium, which is characterized by a spectrum having a low signal-to-noise ratio. 3 2. The Fourier spectrum of a continuous function The basic recourse for the frequency analysis of an aperiodic or transient function is the Fourier integral (Bracewell, 1986). Here we follow the notation of Jenkins and Watts (1968). The spec- trum of an analog function x(t), known for all time t, is given by the continuous Fourier transform (denoted by the operator F)ofx(t), ∞ X(f) ≡F{x(t)} = x(t) e−2πift dt, −∞≤f ≤∞, (1) −∞ which defines the contribution of each frequency f to x(t). The inverse Fourier transform (de- noted by F −1) is defined by the Riemann
Recommended publications
  • Convolution! (CDT-14) Luciano Da Fontoura Costa
    Convolution! (CDT-14) Luciano da Fontoura Costa To cite this version: Luciano da Fontoura Costa. Convolution! (CDT-14). 2019. hal-02334910 HAL Id: hal-02334910 https://hal.archives-ouvertes.fr/hal-02334910 Preprint submitted on 27 Oct 2019 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Convolution! (CDT-14) Luciano da Fontoura Costa [email protected] S~aoCarlos Institute of Physics { DFCM/USP October 22, 2019 Abstract The convolution between two functions, yielding a third function, is a particularly important concept in several areas including physics, engineering, statistics, and mathematics, to name but a few. Yet, it is not often so easy to be conceptually understood, as a consequence of its seemingly intricate definition. In this text, we develop a conceptual framework aimed at hopefully providing a more complete and integrated conceptual understanding of this important operation. In particular, we adopt an alternative graphical interpretation in the time domain, allowing the shift implied in the convolution to proceed over free variable instead of the additive inverse of this variable. In addition, we discuss two possible conceptual interpretations of the convolution of two functions as: (i) the `blending' of these functions, and (ii) as a quantification of `matching' between those functions.
    [Show full text]
  • Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation
    Computers, Materials & Continua Tech Science Press DOI:10.32604/cmc.2021.014752 Article Windowing Techniques, the Welch Method for Improvement of Power Spectrum Estimation Dah-Jing Jwo1, *, Wei-Yeh Chang1 and I-Hua Wu2 1Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University, Keelung, 202-24, Taiwan 2Innovative Navigation Technology Ltd., Kaohsiung, 801, Taiwan *Corresponding Author: Dah-Jing Jwo. Email: [email protected] Received: 01 October 2020; Accepted: 08 November 2020 Abstract: This paper revisits the characteristics of windowing techniques with various window functions involved, and successively investigates spectral leak- age mitigation utilizing the Welch method. The discrete Fourier transform (DFT) is ubiquitous in digital signal processing (DSP) for the spectrum anal- ysis and can be efciently realized by the fast Fourier transform (FFT). The sampling signal will result in distortion and thus may cause unpredictable spectral leakage in discrete spectrum when the DFT is employed. Windowing is implemented by multiplying the input signal with a window function and windowing amplitude modulates the input signal so that the spectral leakage is evened out. Therefore, windowing processing reduces the amplitude of the samples at the beginning and end of the window. In addition to selecting appropriate window functions, a pretreatment method, such as the Welch method, is effective to mitigate the spectral leakage. Due to the noise caused by imperfect, nite data, the noise reduction from Welch’s method is a desired treatment. The nonparametric Welch method is an improvement on the peri- odogram spectrum estimation method where the signal-to-noise ratio (SNR) is high and mitigates noise in the estimated power spectra in exchange for frequency resolution reduction.
    [Show full text]
  • Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage
    E3S Web of Conferences 94, 03001 (2019) https://doi.org/10.1051/e3sconf/20199403001 ISGNSS 2018 Windowing Design and Performance Assessment for Mitigation of Spectrum Leakage Dah-Jing Jwo 1,*, I-Hua Wu 2 and Yi Chang1 1 Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University 2 Pei-Ning Rd., Keelung 202, Taiwan 2 Bison Electronics, Inc., Neihu District, Taipei 114, Taiwan Abstract. This paper investigates the windowing design and performance assessment for mitigation of spectral leakage. A pretreatment method to reduce the spectral leakage is developed. In addition to selecting appropriate window functions, the Welch method is introduced. Windowing is implemented by multiplying the input signal with a windowing function. The periodogram technique based on Welch method is capable of providing good resolution if data length samples are selected optimally. Windowing amplitude modulates the input signal so that the spectral leakage is evened out. Thus, windowing reduces the amplitude of the samples at the beginning and end of the window, altering leakage. The influence of various window functions on the Fourier transform spectrum of the signals was discussed, and the characteristics and functions of various window functions were explained. In addition, we compared the differences in the influence of different data lengths on spectral resolution and noise levels caused by the traditional power spectrum estimation and various window-function-based Welch power spectrum estimations. 1 Introduction knowledge. Due to computer capacity and limitations on processing time, only limited time samples can actually The spectrum analysis [1-4] is a critical concern in many be processed during spectrum analysis and data military and civilian applications, such as performance processing of signals, i.e., the original signals need to be test of electric power systems and communication cut off, which induces leakage errors.
    [Show full text]
  • Fourier Ptychographic Wigner Distribution Deconvolution
    FOURIER PTYCHOGRAPHIC WIGNER DISTRIBUTION DECONVOLUTION Justin Lee1, George Barbastathis2,3 1. Department of Health Sciences and Technology, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 2. Department of Mechanical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 3. Singapore–MIT Alliance for Research and Technology (SMART) Centre Singapore 138602, Singapore Email: [email protected] KEY WORDS: Computational Imaging, phase retrieval, ptychography, Fourier ptychography Ptychography was first proposed as a method for phase retrieval by Hoppe in the late 1960s [1]. A typical ptychographic setup involves generation of illumination diversity via lateral translation of an object relative to its illumination (often called the probe) and measurement of the intensity of the far-field diffraction pattern/Fourier transform of the illuminated object. In the decades following Hoppes initial work, several important extensions to ptychography were developed, including an ability to simultaneously retrieve the amplitude and phase of both the probe and the object via application of an extended ptychographic iterative engine (ePIE) [2] and a non-iterative technique for phase retrieval using ptychographic information called Wigner Distribution Deconvolution (WDD) [3,4]. Most recently, the Fourier dual to ptychography was introduced by the authors in [5] as an experimentally simple technique for phase retrieval using an LED array as the illumination source in a standard microscope setup. In Fourier ptychography, lateral translation of the object relative to the probe occurs in frequency space and is obtained by tilting the angle of illumination (while holding the pupil function of the imaging system constant). In the past few years, several improvements to Fourier ptychography have been introduced, including an ePIE-like technique for simultaneous pupil function and object retrieval [6].
    [Show full text]
  • Effect of Data Gaps: Comparison of Different Spectral Analysis Methods
    Ann. Geophys., 34, 437–449, 2016 www.ann-geophys.net/34/437/2016/ doi:10.5194/angeo-34-437-2016 © Author(s) 2016. CC Attribution 3.0 License. Effect of data gaps: comparison of different spectral analysis methods Costel Munteanu1,2,3, Catalin Negrea1,4,5, Marius Echim1,6, and Kalevi Mursula2 1Institute of Space Science, Magurele, Romania 2Astronomy and Space Physics Research Unit, University of Oulu, Oulu, Finland 3Department of Physics, University of Bucharest, Magurele, Romania 4Cooperative Institute for Research in Environmental Sciences, Univ. of Colorado, Boulder, Colorado, USA 5Space Weather Prediction Center, NOAA, Boulder, Colorado, USA 6Belgian Institute of Space Aeronomy, Brussels, Belgium Correspondence to: Costel Munteanu ([email protected]) Received: 31 August 2015 – Revised: 14 March 2016 – Accepted: 28 March 2016 – Published: 19 April 2016 Abstract. In this paper we investigate quantitatively the ef- 1 Introduction fect of data gaps for four methods of estimating the ampli- tude spectrum of a time series: fast Fourier transform (FFT), Spectral analysis is a widely used tool in data analysis and discrete Fourier transform (DFT), Z transform (ZTR) and processing in most fields of science. The technique became the Lomb–Scargle algorithm (LST). We devise two tests: the very popular with the introduction of the fast Fourier trans- single-large-gap test, which can probe the effect of a sin- form algorithm which allowed for an extremely rapid com- gle data gap of varying size and the multiple-small-gaps test, putation of the Fourier Transform. In the absence of modern used to study the effect of numerous small gaps of variable supercomputers, this was not just useful, but also the only size distributed within the time series.
    [Show full text]
  • Enhanced Three-Dimensional Deconvolution Microscopy Using a Measured Depth-Varying Point-Spread Function
    2622 J. Opt. Soc. Am. A/Vol. 24, No. 9/September 2007 J. W. Shaevitz and D. A. Fletcher Enhanced three-dimensional deconvolution microscopy using a measured depth-varying point-spread function Joshua W. Shaevitz1,3,* and Daniel A. Fletcher2,4 1Department of Integrative Biology, University of California, Berkeley, California 94720, USA 2Department of Bioengineering, University of California, Berkeley, California 94720, USA 3Department of Physics, Princeton University, Princeton, New Jersey 08544, USA 4E-mail: fl[email protected] *Corresponding author: [email protected] Received February 26, 2007; revised April 23, 2007; accepted May 1, 2007; posted May 3, 2007 (Doc. ID 80413); published July 27, 2007 We present a technique to systematically measure the change in the blurring function of an optical microscope with distance between the source and the coverglass (the depth) and demonstrate its utility in three- dimensional (3D) deconvolution. By controlling the axial positions of the microscope stage and an optically trapped bead independently, we can record the 3D blurring function at different depths. We find that the peak intensity collected from a single bead decreases with depth and that the width of the axial, but not the lateral, profile increases with depth. We present simple convolution and deconvolution algorithms that use the full depth-varying point-spread functions and use these to demonstrate a reduction of elongation artifacts in a re- constructed image of a 2 ␮m sphere. © 2007 Optical Society of America OCIS codes: 100.3020, 100.6890, 110.0180, 140.7010, 170.6900, 180.2520. 1. INTRODUCTION mismatch between the objective lens and coupling oil can Three-dimensional (3D) deconvolution microscopy is a be severe [5,6], especially in systems using high- powerful tool for visualizing complex biological struc- numerical apertures.
    [Show full text]
  • Chapter 2: the Fourier Transform Chapter 2
    “EEE305”, “EEE801 Part A”: Digital Signal Processing Chapter 2: The Fourier Transform Chapter 2 The Fourier Transform 2.1 Introduction The sampled Fourier transform of a periodic, discrete-time signal is known as the discrete Fourier transform (DFT). The DFT of such a signal allows an interpretation of the frequency domain descriptions of the given signal. 2.2 Derivation of the DFT The Fourier transform pair of a continuous-time signal is given by Equation 2.1 and Equation 2.2. ∞ − j2πft X ( f ) = ∫ x(t) ⋅e dt (2.1) −∞ ∞ j2πft x(t) = ∫ X ( f ) ⋅e df (2.2) −∞ Now consider the signal x(t) to be periodic i.e. repeating. The Fourier transform pair from above can now be represented by Equation 2.3 and Equation 2.4 respectively. T0 1 − π = ⋅ j2 kf0t X (kf 0 ) ∫ x(t) e dt (2.3) T0 0 ∞ π = ⋅ j 2 kf0t x(t) ∑ X (kf0 ) e (2.4) k=−∞ = 1 where f 0 . Now suppose that the signal x(t) is sampled N times per period, as illustrated in Figure 2.1, with a T0 sampling period of T seconds. The discrete signal can now be represented by Equation 2.5, where δ is the dirac delta impulse function and has a unit area of one. ∞ x[nT ] = ∑ x(t) ⋅δ (t − nT ) (2.5) n=−∞ Now replace x(t) in the Fourier integral of Equation 2.3, with x[nT] of Equation 2.5. The Equation can now be re- written as: ∞ T0 1 − π = ⋅δ − ⋅ jk 2 f0t X (kf0 ) ∑ ∫ x(t) (t nT ) e dt (2.6) T0 n=−∞ 0 1 for t = nT Since the dirac delta function δ ()t − nT = , Equation 2.6 can be modified even further to become 0 otherwise Equation 2.7 below: N −1 1 − π = ⋅ jk 2 f0nT X[kf 0 ] ∑ x[nT] e (2.7) N n=0 University of Newcastle upon Tyne Page 2.1 “EEE305”, “EEE801 Part A”: Digital Signal Processing Chapter 2: The Fourier Transform Sampled sequence δ(t-nT) litude p Am nT Time (t) Figure 2.1: Sampled sequence for the DFT.
    [Show full text]
  • Fourier Analysis
    FOURIER ANALYSIS Lucas Illing 2008 Contents 1 Fourier Series 2 1.1 General Introduction . 2 1.2 Discontinuous Functions . 5 1.3 Complex Fourier Series . 7 2 Fourier Transform 8 2.1 Definition . 8 2.2 The issue of convention . 11 2.3 Convolution Theorem . 12 2.4 Spectral Leakage . 13 3 Discrete Time 17 3.1 Discrete Time Fourier Transform . 17 3.2 Discrete Fourier Transform (and FFT) . 19 4 Executive Summary 20 1 1. Fourier Series 1 Fourier Series 1.1 General Introduction Consider a function f(τ) that is periodic with period T . f(τ + T ) = f(τ) (1) We may always rescale τ to make the function 2π periodic. To do so, define 2π a new independent variable t = T τ, so that f(t + 2π) = f(t) (2) So let us consider the set of all sufficiently nice functions f(t) of a real variable t that are periodic, with period 2π. Since the function is periodic we only need to consider its behavior on one interval of length 2π, e.g. on the interval (−π; π). The idea is to decompose any such function f(t) into an infinite sum, or series, of simpler functions. Following Joseph Fourier (1768-1830) consider the infinite sum of sine and cosine functions 1 a0 X f(t) = + [a cos(nt) + b sin(nt)] (3) 2 n n n=1 where the constant coefficients an and bn are called the Fourier coefficients of f. The first question one would like to answer is how to find those coefficients.
    [Show full text]
  • Sequences Systems Convolution Dirac's Delta Function Fourier
    Sequences Dirac’s Delta Function ( 1 0 x 6= 0 fxngn=−∞ is a discrete sequence. δ(x) = , R 1 δ(x)dx = 1 1 x = 0 −∞ Can derive a sequence from a function x(t) by having xn = n R 1 x(tsn) = x( ). Sampling: f(x)δ(x − a)dx = f(a) fs −∞ P1 P1 Absolute summability: n=−∞ jxnj < 1 Dirac comb: s(t) = ts n=−∞ δ(t − tsn), x^(t) = x(t)s(t) P1 2 Square summability: n=−∞ jxnj < 1 (aka energy signal) Periodic: 9k > 0 : 8n 2 Z : xn = xn+k Fourier Transform ( 0 n < 0 j!n un = is the unit-step sequence Discrete Sejong’s of the form xn = e are eigensequences with 1 n ≥ 0 respect to a system T because for any !, there is some H(!) ( such that T fxng = H(!)fxng 1 n = 0 1 δn = is the impulse sequence Ffg(t)g(!) = G(!) = R g(t)e∓j!tdt 0 n 6= 0 −∞ −1 R 1 ±j!t F fG(!)g(t) = g(t) = −∞ G(!)e dt Systems Linearity: ax(t) + by(t) aX(f) + bY (f) 1 f Time scaling: x(at) jaj X( a ) 1 t A discrete system T transforms a sequence: fyng = T fxng. Frequency scaling: jaj x( a ) X(af) −2πjf∆t Causal systems cannot look into the future: yn = Time shifting: x(t − ∆t) X(f)e f(xn; xn−1; xn−2;:::) 2πj∆ft Frequency shifting: x(t)e X(f − ∆f) Memoryless systems depend only on the current input: yn = R 1 2 R 1 2 Parseval’s theorem: 1 jx(t)j dt = −∞ jX(f)j df f(xn) Continuous convolution: Ff(f ∗ g)(t)g = Fff(t)gFfg(t)g Delay systems shift sequences in time: yn = xn−d Discrete convolution: fxng ∗ fyng = fzng () A system T is time invariant if for all d: fyng = T fxng () X(ej!)Y (ej!) = Z(ej!) fyn−dg = T fxn−dg 0 A system T is linear if for any sequences: T faxn + bxng = 0 Sample Transforms aT fxng + bT fxng Causal linear time-invariant systems are of the form The Fourier transform preserves function even/oddness but flips PN PM real/imaginaryness iff x(t) is odd.
    [Show full text]
  • The Role of the Wigner Distribution Function in Iterative Ptychography T B Edo1
    The role of the Wigner distribution function in iterative ptychography T B Edo1 1 Department of Physics and Astronomy, University of Sheffield, S3 7RH, United Kingdom E-mail: [email protected] Ptychography employs a set of diffraction patterns that capture redundant information about an illuminated specimen as a localized beam is moved over the specimen. The robustness of this method comes from the redundancy of the dataset that in turn depends on the amount of oversampling and the form of the illumination. Although the role of oversampling in ptychography is fairly well understood, the same cannot be said of the illumination structure. This paper provides a vector space model of ptychography that accounts for the illumination structure in a way that highlights the role of the Wigner distribution function in iterative ptychography. Keywords: ptychography; iterative phase retrieval; coherent diffractive imaging. 1. Introduction Iterative phase retrieval (IPR) algorithms provide a very efficient way to recover the specimen function from a set of diffraction patterns in ptychography. The on-going developments behind the execution of these algorithms are pushing the boundaries of the method by extracting more information from the diffraction patterns than was earlier thought possible when iterative ptychography was first proposed [1]. Amongst these are: simultaneous recovery of the illumination (or probe) alongside the specimen [2, 3, 4], super resolved specimen recovery [5, 6], probe position correction [2, 7, 8, 9], three-dimensional specimen recovery [10, 11, 23], ptychography using under sampled diffraction patterns [12, 13, 14] and the use of partial coherent source in diffractive imaging experiments [15, 16, 17, 18].
    [Show full text]
  • Kernel Deconvolution Density Estimation Guillermo Basulto-Elias Iowa State University
    Iowa State University Capstones, Theses and Graduate Theses and Dissertations Dissertations 2016 Kernel deconvolution density estimation Guillermo Basulto-Elias Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/etd Part of the Statistics and Probability Commons Recommended Citation Basulto-Elias, Guillermo, "Kernel deconvolution density estimation" (2016). Graduate Theses and Dissertations. 15874. https://lib.dr.iastate.edu/etd/15874 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Graduate Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Kernel deconvolution density estimation by Guillermo Basulto-Elias A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Major: Statistics Program of Study Committee: Alicia L. Carriquiry, Co-major Professor Kris De Brabanter, Co-major Professor Daniel J. Nordman, Co-major Professor Kenneth J. Koehler Yehua Li Lily Wang Iowa State University Ames, Iowa 2016 Copyright c Guillermo Basulto-Elias, 2016. All rights reserved. ii DEDICATION To Mart´ın.I could not have done without your love and support. iii TABLE OF CONTENTS LIST OF TABLES . vi LIST OF FIGURES . vii ACKNOWLEDGEMENTS . xi CHAPTER 1. INTRODUCTION . 1 CHAPTER 2. A NUMERICAL INVESTIGATION OF BIVARIATE KER- NEL DECONVOLUTION WITH PANEL DATA . 4 Abstract . .4 2.1 Introduction . .4 2.2 Problem description . .7 2.2.1 A simplified estimation framework . .7 2.2.2 Kernel deconvolution estimation with panel data .
    [Show full text]
  • Wavelet Deconvolution in a Periodic Setting
    J. R. Statist. Soc. B (2004) 66, Part 3, pp. 547–573 Wavelet deconvolution in a periodic setting Iain M. Johnstone, Stanford University, USA G´erard Kerkyacharian, Centre National de la Recherche Scientifique and Universit´e de Paris X, France Dominique Picard Centre National de la Recherche Scientifique and Universit´e de Paris VII, France and Marc Raimondo University of Sydney, Australia [Read before The Royal Statistical Society at a meeting organized by the Research Section on ‘Statistical approaches to inverse problems’ on Wednesday, December 10th, 2003, Professor J. T. Kent in the Chair ] Summary. Deconvolution problems are naturally represented in the Fourier domain, whereas thresholding in wavelet bases is known to have broad adaptivity properties. We study a method which combines both fast Fourier and fast wavelet transforms and can recover a blurred function observed in white noise with O{n log.n/2} steps. In the periodic setting, the method applies to most deconvolution problems, including certain ‘boxcar’ kernels, which are important as a model of motion blur, but having poor Fourier characteristics. Asymptotic theory informs the choice of tuning parameters and yields adaptivity properties for the method over a wide class of measures of error and classes of function. The method is tested on simulated light detection and ranging data suggested by underwater remote sensing. Both visual and numerical results show an improvement over competing approaches. Finally, the theory behind our estimation paradigm gives a complete characterization of the ‘maxiset’ of the method: the set of functions where the method attains a near optimal rate of convergence for a variety of Lp loss functions.
    [Show full text]