Uncertainties in interpolated spectral data.

Interpolation is often used to improve the accuracy of integrals over spectral data convolved with various response functions or power distributions. Formulae are developed for propagation of uncertainties through the interpolation process, specifically for Lagrangian interpolation increasing a regular data set by factors of 5 and 2, and for cubic-spline interpolation. The interpolated data are correlated; these correlations must be considered when combining the interpolated values, as in integration. Examples are given using a common spectral integral in photometry.

Correlation coefficients are developed for Lagrangian interpolation where the input data are uncorrelated. It is demonstrated that in practical cases, uncertainties for the integral formed using interpolated data can be reliably estimated using the original data.

Keywords: interpolation; photometry; radiometry; uncertainty.

1. Introduction

Measurements of spectral irradiance, spectral responsivity or spectral reflectance are often made for a limited set of wavelengths and then used to calculate weighted spectral sums for photometry, colorimetry or filter radiometry. It is often necessary to interpolate the spectral data to a finer grid to avoid errors arising from the discrete approximation used to estimate the integral where the weighting function varies strongly between the wavelengths at which the measurements were made (1). Interpolated values are then correlated to the original data and to nearby interpolated points; unless these correlations are taken into account, uncertainty calculations will give misleading results, generally underestimating errors in the spectral sums by significant amounts.

Interpolation may also be required when reference standards are provided for a limited set of wavelengths. Many of the primary reference standards in the national metrology institutes are derived in a functional form, and can be calculated on as fine a wavelength grid as required, along with all the necessary correlations (2). However, calibrations of client lamps or detectors involve a transfer, by comparison, to artifacts that do not necessarily show a spectral variation that is easily modeled, and in the interests of reducing costs may be provided on a limited number of wavelengths. The reference data at different wavelengths may also be correlated. While the primary reference standards are often strongly correlated between wavelengths, the transfer process itself adds uncertainty that is generally random and often reduces the correlations to negligible levels (2)

Where data are available at sufficient wavelengths to avoid errors due to the sum approximation to the integral, it is preferable to interpolate reference tables, such as the photometric response function [V.sub.lambda] and the colonmetric tristumulus response functions (3) to the wavelengths at which measurements are taken, because those tables contain no uncertainty. However, interpolation of measurement data is often applied. One reason is that software for a calculation may require data at a particular interval, another is that instrument measurement programs may provide limited sets of data.

Following a brief description of uncertainty propagation, this paper is divided into two sections. The first covers Lagrangian interpolation, the second cubic spline interpolation; these are the two most commonly used interpolation methods. In each section, interpolation from a data set which is itself correlated is considered. Various simplifications for practical applications are made, and examples are presented. A conclusion is that in practical terms, uncertainty can be accurately derived from the original data set without a complex calculation of correlations.

2. Propagation of Uncertainty

Uncertainty propagation is described in detail in the ISO Guide to the Expression of Uncertainty in Measurement (4). The uncertainty in a quantity y formed by combining measured quantities [x.sub.i] through the relationship y = f([x.sub.1], [x.sub.2], .[x.sub.N]) is given by

[u.sup.2](y) = [summation over (N/i=1)][([partial]f/[partial][x.sub.i]).sup.2][u.sup.2]([x.sub.i]) + [summation over (N/i=1)][summation over (N/j + i=1)][partial]f/[partial][x.sub.i] u ([x.sub.i], [x.sub.j]), (1)

where u([x.sub.i]) is the uncertainty in [x.sub.i] and u([x.sub.i], [x.sub.j]) is the covariance between [x.sub.i] and [x.sub.j]. For uncorrelated input quantities, the covariance between pairs of variables is zero and Eq. (1) reduces to the "sum of squares" commonly applied. The derivatives [partial]f/[partial][x.sub.i] are sensitivity coefficients for the dependence of y on the various measured quantities. Given that [u.sup.2] ([x.sub.i]) = u([x.sub.i], [x.sub.i]), Eq. (1) can be expressed as

[u.sup.2](y) = [f.sup.T.sub.y][U.sub.x][f.sub.y], (2)

where

