News: Imatest 2020.1 (March 2020) Shannon information capacity is now calculated from images of the Siemens star, with much better accuracy than the old slantededge measurements, which have been deprecated and replaced with a new method (convenient, but less accurate than the Siemens Star). Siemens Star measurements are the recommended method for calculating information capacity. The Siemens star method was presented at the Electronic Imaging 2020 conference, and published in the paper, “Measuring camera Shannon information capacity from a Siemens star image“, linked from the Electronic Imaging website. The white paper, described below, is much more readable. (See also the Imatest News Post: Measuring camera Shannon information capacity with a Siemens star image.) 
The white paper, “Camera information capacity: a key performance indicator for Machine Vision and Artificial Intelligence systems“, which briefly introduces information theory, describes the camera information capacity measurement, then shows results (including the effects of artifacts) is now available for download. 
Meaning – Acquiring and framing – Running the Star module – Results – Information capacity plot – Difference plot
3D Surface plot – Equations – Summary – Total information capacity
Nothing like a challenge! There is such a metric for electronic communication channels— one that specifies the maximum amount of information that can be transmitted through a channel without error. The metric includes sharpness (bandwidth) and noise (grain in film). And a camera— or any digital imaging system— is such a channel.
The metric, first published in 1948 by Claude Shannon* of Bell Labs, has become the basis of the electronic communication industry. It is called the Shannon channel capacity or Shannon information transmission capacity C , and has a deceptively simple equation. (See the Wikipedia page on the ShannonHartley theorem for more detail.)
\(\displaystyle C = W \log_2 \left(1+\frac{S}{N}\right) = W \log_2 \left(\frac{S+N}{N}\right)\)
W is the channel bandwidth, which corresponds to image sharpness, S is the signal energy (the square of signal voltage; proportional to MTF^{2} in images), and N is the noise energy (the square of the RMS noise voltage), which corresponds to grain in film. It looks simple enough (only a little more complex than E = mc^{2 }), but it’s not easy to apply.
*Claude Shannon was a genuine genius. The article, 10,000 Hours With Claude Shannon: How A Genius Thinks, Works, and Lives, is a great read. There are also a nice articles in The New Yorker and Scientific American. The 29minute video “Claude Shannon – Father of the Information Age” is of particular interest to me it was produced by the UCSD Center for Memory and Recording Research. which I frequently visited in my previous career.
We will describe how to calculate information capacity from images of the Siemens star, which allows signal and noise to be calculated from the same location. This method is also sensitive to artifacts from demosaicing, clipping, and data compression, resulting in a superior measurement of image quality— better than anything used by the imaging industry until now. Technical details are in the green (“for geeks”) boxes. 
Meaning of Shannon information capacity
(We are preparing a white paper with a brief definition of information. We will add a link here when it’s ready, probably late March 2020.)
In electronic communication channels the information capacity is the maximum amount of information that can pass through a channel without error, i.e., it is a measure of channel “goodness.” The actual amount of information depends on the code— how information is represented. But although coding is integral to data compression (how an image is stored in a file), it is not relevant to digital cameras. What is important is the following hypothesis:
I stress that this statement is a hypothesis— a fancy mathematical term for a conjecture. It agrees with my experience and with numerous measurements, but (as of February 2020) it needs more testing (with a variety of images) before it can be accepted by the industry. Now that information capacity can be conveniently calculated with Imatest, we have an opportunity to learn more about it.
The information capacity, as we mentioned, is a function of both bandwidth W and signaltonoise ratio, S/N.
In texts that introduce the Shannon capacity, bandwidth W is often assumed to be the halfpower frequency, which is closely related to MTF50. Strictly speaking, W log_{2}(1+S/N) is only correct for white noise (which has a flat spectrum) and a simple low pass filter (LPF). But digital cameras have varying amounts of sharpening, and strong sharpening can result in response curves with large peaks that deviate substantially from simple LPF response. For this reason we use the integral form of the Shannon equation: \(\displaystyle C = \int_0^W \log_2 \left( 1 + \frac{S(f)}{N(f)} \right) df = \int_0^W \log_2 \left(\frac{S(f)+N(f)}{N(f)} \right) df \) As explained in the paper, “Measuring camera Shannon Information Capacity with a Siemens Star Image”, we must alter this equation to account for the twodimensional nature of pixels by converting it to a double integral, then to polar form, than back to one dimension. The equations are in the green box, below. 
The beauty of the Siemens Star method is that signal power S(f) and noise power N(f) are calculated from the same location (segments with a range of angles and narrow range of radii, where S and N are subject to the same image processing).
A key challenge in measuring information capacity is how to define mean signal power S. Ideally, the definition should be based on a widelyused test chart. For convenience, the chart should be scaleinvariant (so precise chart magnification does not need to be measured). And as we indicated, signal and noise should be measured at the same location.
For different observers to obtain the same result the chart design and contrast should be standardized.To that end we recommend a sinusoidal Siemens star chart similar to the chart specified in ISO 12233:2014/2017, Annex E. Contrast should be as close as possible to 50:1 (the minimum specified in the standard; close to the maximum achievable with matte media). Higher contrast can make the star image difficult to linearize. Lower contrast is acceptable, but should be reported with the results. The chart should have 144 cycles for high resolution systems, but 72 cycles is sufficient for low resolution systems. The center marker (quadrant pattern), used to center the image for analysis, should have 1/20 the diameter of the star diameter.
Acquiring and framing the image
Acquire a wellexposed image of the Siemens star in even, glarefree light. Exposures should be reasonably consistent when multiple cameras are tested. The mean pixel level of the linearized image inside the star should be in the range of 0.16 to 0.36. (The optimum has yet to be determined.)
The center of the star should be located close to the center of the image to minimize measurement errors caused by optical distortion (if present). For automatic centering to work properly the image should be oriented so the edges in the center market are nearly vertical and horizontal.
The size of the star in the image should be set so the maximum spatial frequency, corresponding to the minimum radius r_{min}, is larger than the Nyquist frequency f_{Nyq}, and, if possible, no larger than 1.3 f_{Nyq}, so sufficient lower frequencies are available for the channel capacity calculation. This means that a 144cycle star with a 1/20 inner marker should have a diameter of 14001750 pixels and a 72cycle star should have a diameter of 700875 pixels. For highquality inkjet printers, the physical diameter of the star should be at least 9 (preferably 12) inches (23 to 30 cm).
Other features may surround the chart, but the average background should be close to neutral gray (18% reflectance) to ensure a good exposure (it is OK to apply exposure compensation if needed). The figure on the right shows a typical star image in a 24megapixel (4000×6000 pixel) camera.
Run the Star module
either in Rescharts (interactive; recommended for getting started) or as a fixed, batchcapable module (Star button on the left of the Imatest main window).
In the Star chart settings window, make sure the Calculate information capacity checkbox (near the bottom of the Settings section) is checked. The SNRI settings will be described later. If other settings are correct, press OK.
When OK is pressed the image will be analyzed. Any of several displays can be selected in Rescharts. The table below shows displays that are only available for information capacity measurements.
Main display  Secondary display  Description 
9. Information capacity, SNRI  SNR (ratio)  SignaltoNoise Ratio (S/N) as a function of spatial frequency for the mean segment and up to 8 individual segments 
SNR (dB)  SNR (dB) as a function of frequency for the mean segment, etc.  
Signal, Noise  Signal, noise, and (S+N)/N (dB) as a function of frequency for the mean segment.  
Signal, 10X Noise  Signal, 10X noise, and (S+N)/N (dB) as a function of frequency for the mean segment. Useful for visualizing low levels of noise  
NEQ  Noise Equivalent Quanta as a function of frequency  
10. Difference image (noiseonly, etc.)  Noiseonly (inputnoiseless)  Display noiseonly (with signal removed). This is a remarkable result — possibly the first time that noise has been measured and visualized in the presence of a signal. 
Loss (inputideal)  Input − Lossless (test chart image). Shows data that has been attenuated. Difficult to interpret.  
Input image  Input image (unmodified)  
Noiseless image  Ideal (noiseless) input image (with MTF loss), derived from S_{ideal}.  
Ideal image (no MTF loss)  “Ideal” image with no MTF loss (represents the original test chart).  
Noiseonly (linear)  Noiseonly linearized. Typically darker than the gammaencoded version.  
Input image (linear)  Input image linearized. Typically darker than the gammaencoded version.  
11. 3D Surface plot  Displays a 3D surface plot of signal as a function of angle (on the chart) and spatial frequency in Cycles/Pixel. Up to 8 chart cycles are shown (more would be cluttered and difficult to interpret. The image can be rotated. Note that the rectangular (angle × frequency) display area is actually pieshaped in the chart. A small plot of signal and noise versus frequency and a summary of results is also shown. 
Results
Theww Rescharts displays are specifically designed for information capacity results: 9. Information capacity, SNRI, 10. Inputnoiseless Diff, etc., and 11. 3D Surface plot. Here are results from Star run in Rescharts for a raw image (converted to TIFF with dcraw using the 24bit sRGB preset; gamma ≅ 2.2) for a high quality 24megapixel Micro FourThirds camera.
Information capacity plot
The plot below shows signal, noise, and (Signal+Noise)/Noise (db) for the 24megapixel Micro FourThirds Sony A6000, set at ISO 400.
Signal, Noise, and Shannon information capacity (3.21 bits/pixel) from a
raw image (converted to TIFF) from a highquality 24megapixel Micro FourThirds camera @ ISO 400.
Difference image plot (Inputnoiseless, etc.)
The noiseonly (inputnoiseless difference) plot is of particular interest because images that allow measurement and visualization of noise measured in the presence of a signal (with the sinusoldal star pattern removed) have not been previously available. Because noise is very low, and hence hard to see, at ISO 400, we illustrate noise at ISO 25600 (the maximum for the Micro FourThirds camera) for both TIFF from raw and JPEG images. The Copy image button on the right copies the image to the clipboard, where it can be pasted into an image editor/viewer or the Image Statistics module for further analysis.
Noiseless image for Micro FourThirds camera, raw/TIFF image, ISO 25600.
The image on the right is an incamera JPEG from the same capture at the above image (ISO 25600). It looks very different from the raw/TIFF image because noise reduction is present. The images below are for raw/TIFF and incamera JPEG images from the same camera acquired at ISO 400. 

3D Surface plot
The 3D surface plot allows you to examine small portions of the image in detail.
3D Surface plot for the highquality 24 Megapixel Micro FourThirds camera analyzed above.
To obtain this display, 3D Surface plot calculation (as well as Calculate information capacity) must be set in the settings window. It shows the signal (for the selected channel) as a function of angle and spatial frequency (in Cycles/Pixel), which is inversely proportional to radius. This plot represents a narrow pieslice of the original image, with angular detail at high spatial frequencies greatly enlarged.
A small plot of MTF and noise as a function of spatial frequency is displayed as well as a summary of key results (information capacity, etc.). This plot was motivated by tests on an iPhone 10, where the image appeared to be saturating at low to middle spatial frequencies, but the degree of saturation was difficult to assess by viewing the image. As we can see on the right, saturation is very strong, apparently as a result of some kind of local tone mapping. It is not evident in MTF curves from the star pattern or from the adjacent slanted edges. The iPhone had some Adobe software installed that allowed both raw (DNG) and JPEG software to be captured. We don’t know if this affected the JPEG processing. The image below shows the response of a TIFF file (converted from a DNG raw image from the same iPhone 10). The response is sinusoidal— wellbehaved with no amplitude visible distortion. The information capacity is nearly identical to the distorted JPEG image, where several things are happening: random noise is zero where the image is saturated, but noise as defined by \(N(\phi) = S(\phi)S_{ideal}(\phi)\) (below), is increased by the amplitude distortion (deviation from the sine function). The insensitivity of information capacity to image processing, observed in other cases, is a remarkable result. By comparison, MTF50 and MTF50P is very much higher in the highlyprocessed JPEG image. 3D Surface plot from iPhone 10 TIFF from raw (DNG), 
3D Surface plot from iPhone 10 JPEG, 3D Surface plot from iPhone 10 JPEG, 
Calculating Shannon capacity with Siemens star images
The pixel levels of most interchangeable images (typically encoded in color spaces such as sRGB or Adobe RGB) are gammaencoded. For these files, pixel level ≅ (sensor illumination)^{1/gamma}, where gamma (typically around 2.2) is the the intended viewing gamma for the color space (display brightness = (pixel level)^{gamma}). To analyze these files they must be linearized by raising the pixel level to the gamma power. RAW files usually don’t need to be linearized (if they were demosaiced without gammaencoding, i.e., gamma = 1).
The image of the n_{total} cycle Siemens star is divided into nr = 32 or 64 radial segments and ns =8 (recommended), 16, or 24 angular segments. Each segment has a period (angular length in radians) P = 2πn_{total}/n_{s} and contains n_{k} =n_{total }/ n_{s}cycles and k_{n} signal points, each at a known angular location φ, in the range {0, P}.
We assume that the ideal signal in the segment has the form
\(\displaystyle S_{ideal}(\phi) = \sum_{j=1}^{2} a_j \cos \bigl(\frac{2 \pi j n_k \phi}{P} \bigr) + b_j \sin \bigl(\frac{2 \pi j n_k \phi}{P} \bigr) \)
a and b are calculated using the Fourier series coefficient equations, derived from the Wikipedia Fourier Series page, Equation 1.
\(\displaystyle a_j = \frac{2}{P}\int_P S(\phi) \cos \bigl(\frac{2 \pi j n_k \phi}{P} \bigr) d\phi;\quad b_j = \frac{2}{P}\int_P S(\phi) \sin \bigl(\frac{2 \pi j n_k \phi}{P} \bigr) d\phi\)
where S(φ) is the measured signal (actually, signal + noise) in the segment. [Note that although this equation is not in the ISO 12233:2017 standard, it fully satisfies the intent of Appendix F, Step 5 (“A sine curve with the expected frequency is fitted into the measured values by minimizing the square error.”)]
Noise is \(\displaystyle N(\phi) = S(\phi)S_{ideal}(\phi)\)
The frequency f in Cycles/Pixel of a segment centered at radius r (in pixels) is \(\displaystyle f = \frac{n_{total}}{2 \pi r}\). An interesting consequence of this equation is that it’s easy to locate the Nyquist frequency (0.5 C/P): \(\displaystyle r = \frac{n_{total}}{\pi}\) = 45.8 pixels for n_{total} = 144 cycles.
A small adjustment (not described here) is made in case f is slightly different from the expected value due to centering errors, optical distortion, or other factors.
Signal power is \(\displaystyle P(f) = \sigma^2(S_{ideal}(f))\). Noise power is \(\displaystyle N(f) = \sigma^2(N)\), where σ^{2} is variance (the square of standard deviation). Note that signal + noise power is \(\displaystyle P(f)+N(f) = \sigma^2(S)\). [Note: From the context of Shannon: “Communication in the presence of noise”, we assume that N(f) is the noise measured in the presence of signal S_{ideal}(f); not narrowband noise of frequency f.]
Transforming Shannon’s equation from onedimension to pixels
The full onedimensional equation for Shannon capacity was presented in Shannon’s second paper in information theory, “Communication in the Presence of Noise,” Proc. IRE, vol. 37, pp. 1021, Jan. 1949, Eq. (32). This equation cannot be used directly because the pixels under consideration are twodimensional.
\(\displaystyle C = \int_0^W \log_2 \left( 1 + \frac{S(f)}{N(f)} \right) df = \int_0^W \log_2 \left(\frac{S(f)+N(f)}{N(f)} \right) df \) [Onedimensional; not used]
This equation has to be converted into two dimensions since pixels (here) have units of area. (They have units of distance for linear measurements like MTF.)
\(\displaystyle C = \int \int_0^W \log_2 \left(\frac{S(f_x,f_y)+N(f_x,f_y)}{N(f_x,f_y)} \right) df_x\: df_y \)
where f_{x} and f_{y} are frequencies in the x and ydirections, respectively. In order to evaluate this integral, we translate x and y into polar coordinates, r and θ.
\(\displaystyle C = \int_0^{2 \pi} \int_0^W \log_2 \left(\frac{S(f_r,f_θ)+N(f_r,f_θ)}{N(f_r,f_θ)} \right) f_r \: df_r\: df_θ \)
Since S and N are only weakly dependent on θ, we can rewrite this equation in onedimension.
\(\displaystyle C = 2 \pi \int_0^W \log_2 \left(\frac{S(f)+N(f)}{N(f)} \right) f \: df \)
Very geeky: The limiting case for Shannon capacity. Suppose you have an 8bit pixel. This corresponds to 256 levels (0255). If you consider the distance of 1 between levels to be the “noise”, then the S/N part of the Shannon equation is log_{2}(1+256^{2}) ≅ 16. The maximum bandwidth where information can be transmitted correctly W— the Nyquist frequency— is 0.5 cycles per pixel. (All signal energy above Nyquist is garbage— disinformation, so to speak.) So C = W log_{2}(1+(S/N)^{2}) = 8 bits per pixel, which is where we started. Sometimes it’s comforting to travel in circles.
Summary
 Shannon information capacity C has long been used as a measure of the goodness of electronic communication channels. It specifies the maximum rate at which data can be transmitted without error if an appropriate code is used (it took nearly a halfcentury to find codes that approached the Shannon capacity). Coding is not an issue with imaging.
 C is ordinarily measured in bits per pixel. The total capacity is \( C_{total} = C \times \text{number of pixels}\).
 The channel must be linearized before C is calculated, i.e., an appropriate gamma correction (signal = pixel levelgamma, where gamma ~= 2) must be applied to obtain correct values of S and N. The value of gamma (close to 2) is determined from runs of any of the Imatest modules that analyze grayscale step charts: Stepchart, Colorcheck., Multicharts, Multitest, SFRplus, or eSFR ISO.
 We hypothesize that C can be used as a figure of merit for evaluating camera quality, especially for machine vision and Artificial Intelligence cameras. (It doesn’t directly translate to consumer camera appearance because they have to be carefully tuned to reach their potential, i.e., to make pleasing images). It provides a fair basis for comparing cameras, especially when used with images converted from raw with minimal processing.
 Imatest calculates the Shannon capacity C for the Y (luminance; 0.212*R + 0.716*G + 0.072*B) channel of digital images, which approximates the eye’s sensitivity. It also calculates C for the individual R, G, and B channels as well as the C_{b} and C_{r} chroma channels (from YC_{b}C_{r}).
 Shannon capacity has not been used to characterize photographic images because it was difficult to calculate and interpret. But now it can be calculated easily, its relationship to photographic image quality is open for study.
 Since C is a new measurement, we will be happy to work with companies or academic institutions who can verify its suitability for Artificial Intelligence systems.
Note: Because the slantededge information capacity measurement used prior to Imatest 2020.1 is inaccurate (correlating poorly with the superior Siemens star measurements) is has been deprecated completely. It has been replaced with a better measurement that is, however, only recommended for use with a Siemens star in the center of the image for calculating the total image information capacity (in units of bits/image). 
Information capacity for the total image
The information capacity we have discussed until now is for a single star, typically located near the center of the image.
The next step after finding the information capacity of a pixel is to find the total capacity, C_{total}, for the camera. Unfortunately, it can’t be reliably obtained by multiplying C by the number of megapixels because lens sharpness (MTF response) tends to be nonuniform, typically decreasing with distance from the image center. To get the total information capacity of the image there are two choices.
 Use a grid of Siemens star charts, similar to the grid illustrated in the ISO 2014/2017 standard. This is the preferred method when available.
\(C_{total} = \text{mean}(C_{star}) \times \text{megapixels}\)  Use a chart with multiple slantededges, preferably one of the four Imatest charts with automatic Region of Interest (ROI) detection. This requires capturing a second image.
Slantededge method— We have developed a convenient method of calculating C_{total} using an array of slantededges. The calculation is based on “starequivalent information capacity,” C_{starEq} — the approximate information capacity that would have been calculated from a Siemens star. The procedure sounds complex, but it is easy to preform and quite effective because C_{starEq} closely tracks the star measurement C_{star} (designated as C, above). It does, however, require that a second image be captured.
This method is only recommended for minimallyprocessed images converted from raw capture because it is strongly affected by sharpening and noise reduction (problematic because sharpness and noise are measured at separate locations for slanted edges). It is not recommended for incamera JPEG images, and it should not be used as a primary information capacity measurement.
To measure C_{starEq}, begin by linearizing the image (if required), then measure using low contrast slantededges. (4:1 contrast ratio, specified in ISO 12233, is recommended.) Be sure the edge is exposed well, i.e., not saturated. Note that the edge contrast doesn’t directly affect the results.
Measure the pixel levels of the light and dark parts of the slantededge region (away from the edge), Px_{lt}_{ }and Px_{dk} , and the noise voltages in the two regions, \(N_{lt} = \sigma(Px_{lt})\)_{ }and \(N_{dk} = \sigma(Px_{dk})\). The mean values of the pixel levels and noise are Px_{mean} and N_{mean}, respectively. Px_{mean} should be close to the mean pixel level of the star measurement. Because C_{starEq} is extremely sensitive to saturation, which reduces the measured values of N_{dk} or N_{lt}, it should not be calculated if Px_{lt} > 0.95_{ }or Px_{dk} < 0.03.
CstarEq is calculated by Imatest slantededge sharpness modules when Secondary readout 1 is set to Secondary readout for star contrast [contrast], where the star contrast (typically 50) should be entered in the [contrast] box. You will need to press the Change readout button in the settings window. The results appear in the Edge & MTF plot.
For assumed equivalent star contrast Contrast_{star}, the modulation (Michelson contrast) is \(\displaystyle modulation = Px_{mean} \frac {Contrast_{star}1}{\sqrt{2}(Contrast_{star}+1)}\) The \(\sqrt{2}\) factor causes modulation to have units of RMS voltage: useful because modulation^{2} has units of mean power. The signal power for calculating C_{starEq} is \(\displaystyle S(f) = (modulation \times MTF(f))^2\) The noise power of the assumed star is \(\displaystyle N(f) = N^2_{mean}\) C_{starEq} is calculated by substituting the above two equations into the Shannon capacity equation. \(\displaystyle C = 2 \pi \int_0^W \log_2 \left(\frac{S(f)+N(f)}{N(f)} \right) f \: df \) The total camera information capacity is \(\displaystyle C_{total} = \frac{C_{star}(\text{center}) \times \text{mean}(C_{starEq})}{C_{starEq}(\text{center})} \times \text{megapixels}\) \(\text{mean}(C_{starEq})\) is the mean of C_{starEq} taken over a large number of slantededges distributed over the image 
.
Links
(Historical) R. Shaw, “The Application of Fourier Techniques and Information Theory to the Assessment of Photographic Image Quality,” Photographic Science and Engineering, Vol. 6, No. 5, Sept.Oct. 1962, pp.281286. Reprinted in “Selected Readings in Image Evaluation,” edited by Rodney Shaw, SPSE (now SPIE), 1976.
C. E. Shannon, “A mathematical theory of communication,” Bell Syst. Tech. J., vol. 27, pp. 379–423, July 1948; vol. 27, pp.
623–656, Oct. 1948.
C. E. Shannon, “Communication in the Presence of Noise”, Proceedings of the I.R.E., January 1949, pp. 1021.
The University of Texas Laboratory for Image & Video Engineering is doing some interesting work on image and video quality assessment. Here are some promising papers. Challenging material!
R. Soundararajan and A.C. Bovik, “Survey of information theory and visual quality assessment (Invited Paper),” Signal, Image, and Video Processing, Special Section on Human Vision and Information Theory , vol. 7, no. 3, pp. 391401, May, 2013.
H. R. Sheikh, and A. C. Bovik, “Image Information and Visual Quality,” IEEE Transactions on Image Processing , vol. 15, no. 2, pp. 430 – 444, February, 2006.
K. Seshadrinathan and A. C. Bovik, “An information theoretic video quality metric based on motion models,” Third International Workshop on Video Processing and Quality Metrics for Consumer Electronics , Scottsdale, Arizona, January, 2007.
H. R. Sheikh and A. C. Bovik, “A Visual Information Fidelity Approach to Video Quality Assessment (Invited Paper),” The First International Workshop on Video Processing and Quality Metrics for Consumer Electronics , Scottsdale, AZ, January, 2005.
H. R. Sheikh and A. C. Bovik, “Image information and visual quality,” Proc. IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ’04) , Montreal, Canada, vol. 3, pp. iii – 709712, May, 2004.
Wikipedia – Shannon Hartley theorem has a frequency dependent form of Shannon’s equation that is applied to the Imatest sine pattern Shannon information capacity calculation. It is modified to a 2D equation, transformed into polar coordinates, then expressed in one dimension to account for the area (not linear) nature of pixels.
\(\displaystyle C=\int_0^B \log_2 \left( 1 + \frac{S(f)}{N(f)} \right) df\)