ASTM E2022-16(2021)
(Practice)Standard Practice for Calculation of Weighting Factors for Tristimulus Integration
Standard Practice for Calculation of Weighting Factors for Tristimulus Integration
SIGNIFICANCE AND USE
5.1 This practice is intended to provide a method that will yield uniformity of calculations used in making, matching, or controlling colors of objects. This uniformity is accomplished by providing a method for calculation of weighting factors for tristimulus integration consistent with the methods utilized to obtain the weighting factors for common illuminant-observer combinations contained in Practice E308.
5.2 This practice should be utilized by persons desiring to calculate a set of weighting factors for tristimulus integration who have custom source, or illuminant spectral power distributions, or custom observer response functions.
SCOPE
1.1 This practice describes the method to be used for calculating tables of weighting factors for tristimulus integration using custom spectral power distributions of illuminants or sources, or custom color-matching functions.
1.2 The values stated in SI units are to be regarded as standard. No other units of measurement are included in this standard.
1.3 This standard does not purport to address all of the safety concerns, if any, associated with its use. It is the responsibility of the user of this standard to establish appropriate safety, health, and environmental practices and determine the applicability of regulatory limitations prior to use.
1.4 This international standard was developed in accordance with internationally recognized principles on standardization established in the Decision on Principles for the Development of International Standards, Guides and Recommendations issued by the World Trade Organization Technical Barriers to Trade (TBT) Committee.
General Information
Relations
Standards Content (Sample)
This international standard was developed in accordance with internationally recognized principles on standardization established in the Decision on Principles for the
Development of International Standards, Guides and Recommendations issued by the World Trade Organization Technical Barriers to Trade (TBT) Committee.
Designation: E2022 − 16 (Reapproved 2021)
Standard Practice for
Calculation of Weighting Factors for Tristimulus Integration
This standard is issued under the fixed designation E2022; the number immediately following the designation indicates the year of
original adoption or, in the case of revision, the year of last revision. A number in parentheses indicates the year of last reapproval. A
superscript epsilon (´) indicates an editorial change since the last revision or reapproval.
1. Scope 3.2 Definitions of Terms Specific to This Standard:
3.2.1 illuminant, n—real or ideal radiant flux, specified by
1.1 This practice describes the method to be used for
its spectral distribution over the wavelengths that, in illuminat-
calculating tables of weighting factors for tristimulus integra-
ing objects, can affect their perceived colors.
tionusingcustomspectralpowerdistributionsofilluminantsor
sources, or custom color-matching functions.
3.2.2 source, n—an object that produces light or other
radiant flux, or the spectral power distribution of that light.
1.2 The values stated in SI units are to be regarded as
3.2.2.1 Discussion—A source is an emitter of visible radia-
standard. No other units of measurement are included in this
tion. An illuminant is a table of agreed spectral power
standard.
distribution that may represent a source; thus, IlluminantAis a
1.3 This standard does not purport to address all of the
standard spectral power distribution and Source A is the
safety concerns, if any, associated with its use. It is the
physical representation of that distribution. Illuminant D65 is a
responsibility of the user of this standard to establish appro-
standard illuminant that represents average north sky daylight
priate safety, health, and environmental practices and deter-
but has no representative source.
mine the applicability of regulatory limitations prior to use.
1.4 This international standard was developed in accor- 3.2.3 spectral power distribution, SPD, S(λ),
n—specificationofanilluminantbythespectralcompositionof
dance with internationally recognized principles on standard-
ization established in the Decision on Principles for the a radiometric quantity, such as radiance or radiant flux, as a
function of wavelength.
Development of International Standards, Guides and Recom-
mendations issued by the World Trade Organization Technical
4. Summary of Practice
Barriers to Trade (TBT) Committee.
4.1 CIE color-matching functions are standardized at 1-nm
2. Referenced Documents
wavelength intervals. Tristimulus integration by multiplication
2.1 ASTM Standards:
of abridged spectral data into sets of weighting factors occurs
E284 Terminology of Appearance
at larger intervals, typically 10-nm; therefore, intermediate
E308 PracticeforComputingtheColorsofObjectsbyUsing
1-nm interval spectral data are missing, but needed.
the CIE System
4.2 Lagrange interpolating coefficients are calculated for the
E2729 Practice for Rectification of Spectrophotometric
Bandpass Differences missing wavelengths. The Lagrange coefficients, when multi-
2.2 CIE Standard: plied into the appropriate measured spectral data, interpolate
CIE Standard S 002 Colorimetric Observers the abridged spectrum to 1-nm interval. The 1-nm interval
spectrum is then multiplied into the CIE 1-nm color-matching
3. Terminology
data, and into the source spectral power distribution. Each
3.1 Definitions—Appearance terms in this practice are in
separate term of this multiplication is collected into a value
accordance with Terminology E284.
associated with a measured spectral wavelength, thus forming
weighting factors for tristimulus integration.
This practice is under the jurisdiction of ASTM Committee E12 on Color and
Appearance and is the direct responsibility of Subcommittee E12.04 on Color and
5. Significance and Use
Appearance Analysis.
Current edition approved June 1, 2021. Published June 2021. Originally
5.1 This practice is intended to provide a method that will
approved in 1999. Last previous edition approved in 2016 as E2022 – 16. DOI:
yield uniformity of calculations used in making, matching, or
10.1520/E2022-16R21.
controlling colors of objects. This uniformity is accomplished
For referenced ASTM standards, visit the ASTM website, www.astm.org, or
contact ASTM Customer Service at service@astm.org. For Annual Book of ASTM
by providing a method for calculation of weighting factors for
Standards volume information, refer to the standard’s Document Summary page on
tristimulus integration consistent with the methods utilized to
the ASTM website.
3 obtain the weighting factors for common illuminant-observer
Available from CIE (International Commission on Illumination), http://
www.cie.co.at or http://www.techstreet.com. combinations contained in Practice E308.
Copyright © ASTM International, 100 Barr Harbor Drive, PO Box C700, West Conshohocken, PA 19428-2959. United States
E2022 − 16 (2021)
FIG. 1 The Values of i in Eq 1 are Plotted Above the Abscissa and the Values of r are Plotted Below for A) the First Measurement
Interval; B) the Intermediate Measurement Intervals; and, C) the Last Measurement Interval Being Interpolated
5.2 This practice should be utilized by persons desiring to ~r!~r 2 2!~r 2 3!
L 5 (3)
calculate a set of weighting factors for tristimulus integration 2
who have custom source, or illuminant spectral power
r 2 1 r r 2 3
~ !~ !~ !
L 5 (4)
distributions, or custom observer response functions.
r 2 1 r 2 2 r
~ !~ !~ !
6. Procedure
L 5 (5)
6.1 Calculation of Lagrange Coeffıcients—Obtain by
calculation, or by table look-up, a set of Lagrange interpolating for the cubic case, and to
coefficients for each of the missing wavelengths.
r 2 1 r 2 2
~ !~ !
L 5 (6)
6.1.1 The coefficients should be quadratic (three-point) in
the first and last missing interval, and cubic (four-point) in all
r r 2 2
~ !~ !
intervals between the first and the last missing interval.
L 5 (7)
6.1.2 Generalized Lagrange Coeffıcients—Lagrange coeffi-
~r 2 1!~r!
cients may be calculated for any interval and number of
L 5 (8)
missing wavelengths by Eq 1:
n
~r 2 r ! for the quadratic case. In each of the above equations, as
i
L ~r! 5 , for j 50,1,…n (1)
j )
r 2 r
~ ! many or as few values of r as required are chosen to generate
i50 ifij
j i
the necessary coefficients.
where:
6.1.3.1 Eq 2-8 are applicable when the spectral data are
n = degree of coefficients being
abridged at 10-nm intervals, and the interpolated interval is
calculated,
regular with respect to the measurement interval, presumably
i and j = indices denoting the location
1-nm.
along the abscissa,
6.1.4 Tables 1 and 2 provide both quadratic and cubic
π = repetitive multiplication of
Lagrange coefficients for 10-nm intervals.
the terms in the numerator
6.2 With the Lagrange coefficients provided, the intermedi-
and the denominator, and
ate missing spectral data may be predicted as follows:
indices of the interpolant, r = chosen on the same scale as
n
the values i and j.
P~λ! 5 L m (9)
( i i
6.1.2.1 Fig. 1 assist the user in selecting the values of i, j,
i50
and r for these calculations.
TABLE 1 The Lagrange Quadratic Interpolation Coefficients
6.1.2.2 Eq1isgeneralandisapplicabletoanymeasurement
Applicable to the First and Last Missing Interval for Calculation
interval or interpolation interval, regular or irregular.
of 10-nm Weighting Factors for Tristimulus Integration
6.1.3 10-nm Lagrange Coeffıcients—Where the measured
Index of Missing
spectral data have a regular or constant interval, the equation
Wavelength L L L
0 1 2
reduces to the following:
1 0.855 0.190 –0.045
2 0.720 0.360 –0.080
~r 2 1!~r 2 2!~r 2 3!
3 0.595 0.510 –0.105
L 5 (2)
4 0.480 0.640 –0.120
5 0.375 0.750 –0.125
6 0.280 0.840 –0.120
Hildebrand, F. B., Introduction to Numerical Analysis, Second Edition, Dover,
7 0.195 0.910 –0.105
New York, 1974, Chapter 3.
8 0.120 0.960 –0.080
Fairman, H. S., “The Calculation of Weight Factors for Tristimulus 9 0.055 0.990 –0.045
Integration,” Color Research and Application, Vol 10, 1985, pp. 199–203.
E2022 − 16 (2021)
TABLE 2 The Lagrange Cubic Interpolation Coefficients
6.5 In the four terms on the right-hand side of this equation,
Applicable to the Interior Missing Intervals for Calculation of
the numerical values of the three factors in the brackets are
10-nm Weighting Factors for Tristimulus Integration
known and should be multiplied into a single coefficient. The
Index of Missing
fourthfactor, m,ineachofthefouradditivetermsisassociated
i
Wavelength L L L L
0 1 2 3
with a different measured wavelength.
1 –0.0285 0.9405 0.1045 –0.0165
2 –0.0480 0.8640 0.2160 –0.0320
6.6 Add all multiplicative coefficients dependent upon each
3 –0.0595 0.7735 0.3315 –0.0455
different measured wavelength into a single coefficient appli-
4 –0.0640 0.6720 0.4480 –0.0560
5 –0.0625 0.5625 0.5625 –0.0625 cable to that wavelength. This results in a single set of
6 –0.0560 0.4480 0.6720 –0.0640
weighting factors that then will contain one value for each
7 –0.0455 0.3315 0.7735 –0.0595
measured wavelength in each of three color-matching func-
8 –0.0320 0.2160 0.8640 –0.0480
9 –0.0165 0.1045 0.9405 –0.0285 tions. The partial contribution to the tristimulus value at
wavelength m is:
@~x¯~λ !S~λ !L !1~x¯~λ !S~λ !L !1… #m 5 wt m (13)
where: 0 0 0 1 1 0 0 0 0
6.7 Normalize the weighting factors by calculating the
P = the value being interpolated at interval λ,
L = the Lagrange coefficients, and following normalizing coefficient:
m = the measured abridged spectral values.
k 5 (14)
Because the measured spectral values are as yet unknown, it
S~λ!y¯~λ!
(
may be best to consider this equation in its expanded form:
where:
P λ 5 L m 1L m 1L m 1L m (10)
~ !
0 0 1 1 2 2 3 3
k = the normalizing coefficient,
6.3 Multiply each P(λ) by the 1-nm interval relative spectral
S(λ) = the power in the 1-nm spectrum, and
power of the source or illuminant being considered.
y(λ) = the CIE Y color-matching function.
6.3.1 It may be necessary to interpolate missing values of
6.8 Multiply the weighting factors by k to normalize the set
the source spectral power distribution S(λ), if the source has
to Y = 100 for the perfect reflecting diffuser.
been measured at other than 1-nm intervals.
6.9 Beginning in January of 2010, rectification of bandpass
6.3.2 Doing so results in the following equation:
differences is no longer accomplished by building the correc-
S λ P λ 5 S λ L m 1S λ L m 1S λ L m 1S λ L m (11)
~ ! ~ ! ~ ! ~ ! ~ ! ~ !
0 0 1 1 2 2 3 3
tion factors into a weight set for tristimulus integration. This is
6.4 Multiply the weighted power at each 1-nm wavelength
because to do so fails to correct the spectrum itself and corrects
by the appropriate custom color-matching function value for
onlythetristimulusvalues.Bandpassrectificationisnowunder
that wavelength. Using the CIE color-matching functions as an
the jurisdiction of Practice E2729.
example, obtain the CIE 1-nm data from CIE Standard S 002,
Colorimetric Observers. Doing so results in the following
7. Precision
equation:
7.1 The precision of the practice is limited only by the
x¯ λ S λ P λ 5 x¯ λ S λ L m 1 x¯ λ S λ L m
~ ! ~ ! ~ ! @ ~ ! ~ ! # @ ~ ! ~ ! #
0 0 1 1 precision of the data provided for the source spectral power
distribution. The CIE color-matching functions are precise to
1@x¯~λ!S~λ!L #m 1@x¯~λ!S~λ!L #m (12)
2 2 3 3
sixdigitsbydefinition.TheLagrangecoefficientsarepreciseto
where:
seven digits.
x¯(λ) = the value of the CIE X color-matching function at
wavelength λ, and the calculations are carried out for
8. Keywords
each of the three CIE color-matching functions, x¯(λ),
8.1 color-matching functions; illuminant; illuminant-
y¯(λ), and z¯(λ).
observer weights; source; tristimulus weighting factors
E2022 − 16 (2021)
APPENDIXES
(Nonmandatory Information)
X1. EXAMPLE OF THE CALCULATIONS
X1.1 TableX1.1givesthespectralpowerdistribution(SPD) same three spectral regions. Tables X1.4-X1.6 illustrate how
of a typical 3-band fluorescent lamp with a correlated color
the same measured data, used to interpolate the missing
temperature of about 3000K. The first step is to multiply each
reflectance data in several different intervals, can be combined
value of the SPD by the appropriate CIE color matching
with the illuminant-color matching function product to form a
function (y¯ in this case), wavelength by wavelength, which is
single weight at a single measurement point. Finally, Table
shown inTable X1.2 for three spectral regions: near 360 nm,
X1.7 shows the resulting weight set for this 3000K source and
560 nm, and 830 nm. Table X1.3 shows a typical interpolation
the 1964 10° color matching functions. Table X1.7 is compat-
of a measured reflectance curve from a 10-nm reported interval
ible with Tables 5 in Practice E308.
to the 1-nm interval that matches the SPD-y¯ product in the
E2022 − 16 (2021)
TABLE X1.1 Spectral Power Distribution of Typical 3-Band Fluorescent Lamp with Correlated Color Temperature of 3000 K (1-nm
measurement interval)
λ SPD λ SPD λ SPD λ SPD λ SPD λ SPD
360 0.004880 450 0.014870 540 0.162400 630 0.111200 720 0.004410 810 0.000000
361 0.004595 451 0.015040 541 0.277600 631 0.102900 721 0.003505 811 0.000000
362 0.004310 452 0.015210 542 0.392800 632 0.094620 722 0.002600 812 0.000000
363 0.020290 453 0.014980 543 0.353900 633 0.062350 723 0.002470 813 0.000000
364 0.036270 454 0.014750 544 0.315100 634 0.030080 724 0.002340 814 0.000000
365 0.047350 455 0.014370 545 0.429800 635 0.027420 725 0.002375 815 0.000000
366 0.058440 456 0.014000 546 0.544600 636 0.024770 726 0.002410 816 0.000000
367 0.031870 457 0.014060 547 0.383500 637 0.023050 727 0.002450 817 0.000000
368 0.005300 458 0.014110 548 0.222500 638 0.021330 728 0.002490 818 0.000000
369 0.004700 459 0.013930 549 0.182100 639 0.020750 729 0.001795 819 0.000000
370 0.004100 460 0.013760 550 0.141700 640 0.020170 730 0.001100 820 0.000000
371 0.003785 461 0.013470 551 0.113500 641 0.019920 731 0.001120 821 0.000000
372 0.003470 462 0.013180 552 0.085290 642 0.019660 732 0.001140 822 0.000000
373 0.003540 463 0.013470 553 0.070050 643 0.019740 733 0.001750 823 0.000000
374 0.003610 464 0.013750 554 0.054810 644 0.019810 734 0.002360 824 0.000000
375 0.003615 465 0.014000 555 0.046030 645 0.019550 735 0.002190 825 0.000000
376 0.003620 466 0.014250 556 0.037250 646 0.019280 736 0.002020 826 0.000000
377 0.004210 467 0.013810 557 0.034310 647 0.019080 737 0.
...
Questions, Comments and Discussion
Ask us and Technical Secretary will try to provide an answer. You can facilitate discussion about the standard in here.