The Impact and Correction of Sensitive Environmental Factors on Spectral Reflectance Measured In Situ

: The spectral reflectance measured in situ is often regarded as the “truth” of objects, which plays an important role in Earth observation applications. However, in situ measurements are influenced by several factors such as atmospheric conditions, illumination and view geometry (I&VG), cloud coverage, and adjacency effects. In order to avoid the influence of these factors, in situ measurements are usually carried out under sunny days and close to noon. However, the impact of I&VG is still present in most cases. At present, people still know little about the influence mechanism of I&VG. Moreover, correcting the impact of I&VG is also a problem that needs to be urgently solved in reflectance spectroscopy. In this work, experiments are carried out using the multi-directional hyperspectral remote sensing simulation facility (MHSRS 2 F), which allows adjustment and control of the I&VG parameters. This paper proposes an uncertainty evaluation model for I&VG and quantifies the uncertainty caused by different I&VG parameters. Then, the sensitivity of reflectance to I&VG at different wavelengths is explored based on uncertainty models. Finally, a correction model for reflectance under different I&VG conditions is proposed. The results reveal that the uncertainty and sensitivity caused by observation height are relatively high, regardless of the surface heterogeneity. It directly affects the size of the field of view and the physicochemical characteristics of the object. For objects that approximate the Lambertian surface, more attention should be paid to the selection and variation of solar and view zenith angles and view azimuth angles. For objects with surface heterogeneity, the selection and variation of solar azimuth angle, view azimuth angle, and solar zenith angle are more crucial. The correction model proposed in this paper has a 41.25% correction effect on different view zenith angles, but the correction effect on other environmental factors is not significant.


Introduction
Field spectroscopy has become an important technology for characterizing the reflectance features of natural surfaces and has received increasing attention in the last two decades [1].The advantage of field spectroscopy is that the portable spectrometers can remain fixed over the object for a longer time, and the shorter path length between the instrument and object reduces the influence of the atmosphere [2].Therefore, in situ measurements are often used as the truth of objects and plays an important role in the applications of Earth observation and environmental monitoring systems [3].The main applications include the following: the vicarious calibration of remote sensors [4]; the correction of atmospheric effects on airborne or spaceborne images [5]; the validation of the quantitative remote sensing data products [6]; the identification and classification of the targets from images [7]; the establishment of models of biophysical attributes, biological processes, and remote sensing feature attributes [8]; the prediction of the best spectral bands, geometric configurations, and the optimum time for a specific task [9]; and the supplementation of the spectral library.
In situ measurements are easily influenced by environmental factors, such as atmospheric conditions, illumination intensity, illumination and view geometry (I&VG), sky cloud coverage, and the adjacency effects, which lead to the phenomenon that the same objects have different spectra [10].In order to avoid the influence of environmental factors, in situ measurements are usually made under sunny days and close to noon.However, the influence of I&VG is still present in most cases.In fact, the prerequisite for achieving ideal results in various applications based on in situ measurement is synchronous observation on the ground and airborne or spaceborne observation.This means that in situ measurements need to ensure that environmental conditions are consistent with the airborne or spaceborne data, especially the angle conditions, field of view, etc. [11,12].However, practice has proven that synchronous observation is difficult to realize.At the same time, it is difficult to maintain the same observation height for in situ measurements.Due to the heterogeneity of most land surfaces, the change in field of view caused by different observation heights inevitably leads to fluctuations in spectral reflectance [13].Moreover, many models are based on the Lambertian surface hypothesis [14].However, most natural features are non-Lambertian, and the changes in spectral reflectance measured in situ caused by illumination conditions, viewing geometry, and field of view will inevitably affect the validation accuracy of a model based on the Lambertian hypothesis.Therefore, this article mainly discusses the influence of I&VG on in situ measurements of objects with different degrees of surface heterogeneity.
At present, the uncertainty of spectral reflectance measured in situ caused by I&VG has been widely studied.There are two common approaches to studying this uncertainty: in the laboratory and in the natural environment.Instruments such as LAGOS, CLabSpeG, and FIGIGO are designed to obtain the bidirectional reflection factor (BRF) of different objects in the laboratory [15][16][17].Meanwhile, a lot of research on the influence of I&VG is conducted outdoors.The BRF varies with solar position, viewing plane and angles, and sensor FOV and distance combinations [18].The measurement uncertainty value of spectral reflectance caused by solar radiation can reach 5% to 10% [10].Gu et al. revealed that around solar noon, the solar zenith angle only changes slightly within 0.2%.When the view zenith angle changes from 0° to 30°, the variation in BRDF can reach 10% [19].Anderson et al. investigated the effect of solar zenith angle on reflectance ranges from 16% to 32% [20].
Despite these advances, many problems have not yet been fully investigated.A challenging problem in uncertainty evaluation based on laboratories is that it is difficult to simulate real field environmental conditions.The main problem of research based on in situ measurements is that I&VG are uncontrollable.This has created an issue where the uncertainty results of measurement data are inconsistent under different I&VG conditions.In addition, it is unfeasible to collect spectral reflectance under all I&VG conditions.This results in the fact that the comprehensive quantitative evaluation of uncertainty results caused by I&VG is rarely performed.Therefore, current applications based on in situ measurements still use a single spectrum as the reference truth, which limits the reliability of subsequent application results [3].
Another problem with these two research methods is that the uncertainty caused by I&VG is formed by the coupling of multiple I&VG parameters.Because the decoupling of multiple I&VG parameters is still a significant challenge, the effect mechanism of any single I&VG parameter is still poorly understood [21].Based on the understanding of the influence of I&VG, correcting the influence of different I&VG conditions on spectral reflectance and normalizing different I&VG conditions are other problems that need to be urgently solved.
To address the above limitations, this study takes rocks as the target object, which can maintain their physicochemical properties for a long time.In this work, five I&VG parameters, namely, solar and view azimuth angles, solar and view zenith angles, and observation height, were experimented using the MHSRS 2 F [22].The MHSRS 2 F can simulate the outdoor environment more realistically and allows adjustment and control of the I&VG parameters.Based on the experimental results, an uncertainty evaluation model for I&VG was constructed.In this work, the uncertainty caused by different I&VG parameters was quantified.Then, the sensitivity of spectral reflectance to I&VG at different wavelengths was explored based on the uncertainty models.On the basis of understanding the influence of I&VG, a correction model for reflectance under different I&VG conditions is proposed.

