Theoretical evaluation of the impact of finite intervals in the measurement of the bidirectional reflectance distribution function

The bidirectional reflectance distribution function (BRDF) is a well-defined quantity that allows the bidirectional reflectance of surfaces to be described. However, the light propagation in a specific direction cannot be perfectly realized in practice, because the physical apertures are not infinitesimal but finite. Consequently, a BRDF measurement cannot be considered fully bidirectional, although the measure coincides with the BRDF within a certain confidence level. To properly understand the impact of the finite apertures on the BRDF measures, the deviation between the real BRDF and the BRDF to be obtained using real experimental conditions with finite apertures was theoretically studied for surfaces with realistic BRDFs. The biconical reflectance factor was used to estimate these “measured BRDFs” in different geometrical conditions, and a family of negative exponential functions was defined to assess the impact on surfaces with different angular scattering distributions. An expression for estimating the relative error from finite apertures is given, which considers the angular variation of the BRDF and the different solid angles involved in the measurement.


Introduction
The complete description of the reflection of the radiant flux by an elementary surface (dA) requires accounting for the partial radiant flux reflected (dU r ) in every direction (r r ), which in turn depends on the incident radiant flux (dU i ) and its direction (r i ). The first complication arises from the fact that the radiant flux within an infinitesimally specified direction is zero. This makes it impossible to define the bidirectional reflectance q(r i ; r r ) as the ratio between the reflected radiant flux and the incident radiant flux for given incidence and reflection directions. Let the directional radiant fluxes be defined as: where L i and L r are the radiances related to the incident and reflected radiant fluxes within the irradiation and collection solid angles dx i and dx r , respectively, and dX i and dX r are their corresponding projected solid angles, calculated as dX = dx cos h, h being the angle of the incidence (i subscript) or reflection (r subscript) directions with respect to the normal to elementary surface, dA. Notice that, whereas dU i only depends on the incidence direction, dU r depends on both incidence and reflection directions, because the radiance obtained at the reflection direction is a function of the incidence direction. In addition, L r also depends on L i , and it is explicitly written in equation (2). The bidirectional reflectance expressed as: does not exist, or, in other words: The bidirectional reflectance would depend on the geometry of the measurement, and not only on surface properties. Fred E. Nicodemus, aware of this inconvenience, proposed an approach to define a bidirectional reflectance quantity based on the following considerations on the involved quantities: 1 1. L r is linear with respect to L i : For given incidence and reflection directions, L r is proportional to L i (L r µ L i ). It is fulfilled as long as the irradiation is not so high to produce nonlinear effects. 2. L r is linear with respect to dX i : For given incidence and reflection directions, L r is proportional to dX i (L r µ dX i ), as long as ''dX i is chosen small enough that there is no significant difference in the radiant absorptance for all rays incident on dA within the solid angle dx i .'' 1 3. q is invariant with respect to L i and dX i : As a consequence of the two previous considerations, for given incidence and reflection directions, q, as defined in equation (3), is invariant with respect to L i and dX i . Since L r = f r L i dX i (from i and ii, where f r is a constant for a given geometry), then q = f r dX r . 4. q is linear with respect to dX r : For given incidence and reflection directions, q is proportional to dX r (q µ dX r ), which ''holds only for values of dX r small enough so that the partial reflectance does not change significantly with direction for rays within the solid angle dx r .'' 1 These are the foundations of the quantity proposed by Nicodemus to measure the bidirectional reflectance. He correctly reasoned that, bearing in mind these considerations, the quantity representing the bidirectional reflecting properties of a surface, independent of the geometry of the solid angles used in the measurement, could be defined as the proportionality constant between q and dX r : Nicodemus preferred to use the differential term dL r instead of L r because ''other incident radiation (from other sources in other directions) may also be diffusely reflected (scattered) into the same reflected directions along with that from the source under consideration.'' 1 What Nicodemus actually proposed was to avoid the dependence of the bidirectional reflectance q on the reflection solid angle dX r by normalizing q by dX r . This new quantity, denoted by f r in equation (5), is the well-known bidirectional reflectance distribution function (BRDF), and it is widely used nowadays in diverse fields such as remote sensing, radiometry, or 3D-rendering as the scattering function. The Optics Classification and Indexing Scheme (OCIS) has reserved a code for it, which proves its general acceptance. The monograph on reflectance published by Nicodemus et al. in 1977, 2 where the concept of BRDF is thoroughly explained, has been cited more than 1600 times by scientific publications in this century.
However, there are some aspects of the BRDF that are difficult to understand. Perhaps the most obvious is why the BRDF can reach any positive value, even an infinite value for perfectly specular reflectors. The reason is that f r is defined as a division by a solid angle element [equation (5)]. The magnitude of these solid angle elements must be small enough to contain only directions at which the involved quantities (radiance or bidirectional reflectance) do not change significantly. For glossy surfaces, for which the reflectance changes rapidly approaching the specular peak, dx r needs to be very small in order to meet this requirement, resulting in very high BRDF values near the specular peak [equation (5)]. This partial reflectance is a Dirac delta in the case of perfectly specular reflectors, so the BRDF has no defined value in specular directions.
According to the terminology introduced by Judd 3 and later adopted by Nicodemus 2 , three adjectives are used to describe the reflectance under general geometrical conditions of irradiation and collection: directional (flux confined to an infinitesimal solid angle), conical (flux distributed over a finite solid angle), and hemispherical (flux distributed over the full hemisphere). These forms of reflectances are obtained by integration of f r (r i ;r r ).
Biconical reflectance, with both irradiation and collection in conical geometry, can be regarded as the general case, and the other reflectance coefficients are obtained by using extreme values of either irradiation or collection solid angles. It is calculated as the integrated value of the function f r (r i ;r r ) over the specified solid angles 2 , x i and x r : where r i0 and r r0 are the central incidence and collection directions of the solid angles. Similarly, the biconical reflectance factor (ratio between the reflected radiant flux and the ideally reflected radiant flux by a perfect reflecting diffuser, PRD) is given by 2 : In practice, x i and x r are not completely zero, and a biconical quantity, instead of a bidirectional one, is measured. Therefore, these concepts will be used in the next section to establish the relation between the BRDF and the quantity actually measured. Some of the assumptions in the BRDF definition, as described above (''dX i chosen small enough that there is no significant difference in the radiant absorptance for all rays incident on dA within the solid angle dx i ,'' or ''dX r small enough so that the partial reflectance does not change significantly with direction for rays within the solid angle dx r ''), need to be considered in the measurement of the BRDF, in order to avoid systematic errors or to account for them in the uncertainty budget. In practice, since physical irradiation and collection apertures are finite (neither infinitely big nor infinitesimal), and the measurement area on the surface has a certain size, it is irradiated with radiant fluxes from different directions, and the collected radiant flux comes from different directions too. To understand properly the limitations of a BRDF measurement, the deviation between the real BRDF and the BRDF to be obtained using real experimental conditions with finite apertures needs to be assessed.
The consequence of using a finite measurement area and finite irradiation and collection solid angles is that the BRDF is evaluated as a weighted average over a set of pairs of irradiation and collection directions, which can be modeled with the expression in equation (7). As long as this weighted average coincides with the real BRDF with a negligible deviation, the measurement conditions can be considered adequate.
The following simple examples illustrate the impact of using finite interval. Four different BRDF angular distributions normalized to the value at h asp = 0°are shown in Fig. 1, where h asp denotes the aspecular angle, which is the angular deviation between collec-tion and specular directions. For each distribution, the average BRDF at h asp = 0°was estimated by integration, with a collection solid angle subtending an angle of j = 5°, much larger than usual to show better the effect in the graphs. For simplicity, the sizes of the measurement area and the irradiation solid angle were not considered in these illustrating examples.
The upper left plot shows a BRDF distribution without curvature. In this case, the average BRDF coincides with the value of the BRDF at h asp = 0°. For this distribution, the chosen experimental condition regarding the solid angle does not limit the accuracy of the measurement. This would not be true for the distributions shown in the upper right plot and the lower left plot, which are very representative of the BRDFs of matte materials at large incidence angles. 4 In both cases, there is a curvature in the distribution, and a slight deviation is observed between the average BRDF and the actual value of the BRDF at h asp = 0°. This deviation is almost three times larger for the distribution in the upper right plot, which can be explained by a steeper slope of the distribution. For both cases, the deviation is always smaller with a smaller collection angle. Finally, the distribution shown in the lower right plot is not monotonic and has a sharp peak. This example may represent specular reflection. In this case, the deviation is quite large, and the measurement would be meaningless. Again, the devi-°°°°°F ig. 1: Four representative BRDF distributions and the impact of a finite solid angle when measuring these distributions. The solid line represents the BRDF distribution as a function of the aspecular angle. The circle is the average BRDF over a ± 5 degree window ation would always be lower using a smaller collection solid angle.
From this discussion, it can be asserted that the usage of finite apertures limits the BRDF measurements when its distribution has a curvature in the measuring angular range, being more critical when the relative variation of the distribution is higher. The maximum finite sizes of involved intervals that would give an acceptable error are what this article intends to provide. Specular reflection is an important case, but it is not the only one. Any reflectance highly nonlinear with respect to the geometry must be evaluated.
It is an important metrological issue. The International Commission on Illumination (CIE) has recently created a technical committee ''Recommendation on the geometrical parameters for the measurement of the bidirectional reflectance distribution function (BRDF)'' to ''provide geometrical recommendations for the BRDF measurement according to the type of sample under investigation, in order to allow better comparison between the different instruments, to improve the traceability of the measurements, and to help the user to choose the right angular configuration.'' A description of the metrological issues to be considered in the BRDF measurement around the specular peak was published, 5 and it is a very good reference to understand the problem in practice.
This work examines the impact of using finite intervals in the measurement of BRDF, for BRDFs with distributions with different relative variations, and it is proposed as a practical approach to evaluate the limitations of instruments designed for the measurement of the BRDF, and the errors produced by these limitations. Only the impact due to the geometry is taken into account in this article. Out of scope are other effects from using finite areas, such as nonuniformity, which need to be separately considered.

