DICOM PS3.14 2023e - Grayscale Standard Display Function

## C.2 Methodology

Step (1)

The Characteristic Curve of the test Display System should be determined with as many measurements as practical (see Section D.1, Section D.2, and Section D.3). Using the Grayscale Standard Display Function, the fractional number of JNDs are calculated for each Luminance interval between equally spaced P-Value steps. The JNDs/Luminance interval may be calculated directly, or iteratively. For example, if only a few JNDs belong to every Luminance interval, a linear interpolation may be performed. After transformation of the grayscale response of the Display System, the Luminance Levels for every P-Value are Li and the corresponding Standard Luminance Levels are Lj; dj specifies the JNDs /Luminance-Interval on the Grayscale Standard Display Function for the given number of P-Values. Then, the JNDs/Luminance interval for the transformed Display Function are

Equation C-1.

r = dj(Li+1 - Li)(Lj+1 + Lj) / ((Li+1 + Li)(Lj+1 - Lj))

Additionally, an iterative method can be used to calculate the number of JNDs per Luminance interval, requiring only the Grayscale Standard Display Function that defines a JND step in Luminance given a Luminance value. This is done by simply counting the number of complete JND steps in the Luminance interval, and then the remaining fractional step. Start at the Luminance low end of the interval, and calculate from the Grayscale Standard Display Function the Luminance step required for one JND step. Then continue stepping from the low Luminance value to the high Luminance value in single JND steps, until the Luminance value of the upper end of the Luminance Range is passed. Calculate the fraction portion of one JND that this last step represents. the total number of completed integer JND steps plus the fractional portion of the last uncompleted step is the fractional number of JND steps in the Luminance interval.

Plot the number of JNDs per Luminance interval (vertical axis) versus the index of the Luminance interval (horizontal axis). This curve is referred to as the Luminance intervals vs JNDscurve. An example of a plot of Luminance intervals vs JNDs is shown in figure C-1. The plot is matched very well by a horizontal line when a linear regression is applied. Figure C-1. Illustration for the LUM and FIT conformance measures

The JNDs/Luminance interval data are evaluated by two statistical measures [C4]. The first assesses the global match of the test Display Function with the Grayscale Standard Display Function. The second measure locally analyses the approximation of the Grayscale Standard Display Function to the test Display Function.

Step (2)

Two related measures of a regression analysis are applied after normal multiple linear regression assumptions are verified for the data [C3].The first measure, named the FITtest, attempts to match the Luminance-Intervals-vs-JNDs curve of the test Luminance distribution with different order polynomial fits. The Grayscale Standard Display Function is characterized by exactly one JND per Luminance interval over the entire Luminance Range. Therefore, ideally, the data of JNDs/Luminance intervals vs index of the Luminance interval are best fit by a horizontal line of a constant number of JNDs/Luminance interval, indicating that both the local and global means of JNDs/Luminance interval are constant over the given Luminance Range. If the curve is better matched by a higher-order curve, the distribution is not closely approximating the Grayscale Standard Display Function. The regression analysis should test comparisons through third-order curves.

The second measure, the Luminance uniformity metric (LUM), analyzes whether the size of Luminance steps are uniform in perceptual size (i.e., JNDs) across the Luminance Range. This is measured by the Root Mean Square Error (RMSE) of the curve fit by a horizontal line of the JNDs/Luminance interval. The smaller the RMSE of the JNDs/Luminance interval, the more closely the test Display Function approximates the Grayscale Standard Display Function on a microscopic scale.

Both the FIT and LUM measures can be conveniently calculated on standard statistical packages.

Assuming the test Luminance distribution passes the FIT test, then the measure of quality of the distribution is determined by the single quantitative measurement (LUM) of the standard deviation of the JNDs/Luminance interval from their mean. Clinical practice is expected to determine the tolerances for the FIT and LUM values.

An important factor in reaching a close approximation of a test Display Function to the Grayscale Standard Display Function is the number of discrete output levels of the Display System. For instance, the LUM measure can be improved by using only a subset of the available DDLs while maintaining the full available output digitization resolution at the cost of decreasing contrast resolution.

While the LUM is influenced by the choice of the number of discrete output gray levels in the Grayscale Standard Display Function, the appropriate number of output levels is determined by the clinical application, including possible gray scale image processing that may occur independently of the Grayscale Standard Display Function standardization. Thus, PS3.14 does not prescribe a certain number of gray levels of output. However, in general, the larger the number of distinguishable gray levels available, the higher the possible image quality because the contrast resolution is increased. It is recommended that the number of necessary output driving levels for the transformed Display Function be determined prior to standardization of the Display System (based on clinical applications of the Display System), so that this information can be used when calculating the transformation in order to avoid using gray scale distributions with fewer output levels than needed.

DICOM PS3.14 2023e - Grayscale Standard Display Function