f = [([partial]f/[partial][x.sub.1][partial]f/[partial][x.sub.2]..[partia l]f/[partial[x.sub.n]).sup.T] (3)

is a column vector of sensitivity coefficients (T indicates the transpose) and

[U.sub.x] = [u ([x.sub.i], [x.sub.j])] (4)

is the Nx N uncertainty matrix.

Interpolation of spectral data is generally performed to produce a set of values ([p.sub.k], [x.sub.k]) from set ([y.sub.i], [x.sub.i]) where [x.sub.i] is the independent variable (generally wavelength in radiometry). Quantities in the set ([P.sub.k], [x.sub.k]) that depend on the same ([y.sub.i], [x.sub.i]) are correlated through this common dependence. The covariance between two values [P.sub.k] and [P.sub.m] (and, when k = m, the square of the uncertainty) is given by

U ([p.sub.k], [P.sub.m]) = [summation over (N/i=1)[summation over (N/(j=1)][partial][p.sub.k]/[partial][y.sub.i][partial][p.sub.m]/[par tial][y.sub.j] u ([y.sub.i],[y.sub.j]). (5)

In matrix form, this is simply expressed as

U([p.sub.k], [p.sub.m]) = [f.sub.Pk.sup.T][U.sub.y][f.sub.Pm] (6)

It is sometimes convenient to use correlation coefficients rather than covariances, defined as

r ([p.sub.k],[P.sub.m] = u([P.sub.k], [P.sub.m])/[square root of (u([P.sub.k], [P.sub.k]) u ([P.sub.m], [P.sub.m))]. (7)

A matrix of correlation coefficients is square, symmetric about the diagonal and has value 1 in the diagonal elements.

3. Lagrangian Interpolation

We have a tabulated function [y.sub.i] at values [x.sub.i]:

[y.sub.i] = y([x.sub.i]), i = 1....N. (8)

For x in the range [x.sub.m] to [x.sub.m+n], the formula for Lagrangian interpolation is (5)

P(x) = [summation over (n/j=1)][[PI].sup.n.sub.i=j=1](x - [x.sub.m+i-1])/[[PI].sup.n.sub.i[not equal to]j=1]([x.sub.m+j-1] - [X.sub.m+i-1]) [y.sub.m+j-1)](9)

= [summation (n/j=1)][W.sub.j][y.sub.m+j-1]

Equation (9) represents an (n-l)th order polynomial fitted through the n original values. Interpolated data are formed as a linear combination of nearby existing data. Sensitivity coefficients for the dependence of the interpolated data on the input data are simply the weights [W.sub.j] in Eq. (9). Covariances between the output values, that is, the original input values plus the interpolated values, take several forms. Correlations present in the original values remain. The values newly formed by interpolation are correlated to both the input values forming them (and through them to the remaining input values if correlations are present in the original set) and to any new values formed from common input values. All of these covariances, including the uncertainty of the interpolated values, can be calculated through Eq. (6).

In many instances, the input data [y.sub.i] are present at regular values of [x.sub.i] and further, the output data are also required at regular intervals. Interpolation is usually then performed by running Eq. (9) through the data and forming a new value or values in the center of the range. The multipliers convolving the data are determined by the interval spacings only and can be calculated prior to the interpolation. This was clearly demonstrated by Savitsky and Golay (6) in the development of their algorithms for smoothing and differentiation of spectral data where polynomial expressions of various order are fitted to regularly-spaced data; in those cases, the coefficients are determined as fixed linear combinations of the input dependent data. The arguments presented here for Lagrangian interpolation can easily be extended to cover smoothing of data using the Savitsky-Golay routines. Two-point Lagrangian interpolation forming a new value in the center of the existing values has weights ([w.sub.1], [w.sub.2])=( 1/2, 1/2), equivalent to a linear interpolation.

Propagation of uncertainty through two common examples of Lagrangian interpolation used in photometry are now discussed. In both of these we consider the calculation of illuminance response to CIE illuminant D65, a tabulated distribution carrying no uncertainty, for a photometer whose spectral response is a close approximation to [V.sub.[lambda]], measured at different spectral intervals and where the measured response values are uncorrelated. These distributions are shown in Fig. 1 for a wavelength interval of 5 nm. As the response function tapers to zero at each end, the luminance response is given as

[R.sub.v] = [DELTA][lambda][summation over (N/(i=1))][R.sub.i][E.sub.D65,i], (10)

where [R.sub.i] is the photometer response and [E.sub.D65,i] is the illuminant value at the ith wavelength, respectively, and [DELTA][lambda] is the 5 urn wavelength separation between the values.

For uncorrelated spectral response values, the uncertainty in [R.sub.v] is given by

[u.sup.2]([R.sub.v]) = [([DELTA][lambda]).sup.2][summation over (i=1/N)][E.sup.2.sub.D65,i][u.sup.2]([R.sub.i]). (11)

For values as tabulated by CIE (3), the value of [R.sub.v] is 10567.41, with an uncertainty 19.60 (relative uncertainty 0.1855%) if the responsivity values have a relative uncertainty of 1% and are uncorrelated. (Note that extra significant figures for the uncertainty are presented above those of normal practice for the purpose of comparison.)

3.1 Photometer Measured at 10 nm Intervals, Interpolated to 5 nm

The spectral integral Eq. (10) for input values on a 10 nm grid evaluates to 10568.18, a small change compared to the 5 nm data due to the discrete approximation to the integral; the relative uncertainty for uncorrelated spectral response values with a relative uncertainty of 1% becomes 0.2623%, or the expected [square root of (2)] increase compared to the response measured at 5 nm intervals. We wish to interpolate the photometer response with a four point Lagrangian function to data on a 5 nm interval. The weights for Eq. (9) are then

[w.sub.1],[w.sub.2],[w.sub.3],[w.sub.4] = -1/16,9/16,9/16,-1/16 (12)

and, as the input data are uncorrelated, the uncertainty for an interpolated value is given by

[u.sub.2](p) = [summation over (4/j=l))][w.sup.2.sub.j][u.sup.2]([y.sub.m+j-l]). (13)

In any given interval spanning four input values, the uncertainty value u([y.sub.i]) is approximately constant, and the interpolated values have an uncertainty approx. 80% of the original values in that range. If we ignore the correlations that have been introduced, the relative uncertainty of the integral evaluates as 0.17%, low compared to the original value and clearly incorrect as it would imply that it could be reduced to zero by repeated interpolation.

A full correlation matrix for the output values is formed as follows. First the interpolation is performed (using linear interpolation in the first and last intervals of the input values). From N original values, we now have (2N-1) values

[y.sub.i] = [y.sub.original,j], iodd,j = (i+1)/2

[y.sub.i] = [w.sub.i][y.sub.i-3] + [w.sub.2][y.sub.i-1] + [w.sub.3][y.sub.i-1][w.sub.4][y.sub.i-3], (14)

4<i<2N-4, i even.

Uncertainties for these values are known for the original data and for the interpolated data from Eq. (13) and these are used to populate the diagonal elements of the uncertainty matrix [U.sub.y]. We then have to populate the elements to the right of the diagonal only before filling to the left of the diagonal by symmetry. As the original input values (now at i odd) are uncorrelated, we have for all i, from Eq. (5),

U([y.sub.i],[y.sub.i+1]) = [w.sub.3][u.sup.2]([y.sub.i+1])

U([y.sub.i],[y.sub.i+3]) = [w.sub.4][u.sup.2]([y.sub.i+3]) (15)

For i even (interpolated values),

[y.sub.i+2] = [w.sub.1][y.sub.i-1] + [w.sub.2][y.sub.i+1] + [w.sub.3][y.sub.i+3] + [w.sub.4][y.sub.i+5]

[y.sub.i+4] = [w.sub.1][y.sub.i+1] + [w.sub.2][y.sub.i+3] + [w.sub.3][y.sub.i+5] + [w.sub.4][y.sub.i+7], (16)

[y.sub.i+6] = [w.sub.1][y.sub.i+3] + [w.sub.2][y.sub.i+5] + [w.sub.3][y.sub.i+7] + [w.sub.4][y.sub.i+9]

and from Eq. (5) we have

u(y.sub.i], [y.sub.i+2]) = [w.sub.1][w.sub.2][u.sup.2] ([y.sub.i-1]) + [w.sub.2][w.sub.3][u.sup.2] ([y.sub.i+1]) + [w.sub.3][w.sub.4][u.sup.2] ([y.sub.i+3])

u(y.sub.i], [y.sub.i+4]) = [w.sub.1][w.sub.3][u.sup.2] ([y.sub.i+1]) + [w.sub.2][w.sub.4][u.sup.2] ([y.sub.i+3])

u(y.sub.i],[y.sub.i+6]) = [w.sub.1][w.sub.4][u.sup.2] ([y.sub.i+3]) (17)