Estimation of the measured BRDF
The value obtained when attempting to measure the BRDF of a sample at a given geometry is denoted by ''measured BRDF'' here. Any instrument designed for this task has irradiation and collection apertures with finite areas, and the radiant flux reflected by the sample and collected by the instrument receiver comes from a finite area: the ''measurement area.'' These areas are graphically depicted in Fig. 2. Four angles are subtended by the irradiation, collection, and measurement areas (j i , j r , j ma,i , and j ma,r ), and all of them affect the measurement. Notice that, whereas j i and j r are constant for any pair of incidence and collection directions (h i0 and h r0 ) (right cone), j ma,i and j ma,r can take slightly different values depending on the illumination and collection directions and the selected point at the perimeter of the measurement area (oblique cone). The values given here for j ma,i and j ma,r are for the geometry with h i0 and h r0 = 0°, for which these values are independent of the selected point at the perimeter of the measurement area.
Only the pair of axial directions of the solid angles in Fig. 2 reproduces exactly the bidirectional geometry at which the BRDF is intended to be measured. The inclusion of the other pairs of directions has an impact on the measurement which depends on the sensitivity of the BRDF to geometry and to the size of the solid angles. This impact is what we want to examine here.
The measured BRDF can be estimated if the scattering distribution function and the finite intervals involved in the measurement are known. The measured BRDF, f r,m , is, by definition, the biconical reflectance factor given in equation (7) divided by p. Then, the measured BRDF (for incidence and collection directions r i0 and r r0 ) can be expressed as: Notice that this expression does not include explicitly the measurement area. The BRDF was defined by Nicodemus et al. 2 as the integration of a more complex scattering distribution function (S) over an elementary area (A i ): where x i and x r are two-dimensional vectors representing the positions of the points at the surface of the incident radiant flux and the exiting radiant flux, respectively. The value of the quantity S (known as bidirectional scattering-surface reflectance distribution function or BSSRDF) depends on r i , x i , r r , and x r , and it is much more complex to measure than the BRDF. For this analysis, it was assumed that the elementary area (A i ) is large enough to be considered as a measurement area and that its reflectance properties are highly uniform. In these conditions, the value of S depends much more on r i and r r than on x i and x r , and the latter dependence can be neglected. Notice that indistinguishable but separated points on the measurement The expression of f r (r i ;r r ) [equation (9)] is inserted in equation (8), obtaining: where the dependence on x i and x r is already neglected. It can be given as a function of the involved subtended angles using McCamy's notation for right circular cones: 2,6 and a h n ð Þ ¼ cos À1 cos j n À cos h n cos h n0 sin h n sin h n0 ð13Þ with n = i or r. R A is the radius of the measurement area, r and c are, respectively, the radial and polar coordinates used to integrate over the measurement area, j i and j r are the angles subtended by the irradiation and collection apertures (see Fig. 2), and h and u, respectively, refer to the polar and azimuth spherical coordinates defining corresponding directions denoted as r. Notice that the angles subtended by the measurement area j ma,i and j ma,r (see Fig. 2) are used in the area integral R 2p 0 dc