Experimental Setup
The experimental instruments included the Multi-Directional Hyperspectral Remote Sensing Simulation Facility (MHSRS 2 F) and the Analytical Spectroscopy DeviceFieldSpec Pro (ASD) spectroradiometer.The MHSRS 2 F (Figure 1) consists of three parts: a solar simulator, a skylight simulator, and a multi-angle observation simulator [22].This is specifically designed to obtain the same value as the reflectance factor (RF) acquired in situ by simulating the illumination conditions of sunlight and skylight.After statistical analysis, the spectral correlation coefficients of the simulated solar spectrum and skylight spectrum with the standard spectrum can reach 0.928 and 0.912, respectively.Appendices A and B provide field environmental simulation experiments and instrument stability tests based on the MHSRS 2 F. The above experiments prove that the MHSRS 2 F can stably and accurately simulate the field environmental conditions in the laboratory to collect the hemispherical conical reflectance factor (HCRF).In addition, the solar and view azimuth angles can be adjusted from 0° to 360°.Solar zenith angles range from 15° to 60°, while view zenith angles span from 0° to 90°.All the spectral reflectance data described in this paper were obtained by using the ASD spectroradiometer with a 5° fore optic lens, which covers the 350-2500 nm spectral range.In this experiment, the Spectralon TM White Diffuse Reflectance Standard (99%) was used as a reference panel to acquire the RF.

Samples, Sampling Sites, and In Situ Measurements
Ten rock samples collected from the Gobi Desert Mining Area in Hami, Western China, were selected for the experiment (Figure 2).As shown in Figure 3, the collected rock samples encompassed a variety of types, including common sulfides, halides, oxides, carbonates, and silicates.These samples exhibited diverse properties such as different colors, roughness, crystal structures, grain sizes, and gloss levels.The surface area of rock samples was sufficient to fill the field of view of the instrument.Since this research focuses on the impact of environmental conditions, the original-state rocks were used in this experiment.The spectral reflectance of corresponding rocks was measured at the place where the rock samples were collected, and then the environmental conditions were recorded (Table 1).At the position of each rock sample, 7 measurements were made on the panel and the rock surface, and the measurement sequence was panel, 5 times rock surface, and panel.Each individual measurement consisted of 5 ASD scans.Therefore, a total of 35 in situ measured spectral reflectance were generated for each rock sample.

Simulation Experiments of Field Environmental Factors
The main purpose of this experiment is to explore the influence of I&VG within the environmental factors on measured reflectance spectra.For near-surface observations, the rock surface is a three-dimensional rather than a smooth plane, and the bidirectional reflection characteristics of rocks formed by different observation and illumination azimuth angles are different.Therefore, this experiment selected the use of absolute azimuth angles instead of relative azimuth angles.In addition, in the process of in situ measurements, the variation in the instrument field of view is usually caused by differences in observation heights.Therefore, this experiment selected five parameters of I&VG, namely solar zenith angle, solar azimuth angle, view zenith angle, view azimuth angle, and observation height, as the research objects.
Based on the latitude of the Mining Area, the annual variation range of solar zenith angle in this region is around 19.2-90°.And the in situ measurement is usually carried out at noon.Therefore, the variation range of solar zenith angle is from 20° to 60°.During in situ measurement, different operators may use different measurement methods, which may lead to differences in observation height.When an adult stands and holds the spectroradiometer probe, the distance between the probe and the target object is usually around 100 cm.Simultaneously, when the observation height is less than 10 cm, it is easy for the shadow caused by the probe to appear on the surface of the measured object.Therefore, the variation range considered in this simulation experiment of observation height ranges from 10 cm to 100 cm.The HCRF of rock samples with different environments is collected by adjusting the sunlight simulator, the skylight simulator, and the multi-angle observation simulator in the MHSRS 2 F. The specific experimental protocol is shown in Table 2 and Figure 4.