For completeness, similar expressions were applied for the values formed by linear interpolation in the first and last intervals. Uncertainty calculated for the integral including all these correlations was then exactly that calculated with the original data for the 10 nm grid.

One further simplification can be made for uncorrelated input values. In practical terms, nearby values have the same uncertainty. Hence the sets of Eqs. (14) to (17) can be reduced to a matrix of correlation coefficients determined purely from the weights,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where the first row corresponds to an interpolated value. Uncertainty for the integral using the interpolated data set is then found by modifying the sensitivity column vector to include the uncertainty at each value,

f = ([partial]f/[parital[y.sub.i] u([y.sub.i])), (19)

and then performing the matrix multiplication Eq. (2). A negligible change relative to the true values is due to averaging through regions where the response is changing rapidly; while the relative uncertainty in these regions is constant, the absolute value is not.

Table 1 shows uncertainties calculated for all the options discussed in this section. It can be seen that proper accounting for the correlations introduced by the interpolation reproduces the uncertainty calculated for the input values alone, and that using the correlation coefficients provides a practical calculation of the uncertainty matrix.

3.2 Photometer Measured at 5 nm Intervals, Interpolated to 1 nm

After re-arranging the N input values to their new positions, we add four values between input values [y.sub.i], and [y.sub.i+5] that can be represented as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Fig. 2 shows the relative uncertainty of the central data set where a photometer with a [V.sub.[lambda]] response is interpolated from 5 nm to 1 nm and the input values are assumed uncorrelated with a relative uncertainty of 1 %. Uncertainties of the interpolated values are lower than those of the input values. If we ignore the correlations introduced by the interpolation, the relative uncertainty of the integral with the D65 illuminant (interpolated to a I nm grid with a cubic-spine routine), reduces to 0.074 %, a value too low by the order of [square root of (5)]. The integral itself has the same value as that shown in Table I for the 5 nm grid.

For a four point Lagrange interpolation adding four values between each of the input values, correlations in the output set extend over the 19 values following an input value. The matrix of correlation coefficients is shown as the transpose relative to Eq. (18) in the interests of printing; that is, it is equivalent to filling to the lower left of the diagonal of the correlation matrix prior to filling the upper right by symmetry. Beginning at a column corresponding to an input value, the matrix of correlation coefficients is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

A number of these, for the values furthest from the input set, are negligible. The relative uncertainty of the integral for the interpolated data set calculated with these correlation coefficients is 0.197 %, equivalent in practical terms to that calculated using the original 5 nm data set and shown in Table 1.

4. Cubic-Spline Uncertainty Propagation

For our set of data Eq. (8), cubic-spline interpolation (7) calculates a value y at x in the interval [x.sub.i] to [x.sub.i+1] as

y = A[y.sub.i]+B[y.sub.i+1] + [Cy.sub.i] + D[y".sub.i+1] (22)

where

A = [x.sub.i+1] - x/[x.sub.i+1] - [x.sub.i], (23)

B = x - [x.sub.i]/[x.sub.i+1] - [x.sub.i], (24)

C = 1/6([A.sup.3] - A)[([x.sub.i+1] - [x.sub.i]).sup.2], (25)

D = 1/6([B.sup.3] - B)[([x.sub.i+1] - [x.sub.i]).sup.2], (26)

The first two terms of Eq. (22) represent simple linear interpolation. Including the second derivatives y" yields a function that has first and second derivatives continuous at the boundaries between intervals.

The second-derivatives are unknown. The relation between them is given by

[x.sub.i] - [x.sub.i-1]/6 [y".sub.i+1] + [x.sub.i+1] - [x.sub.i-1]/3[y".sub.i] + [x.sub.i+1] - [x.sub.i]/6 [y"i+1] =

[y.sub.i+1] - [y.sub.i]/[x.sub.i+1] - [x.sub.i] - [y.sub.i] - [y.sub.i-1]/[x.sub.i] - [x.sub.i+1] (27)

which is a system of N-2 equations in the N unknowns [y".sub.i]. The natural cubic-spline, which is commonly used, sets

[y".sub.i] = [y".sub.N] = 0 (28)

and solves for the remaining terms (7). We are interested in using the cubic-spline interpolation on spectral data of known uncertainties, including the possibility of correlations, where the interpolated data may then be combined in various ways, so that not only the uncertainties in the interpolated data but also the correlations present are important in propagating uncertainties in the combinations.

The (N-2) values of [y".sub.i] depend on each of the input values, i.e., are correlated to each input value. Then even for uncorrelated input data, the output data are correlated over the whole set of interpolated values. We wish to calculate the covariance u ([y.sub.m] [y.sub.m]) between two interpolated values

[y.sub.m] = [A.sub.m][y.sub.i] + [B.sub.m][y.sub.i+1] + [C.sub.m][y".sub.i] + [D.sub.m][y".sub.+i+1],

[y.sub.n] = [A.sub.n][y.sub.j] + [B.sub.n][y.sub.j+1] + [C.sub.n][y".sub.j] + [D.sub.n][y".sub.+j+1], (29)

where the [y.sub.n] and [y.sub.m] values may be in the same or different intervals denoted by i,j. The uncertainty in [y.sub.n] is given by

[u.sup.2]([y.sub.n]) = u([y.sub.n],[y.sub.n]). (30)

The covariance between (and uncertainty of) the input values is known, carried in the matrix

[U.sup.y][u[y.sub.i], [y.sub.j])]. (31)

To propagate uncertainties through Eq. (29) we need the covariance between the second-derivative values and the input values, and between the second-derivatives themselves. These in turn require the sensitivity coefficients for the dependence of the second-derivative on the input values, [partial][y".sub.k]/[partial][y.sub.i]

The set of N-2 equations Eq. (27) can be written in matrix form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where

[h.sub.i,i-1] = [x.sub.i] - [x.sub.i-1]/6

[h.sub.i,i] = [x.sub.i+2] - [x.sub.i]/3

[h.sub.i,i+1] = [x.sub.i+2] - [x.sub.i+1]/6 (33)

is a tri-diagonal matrix with remaining terms zero. Multiplying both sides of Eq. (32) on the left by the inverse matrix [(h.sub.i,j]).sup.-1] yields the required second-derivatives. By selecting the ith row of this multiplication we have the relationship

[y.sub.i+1]= [summation over (N-2/j=1)] [h.sup.-1.sub.i,j]([y.sub.j+2] - [y.sub.j+1]/[x.sub.j+2]-[x.sub.j+1]-[y.sub.j+1]-[y.sub.j]/[x.sub.j+1] -[x.sub.j]) (34)

which when rearranged to

[y.sub.i+1] = [summation over (N/j=3)][h.sup.-1.sub.i,j-2][y.sub.j+2]/[x.sub.j+2] - [x.sub.j+1]

- [summation over (N-1/j=2)][h.sup.-1.sub.i,j-1](1/[x.sub.j+2]-[x.sub.j+1])+(1/[x.sub.j +1]-[x.sub.j])[y.sub.j]

+[summation over (N-2/j=1)][h.sup.-1.sub.i,j] [y.sub.j]/[(x.sub.j+1] - [x.sub.j] (35)

yields

[partial][y.sub.i+1]/[partial][y.sub.j] = [h.sup.-1.sub.i,j-2]/[x.sub.j] - [x.sub.j-1])

-[h.sup.-1.sub.i,j-1](1/[x.sub.j+1]-[x.sub.j] + 1/[x.sub.j] - [x.sub.j-1])

+ [h.sup.-1.sub.i,j]/[x.sub.j+1]-[x.sub.j] (36)

These are augmented by

[partial][y.sub.k]/[partial][y.sub.i] = 0, k =1, N. (37)

If the input [y.sub.i] values are uncorrelated, the required covariances are given by

U([y.sub.i],[y.sub.j]) = [summation over (N/k=1)][partial][y.sub.i]/[partial][y.sub.k] [partial][y.sub.j]/[partial][y.sub.k] [u.sup.2]([y.sub.k]) (38)

and

Note that the N x N matrix represented by Eq. (39) is not symmetric. For correlated inputs, the sums required are the matrix products

U([y".sub.i],[y.sub.j]) = [summation over (N/k=1)][partial][y".sub.i]/[partial][y.sub.k] [partial][y.sub.j]/[partial][y.sub.k][u.sup.2]([y.sub.k]) = [partial][y".sub.i]/[partial][y.sub.k][u.sup.2]([y.sub.k]) (39)

u(y".sub.i],[y".sub.j]) = [f.sup.[tau].sub.i] [U.sub.u][f.sub.j]

u(y".sub.i],[y".sub.j]) = [f.sup.[tau].sub.i] [U.sub.u][f.sub.j] (40)

where

[f.sub.i] = ([partial][y".sub.i]/[partial][y.sub.k]), k = 1..N (41)

is a column-vector of sensitivity coefficients of the second-derivatives vs the input values and [g.sub.j] is a column vector of length N with i in the jth row, 0 elsewhere.

Covariances and uncertainties for the interpolated values are then found by recognizing that the values of [y.sub.m] and [y.sub.n] in Eq. (29) are represented by

[y.sub.m] = [f.sub.m] [([y.sub.i] [y.sub.i+1] [y".sub.i][y".sub.i+1][y.sub.j][y.sub.j+1][y".sub.j][y".sub.j+1]).sup .T]

[y.sub.n] = [f.sub.n] [([y.sub.i] [y.sub.i+1] [y".sub.i][y".sub.i+1][y.sub.j][y.sub.j+1][y".sub.j][y".sub.j+1]).sup .T] (42)

where

[f.sub.m] = ([A.sub.i] [B.sub.i] [C.sub.i] [D.sub.i] 0 0 0 0)

[f.sub.n] = (0 0 0 0 [A.sub.j] [B.sub.j] [C.sub.j] [D.sub.j]) (43)

are sensitivity vectors tar [y.sub.m] and [y.sub.n] against the eighty variables in the vectors shown in Eq. (42). Covariances between the interpolated values are then given by

u([y.sub.m],[y.sub.n]) = [f.sup.T.sub.m][Uf.sub.n], (44)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

with the lower half symmetric about the diagonal. Only half of this matrix is required. The two right quadrants can be separately multiplied on the right by the column vector [([A.sub.j] [B.sub.j] [C.sub.j] [D.sub.j]).sup.T] and the two results combined into a single eight-element column vector for the final multiplication.

4.1 Cubic-Spline Interpolation Examples

Consider again the photometer response curve of Fig. 1, to be integrated over the wavelength range from 360 am to 830 nm (effectively the photometric response to an equal-energy source). The function was interpolated over the same input range (one less value) but shifted by 2.5 nm and the integral recalculated. Similarly, the function was interpolated to 1 nm intervals and the integral recalculated. Table 2 shows the results for these interpolations, where the uncertainty in the integral was calculated based on i % uncertainty in the input values, uncorrelated between values. The consequence of ignoring the correlation between the interpolated values is also shown. Correlations between distant points, introduced through the dependence of the second-derivatives, were negligible (largely because the response curve is relatively smooth), but strong correlations were found between near-neighbours.

Figure 3 shows propagated uncertainties for the interpolation shifting the input by 2.5 nm; for an interpolation to the mid-point, we would expect the interpolated value for a smooth function to be near the mean of the two interval boundaries with a propagated uncertainty 1/[square root of (2)] = 0.71 of the input (but of course correlated to adjacent values).

Figure 4 shows the variation in uncertainty for values interpolated at different positions within the 5 nm interval. This is a practical concern where a wavelength offset may be present in measurements, although in general it is a better practice to retain the wavelength values of the measured points and interpolate weighting functions such as the illuminant or the colorimetric response functions as these carry no uncertainty. At the input points, uncertainties equal that of the input; at the midpoints they fall to 1/[square root of (2)] of the input points. Again this is expected as the input data are smooth, and cubic-spline results are not much different from linear interpolation. The cubic-spline reproduces the input ordinate values for abscissa values equal to an input value; for these points, we expect the propagated uncertainty to be that of the input and correlations between similar such values to match that of the input matrix [U.sub.y]. These conditions can be used to test the coding.

Figure 5 shows [V.sub.[lambda]] interpolated from a 20 nm grid to 2 nm. Where the input curve is changing rapidly relative to the magnitude of the data, linear interpolation would be discontinuous and for these regions, the relative uncertainty of the interpolated values rises above that of the 1 % assumed for the input values. This is shown more strongly in Fig. 6 for input data spaced at 40 nm, interpolated to 5 rim. Here the interpolation does not provide a good representation of VA (as shown in Fig. 7); where the interpolation is poor, the uncertainties for the interpolated values rise above those of the input.

5. Conclusion

Interpolation of spectral data is a common occurrence in radiometric and photometric measurements. Those data are often then combined in forming integral values such as photometric or colorimetric responses or filter radiometer responses. Interpolation is particularly important when a relatively smooth curve available only on a wide spectral spacing may need to be convolved with a more-rapidly changing curve and then integrated. Uncertainties in the interpolated values will generally be smaller than those of the input data, although this is not always true in the case of cubic-spline interpolation. Uncertainties in combinations calculated using the interpolated data set then must include correlations introduced by the interpolation Ignoring the correlation will lead to significant underestimation of uncertainties. The calculations in this paper for both Lagrange interpolation and for cubicspline interpolation show that the uncertainty can be reliably estimated, in practical terms, by propagating the uncertai nty through the combination using only the original set of data.

[FIGURE 1 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

Table 1

Integral value and relative uncertainty for various calculation options;
four point Lagrange interpolation adding one value in each range. Input
values uncorrelated, 1 % relative uncertainty.

         Method              Integral  Relative uncertainty

Original data on 5 nm grid   10567.4          0.185
Original data on 10 nm grid  10568.2          0.262
Interpolated data; ignore    10567.5          0.168
correlations
Interpolated data; with      10567.5          0.262
correlations
Interpolated data; use       10567.5          0.266
correlation coefficients

Table 2

Integral of [V.sub.[lambda]] and its uncertainty for input values on a 5
nm grid and 1 % relative uncertainty, uncorrelated

                     Integral  Uncertainty  Uncertainty
                                              ignoring
                                            correlations

Original data        106.8561    0.19647
Shift 2.5 nm         106.8559    0.19647       0.1459
Interpolate to 1 nm  106.8559    0.19647       0.0746

Accepted: December, 19, 2002

6. References

(1.) ASTM E 308-9, Standard practice for computing the colors of objects by using the CIE system, American Society for Testing and Materials, PA (2000).

(2.) J. L. Gardner, Correlations in primary spectral standards, Metrologia, to be published (NEWRAD 2001 proceedings).

(3.) CIE 15.2, Colorimetry, International Commission on Illumination, Vienna (1986).

(4.) Guide to the expression of uncertainty of uncertainty in measurement, International Organization for Standardization, Geneva (1993).

(5.) Numerical Recipes in C: the art of scientific computing, Cambridge Press (1992) P. 108.

(6.) A. Savitsky and M. J. E.Golay, Smoothing and differentiation of data by simplified least squares procedures, Anal. Chem. 36, 1627-1639 (1964).

(7.) Numerical Recipes in C: the art of scientific computing, Cambridge Press (1992) p. 113.

About the author: James L. Gardner is a research fellow with the CSIRO National Measurement Laboratory, Sydney, Australia and was a Guest Researcher at NIST in the Optical Technology Division of the NIST Physics Laboratory The National Institute of Standards and Technology is an agency of the Technology Administration, US. Department of Commerce.