Evaluation of the measured BRDF for different conditions
In order to evaluate the impact of finite intervals in the BRDF measurement for non-Lambertian surfaces, the relative deviation of the measured BRDF [estimated by equation (11)] from the actual BRDF (f r ) was examined for different conditions. For this purpose, a family of bidirectional scattering-surface reflectance distribution functions (S) was defined, with a very simple dependence on r i and r r . This family provides a wide range of curvatures and includes angular scattering distribution functions from Lambertian to highly specular with high second derivatives (large curva-tures). These functions are therefore adequate to study the measured BRDF around specular peaks, which is an especially challenging case. 5 The family of curves that represents S can be defined as: where h asp is the aspecular angle, or angular deviation between collection and specular directions [cos À1 (cos h i0 cos h r0 + sin h i0 sin h r0 sin u r0 ), always positive], and b is a free parameter. Notice that the larger the b, the larger is the curvature too, and that  Fig. 3. In the literature 5 , it is claimed that glossy surfaces can be encountered with very narrow BRDF, with typical full width at half maximum (FWHM) inferior to 0.5°, which correspond to b values larger than 3. However, until now, conventional goniospectrophotometers 7-12 do not allow those peaks to be resolved with metrological accuracy. Only special conoscopic systems 5 have been proved adequate to achieve that angular resolution, and they require a separate discussion. We think that the curvatures given by the selected functions are enough to model most of the experimental conditions. In the case of highly glossy samples, there may be a minimum aspecular angle from which the results of this work cannot be applied.
To understand the impact of finite angles, these representative scattering functions were used along with different realistic geometrical conditions, with combinations of the subtended angles j i , j r , j ma,i , and j ma,r (as defined in Fig. 2). To reduce the analysis, a negligible value of j i was used (0.0006°, very small irradiation aperture), which means that the effect produced by the values of j r prevails over the effect produced by j i . There is reciprocity between j i and j r since, according to the Helmholtz reciprocity principle, the result from the analysis is equal to the inverted case at which apertures and distances are interchanged between source and collector, as long as the measurement is for unpolarized reflectance 13 . Values of j r , j ma,i , and j ma,r were selected as 0.11°, 0.23°, 0.49°, or 0.92°, with the additional condition that j ma,i is always equal to j ma,r , equivalent to considering that the selected distances from source and receiver to the measurement area are identical. The selected j values are realistic. If source and detector were located at a distance of 1 m from the sample, the diameter of their areas (assumed circular) would be 3.8 mm for the minimum (0.11°) and 32.1 mm for the maximum (0.92°). It must be reminded here that the solid angles subtended by the measurement area are projected under different incidence and collection directions and that the values of j ma,i and j ma,r reported in this work are those at normal irradiation or collection directions.
The measured BRDF was calculated using equation (11) for the four functions shown in Fig. 3, and for the different geometries, as described in the above paragraph.
The complete evaluation was carried out at two irradiation angles (h i0 = 20°and h i0 = 70°), and for those aspecular angles for which the function is 90% of its maximum value (h asp,90 ). At these aspecular angles, the curvature of the reflectance is 90% of the maximum curvature observed in the sample. The deviations, calculated as the relative difference between the actual BRDF and the estimated measured BRDF, are shown in Fig. 4. Each subplot corresponds to a different surface, with a BRDF calculated at different values of the parameter b, and, consequently, described with a different value of h asp,90 . It is observed that this deviation is always larger for lower irradiation angles, for larger j values, and for larger b values, if the other involved variables are fixed. The awkward behavior of the two lower lines in the lower left plot is precision artifacts from the numerical calculation.
As expected, the larger the b (related to sample's scattering curvature) and j values (subtended angles) are, the larger the deviation is. The relative deviation obtained from the evaluation is similar for the two selected irradiation angles, but slightly lower for the largest one. (Note that different irradiation angles are represented with lines of different thicknesses in Fig. 4.) The reason for this difference is that at larger angles the incidence angles h i on the measurement area are more closely distributed around the irradiation angle h i0 , since the solid angle subtended by this area is projected.
This relative deviation is represented as a function of j r , and every curve corresponds to a different value of j ma,i , which is equal to j ma,r in this case. The curves are constant when the effect of j ma,r dominates that of j r .

General rule to estimate error
The error for any condition can be estimated as the difference between f r and the result from equation (11). However, it may be of interest for the experimentalist to have a general rule or metric to estimate errors. The described methodology was used to obtain the following general rule to calculate the error introduced by finite intervals when measuring BRDFs varying with the incidence and collection directions within a small angular range. Equations (11)-(13) were used to estimate the ''measured'' BRDF under different measuring conditions of surfaces with the scattering distribution functions described by equation (14).
An estimation of a ''measured'' BRDF is defined by b, S/S max [or S (b, h asp )/S(b,0)], h i0 , j i , j r , j ma,i , and j ma,r . Notice that h asp can be determined from b and S/ S max (h asp = -ln(S/S max )/b) and that h r0 can be obtained from the values of h asp and h i0 . In order to assess the impact of the different variables on the error due to finite intervals, two sets of estimations were defined, both of them for b = 1/2 degrees À1 , b = 1/4 degrees À1 , b = 1/8 degrees À1 , and b = 1/16 degrees À1 . The motivation of this strategy was to make the analysis more efficient in terms of computational time, while providing the most valuable information. Note that each estimation is the sum of around 50 millions terms used to calculate the integral in equation (11) with a reduced level of noise and that the calculation of the second set of estimations described below took 80 h.
The first set of estimations assessed the sensitivity of the error to S/S max , which is closely related to the curvature of S. This value was selected from 0.25 to 0.95, with steps of 0.05. The variable h i0 and the relation between the variables j i , j r , j ma,i , and j ma,r were fixed (h i0 = 5°, and j i : j r : j ma,i : j ma,r = 0.25:1:0.5:1), as described in the second column of Table 1. Values of j r = 0.23°, 0.46º, 0.92°, and 1.83°w ere selected. The second set of estimations was defined to assess the variation of the error at different values of h i0 and to confirm the relationships observed after examining the estimations of the first set. In this set (see third column of Table 1), fewer values of S/S max were used (0.3, 0.6, and 0.9). h i0 ran from 15°to 75°, with steps of 15°, and estimations were done for all combinations of the values 0.11°, 0.34°, and 0.92°for j i , j r , j ma,i , and j ma,r .
Insight into the results of the calculations and the steps followed in the analysis is given in the two plots shown in Fig. 5. In the upper plot, the relative errors of the two sets of estimations are represented as a function of S. The data of the second set of estimations are arranged in nine vertical lines, which correspond with the nine S values that were obtained for this set [see attainable values for S = (S/S max ) 9 b]. The data of the first set are scattered uniformly over the entire S range, but are somewhat arranged too on four curves. The estimated errors are related to the value of S because they are related to curvature and curvature is higher for higher S. But it is not the only variable that affects the error, and that is the reason for the four observed curves and the data scattering.
In the lower plot, the estimated errors of the two sets are again represented, but in this case as a function of S 9 F j , where F j is a factor including the finite intervals. It is defined as: with the j values expressed in degrees. There is a clear relation between S 9 F j and the relative error due to finite intervals combined with the curvature growth with S value. By using the 1860 estimations, it was possible to obtain a function expressing the dependence of the 95th percentile of the estimated error on S 9 F j . After trying different types of functions, it was found that an exponential function was the most adequate to relate estimated errors to S 9 F j . The values of S used in the fitting were calculated as the product of the parameters (S/S max ) and b (see Table 1). The equation resulting from the fitting is given below in equation (18). This function is shown in the lower plot of Fig. 5 as a continuous line. This function allows the 95th percentile of the relative error [d 95 (f r )] to be estimated from the j values and the curvature of S, and it can be used as metric for this kind of error. The 95th percentile of error function was obtained using a specific family of scattering distribution functions S. Since it is not obvious that this result can be applied to other functions, it must be clarified here how to apply this function to practical cases. Notice that the function in equation (14) is defined in such a way that the larger its value, the larger its curvature, and that is why d 95 (f r ) is proportional to S (b, h asp ) for this family of functions for a given instrument with constant F j . This is not necessarily the case for other functions, but, within a small range of an angular variable h, an exponential function: as in equation (14), can model any reflectance data with monotonic concave curvature, where the parameter A can be regarded as a scale constant which does not affect the estimated error, and b and h 0 are other free parameters. For monotonic reflectance data with convex curvature, the convex function (h m being the measurement angle) can be used as model, which is symmetric to the previous one with respect to S = S(h m ) and produces the same value of d 95 (f r ).
Once the parameters b and h 0 are obtained, d 95 (f r ) can be estimated for any surface using the function for the 95th percentile of error found above, as: This expression is obtained from the fit shown in the lower plot in Fig. 5. The value of b must be given in degree À1 , and the values of the different js, and any angular variable, in degrees.

Discussion
Equations (16)-(18) developed from the simulations described above can be used to evaluate the value of the 95th percentile of the relative finite-intervalsrelated error for a particular experimental setup. Equation (18) is composed of two important factors: b Â e ÀbÂ hÀh 0 j j and F j which are, respectively, related to the scattering distribution function of the surface and the size of the finite intervals. Any of these two factors can make the error negligible if it is small enough, no matter the value of the other factor. For instance, a small b, which makes the first factor very small, represents a Lambertian surface and produces a negligible error regardless of the size of the finite intervals (the values of j). On the other hand, making j values small enough, the committed error can be made negligible even in the extreme case of very large values of b Â e ÀbÂ hÀh 0 j j , at geometries close to specular peaks.
In general, keeping similar values for j i and j ma,i (and for j r and j ma,r ) is optimal in a goniospectrophotometer, since a lower value of only one of these variables does not reduce significantly the relative error due to finite intervals, and reduces the available light, which increases other uncertainty sources.
Equation (18) supports design or characterization of goniospectrophotometers or multiangle spectrophotometers, since it allows the applicability range of the instrument to be specified for a target uncertainty. For instruments with a fixed value of F j , as it is the case of commercially available portable multiangle spectrophotometers, it may help to estimate for different types of samples the finite intervals uncertainty at their measurement geometries, which usually have different aspecular angles. Aspecular angles would be used as |h À h 0 | value in equation (18), and different values for b would be used to illustrate different types of samples, from very diffuse to very glossy. Instruments allowing the variation of the size of apertures, and consequently F j , can be adapted to the specific type of sample to be measured. For these instruments, the necessary F j is worked out from equation (18), using a value of b which, based on the experience, corresponds with the observed glossiness of the sample, and, as d 95 , the aimed maximum relative uncertainty. Those instruments designed to measure at the specular peak, as glossimeters, introduce an error which depends on their apertures. In such cases, that error cannot be estimated by equation (18), but by comparing the result from equation (11) with the BRDF. To illustrate the error for different types of samples, the BRDF can be modeled with the family of curves in equation (14). In general, equation (11) is recommended to simulate complex conditions for which the derived equation (18) is not suitable.
To estimate experimentally the relative error for a given specimen, the values of b and h 0 must be known for every measurement geometry. If the expected error due to finite intervals is relatively low, the ''measured'' BRDF itself can be used to estimate b and h 0 . However, if a high error is expected in the measurement, as for low aspecular angles in highly glossy specimens, some deconvolution algorithm needs to be applied to estimate b and h 0 from the deconvolved BRDF. 14,15 In these cases, equation (18) is recommended to be used to characterize the limitations of the instrument to measure BRDFs with high angular dependence, and to redesign it if required.
It must be noticed that there are limitations in the applicability of equation (18), since it was derived for a range of curvature and j values. The maximum curvature used in the derivation was slightly larger than 0.1, so that it can be recommended to use equation (18) only when it is expected that: curvature : with h 0 and b obtained from equations (16) or (17).
Regarding j values, the use of equation (18) is restricted to values below 1º.

Conclusions
The impact of using finite intervals (for the measurement area, irradiation and collection apertures, and solid angles) has been examined for the measurement of BRDF. A theoretical approach has been proposed to evaluate the limitations of experimental systems for the measurement of the BRDF and to account for the relative error exclusively due to finite intervals. This error, calculated as the relative difference between a biconical integration of the BRDF and the actual BRDF, is composed of two main factors, which are related to the scattering distribution function of the surface, and the geometries of the incident and collected radiant fluxes. Any of these factors can reduce to negligible the total error if it is set to a value small enough. The general rule for experimentalists provided here should be useful for the design of instruments for measuring BRDFs very dependent on incidence and collection geometries, such as those intended to measure gloss or iridescence.