Point spread function
This article needs additional citations for verification. (September 2010) |
The point spread function (PSF) describes the response of a focused optical imaging system to a
The degree of spreading (blurring) in the image of a point object for an imaging system is a measure of the quality of the imaging system. In
Introduction
By virtue of the linearity property of optical non-coherent imaging systems, i.e.,
- Image(Object1 + Object2) = Image(Object1) + Image(Object2)
the image of an object in a microscope or telescope as a non-coherent imaging system can be computed by expressing the object-plane field as a weighted sum of 2D impulse functions, and then expressing the image plane field as a weighted sum of the images of these impulse functions. This is known as the superposition principle, valid for
When the object is divided into discrete point objects of varying intensity, the image is computed as a sum of the PSF of each point. As the PSF is typically determined entirely by the imaging system (that is, microscope or telescope), the entire image can be described by knowing the optical properties of the system. This imaging process is usually formulated by a convolution equation. In microscope image processing and astronomy, knowing the PSF of the measuring device is very important for restoring the (original) object with deconvolution. For the case of laser beams, the PSF can be mathematically modeled using the concepts of Gaussian beams.[3] For instance, deconvolution of the mathematically modeled PSF and the image, improves visibility of features and removes imaging noise.[2]
Theory
The point spread function may be independent of position in the object plane, in which case it is called shift invariant. In addition, if there is no distortion in the system, the image plane coordinates are linearly related to the object plane coordinates via the magnification M as:
- .
If the imaging system produces an inverted image, we may simply regard the image plane coordinate axes as being reversed from the object plane axes. With these two assumptions, i.e., that the PSF is shift-invariant and that there is no distortion, calculating the image plane convolution integral is a straightforward process.
Mathematically, we may represent the object plane field as:
i.e., as a sum over weighted impulse functions, although this is also really just stating the sifting property of 2D delta functions (discussed further below). Rewriting the object transmittance function in the form above allows us to calculate the image plane field as the superposition of the images of each of the individual impulse functions, i.e., as a superposition over weighted point spread functions in the image plane using the same weighting function as in the object plane, i.e., . Mathematically, the image is expressed as:
in which is the image of the impulse function .
The 2D impulse function may be regarded as the limit (as side dimension w tends to zero) of the "square post" function, shown in the figure below.
We imagine the object plane as being decomposed into square areas such as this, with each having its own associated square post function. If the height, h, of the post is maintained at 1/w2, then as the side dimension w tends to zero, the height, h, tends to infinity in such a way that the volume (integral) remains constant at 1. This gives the 2D impulse the sifting property (which is implied in the equation above), which says that when the 2D impulse function, δ(x − u,y − v), is integrated against any other continuous function, f(u,v), it "sifts out" the value of f at the location of the impulse, i.e., at the point (x,y).
The concept of a perfect point source object is central to the idea of PSF. However, there is no such thing in nature as a perfect mathematical point source radiator; the concept is completely non-physical and is rather a mathematical construct used to model and understand optical imaging systems. The utility of the point source concept comes from the fact that a point source in the 2D object plane can only radiate a perfect uniform-amplitude, spherical wave — a wave having perfectly spherical, outward travelling phase fronts with uniform intensity everywhere on the spheres (see
The quadratic
Therefore, the converging (partial) spherical wave shown in the figure above produces an
The figure above illustrates the truncation of the incident spherical wave by the lens. In order to measure the point spread function — or impulse response function — of the lens, a perfect point source that radiates a perfect spherical wave in all directions of space is not needed. This is because the lens has only a finite (angular) bandwidth, or finite intercept angle. Therefore, any angular bandwidth contained in the source, which extends past the edge angle of the lens (i.e., lies outside the bandwidth of the system), is essentially wasted source bandwidth because the lens can't intercept it in order to process it. As a result, a perfect point source is not required in order to measure a perfect point spread function. All we need is a light source which has at least as much angular bandwidth as the lens being tested (and of course, is uniform over that angular sector). In other words, we only require a point source which is produced by a convergent (uniform) spherical wave whose half angle is greater than the edge angle of the lens.
Due to intrinsic limited resolution of the imaging systems, measured PSFs are not free of uncertainty.[4] In imaging, it is desired to suppress the side-lobes of the imaging beam by apodization techniques. In the case of transmission imaging systems with Gaussian beam distribution, the PSF is modeled by the following equation:[5]
where k-factor depends on the truncation ratio and level of the irradiance, NA is numerical aperture, c is the speed of light, f is the photon frequency of the imaging beam, Ir is the intensity of reference beam, a is an adjustment factor and is the radial position from the center of the beam on the corresponding z-plane.
History and methods
The diffraction theory of point spread functions was first studied by
Applications
Microscopy
In microscopy, experimental determination of PSF requires sub-resolution (point-like) radiating sources.
Theoretical models as described above, on the other hand, allow the detailed calculation of the PSF for various imaging conditions. The most compactAstronomy
In
For
For ground-based optical telescopes, atmospheric turbulence (known as astronomical seeing) dominates the contribution to the PSF. In high-resolution ground-based imaging, the PSF is often found to vary with position in the image (an effect called anisoplanatism). In ground-based adaptive optics systems, the PSF is a combination of the aperture of the system with residual uncorrected atmospheric terms.[8]
Lithography
The PSF is also a fundamental limit to the conventional focused imaging of a hole,[9] with the minimum printed size being in the range of 0.6-0.7 wavelength/NA, with NA being the numerical aperture of the imaging system.[10][11] For example, in the case of an EUV system with wavelength of 13.5 nm and NA=0.33, the minimum individual hole size that can be imaged is in the range of 25-29 nm. A phase-shift mask has 180-degree phase edges which allow finer resolution.[9]
Ophthalmology
Point spread functions have recently become a useful diagnostic tool in clinical
See also
- Airy disc
- Circle of confusion, for the closely related topic in general photography.
- Deconvolution
- Encircled energy
- Impulse response function
- Microscope
- Microsphere
- PSF Lab
References
- ISBN 978-0-08-055768-7.
- ^ S2CID 114994724.
- S2CID 124315172.
- .
- S2CID 11781848.
- ^ Light transmitted through minute holes in a thin layer of silver vacuum or chemically deposited on a slide or cover-slip have also been used, as they are bright and do not photo-bleach.
S. Courty; C. Bouzigues; C. Luccardini; M-V Ehrensperger; S. Bonneau & M. Dahan (2006). "Tracking individual proteins in living cells using single quantum dot imaging". In James Inglese (ed.). Methods in Enzymology: Measuring biological responses with automated microscopy, Volume 414. Academic Press. pp. 223–224. ISBN 978-0-12-182819-6.
- ^
P. J. Shaw & D. J. Rawlins (August 1991). "The point-spread function of a confocal microscope: its measurement and use in deconvolution of 3-D data". Journal of Microscopy. 163 (2): 151–165. S2CID 95121909.
- ^ "POINT SPREAD FUNCTION (PSF)". www.telescope-optics.net. Retrieved 2017-12-30.
- ^ a b The Natural Resolution
- ^ Principles and Practice of Light Microscopy
- ^ Corner Rounding and Line-end Shortening
- S2CID 21971504.
- Hagai Kirshner; François Aguet; Daniel Sage; Michael Unser (2013). "3-D PSF Fitting for Fluorescence Microscopy: Implementation and Localization Application" (PDF). S2CID 5318333.
- Rachel Noek; Caleb Knoernschild; Justin Migacz; Taehyun Kim; Peter Maunz; True Merrill; Harley Hayden; C.S. Pai; Jungsang Kim (2010). "Multi-scale Optics for Enhanced Light Collection from a Point Source" (PDF). S2CID 6838852.