Documentation – Previous v2021.1

Diffraction and Optimum Aperture

Lens aberrations and diffraction are two basic factors that limit lens sharpness. Details regarding these basis factors are provided in the following sections.

Lens (Optical) Aberrations 

Imperfections in optical systems arise from a number of causes that include different bending of light at different wavelengths, the inability of spherical surfaces (or even ashperes) to provide clear images over large fields of view, changes in focus for light rays that don’t pass through the center of the lens, and many more (i.e., coma, stigmatism, spherical aberration, and chromatic aberration (longitudinal and lateral)). Aberration correction is the primary purpose of sophisticated lens design and manufacturing, and it is what distinguishes an excellent from a mediocre optical design.

Optical aberrations (sometimes called Seidel Aberrations after a low-order mathematical model) can be broken down into five or six individual components: . Imatest cannot separate them— that requires an (expensive) optical lens tester.

Lens aberrations tend to increase (become more difficult to correct) at large apertures (small f-numbers) and vary greatly for different lenses, even among different samples of the same lens; i.e., the quality control of mass-produced lenses does not always produce lenses of identical quality

Diffraction 

Diffraction is a fundamental physical property that blurs images. It is caused by the bending of light waves near boundaries. The smaller the aperture (the larger the f-number), the worse the diffraction blur. Since diffraction is a fundamental physical effect, it is the same for all lenses. Lens performance does not vary significantly at small apertures (large f-numbers). The f-number is equal to the focal length divided by its aperture diameter. In the classic f-stop sequence

{ 1 1.4 2 2.8 4 5.6 8 11 16 22 32 45 64 …}

each stop admits half the light of the previous stop, while the f-stop number is multiplied by the square root of 2 (1.414). When a photographer says, “I increased the exposure by one f-stop,” then the sequence is decreased by one step; e.g., the aperture changes from f/8 to f/5.6.

Aberrations tend to get worse at large apertures (small f-numbers), while diffraction gets worse at small apertures (large f-numbers).

As a result, lenses tend to have an optimum aperture where they are sharpest, typically around 2 to 3 f-stops below a lens’ maximum aperture (often less for premium lenses). The optimum is fairly broad.

 
You can find the optimum aperture by running a batch of images (
Imatest SFRplus or eSFR ISO recommended) taken at different apertures, then entering the combined output (a CSV file) into Batchview (Figure 1).

Figure 1. Set of images taken at f/4.5 to f/22

In Figure 1, bars show MTF50 in line widths per picture height (LW/PH) for the weighted mean (black), center area (red), part-way area (green), and corners (blue). The procedure is described in detail here. For the lens results in Figure 1, the optimum aperture is 2 to 3 stops from the maximum and edge sharpness is unimpressive.

In Figure 2, diffraction-limited MTF is displayed as a pale brown dotted curve in the MTF figures produced by SFRSFRplus, and eSFR ISO when the pixel spacing (usually in microns) has been manually entered in the appropriate dialog box.

Figure 2. MTF plots showing diffraction-limited MTF (as a pale brown dotted line)

In this figure, the curve on the left is for the Canon EOS-40D (5.7-micron pixel spacing) with the lens set at f/22, which is a small aperture that should only be used when large depth of field is required and sharpness can be sacrificed.

Diffraction causes a point of light to spread out, forming a circular pattern called an Airy disk. We’ll leave the full equation for Wikipedia (it’s complicated; it involves Bessel functions). But the radius of the disk (center to first null) is quite simple.

\(\displaystyle r_{airy} = 1.22 \ \lambda N\)  where λ is the wavelength and N is the F-number (f-stop) of the lens.

The equation for diffraction-limited MTF can be found in Diffraction Modulation Transfer Function from the SPIE OPTIPEDIADavid Jacobson’s Lens Tutorial, and normankoren.comThe diffraction cutoff frequency is 

\(\displaystyle f_{\mathit{cutoff}} = \frac{1}{\lambda N}\)

CIE 1931 chromaticity diagram showing colors & wavelengths

λ is typically 0.555 microns (0.000555 mm) for visible light (yellow-green; the approximate wavelength where the eye is most sensitive), but it can be changed for cameras with different spectral response (such as Infrared). 

Let

\(\displaystyle s=\frac{f}{f_{\mathit{cutoff}}} = f \lambda N\)

The diffraction-limited MTF is

\(\displaystyle \mathit{MTF}_{\mathit{diffraction}-\mathit{ltd}}(s) = \frac{2}{\pi}\bigl(\arccos (s) -s \sqrt{1-s^2} \bigr) \text{     for } s < 1\)

\(\displaystyle \mathit{MTF}_{\mathit{diffraction}-\mathit{ltd}}(s) = 0 \text{     for } s \geq 1 \)

Figure 2 is also shown with a Data Cursor Datatip, which allows you to examine plot or image pixel values. The Data Cursor Datatip is available in all figures and interactive GUIs in Imatest 3.6+.

Lens MTF response can never exceed the diffraction-limited response, but system MTF response often exceeds it at medium spatial frequencies as a result of sharpening, which is (and should be) present in most digital imaging systems.

In addition to lens response, system MTF response is affected by the sensor (which has a null at 1 cycle/pixel), the anti-aliasing filter (designed to suppress energy above 0.5 cycles/pixel), and signal processing (which can be very complex and can be different in different regions of an image).

Diffraction, pixel size, and test chart requirements

Imatest MTF tests are performed by photographing and analyzing test charts— most often slanted-edge charts. We have made a strong effort to ensure that chart quality is sufficient to give accurate results. We recently (August 2021) realized that because of the diffraction limit, chart size requirements can be relaxed somewhat for cameras with tiny pixels, roughly defined as under 1.2 microns. 