Surface Heterogeneity Experiment
In order to explore the influence of different I&VG parameters on the measured reflectance spectra with different degrees of surface heterogeneity, the goal of this experiment is to quantify the surface heterogeneity in 10 rock samples.Under the same I&VG conditions, 45 different observation positions of rock samples were selected along the X and Y axes at 2 cm intervals, by using an electronic control translation stage.Since the observation height of the instrument in this experiment is 100 cm and the fore optic lens is 5°, the surface area of the rock sample evaluated in this experiment is around 20 cm × 12 cm (Figure 5).The experimental results show that the root-mean-square error (RMSE) caused by the heterogeneity of 10 rock surfaces is around 0.42-1.62%.

Field Data Uncertainty Propagation
For this research based on the MHSRS 2 F, the simulated in situ measurement values can be decomposed into four independent parts: the part related to geometric variables, the part related to sample surface heterogeneity, the part related to ASD spectroradiometer, and the part related to radiation intensity.Therefore, the simulated in situ measurement model can be expressed as the sum of four independent functions by using an additive model: where ) , , , , ( According to the measurement model (Equation ( 1)), a complete evaluation model of spectral uncertainty for in situ measurements caused by I&VG under the MHSRS 2 F is proposed.Since the measurement model is an implicit model, quantifying the sensitivity coefficient is difficult.Therefore, this section assumes that the uncertainty of each term has a consistent influence on the sensitivity of the measurement uncertainty, and the sensitivity coefficient of each term is ignored in the uncertainty calculation.The total uncertainty caused by I&VG total  can be given by the following expression: where total  refers to the mean of the standard deviation (SD) of the reflectance spectra collected under different I&VG parameters that is simulated by the MHSRS 2 F; geometry Therefore, the uncertainty of the objects caused by a single I&VG parameter geometry  can be expressed as The uncertainty of each band for the reference panel caused by five parameters of I&VG is calculated based on Formula (3) (Figure 6).Since these five parameters are all geometric variables, their changes in different bands are consistent.In addition, according to the stability experiment in Appendix B, the SD caused by the stability of solar radiation intensity, scattered radiation intensity, and instruments during the long-term experiment ranges from 0.028% to 0.059%.If their influence is ignored, the uncertainty of each band caused by any parameter of I&VG should be the same.However, Figure 6 shows that the uncertainty greatly fluctuates in the range of 350 nm-400 nm, NIR and SWIR regions.This is because the noise of the two sensors in the SWIR of the ASD spectroradiometer is greater than that of the sensors in the VNIR, and using the average noise of all bands as ASD  , the uncertainty caused by noise in SWIR regions has not been eliminated [23].The difference between the average uncertainty of the reference panel in all wavebands and in the range of 400 nm-900 nm can be used to represent the portion of the uncertainty caused by the noise of the ASD spectroradiometer in the SWIR region ASD_SWIR  that exceeds the average noise value of all bands ASD  .Therefore, this paper uses the average uncertainty of the reference panel in the 400 nm-900 nm range to characterize the uncertainty caused by I&VG, which can be considered as noise having no additional impact in the SWIR range.Since the reference panel is considered as a homogeneous surface, the uncertainty of the reference panel caused by a single I&VG anel geometry_p  can be written as where n is the number of I&VG parameters; i  refers to the uncertainty caused by the I&VG parameters i ; i x refers to the average value of repeated measurements under is the correlation coefficient between two I&VG parameters; and k is the coverage factor, which is taken as 2 in this work.

Sobol's Global Sensitivity Analysis
Since the in situ measurement model is an implicit function, it is difficult to quantify the sensitivity coefficient.Therefore, the uncertainty calculation in Section 3.1 ignores the sensitivity coefficients of each term.The innovation of this paper is the use of the spectral reflectance collected by the MHSRS 2 F to replace the results of the in situ measurement model.According to the variation range of the five I&VG parameters in Table 2, the combination of five I&VG parameters is randomly selected using the random sampling method.The combination of five I&VG parameters obtained through random sampling is set in the MHSRS 2 F, and then the HCRF of the measured object is collected.A total of 96 experiments were conducted in this work, and a total of 96 sets of experimental data were obtained for subsequent sensitivity analysis.
Sensitivity analysis (SA) of the model output can be used to determine how the model (numerical value or otherwise) depends on its input factors.I.M. proposed Sobol's Global Sensitivity Analysis method [24].Global Sensitivity Analysis can not only consider the influence of multiple parameter changes on the output of the model at the same time, but also quantify the influence of the interaction between parameters on the output of the model [25].The advantage of this method is that the calculation form is simple and there is no requirement for the monotonicity, linearity, and distribution characteristics of the model.Sobol's method is to decompose the total variance of the output result of the model into the first-order variance of the independent parameters and the higher-order variance of the parameter combination.Then, the ratio of the first-order variance and the higherorder variance to the total variance is calculated to determine the influence of the parameters or parameter interaction on the output of the model.Therefore, the first-order sensitivity index for parameter i x and the total sensitivity index for parameter i x can be de- fined as  x .

Field Data Correction Model
Since reflectance cannot be directly measured, the infinitesimal elements of the solid angle do not include measurable amounts of radiant flux [26].Therefore, the reflectance factor (RF) is proposed as the equivalent coefficient used in practice.It is the ratio of the radiant flux reflected by a surface to that reflected by an ideal diffuse Lambertian standard surface, irradiated under the same illumination and observation conditions [27].Since a perfectly reflecting panel does not exist in practice, the spectral reflectance of the reference panel needs to be corrected [14]: where ) , , , , ( refer to the radiance of the target objects and reference panel under the same illumination and observation conditions, respectively; and K refers to the panel correction factor.
With the change in I&VG, the radiance of the object and the reference panel will also change at the same time.Because the variation in object radiance with I&VG is related to its physicochemical properties, it is difficult to find the universal law of variation of object radiance with I&VG.However, because the reference panel used in spectrum collection is unchanged, the law of variation of the reference panel radiance with the I&VG can be determined through multiple measurements based on the MHSRS 2 F.
Therefore, a correction model of spectral reflectance under different I&VG parameters is proposed in this paper.The core concept is to correct the influence of different I&VG parameters on the reference panel, to eliminate the partial of influence of I&VG on the spectral reflectance of the object:  As for the correction work of this article, firstly, the spectral reflectance of two groups of rock samples under different I&VG conditions were collected in the field as experimental data before and after correction, and the corresponding I&VG conditions were recorded.Then, the same I&VG conditions were set in the MHSRS 2 F and the radiance of the reference panel was collected.And based on the method of locally weighted regression, a correction factor ( C ) under the corresponding I&VG conditions is obtained.Finally, mul- tiplying the RF of the object before correction by the corresponding correction factor ( C ),the RF after correction can be obtained.RMSE is used as an indicator to evaluate correction ability in this work.

Reference Panel
The reference panel used is the Spectralon TM White Diffuse Reflectance Standard (99%), which is considered as an ideal diffuse reflectance reference surface.Therefore, this article uses the variation in the spectral reflectance for the reference panel under different I&VG conditions to characterize the influence of different I&VG conditions on the spectral reflectance.The uncertainty in the reference panel caused by a single I&VG parameter and the coupling of multiple parameters are calculated and shown in Table 3.A box-plot is used to compare the magnitude of uncertainty (Figure 7).The comparison results in Figure 7 reveal that the uncertainty of reflectance caused by the observation height is the highest, because it involves the area of the measured object detected using the spectroradiometer.Simultaneously, the uncertainty caused by the solar zenith angle is the second highest, because changes in the solar zenith angle directly affect the magnitude of solar irradiance.In addition, the uncertainty caused by the solar azimuth angle is the lowest.Another noteworthy phenomenon is that the uncertainty caused by the azimuth angle is significantly lower than that of zenith angle.The main reason may be that the surface of the reference panel is a Lambertian and the different view and solar azimuth angles have relatively little influence on spectral reflectance.The above findings indicate that attention should be paid to selecting and unifying the most suitable observation height, viewing direction, and time of spectrum acquisition during field operations.

Rock Samples
In this section, the variation in reflectance of rock samples under different I&VG conditions is used to explore the uncertainty caused by I&VG at different degrees of surface heterogeneity.The uncertainty range caused by a single I&VG parameter and the coupling of multiple parameters for 10 rock samples are presented in Table 4.A dumbbell plot is used to represent the distribution of the uncertainty range (Figure 8).The dumbbell plot reflects the consistency of uncertainty caused by I&VG on 10 rock samples.When the uncertainty range is small, it indicates that the uncertainties caused by I&VG on 10 rock samples are either very large or very small.On the contrary, a large uncertainty range indicates that I&VG has higher uncertainty for some rock samples, while it has lower uncertainty for others.Table 4. Summary of uncertainty results of the spectral reflectance for rock sample, which are caused by five I&VG parameters and the coupling of five parameters (k = 2).The uncertainty range refers to the calculations given in Equations ( 5) and (6).

Source of Uncertainty
Uncertainty Range (%)  According to Figure 8, it is apparent that the uncertainty range and value of rock samples caused by the solar zenith angle is the smallest.This indicates that for most rock samples, the uncertainty caused by the solar zenith angle is the smallest.Another finding is that the uncertainty range caused by observation height is relatively small, but the uncertainty value is relatively high.This shows that for most rock samples, the uncertainty caused by the observation height is relatively high.In addition, the uncertainty range caused by the solar and view azimuth angles is the largest.This indicates that for some rock samples, the uncertainty caused by azimuth angle is high, while for other rock samples, the uncertainty caused by azimuth angle is very low.And the uncertainty caused by azimuth angle depends more on the surface heterogeneity of the rock samples.From the above analysis, it can be seen that when there is heterogeneity on the surface, it is necessary to reasonably select and determine the observation height, the direction of sunlight incidence, and the viewing direction.

Comparison of Reference Panel and Rock Samples
According to the surface heterogeneity simulation experiment (see Section 2.4), two rock samples with the highest and lowest surface heterogeneity are selected in this section for comparison with the reference panel (DTS21-1 and DTS21-6 seen in Figure 3).The comparison results of the uncertainty caused by single I&VG parameters are shown in Figure 9. Figure 9 shows that the influence of coupling I&VG with the surface heterogeneity of rock samples is much higher than that caused by I&VG alone.Regardless of the surface heterogeneity, the uncertainty caused by the observation height is relatively high.
Another finding is that the uncertainty of the reference panel caused by the solar zenith angle is the second highest.However, when there is heterogeneity on the surface, the uncertainty caused by the solar zenith angle becomes the lowest.This indicates that the uncertainty caused by the coupling of the solar zenith angle with the surface heterogeneity of rocks is smaller than the influence of other I&VG parameters coupling with the surface heterogeneity.The influence of heterogeneity on the solar zenith angle is less than that on other parameters.In addition, with the increase in the heterogeneity of the rock surface, the uncertainty caused by the solar and view azimuth angles is significantly greater than that caused by the zenith angle.Taken together, the uncertainty caused by different I&VG parameters is closely related to the surface heterogeneity of the object.Regardless of the surface heterogeneity, the selection of the observation height is very important for spectral reflectance.When the surface of the object is approximately Lambertian, the influence of zenith angle is greater than that of the azimuth angle.This means that more attention needs to be paid to the selection of operation time, as well as the unified regulations of the view zenith angle.When the surface heterogeneity increases to a certain extent, the influence of the azimuth angle is greater than that of the zenith angle.It is necessary to pay more attention to the selection and unified regulations of the incident direction of sunlight and the viewing direction.

First-Order Sensitivity Analysis
In the previous chapter, the uncertainty calculation caused by the coupling of multiple I&VG parameters assumes that the sensitivity of the uncertainty caused by each I&VG parameter on the spectral reflectance uncertainty is consistent.The purpose of this section is to explore the sensitivity of reflectance to each I&VG parameter.The first-order sensitivity coefficient represents the sensitivity of the reflectance data to a single I&VG parameter.Sensitivity analysis still selects the reference panel and the two rock samples with the minimum and maximum heterogeneity as the objects.Since the sensitivity of different bands to I&VG parameters is inconsistent, this work separately counted the most sensitive I&VG parameters for 2150 bands (350 nm-2500 nm) and used a histogram to characterize the number of bands most sensitive to these five I&VG parameters (Figure 10).The comparison results characterize the sensitivity of spectral reflectance to I&VG and the impact of heterogeneity on the first-order sensitivity.In Figure 10, it can be seen that for the reference panel, the differences between the columns representing the five I&VG parameters are relatively small, which indicates that the sensitivity of each band to I&VG parameters is different.Comparing the histograms of the three measured objects, an interesting phenomenon is that when there is heterogeneity on the surface, the sensitivity of different bands to I&VG begins to show some consistency.The number of bands that are most sensitive to the observation height increases significantly, from minimum to maximum.However, the number of bands that are most sensitive to the view azimuth angle initially changes from maximum to minimum.This reveals that when I&VG parameters are coupled to surface heterogeneity, the reflectance of most bands is more sensitive to observation height than other parameters, but the sensitivity to view azimuth angle is much less than those of the others.
The comparison of the two rock samples in Figure 10 shows that for non-Lambertian objects, the reflectance data are more sensitive to zenith angle than to azimuth angle.When the surface heterogeneity is relatively high, the reflectance is more sensitive to two angles in the direction of sunlight illumination than to two angles in the viewing direction.Another finding is that relatively low surface heterogeneity is more sensitive to the view zenith angle, while relatively high surface heterogeneity is more sensitive to the solar zenith angle.

Total Sensitivity Analysis
Compared to the first-order sensitivity coefficient, the interaction between I&VG parameters is also considered in the total sensitivity coefficient.In this section, the number of bands with the highest total sensitivity coefficients is calculated for the five I&VG parameters for each band of the reference panel and two rock samples (Figure 11).The purpose is to investigate the sensitivity of reflectance at different bands to I&VG after considering the interaction.Figures 10 and 11 show that, after considering the interaction between different I&VG parameters, most of the bands of two rock samples have changed from being more sensitive to the observation height and zenith angle to being more sensitive to the solar azimuth angle.Since the total sensitivity of each band for the reference panel still has significant differences, it is difficult to analyze the influence of the interactions of I&VG parameters on the sensitivity.Therefore, this section uses a bump chart to characterize the sequence changes in the number of bands that are sensitive to I&VG parameters for the first-order and the total sensitivity coefficients of the reference panel (Figure 12).
Figure 12 shows that the reference panel is the least sensitive to observation height and solar azimuth angle, regardless of whether the interaction between I&VG parameters is considered.This reveals that when the surface of the measured object is close to Lambertian, the slight changes in the observation height and solar azimuth angle have little influence on the reflectance data.Another novel finding is that, for the total sensitivity coefficient, the number of bands sensitive to the zenith angle of the reference panel is greater than that sensitive to the azimuth angle.These results are consistent with the previous uncertainty results of the reference panel caused by I&VG: the uncertainty caused by the zenith angle is greater than by the azimuth angle.The implication of these findings is that the impact of solar and view zenith angles on reflectance is greater than that of solar and view azimuth angles.

Comparison of First-Order Sensitivity and Uncertainty
The difference between sensitivity and uncertainty is that uncertainty describes the degree of dispersion for the reflectivity of the object within a certain range of I&VG conditions, while sensitivity represents the contribution of various I&VG parameters to the uncertainty of the reflectivity of object.Because the first-order sensitivity coefficient represents the sensitivity to a single I&VG parameter, this section compares the first-order sensitivity coefficient of the reference panel and rock samples with the uncertainty caused by a single I&VG parameter.This section uses a bump chart to characterize the difference between the first-order sensitivity coefficient and the uncertainty caused by a single parameter (Figure 13).
An important finding in Figure 13a is that the uncertainty of the reference panel caused by the observation height is highest, but the reflectance is the least sensitive to the observation height.This discovery reveals that there is a relatively large difference in reflectance collected at observation heights of 20 cm and 100 cm, but the differences between observation heights of several centimeters are small.This means that when objects are close to Lambertian, it is necessary to uniformly specify the observation height, and the influence of the error of several centimeters on the reflectance can be negligible.A further novel finding is that the uncertainty and sensitivity of the solar zenith angle are relatively large.This experimental result definitely supports the need to pay attention to the selection of operation time and changes in the solar zenith angle during in situ measurement.In addition, the uncertainty caused by the view azimuth angle is relatively low, but the sensitivity is highest.This result indicates that when the surface of the object is close to Lambertian, it is necessary to pay attention to changes in the view azimuth angle during the operation process, and slight changes may lead to changes in reflectance.Comparing these two bump charts in Figure 13b,c, a consistent phenomenon is that the uncertainty and sensitivity caused by the observation height are high regardless of heterogeneity.This reveals that it is necessary to reasonably set a unified observation height and pay attention to slight changes during the spectral acquisition process.It is also worth noting that the uncertainty caused by the solar zenith angle is relatively low, but the sensitivity is relatively high.This indicates that when collecting spectra in situ, it is necessary to pay attention to the instantaneous changes in the solar zenith angle.In addition, the uncertainty caused by the view azimuth angle is relatively high but the sensitivity is the lowest, which means that it is important to carefully consider the viewing direction, but slight changes in the view azimuth angle have little effect on reflectance.

Correction of Field Data under Different I&VG Conditions
This section conducts reflectance correction work under different I&VG conditions for the 10 rock samples.Figure 14a,b give the correction results for two rock samples from view zenith angles of 40° to 20°.The results of the mean value of the RMSE for 10 rock samples at various wavelengths show that this model has a correction ability of 41.25% for reflectance under different view zenith angles.Moreover, the applicability of this correction model is also tested on soil samples.Figure 14c shows the correction results of a soil sample, and the correction ability can reach 33.9% under different view zenith angles.These results show that this correction model has a good correction effect for reflectance under different view zenith angles.Compared to the view zenith angle, the correction result of this model for other I&VG parameters is not ideal.The main reason can be attributed to a significant negative linear correlation between the radiance data of the reference panel under different view zenith angles (Figure 15).The average determination coefficient of each band obtained through linear fitting is about 0.9516.However, for the other four I&VG parameters, the goodness of fit of the reference panel radiance under different I&VG parameters is poor or there is no significant trend.Moreover, it can be seen from the red fitted curve in Figure 16c that there is no significant trend in the radiance of different bands under different solar zenith angles.Previous studies have shown that there is only a small correlation between the SD of the reference panel radiance and the solar zenith angle [13], hence the difficulty in correcting the impact of solar zenith angles through the correction factors obtained through curve fitting.
Figure 16d shows the radiance distribution of the reference panel at different observation heights.From the fitting results, it is obviously better than the goodness of fit of azimuth angle.However, the correction model proposed in this paper has little effect on correcting observation heights.One possible explanation for this is that there is a very small difference in the radiance of the reference panel at different observation heights.Therefore, the correction factor is very close to 1, and the correction results for different observation heights through Equation ( 11) are very small and difficult to observe.

Conclusions
In this study, the uncertainty caused by a single I&VG parameter and multiple parameters is quantified (Tables 3 and 4).The quantification of these uncertainty ranges can provide reference for optimization and improvement of the application of a single truth value based on the object.Uncertainty describes the dispersion of the reflectivity within a certain range of I&VG conditions, while sensitivity represents the influence of slight changes in various I&VG parameters.Therefore, the results of uncertainty and sensitivity analysis can provide the following suggestions and guidance for the optimization design of experimental plans for in situ measurement: 1.The uncertainty and sensitivity caused by different I&VG conditions are closely related to the surface heterogeneity of the object.The greater the surface heterogeneity, the greater the uncertainty.2. Regardless of surface heterogeneity, the uncertainty and sensitivity caused by observation height are greater than those caused by I&VG.Because the observation height directly affects the size of the field of view and the physicochemical characteristics of the measured object within the field of view, the selection of observation height and the avoidance of changes during the experimental process are crucial.And the scale effect can be a noteworthy issue in the future.3.For approximate Lambertian objects, the results of uncertainty and sensitivity are relatively consistent.The uncertainty and sensitivity caused by the solar and view zenith angles are relatively high.This indicates that the selection and variation of the zenith angle are crucial.4. When there is surface heterogeneity on the measured object, the uncertainty caused by the solar and view azimuth angles is relatively high, but it is more sensitive to the solar azimuth angle and solar zenith angle.This indicates that the selection of the solar azimuth angle and the avoidance of changes during the experimental process are crucial.Additionally, more attention should be paid to the changes in the solar zenith angle and the selection of view azimuth angle.
The correction method for reflectance data under different I&VG conditions proposed in this study has been successfully applied to correct the view zenith angle and has achieved good results, with a correction ability of 41.25%.However, the correction effect for the other four I&VG parameters is not ideal.Therefore, further exploration of correction models for these I&VG parameters will be required in the future.
Notwithstanding the relatively limited rock samples, this work provides valuable insights into the influence of different I&VG parameters on the spectral reflectance.One of the limitations of this study is the unified use of surface heterogeneity to represent differences in physicochemical properties between objects.The experimental results show that the uncertainty caused by the coupling of surface heterogeneity with I&VG parameters is much greater than that caused by I&VG alone.Therefore, on the basis of understanding the impact of I&VG on spectral reflectance, exploring the impact of coupling different physicochemical properties with I&VG will be a fruitful area for further research.This is helpful for us to better understand the influence mechanism of I&VG for the spectral reflectance of objects.In addition, spectral characteristic parameters are also selected as an index for evaluating spectral consistency in this experiment.And a SpectralonTM Color Diffuse Reflectance Standards panel with absorption characteristics is selected as the measurement object (Figure A1).The HCRF of the reference panel is collected under the same environmental conditions in the field and the MHSRS 2 F. Three commonly used spectral characteristic parameters, absorption position, absorption depth, and spectral absorption index (SAI), are used to measure the simulation degree of the MHSRS 2 F to the field environmental conditions.The calculation results indicate that the maximum difference in absorption peak is 2 bands, and the maximum differences in absorption depth and SAI are 0.0159 and 0.01486, respectively (see in Table A2).The above experimental results have verified that the MHSRS 2 F has a good ability to simulate field environmental conditions.

Appendix B. Instrument and Illumination Stability Experiment
The main focus of this experiment is to verify that the MHSRS 2 F can maintain the stability of illumination conditions for a long time.According to the experimental scheme in Table A3, the spectra of rock samples are collected every 2 min for 1 h.The environmental conditions of the MHSRS 2 F are kept unchanged during the experiment.This study uses the root-mean-square error (RMSE) to characterize the reproducibility of the MHSRS 2 F. It turns out that the reproducibility of the MHSRS 2 F is between 0.028% and 0.059%.Milton et al. found that under temperate-latitude clear-sky conditions, the reproducibility of continuous-reflectance measurement using the single spectrometer approach is around 7%, and the reproducibility of the dual spectrometer approach is about 3% [28].Comparison shows that the reproducibility of the single spectrometer approach based on the MHSRS 2 F is significantly better than the reproducibility of the single spectrometer and dual spectrometer approach in the field.
Table A3.Three groups of stability experiments were conducted, and the environmental condition parameters set for the experiment were as follows.

Experiment
No.

Sample
No.

Figure 1 .
Figure 1.Actual picture of the components and structure of MHSRS 2 F.

Figure 2 .
Figure 2. The 10 rock samples were collected from 6 sampling sites in the Gobi Desert Mining Area in the west of Hami City, Xinjiang Uygur Autonomous Region, China.The 6 sampling sites are represented by colored five-pointed stars.

Figure 3 .
Figure 3.The 10 original-state rock samples used in this experiment.

Figure 4 .
Figure 4. Adjustment trajectories of different I&VG parameters in the MHSRS 2 F.

Figure 5 .
Figure 5.The observation range for heterogeneity evaluation of the sample is shown in the red box (20 cm × 12 cm).The red dots indicate the 45 observation positions at intervals of 2 cm.
is a function related to five geometric parameters of I&VG within environmental factors; H refers to the observation height; i  and v  are the solar and view zenith angles, respectively; i  and v  are the solar and view azimuth angles, respectively; ity heterogene f is a function related to the surface heterogeneity of samples; ASD f is a function related to the ASD spectroradiometer noise; and on illuminati f is related to the stability of solar radiation intensity and scattered radiation intensity simulated by the MHSRS 2 F.
by the I&VG and surface heterogeneity of the object, respectively; and ASD  and on illuminati  are uncertainties caused by the ASD spectroradiometer noise and the stability of solar radiation intensity and scattered radiation intensity simulated by the MHSRS 2 F, respectively.

Figure 6 .
Figure 6.Scatter charts with different colors and shapes indicate the uncertainty of the reference panel at different wavelengths with different I&VG parameters.Solid lines of different colors represent the average uncertainty of different wavelengths.
mean of the SD of the reference panel in the range of 400 nm-900 nm, and ASD  refers to the uncertainty caused by the ASD spectroradiometer noise in all bands.The uncertainty of the rock samples caused by a single I&VG ample geometry_s mean of the SD of the reference panel in all bands, and ASD_SWIR  refers to the portion of the uncertainty caused by ASD spectroradiometer noise in the SWIR region that exceeds the average noise value of all bands ASD  .For uncertainty caused by the coupling of multiple I&VG parameters, geometry  is given by the following expression, which is equal to geometry  in Equation (2): where

S
refers to the first-order sensitivity coefficient, which measures the main effects of i x on the output result total sensitivity coefficient, which measures the main effects of i x and the effects of the interaction between i x and other parameters; andi x V ~is the sum of the variance caused by all parameters except i

__
is the RF of the target object under a certain I&VG condition before correction;after T R _ represents the RF of the target object after correction to the desired I&VG conditions; refer to the radiance of the target object and panel before correction, respectively; and as a correction factor ( C ) related to the I&VG conditions.

Figure 7 .
Figure 7.The box-plot shows the uncertainty distribution of the reference panel within the 400 nm-900 nm range, and the red cross represents the average value.

Figure 8 .
Figure 8.The dumbbell plot shows the uncertainty range of 10 rock samples at different single I&VG parameters and coupling of multiple parameters.

Figure 9 .
Figure 9.Comparison of the uncertainty caused by five I&VG parameters for the reference panel and two rock samples (k = 2)."Rock_low_heterogeneity" and "Rock_high_heterogeneity" represent the two rock samples DTS21-6 and DTS21-1 with the lowest and highest heterogeneity, respectively.

Figure 10 .
Figure 10.Comparison of the number of bands with the highest first-order sensitivity coefficient to five I&VG parameters for each band of the reference panel and two rock samples."Rock_low_heterogeneity" and "Rock_high_heterogeneity" represent the two rock samples DTS21-6 and DTS21-1 with the lowest and highest heterogeneity, respectively.

Figure 11 .
Figure 11.Comparison of the number of bands with the highest total sensitivity coefficient to five I&VGs for each band of the reference panel and two rock samples."Rock_low_heterogeneity" and "Rock_high_heterogeneity" represent the two rock samples DTS21-6 and DTS21-1 with the lowest and highest heterogeneity, respectively.

Figure 12 .
Figure 12.Comparison of sequential changes in the first-order and total sensitivity coefficients of the reference panel for the five I&VG parameters.The sequence from 1 to 5 ranges from the most sensitive to the least sensitive.

Figure 13 .
Figure 13.Comparison between first-order sensitivity coefficients and uncertainties caused by a single I&VG parameter for reference panel and rock samples.(a) Reference panel, (b) rock sample (DTS21-6) with lowest heterogeneity, (c) rock sample (DTS21-1) with highest heterogeneity.The sequence from 1 to 5 ranges from the most sensitive to the least sensitive.

Figure 14 .
Figure 14.Correction results of spectral reflectance from view zenith angles of 40° to 20°.(a) Rock sample DTS21-1, (b) rock sample DTS21-6, (c) soil sample.These spectra are consistent in other I&VG conditions except for the view zenith angle.

Figure 15 .
Figure 15.Reflected radiance of reference panel at different wavelengths as a function of view zenith angle with other I&VG conditions unchanged.The blue dot represents the observation result, the red line represents a fit to the observed results.

Figure
Figure16a,b show the variation in the radiance of the reference panel under different view and solar azimuth angles.It is worth noting that there is a significant nonlinear correlation between the bands.The correction factors obtained through locally weighted regression do not provide satisfactory correction results for correcting different view and solar azimuth angles.This result can be explained by the fact that most of the blue dots, which represent the measured values, do not fall on the red fitted curve.Moreover, it can be seen from the red fitted curve in Figure16cthat there is no significant trend in the radiance of different bands under different solar zenith angles.Previous studies have shown that there is only a small correlation between the SD of the reference panel radiance and the solar zenith angle[13], hence the difficulty in correcting the impact of solar zenith angles through the correction factors obtained through curve fitting.Figure16dshows the radiance distribution of the reference panel at different observation heights.From the fitting results, it is obviously better than the goodness of fit of azimuth angle.However, the correction model proposed in this paper has little effect on correcting observation heights.One possible explanation for this is that there is a very small difference in the radiance of the reference panel at different observation heights.Therefore, the correction factor is very close to 1, and the correction results for different observation heights through Equation (11) are very small and difficult to observe.

Figure 16 .
Figure 16.Reflected radiance of reference panel at different wavelengths as a function of different I&VG parameters: (a) view azimuth angle, (b) solar azimuth angle, (c) solar zenith solar, (d) observation height.Other I&VG parameter conditions remain unchanged.The blue dot represents the observation result, the red line represents a fit to the observed results.
Author Contributions: Conceptualization, Z.W., G.J. and H.Z.; Methodology, Z.W. and G.J.; Writing-Original Draft, Z.W.; Writing-Review and Editing, G.J. and J.T.; Data Curation, Z.W.; Investigation, S.J.; Resources, S.L.; Visualization, Y.L.All authors have read and agreed to the published version of the manuscript.Funding: This work is supported by the National Natural Science Foundation of China (61675012).

Figure A1 .
Figure A1.With other environmental conditions unchanged, the HCRF of the reference panel at four solar zenith angles are collected in the field and MHSRS 2 F. The yellow dashed line shows the three relatively obvious absorption characteristics used for comparison in this section.

Table 1 .
Date, time, location, and environmental conditions for in situ measurements.

Table 2 .
Setting of each parameter of I&VG in the simulation experiment of field environmental factors.

Table 3 .
(6)mary of uncertainty results for reference panel, which are caused by five I&VG parameters and the coupling of five parameters (k = 2).The uncertainty caused by a single I&VG parameter is calculated using Equation (4).The combined standard uncertainty refers to the calculations given in Equation(6).

Table A2 .
The three popular spectral characteristic parameters used at 1266 nm, 1499 nm, and 1757 nm are compared between the HCRF of the Diffuse Reflectance Standards panel collected based on MHSRS 2 F and that collected in the field.