Tutorial Solutions 3 Point Spread Functions

Total Page:16

File Type:pdf, Size:1020Kb

Tutorial Solutions 3 Point Spread Functions Tutorial Solutions 3 Point Spread Functions This set of questions look in detail at the point spread function of lenses and its implications. The details of problems 6 and 7 are beyond the scope of this course but the results are very important since they put produce real practical numbers. The computer simulation programs in problem 8 are strongly recommended as they, hopefully, give insight into what a PSF looks like. dx 3.1 Bessel Functions Given that the expansion of Bessel functions of integer order is, ∞ k n ·2k µ ´ = µ ´ 1 x 2 µ= Jn ´x ∑ · · µ k!Γ´n k 1 k =0 µ ´ µ write out the first three terms of the expansion for J0 ´x and J1 x . Use this expansion to determine the value of µ J1 ´x when x = 0 x µ ´ µ ´ µ= = ! Plot the functions J0 ´x ,J1 x and J1 x x for x 10 10. (Maple would be a good idea). Solution Part a: Firstly note the identity, that for n an integer Γ · µ= ´n 1 n! so for we can write the expansion as: ∞ k n ·2k µ ´ = µ ´ 1 x 2 µ= Jn ´x ∑ · µ k! ´n k ! k =0 So the first three terms are x2 x4 µ = · ::: J ´x 1 0 4 64 x x3 x5 µ = · J ´x 1 2 16 192 So from these expansions, we get that 2 4 µ J1 ´x 1 x x · x = 2 16 192 Department of Physics and Astronomy Revised: August 2000 so that when x = 0, µ J1 ´0 1 0 = 2 Part b: the plots are: µ ! Plot of J0 ´x from 10 10: 1 besj0(x) 0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -10 -8 -6 -4 -2 0 2 4 6 8 10 µ ! Plot of J1 ´x from 10 10: 0.6 besj1(x) 0.4 0.2 0 -0.2 -0.4 -0.6 -10 -8 -6 -4 -2 0 2 4 6 8 10 µ= ! Plot of J1 ´x x from 10 10: Department of Physics and Astronomy Revised: August 2000 0.5 besj1(x)/x 0.4 0.3 0.2 0.1 0 -0.1 -10 -8 -6 -4 -2 0 2 4 6 8 10 3.2 PSF of Simple Lens Given that the Intensity PSF of a lens is given by the scaled Power Spectrum of the pupil function derive an expression for the PSF of a circular lens at geometric focus. This is the full solution to the round lens outlined in lectures. Solution ; µ δ ∞ For a lens with pupil function p´x y then if the object is a -function at , the amplitude in the back focal plane is: ZZ κ ´ ; µ= ´ ; µ ´ · µ u x y Bˆ p s t exp ı xs yt dsdt 2 0 f where κ 2 2 ´ · µ Bˆ = B exp ı x y 0 0 2 f the PSF is the intensity in the back focal plane, given by 2 ; µ=j ´ ; j g´x y u2 x y So for a simple round lens of radius a,wehave 2 2 2 ; µ = · p´x y 1forx y a = 0else If we shift to Polar Coordinates, with ρ θ = ρ θ s = cos & t sin we get that Z Z a 2π κ ; µ= ´ ρ θ · µρ θ ρ ρ θ u2 ´x y Bˆ0 exp ı x cos y sin d d 0 0 f Department of Physics and Astronomy Revised: August 2000 ; µ ´ ; µ We have that p´s t is circularly symmetric, so u2 x y must also be circularly symmetric. Note: if we rotate a circular lens we do not do anything to the PSF! As a consequent, we can ; µ = calculate u2 ´x y along one radial line. Select the line y 0, so we get, Z Z a 2π κ ; µ= ρ ρ θ u2 ´x 0 Bˆ0 exp ı xρcosθ d d 0 0 f We now use the standard result that Z 2π θµ θ = π ´ µ exp ´ırcos d 2 J0 r 0 so if we let κ r = xρ f we get that Z a κ ρ ; µ= π ρ ρ u2 ´x 0 2 Bˆ0 J0 x d 0 f now noting the second standard result that d µ= ´ ´ µµ rJ ´r rJ r 0 dr 1 so by integration we have Z r µ= ´ µ rJ1 ´r J0 t tdt 0 If we substitute, κ α = xρ f we can write κ Z xa 2 f f ; µ= π ´αµ α α u2 ´x 0 2 Bˆ0 J0 d 0 κx This is now in the right form to be integrated, to get κax J1 2 f ; µ= π ´ ˆ u2 x 0 2 B0a κax f ; µ= The PSF is usually normalised so that u2 ´0 0 1 and the phase term associated with Bˆ0 ignored, so that κax 2J1 f ; µ= u2 ´x 0 κax f µ= = = Note the factor of 2, since J1 ´0 0 1 2 from previous question. The system is rotational symmetry, so in two dimensions we have that κar 2J1 f ; µ= u2 ´x y κar f Department of Physics and Astronomy Revised: August 2000 2 2 2 · where r = x y . The intensity PSF is then ¬ ¬ 2 ¬ ¬ κar ¬ ¬ 2J1 f ¬ ¬ ; µ= g´x y κ ¬ ¬ ar ¬ ¬ f which is the expressions given in lectures and in Physics 3 Optics course. (See next question for plots of this function). dx 3.3 Annular Aperture Many telescopes, and mirror objectives have a central stop, giving an annular aperture. Extend the above derivation to give an expression for the PSF of annular lens. Hint: The method is almost identical to the full aperture case except for the limits of integration. For the special case of the central stop being half the diameter of the aperture find the location of the first zero, and compare this with full aperture case. Plot the cross-section of the intensity PSF of both apertures, and comment on the result. Hint: You will been a numerical solution to obtain the location of the first zero, and use Maple to produce the plots. Solution The annular aperture or outer radius a and inner b is described as 2 2 2 2 ; µ = · p´x y 1forb x y a = 0else being of shape a b The amplitude in the back focal plane is given by the scaled Fourier Transform of the pupil function. The main part of the integration is identical to the previous question for the open circular aperture except that the radial integration is between b ! a. The amplitude in tha back focal plane is now Z a κ ; µ= π ρ ρ u2 ´x 0 2 Bˆ0 J0 xρ d b f This can separated and written as, Z Z a κ 0 κ ρ ρ ; µ= π ρ ρ · ρ ρ u2 ´x 0 2 Bˆ0 J0 x d J0 x d 0 f b f Department of Physics and Astronomy Revised: August 2000 µ ´ µ We now note that J1 ´x is even so that xJ1 x is odd. Using this we can rearrange the integral to give Z Z a κ b κ ; µ= π ρ ρ ρ ρ u2 ´x 0 2 Bˆ0 J0 xρ d J0 xρ d 0 f 0 f Both terms here have the same form, that being identical to the form for the normal circular aperture. This can then be integrated to give, ¿ ¾ κax κbx J1 J1 2 f 2 f 4 5 ; µ= π ´ ˆ u2 x 0 2 B0 a κax b κbx f f ; µ= so again, we “spin” about the optical axis, normalise to get u2 ´0 0 1 and ignore the phase term associated with Bˆ0 to get ¾ κar κbr ¿ J1 J1 2 2 f 2 f 4 5 ; µ= u2 ´x y 2 2 a κar b κbr µ ´a b f f and then the intensity PSF is given by 2 ; µ=j ´ ; µj g´x y u2 x y Part b: Look at locations of first zeros. For the full apetrure b = 0, we get the usual case ¬ ¬ 2 ¬ ¬ ¬ ¬ κ ar 2 ¬ ¬ ¬ J1 ¬ αµ f J1 ´ ¬ ¬ ¬ ¬ ; µ= = g´x y κ ¬ ¬ ¬ ar ¬ α ¬ ¬ f κ = α = : π where α = ar f The first zero occurs at 0 1 22 , so, as covered in lectures, the first zero is at λ 0:61 f r = 0 a = Look at the special case of b = a 2. By substitution we get that 2 2 αµ ´α= µ 4 2 J1 ´ a J1 2 ; µ= g´x y 2 2 2 a = µ α α= ´a a 4 4 2 which can we rewritten as: 64 1 1 2 ; µ= ´αµ ´α= µ g´x y J J 2 9 α 1 2 1 The first zero is thus when 1 αµ= ´α= µ J ´ J 2 1 2 1 This has a numerical solution of α : < : π 0 3 144 1 22 so the first zero occurs at λ 0:5 f r 0 a which is narrower then the full aperture case, so the PSF of the annular aperture in narrower, so the resolution is increased.Thisisnot what you would expect! The radial, normalsied, plot is Department of Physics and Astronomy Revised: August 2000 1 Annular(x) 0.9 Full(x) 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 -10 -8 -6 -4 -2 0 2 4 6 8 10 The full aperture is the wider plot with the lower secondard maximas, while the annular aperture has a narrower central peak but higher secondary maximas.
Recommended publications
  • Journal of Interdisciplinary Science Topics Determining Distance Or
    Journal of Interdisciplinary Science Topics Determining Distance or Object Size from a Photograph Chuqiao Huang The Centre for Interdisciplinary Science, University of Leicester 19/02/2014 Abstract The purpose of this paper is to create an equation relating the distance and width of an object in a photograph to a constant when the conditions under which the photograph was taken are known. These conditions are the sensor size of the camera and the resolution and focal length of the picture. The paper will highlight when such an equation would be an accurate prediction of reality, and will conclude with an example calculation. Introduction may be expressed as a tan (opposite/adjacent) In addition to displaying visual information relationship: such as colour, texture, and shape, a photograph (1.0) may provide size and distance information on (1.1) subjects. The purpose of this paper is to develop such an equation from basic trigonometric principles when given certain parameters about the photo. Picture Angular Resolution With the image dimensions, we can make an Picture Field of View equation for one-dimensional angular resolution. First, we will determine an equation for the Since l2 describes sensor width and not height, we horizontal field of view of a photograph when given will calculate horizontal instead of vertical angular the focal length and width of the sensor. resolution: Consider the top down view of the camera (2.0) system below (figure 1), consisting of an object (left), lens (centre), and sensor (right). Let d be the (2.1) 1 distance of object to the lens, d2 be the focal length, Here, αhPixel is the average horizontal field of l1 be the width of the object, l2 be the width of the view for a single pixel, and nhPixel is the number of sensor, and α be the horizontal field of view of the pixels that make up a row of the picture.
    [Show full text]
  • Sub-Airy Disk Angular Resolution with High Dynamic Range in the Near-Infrared A
    EPJ Web of Conferences 16, 03002 (2011) DOI: 10.1051/epjconf/20111603002 C Owned by the authors, published by EDP Sciences, 2011 Sub-Airy disk angular resolution with high dynamic range in the near-infrared A. Richichi European Southern Observatory,Karl-Schwarzschildstr. 2, 85748 Garching, Germany Abstract. Lunar occultations (LO) are a simple and effective high angular resolution method, with minimum requirements in instrumentation and telescope time. They rely on the analysis of the diffraction fringes created by the lunar limb. The diffraction phenomen occurs in space, and as a result LO are highly insensitive to most of the degrading effects that limit the performance of traditional single telescope and long-baseline interferometric techniques used for direct detection of faint, close companions to bright stars. We present very recent results obtained with the technique of lunar occultations in the near-IR, showing the detection of companions with very high dynamic range as close as few milliarcseconds to the primary star. We discuss the potential improvements that could be made, to increase further the current performance. Of course, LO are fixed-time events applicable only to sources which happen to lie on the Moon’s apparent orbit. However, with the continuously increasing numbers of potential exoplanets and brown dwarfs beign discovered, the frequency of such events is not negligible. I will list some of the most favorable potential LO in the near future, to be observed from major observatories. 1. THE METHOD The geometry of a lunar occultation (LO) event is sketched in Figure 1. The lunar limb acts as a straight diffracting edge, moving across the source with an angular speed that is the product of the lunar motion vector VM and the cosine of the contact angle CA.
    [Show full text]
  • The Most Important Equation in Astronomy! 50
    The Most Important Equation in Astronomy! 50 There are many equations that astronomers use L to describe the physical world, but none is more R 1.22 important and fundamental to the research that we = conduct than the one to the left! You cannot design a D telescope, or a satellite sensor, without paying attention to the relationship that it describes. In optics, the best focused spot of light that a perfect lens with a circular aperture can make, limited by the diffraction of light. The diffraction pattern has a bright region in the center called the Airy Disk. The diameter of the Airy Disk is related to the wavelength of the illuminating light, L, and the size of the circular aperture (mirror, lens), given by D. When L and D are expressed in the same units (e.g. centimeters, meters), R will be in units of angular measure called radians ( 1 radian = 57.3 degrees). You cannot see details with your eye, with a camera, or with a telescope, that are smaller than the Airy Disk size for your particular optical system. The formula also says that larger telescopes (making D bigger) allow you to see much finer details. For example, compare the top image of the Apollo-15 landing area taken by the Japanese Kaguya Satellite (10 meters/pixel at 100 km orbit elevation: aperture = about 15cm ) with the lower image taken by the LRO satellite (0.5 meters/pixel at a 50km orbit elevation: aperture = ). The Apollo-15 Lunar Module (LM) can be seen by its 'horizontal shadow' near the center of the image.
    [Show full text]
  • Understanding Resolution Diffraction and the Airy Disk, Dawes Limit & Rayleigh Criterion by Ed Zarenski [email protected]
    Understanding Resolution Diffraction and The Airy Disk, Dawes Limit & Rayleigh Criterion by Ed Zarenski [email protected] These explanations of terms are based on my understanding and application of published data and measurement criteria with specific notation from the credited sources. Noted paragraphs are not necessarily quoted but may be summarized directly from the source stated. All information if not noted from a specific source is mine. I have attempted to be as clear and accurate as possible in my presentation of all the data and applications put forth here. Although this article utilizes much information from the sources noted, it represents my opinion and my understanding of optical theory. Any errors are mine. Comments and discussion are welcome. Clear Skies, and if not, Cloudy Nights. EdZ November 2003, Introduction: Diffraction and Resolution Common Diffraction Limit Criteria The Airy Disk Understanding Rayleigh and Dawes Limits Seeing a Black Space Between Components Affects of Magnitude and Color Affects of Central Obstruction Affect of Exit Pupil on Resolution Resolving Power Resolving Power in Extended Objects Extended Object Resolution Criteria Diffraction Fringes Interfere with Resolution Magnification is Necessary to See Acuity Determines Magnification Summary / Conclusions Credits 1 Introduction: Diffraction and Resolution All lenses or mirrors cause diffraction of light. Assuming a circular aperture, the image of a point source formed by the lens shows a small disk of light surrounded by a number of alternating dark and bright rings. This is known as the diffraction pattern or the Airy pattern. At the center of this pattern is the Airy disk. As the diameter of the aperture increases, the size of the Airy disk decreases.
    [Show full text]
  • Single‑Molecule‑Based Super‑Resolution Imaging
    Histochem Cell Biol (2014) 141:577–585 DOI 10.1007/s00418-014-1186-1 REVIEW The changing point‑spread function: single‑molecule‑based super‑resolution imaging Mathew H. Horrocks · Matthieu Palayret · David Klenerman · Steven F. Lee Accepted: 20 January 2014 / Published online: 11 February 2014 © Springer-Verlag Berlin Heidelberg 2014 Abstract Over the past decade, many techniques for limit, gaining over two orders of magnitude in precision imaging systems at a resolution greater than the diffraction (Szymborska et al. 2013), allowing direct observation of limit have been developed. These methods have allowed processes at spatial scales much more compatible with systems previously inaccessible to fluorescence micros- the regime that biomolecular interactions take place on. copy to be studied and biological problems to be solved in Radically, different approaches have so far been proposed, the condensed phase. This brief review explains the basic including limiting the illumination of the sample to regions principles of super-resolution imaging in both two and smaller than the diffraction limit (targeted switching and three dimensions, summarizes recent developments, and readout) or stochastically separating single fluorophores in gives examples of how these techniques have been used to time to gain resolution in space (stochastic switching and study complex biological systems. readout). The latter also described as follows: Single-mol- ecule active control microscopy (SMACM), or single-mol- Keywords Single-molecule microscopy · Super- ecule localization microscopy (SMLM), allows imaging of resolution imaging · PALM/(d)STORM imaging · single molecules which cannot only be precisely localized, Localization microscopy but also followed through time and quantified. This brief review will focus on this “pointillism-based” SR imaging and its application to biological imaging in both two and Fluorescence microscopy allows users to dynamically three dimensions.
    [Show full text]
  • Mapping the PSF Across Adaptive Optics Images
    Mapping the PSF across Adaptive Optics images Laura Schreiber Osservatorio Astronomico di Bologna Email: [email protected] Abstract • Adaptive Optics (AO) has become a key technology for all the main existing telescopes (VLT, Keck, Gemini, Subaru, LBT..) and is considered a kind of enabling technology for future giant telescopes (E-ELT, TMT, GMT). • AO increases the energy concentration of the Point Spread Function (PSF) almost reaching the resolution imposed by the diffraction limit, but the PSF itself is characterized by complex shape, no longer easily representable with an analytical model, and by sometimes significant spatial variation across the image, depending on the AO flavour and configuration. • The aim of this lesson is to describe the AO PSF characteristics and variation in order to provide (together with some AO tips) basic elements that could be useful for AO images data reduction. Erice School 2015: Science and Technology with E-ELT What’s PSF • ‘The Point Spread Function (PSF) describes the response of an imaging system to a point source’ • Circular aperture of diameter D at a wavelenght λ (no aberrations) Airy diffraction disk 2 퐼휃 = 퐼0 퐽1(푥)/푥 Where 퐽1(푥) represents the Bessel function of order 1 푥 = 휋 퐷 휆 푠푛휗 휗 is the angular radius from the aperture center First goes to 0 when 휗 ~ 1.22 휆 퐷 Erice School 2015: Science and Technology with E-ELT Imaging of a point source through a general aperture Consider a plane wave propagating in the z direction and illuminating an aperture. The element ds = dudv becomes the a sourse of a secondary spherical wave.
    [Show full text]
  • Computation and Validation of Two-Dimensional PSF Simulation
    Computation and validation of two-dimensional PSF simulation based on physical optics K. Tayabaly1,2, D. Spiga1, G. Sironi1, R.Canestrari1, M.Lavagna2, G. Pareschi1 1 INAF/Brera Astronomical Observatory, Via Bianchi 46, 23807 Merate, Italy 2 Politecnico di Milano, Via La Masa 1, 20156 Milano, Italy ABSTRACT The Point Spread Function (PSF) is a key figure of merit for specifying the angular resolution of optical systems and, as the demand for higher and higher angular resolution increases, the problem of surface finishing must be taken seriously even in optical telescopes. From the optical design of the instrument, reliable ray-tracing routines allow computing and display of the PSF based on geometrical optics. However, such an approach does not directly account for the scattering caused by surface microroughness, which is interferential in nature. Although the scattering effect can be separately modeled, its inclusion in the ray-tracing routine requires assumptions that are difficult to verify. In that context, a purely physical optics approach is more appropriate as it remains valid regardless of the shape and size of the defects appearing on the optical surface. Such a computation, when performed in two-dimensional consideration, is memory and time consuming because it requires one to process a surface map with a few micron resolution, and the situation becomes even more complicated in case of optical systems characterized by more than one reflection. Fortunately, the computation is significantly simplified in far-field configuration, since the computation involves only a sequence of Fourier Transforms. In this paper, we provide validation of the PSF simulation with Physical Optics approach through comparison with real PSF measurement data in the case of ASTRI-SST M1 hexagonal segments.
    [Show full text]
  • Topic 3: Operation of Simple Lens
    V N I E R U S E I T H Y Modern Optics T O H F G E R D I N B U Topic 3: Operation of Simple Lens Aim: Covers imaging of simple lens using Fresnel Diffraction, resolu- tion limits and basics of aberrations theory. Contents: 1. Phase and Pupil Functions of a lens 2. Image of Axial Point 3. Example of Round Lens 4. Diffraction limit of lens 5. Defocus 6. The Strehl Limit 7. Other Aberrations PTIC D O S G IE R L O P U P P A D E S C P I A S Properties of a Lens -1- Autumn Term R Y TM H ENT of P V N I E R U S E I T H Y Modern Optics T O H F G E R D I N B U Ray Model Simple Ray Optics gives f Image Object u v Imaging properties of 1 1 1 + = u v f The focal length is given by 1 1 1 = (n − 1) + f R1 R2 For Infinite object Phase Shift Ray Optics gives Delta Fn f Lens introduces a path length difference, or PHASE SHIFT. PTIC D O S G IE R L O P U P P A D E S C P I A S Properties of a Lens -2- Autumn Term R Y TM H ENT of P V N I E R U S E I T H Y Modern Optics T O H F G E R D I N B U Phase Function of a Lens δ1 δ2 h R2 R1 n P0 P ∆ 1 With NO lens, Phase Shift between , P0 ! P1 is 2p F = kD where k = l with lens in place, at distance h from optical, F = k0d1 + d2 +n(D − d1 − d2)1 Air Glass @ A which can be arranged to|giv{ze } | {z } F = knD − k(n − 1)(d1 + d2) where d1 and d2 depend on h, the ray height.
    [Show full text]
  • Modern Astronomical Optics 1
    Modern Astronomical Optics 1. Fundamental of Astronomical Imaging Systems OUTLINE: A few key fundamental concepts used in this course: Light detection: Photon noise Diffraction: Diffraction by an aperture, diffraction limit Spatial sampling Earth's atmosphere: every ground-based telescope's first optical element Effects for imaging (transmission, emission, distortion and scattering) and quick overview of impact on optical design of telescopes and instruments Geometrical optics: Pupil and focal plane, Lagrange invariant Astronomical measurements & important characteristics of astronomical imaging systems: Collecting area and throughput (sensitivity) flux units in astronomy Angular resolution Field of View (FOV) Time domain astronomy Spectral resolution Polarimetric measurement Astrometry Light detection: Photon noise Poisson noise Photon detection of a source of constant flux F. Mean # of photon in a unit dt = F dt. Probability to detect a photon in a unit of time is independent of when last photon was detected→ photon arrival times follows Poisson distribution Probability of detecting n photon given expected number of detection x (= F dt): f(n,x) = xne-x/(n!) x = mean value of f = variance of f Signal to noise ration (SNR) and measurement uncertainties SNR is a measure of how good a detection is, and can be converted into probability of detection, degree of confidence Signal = # of photon detected Noise (std deviation) = Poisson noise + additional instrumental noises (+ noise(s) due to unknown nature of object observed) Simplest case (often
    [Show full text]
  • Beat the Rayleigh Limit: OAM Based Super-Resolution Diffraction Tomography
    Beat the Rayleigh limit: OAM based super-resolution diffraction tomography Lianlin Li1 and Fang Li2 1Center of Advanced Electromagnetic Imaging, Department of Electronics, Peking University, Beijing 100871, China 2 Institute of Electronics, CAS, Beijing 100190, China This letter is the first to report that a super-resolution imaging beyond the Rayleigh limit can be achieved by using classical diffraction tomography (DT) extended with orbital angular momentum (OAM), termed as OAM based diffraction tomography (OAM-DT). It is well accepted that the orbital angular momentum (OAM) provides additional electromagnetic degrees of freedom. This concept has been widely applied in science and technology. In this letter we revisit the DT problem extended with OAM, demonstrate theoretically and numerically that there is no physical limits on imaging resolution in principle by the OAM-DT. This super-resolution OAM-DT imaging paradigm, has no requirement for evanescent fields, subtle focusing lens, complicated post- processing, etc., thus provides a new approach to realize the wavefield imaging of universal objects with sub-wavelength resolution. In 1879 Lord Rayleigh formulated a criterion that a conventional imaging system cannot achieve a resolution beyond that of the diffraction limit, i.e., the Rayleigh limit. Such a “physical barrier” characters the minimum separation of two adjacent objects that an imaging system can resolve. In the past decade numerous efforts have been made to achieve imaging resolution beyond that of Rayleigh limit. Among various proposals dedicated to surpass this “barrier”, a representative example is the well-known technique of near-field imaging, an essential sub-wavelength imaging technology. Near-field imaging relies on the fact that the evanescent component carries fine details of the electromagnetic field distribution at the immediate vicinity of probed objects [1].
    [Show full text]
  • Diffraction Optics
    ECE 425 CLASS NOTES – 2000 DIFFRACTION OPTICS Physical basis • Considers the wave nature of light, unlike geometrical optics • Optical system apertures limit the extent of the wavefronts • Even a “perfect” system, from a geometrical optics viewpoint, will not form a point image of a point source • Such a system is called diffraction-limited Diffraction as a linear system • Without proof here, we state that the impulse response of diffraction is the Fourier transform (squared) of the exit pupil of the optical system (see Gaskill for derivation) 207 DR. ROBERT A. SCHOWENGERDT [email protected] 520 621-2706 (voice), 520 621-8076 (fax) ECE 425 CLASS NOTES – 2000 • impulse response called the Point Spread Function (PSF) • image of a point source • The transfer function of diffraction is the Fourier transform of the PSF • called the Optical Transfer Function (OTF) Diffraction-Limited PSF • Incoherent light, circular aperture ()2 J1 r' PSF() r' = 2-------------- r' where J1 is the Bessel function of the first kind and the normalized radius r’ is given by, 208 DR. ROBERT A. SCHOWENGERDT [email protected] 520 621-2706 (voice), 520 621-8076 (fax) ECE 425 CLASS NOTES – 2000 D = aperture diameter πD πr f = focal length r' ==------- r -------- where λf λN N = f-number λ = wavelength of light λf • The first zero occurs at r = 1.22 ------ = 1.22λN , or r'= 1.22π D Compare to the course notes on 2-D Fourier transforms 209 DR. ROBERT A. SCHOWENGERDT [email protected] 520 621-2706 (voice), 520 621-8076 (fax) ECE 425 CLASS NOTES – 2000 2-D view (contrast-enhanced) “Airy pattern” central bright region, to first- zero ring, is called the “Airy disk” 1 0.8 0.6 PSF 0.4 radial profile 0.2 0 0246810 normalized radius r' 210 DR.
    [Show full text]
  • Improved Single-Molecule Localization Precision in Astigmatism-Based 3D Superresolution
    bioRxiv preprint doi: https://doi.org/10.1101/304816; this version posted April 19, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. Improved single-molecule localization precision in astigmatism-based 3D superresolution imaging using weighted likelihood estimation Christopher H. Bohrer1Ϯ, Xinxing Yang1Ϯ, Zhixin Lyu1, Shih-Chin Wang1, Jie Xiao1* 1Department of Biophysics and Biophysical Chemistry, Johns Hopkins School of Medicine, Baltimore, Maryland, 21205, USA. Ϯ These authors contributed equally to this work. *Correspondence and requests for materials should be addressed J.X. (email: [email protected]). Abstract: Astigmatism-based superresolution microscopy is widely used to determine the position of individual fluorescent emitters in three-dimensions (3D) with subdiffraction-limited resolutions. This point spread function (PSF) engineering technique utilizes a cylindrical lens to modify the shape of the PSF and break its symmetry above and below the focal plane. The resulting asymmetric PSFs at different z-positions for single emitters are fit with an elliptical 2D-Gaussian function to extract the widths along two principle x- and y-axes, which are then compared with a pre-measured calibration function to determine its z-position. While conceptually simple and easy to implement, in practice, distorted PSFs due to an imperfect optical system often compromise the localization precision; and it is laborious to optimize a multi-purpose optical system. Here we present a methodology that is independent of obtaining a perfect PSF and enhances the localization precision along the z-axis.
    [Show full text]