In developing the procedure for dealing with diffraction, we note that the spatial frequencies where diffraction-limited MTF equals specified values (fnn for MTF = nn%) are a simple fractions of \(f_{\mathit{cutoff}} =1/(\lambda N)\)

\( f_{90} = 0.0786 /(\lambda N); \ \ \ \ \ f_{70} = 0.237/(\lambda N); \ \ \ \ \ f_{50} = 0.404/(\lambda N); \ \ \ \ \ f_{30} = 0.585 /(\lambda N)\)

For pixel size p,  \(f_{Nyq} = 1/(2p)\)

If \(f_{30} < f_{Nyq}\), the best possible sharpness of the lens is less than we normally assume. This allows Sensor height in pixels (the left y-axis in the figures in Test chart suitability for MTF measurements) to be multiplied by \(f_{30} / f_{Nyq}\), which reduces the recommended chart height.

The full procedure explained and illustrated for a 108MP sensor with really tiny 0.7 μm pixels and an f/2 aperture, in Test chart suitability for MTF measurements.

Diffraction, pixel response limits, and Q 

The question naturally arises, “When is a camera’s resolution limited by the optics, and when it it limited by the sensor pixel size?” In other words,

What is the boundary between diffraction-limited and sensor pixel size-limited optical systems?

The question can be addressed with the concept of Q, presented in “Modeling the Imaging Chain of Digital Cameras” by Robert Feite, SPIE Digital Library, 8.3.,    (Special thanks to J. Gordon Arkenberg for pointing this out to me.)

Essentially, Q is the ratio between the sensor cutoff frequency frequency \(f_\mathit{sensor-cutoff} = 1/p\), which is twice the Nyquist frequency fNyq, and the diffraction-limited cutoff frequency \(f_{\mathit{cutoff}} =1/(\lambda N)\).  

\(\displaystyle Q = \frac{f_\mathit{sensor-cutoff}}{2 f_{Nyq}} = \frac{\lambda (f/\#)}{p} = \frac{\lambda N}{p} = 2 \lambda N f_{Nyq} \)

When Q = 2, \(f_{\mathit{cutoff}} =1/(\lambda N) = 1/(2p) = f_{Nyq}\), which is the highest frequency where image energy can be reproduced correctly, without aliasing. 

To quote Feite, “Cameras with Q < 2 have a higher MTF but introduce aliasing, whereas cameras with Q > 2 have a lower MTF and produce blurrier images, but have no aliasing.” So for Feite, Q = 2 is the dividing line. In Imatest’s (mostly NLK’s) experience, a small amount of aliasing improves sharpness, but does not introduce unacceptable visual artifacts, so a lower Q is acceptable. I consider a system with MTF @ fNyq = 0.3 (30%) to have acceptable visual artifacts but potentially excellent perceptual sharpness. At this level,

\(\displaystyle  f_{30} = 0.585/(\lambda N) = f_{Nyq} = 1/(2 p)\) ;   hence

\(\lambda N / p = Q = 2 \times 0.585 = 1.17\)

This is significantly lower than Q = 2 , which ensures zero aliasing at the expense of perceptual sharpness. Q should not be thought of as a quality factor.

When the imaging system has no aliasing, i.e., the diffraction limit ≤ the Nyquist frequency (Q = 2),
the radius of the Airy disk will be nearly five times larger than the pixel spacing.
The image will not appear to be sharp at the pixel level.

Note that the Airy disk diameter is \(d_{airy} = 2 \ r_{airy} = 2.44 \ \lambda N = 2.44\ Q p\).  At Q = 2 (no aliasing),  dairy is 4.88 pixels wide. At Q = 1 (where the diffraction and sensor cutoff frequencies are identical), dairy is still 2.44 pixels.

Q is closely related to Depth of Field (DoF), which will not be discussed in detail on this page because it involves lens focal length and camera Field of View (which is used in defining the FoV limits). Briefly, large DoF requires small apertures with high values of Q, which may limit the sharpness.

While Q is not of direct concern to Imatest testing, it affects lens and system design and the interpretation of Imatest results. Examples:

  • Cameras used in flatbed scanners, bar code readers, and certain medical devices (e.g., endoscopes) may use fixed-focus lenses and require large depth of field. They often have built-in light sources. The small apertures (large f/numbers) allow relatively little light to enter the lens, but the built-in light sources provide adequate illumination. These systems tend to have high values of Q. Their expected sharpness (when they’re in focus) can be considerably lower than the pictorial cameras described below.
  • Cinema cameras often require very narrow DoF (to keep backgrounds uncluttered). This means large aperture lenses (set to large apertures) that typically have small values of Q.
  • Most pictorial cameras (consumer cameras) have lenses than can be focused (either auto or manual focus), and are expected to produce sharp images. These include small cameras with large aperture prime lenses (typically around f/2) as well as larger cameras (DSLR and mirrorless) with adjustable apertures. These cameras have values of Q well under 2; often under 1. 

We have been looking at a 108 megapixel camera with 0.7 micron pixel size and an f/2 aperture, which admits a large amount of light. F/2 apertures are relatively common in small (mobile) cameras, and are often diffraction-limited. Increasing the aperture (decreasing the f/number), makes it difficult to achieve diffraction-limited performance because lens aberrations start to dominate. This camera has Q = 0.555 × 2 / 0.7 = 1.59.

Lenses in cameras with large pixels and sensors— mirrorless and DSLR cameras with 1 inch or larger sensors— are rarely diffraction-limited at their maximum aperture. The have low values of Q, normally well under 1 at their maximum aperture, but due to aberrations at large apertures they have have an “optimum aperture”, usually about 2 stops (twice the f/number) below the maximum, after which they tend to slowly approach diffraction-limited performance.