Point spread function

Image formation in a confocal microscope: central longitudinal (XZ) slice. The 3D acquired distribution arises from the convolution of the real light sources with the PSF. A point source as imaged by a system with negative (top), zero (center), and positive (bottom) spherical aberration. Images to the left are defocused toward the inside, images on the right toward the outside.

The point spread function (PSF) describes the response of a focused optical imaging system to an idealized point source of light. In casual terms, for a given camera, it is the blurry blob image captured from pointing that camera at a single speck of light.

More technically, a PSF is a form of impulse response function (IRF) for a focused optical imaging system, in spatial terms (as opposed to temporal terms). In functional terms, it is the spatial domain version (i.e., the inverse Fourier transform) of the optical transfer function (OTF) of an imaging system. It is a useful concept in Fourier optics, astronomical imaging, medical imaging, electron microscopy and other imaging techniques such as 3D microscopy (like in confocal laser scanning microscopy) and fluorescence microscopy.

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 non-coherent imaging systems, such as fluorescent microscopes, telescopes or optical microscopes, the image formation process is linear in the image intensity and described by a linear system theory. This means that when two objects A and B are imaged simultaneously by a non-coherent imaging system, the resulting image is equal to the sum of the independently imaged objects. In other words: the imaging of A is unaffected by the imaging of B and vice versa, owing to the non-interacting property of photons. In space-invariant systems, i.e. those in which the PSF is the same everywhere in the imaging space, the image of a complex object is then the convolution of that object and the PSF. The PSF can be derived from diffraction integrals.

O(xo,yo)=∬O(u,v) δ(xo−u,yo−v) dudv

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., PSF(xi/M−u,yi/M−v) is the image of the impulse function (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 Huygens–Fresnel principle). Such a source of uniform spherical waves is shown in the figure below. We also note that a perfect point source radiator will not only radiate a uniform spectrum of propagating plane waves, but a uniform spectrum of exponentially decaying (evanescent) waves as well, and it is these which are responsible for resolution finer than one wavelength (see Fourier optics). This follows from the following Fourier transform expression for a 2D impulse function,

J1(x) is the first-order Bessel function of the first kind. That is, a uniformly-illuminated circular aperture that passes a converging uniform spherical wave yields an Airy disk image at the focal plane. A graph of a sample Airy disk is shown in the adjoining figure.

Airy disk

Therefore, the converging (partial) spherical wave shown in the figure above produces an Airy disc in the image plane. The argument of the function 1</sub>(''x'')/''x''"}},"i":0}}]}'>J1(x)/x is important, because this determines the scaling of the Airy disc (in other words, how big the disc is in the image plane). If Θmax is the maximum angle that the converging waves make with the lens axis, r is radial distance in the image plane, and wavenumber k = 2π/λ where λ = wavelength, then the argument of the function is: max</sub>)"}},"i":0}}]}'>kr tan(Θmax). If Θmax is small (only a small portion of the converging spherical wave is available to form the image), then radial distance, r, has to be very large before the total argument of the function moves away from the central spot. In other words, if Θmax is small, the Airy disc is large (which is just another statement of Heisenberg's uncertainty principle for Fourier Transform pairs, namely that small extent in one domain corresponds to wide extent in the other domain, and the two are related via the space-bandwidth product). By virtue of this, high magnification systems, which typically have small values of Θmax (by the Abbe sine condition), can have more blur in the image, owing to the broader PSF. The size of the PSF is proportional to the magnification, so that the blur is no worse in a relative sense, but it is definitely worse in an absolute sense.

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. 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:

Ahi, Kiarash; Anwar, Mehdi (May 26, 2016). Anwar, Mehdi F; Crowe, Thomas W; Manzur, Tariq (eds.). "Modeling of terahertz images based on x-ray images: a novel approach for verification of terahertz images and identification of objects with fine details beyond terahertz resolution". Proc. SPIE 9856, Terahertz Physics, Devices, and Systems X: Advanced Applications in Industry and Defense, 98560N. Terahertz Physics, Devices, and Systems X: Advanced Applications in Industry and Defense. 9856: 985610. Bibcode:2016SPIE.9856E..10A. doi:10.1117/12.2228685. S2CID 124315172.

↑ 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.

↑ "QLP Data Release Notes 004: TESS-Gaia Light Curve Photometry Implementation". arxiv.org. Retrieved 2026-03-28.

↑ Hagai Kirshner; François Aguet; Daniel Sage; Michael Unser (2013). "3-D PSF Fitting for Fluorescence Microscopy: Implementation and Localization Application" (PDF). Journal of Microscopy. 249 (January 2013): 13–25. doi:10.1111/j.1365-2818.2012.03675.x. PMID 23126323. S2CID 5318333.