Abstract

In the spectrum measurement experiment, the roughness of the object surface is an essential factor that cannot be ignored. In this experiment, a group of mixed pixel samples with different mixing ratios were designed, and these samples were printed on four kinds of papers with different roughness. The spectral characteristics of mixed pixels with different roughness are quantitatively analyzed by using the measured spectral data. The linear spectral mixture model is used for spectral decomposition, and the effect of roughness on the unmixing precision of mixed pixels was studied. The surface roughness will affect the reflectivity of the mixed pixel. Specifically, the higher the roughness is, the higher the reflectivity of the sample is. This phenomenon is more noticeable when the proportion of white endmember (PWE) is large, and as the white area ratio decreases, the reflectance difference gradually decreases. When the surface roughness of the sample is less than 3.339 μm, the spectral decomposition is performed using a linear spectral mixing model in the visible light band. The average error of the unmixing is less than 0.53%, which is lower than the conventional standard spectral measurement error. In other words, when the surface roughness of the sample is controlled within a specific range, the effect of roughness on the unmixing accuracy of the mixed pixels is small, and this effect can be almost ignored. Multiple scattering within the pixels is the key to model selection and unmixing accuracy, when using the ASD FieldSpec3 spectrometer to perform spectral reflectance measurement and linear spectral unmixing experiments. If the surface roughness of the sample to be measured is less than the maximum wavelength of the spectrometer, the experimental results believe that the photon energy is mainly mirror reflection on the surface of the object and diffuse reflection. At this time, it is still a better choice to use a linear spectral mixing model to decompose the mixed pixels.

1. Introduction

The information recorded by each pixel in the remote sensing image is the sum of the radiant energy of all target scenes in the ground range corresponding to the instantaneous field of view angle of the detection unit; such a pixel is called a mixed pixel [1, 2]. Whether it is aerial remote sensing or aerospace remote sensing, the spectral information obtained by the sensor is a mixed spectrum containing multiple target scenes, so almost all pixels in remote sensing observations are mixed pixels [3]. There are three main reasons for mixed pixels. First is the composition of the detection target itself: the influence of geometric structure, position distribution, and other factors, resulting in the generation of mixed pixels similar to a specific type of scene [4]. Second reason is due to the different observation scale of the sensor, resulting in different spatial resolution of the image [5]. Third, this situation is caused by the background influence of the detection target, such as mixed pixels caused by target shadows, ground fluctuations, and atmospheric scattering [6]. Because the above reasons are unavoidable in actual remote sensing observations, these factors are indispensable when conducting mixed pixel spectral measurement experiments should all be considered.

The surface of most objects in nature is a rough surface relative to visible wavelengths. The surface of most objects is a rough surface that radiates relatively to visible wavelengths [7]. When the radiant brightness of an object’s surface does not change with the angle of view, the object is called a Lambert reflector [8]. However, the surface of the object has anisotropic characteristics for the scattering of solar short-wave radiation, and almost all objects have non-Lambertian features. The reflection of non-Lambert reflector surfaces is between diffuse reflection and specular reflection. There are multiangle and multidirectional spatial distribution characteristics [9, 10]. The surface roughness of objects is one of the main factors affecting emission, scattering, and radiation, which is essential for the research and application of passive microwave remote sensing. In the visible light and near-infrared band, roughness is one of the essential structural parameters on the surface of the scene, which directly affects the distribution characteristics of the target’s retro reflection [11]. The phenomenon of visible light reflected by the surface of the object has two types, specular reflection and diffuse reflection, which is usually related to the roughness of the surface. When the roughness of the surface is less than the wavelength of the human wave, specular reflection is formed. At this time, the reflection angle is equal to the human angle. The highly polished metal plate is secularly reflected [1214]. Diffuse reflection is formed when the uneven surface size is larger than the wavelength of the incident wave, and the reflection line is very irregular. For example, the heating furnace tube, the refractory brick surface, and so on are mostly diffuse reflection [15, 16]. The spectral reflectance of a mixed pixel not only depends on the internal determining factors, such as endmember composition and composition structure, but also is influenced by external observation conditions factors such as observation azimuth, observation distance, atmospheric pressure, temperature, humidity, and meteorological conditions. In addition to the above factors, the influence of the surface roughness of the mixed pixels cannot be neglected. Regarding the effect of degree, especially when the roughness of the sample surface is substantially different, roughness is an essential factor that affects the spectral characteristics and the effect of spectral unmixing. Xie et al. used the atomic force scanning probe microscope to obtain the surface roughness data of granite. Taking into account the incident shielding and multiple scattering effects in the surface scattering process, the bidirectional reflection distribution function of the actual surface was simulated [17]. Based on the field-measured multiangle hyperspectral reflection information, Lu and Zhao analyzed the influence of surface roughness variation on the reflection information distribution of sandy land and used the model parameters to invert the roughness characteristics [18]. Kang et al. use multiangle observation information to soil surface roughness, exploring a new method for estimating surface parameters by combining optical and microwave remote sensing data. In order to simulate the distribution characteristics of the surface reflection of soil surface pixels under different roughness conditions, a geometric optical model of bare soil surface and reflection is established [19].

In terms of mixed pixel spectral unmixing, researchers have proposed a large number of mixed pixel decomposition models, mainly based on linear spectral mixing models and nonlinear spectral mixing models [20]. Among them, the linear spectral mixing model regards the mixed pixel spectrum. It is an area-weighted linear combination of component spectra within a pixel. The model believes that light does not interact between pixel components. The mixed spectrum is formed after the spectrum of each pixel component enters the sensor. A large number of laboratory and field experiments show that, in the case of uniform illumination and a smooth surface, the signal received by the detector is mainly a single scattering, and the linear spectral mixture model has a better decomposition effect. Because the linear spectral mixing model is more close to remote sensing transmission theory and convenient and practical, it has been widely used in the study of mixed pixel decomposition [21, 22]. However, it is affected by factors such as atmospheric scattering, topographic fluctuations, and target surface roughness. In these cases, the multiple scattering effects between the internal components of the pixel are enhanced. At this time, the nonlinear spectral mixing model is more in line with the actual situation, so it is particularly important to choose different models according to the actual situation [23]. In terms of model decomposition accuracy, Wang and others compared the effects of linear spectral mixing model (LSMM) and Constraint Linear Spectral Mixing Model (CLSMM) in spectral mixing decomposition. The results show that whether using all spectral band data or using the band-selected data, LSMM is in the mixed spectrum. The decomposition effect is better than CLSMM [24]. Luo and Gong simulated the mixed spectrum of soil and vegetation in the Poyang Lake with a mixed model of linear and nonlinear spectra, and the results show that the effect of the linear model is better than that of the nonlinear model. Each of the nonlinear models has its applicable conditions and ranges [25]. Chen and Wang found that atmospheric conditions and terrain factors will affect the decomposition accuracy of the linear spectral mixing model. In small and medium spatial scales with terrain fluctuations, the influence of terrain factors is more significant [26]. Yang et al. measured the laboratory spectra of multiple sets of samples with different materials and component ratios and found that the decomposition accuracy of the linear mixing model is lower when the mixed samples have similar mixing ratios [5]. Similar studies show that the characteristics of spectral reflection of mixed pixels are related to the topological structure of the target and the incident zenith angle. Studies have shown that the more scattered the targets are, the smaller the reflectivity will be, and the larger the incident zenith angle will be, the smaller the reflectivity will be. So considering these two factors will improve the accuracy of mixed pixel decomposition [27].

The spectral characteristics of mixed pixels are the result of a variety of internal and external factors. At present, there are many kinds of research on the factors affecting the decomposition of mixed pixels. Most of them focus on the influence of factors such as structural differences within the pixels, atmospheric scattering, topographic fluctuations, and target roughness on the decomposition accuracy. Nevertheless, the effects of target surface roughness on spectral reflectance and spectral unmixing are often ignored. In order to reveal the effect of roughness on spectral unmixing, other factors need to be kept constant during the design of the experiment, and only one factor of roughness is allowed to change. The effect of roughness on the spectral unmixing of mixed pixels was studied. By changing the type of printing paper, four sets of samples with different roughness were designed. Comparing the results of different unmixing with the same roughness and different roughness reveals the effect of surface roughness on the accuracy of spectral unmixing of mixed pixels.

2. Materials and Methods

2.1. Sample Design and Spectral Measurement

The ASD FieldSpec3 was used in this experiment, and the probe is equipped with a halogen light source. The inner diameter of the detector is 21 mm, and the outer diameter is 27 mm. Besides, the front field of view angle of the probe is 25°, the waveband range is 350–2500 nm, and the resolution of the measured spectral data is 1 nm [28]. According to the size of the probe, the design sample diameter is 21 mm, the inner diameter of the black ring is 21 mm, and the outer diameter of the black ring is 27 mm. The sample mixing design method is shown in Figure 1. It is ensured that the probe is aligned with the sample during the spectral measurement. Samples of mixed pixels are printed on four kinds of photo papers: glossy paper, matte paper, frosted paper, and satin paper. The specifications of the four types of photo paper are standard A4 size (210  297 mm), and the weight of each sheet is 13.5 g.

In the laboratory, the endmember spectrum and the mixed pixel spectrum are measured separately. When measuring the spectrum of each sample, the probe is vertically close to the sample and horizontal. Each sample needs to be measured four times at a 90-degree interval; then the spectral processing tool was used to smooth the data. After preprocessing, such as spike correction, the arithmetic average is taken as the original spectral reflectance data for each sample. This experimental study involves spectral data in the visible range. For the impact of the infrared band, the data are resampled from 380 to 780 nm at 10 nm intervals. The purpose of resampling the spectral data is to reduce the amount of data and eliminate the data redundancy of adjacent bands by compressing the bands. Doing so will not reduce the method’s ability to distinguish the roughness characteristics of mixed pixels. The design diagram of the spectral measurement experiment is shown in Figure 2.

2.2. Linear Spectral Mixing Model

In recent years, the linear spectral mixing model has been widely used in the field of mixed pixel decomposition because it is closer to the remote sensing transmission theory and is convenient and practical. The linear spectral mixing model is to treat the mixed pixel spectrum as the area-weighted linearity of the component spectrum within the pixel. In addition, the model shows that light does not interact between pixel components, and the mixed spectrum is formed after the spectrum of each pixel component enters the sensor. The model assumes that the spectral reflectance of the pixel is the spectral reflection of each endmember inside the pixel. The ratio is a linear combination of weighting coefficients based on its percentage of the pixel area, which is the area-weighted average of the endmember spectrum (equation (1)). Based on LSMM, constraints are added (equation (2)) to form a constrained linear spectral mixture model [23].

RMSE represents the square root of the square sum of the deviation between the observed value and the true value and the ratio of the number of observations. It is used to measure the deviation between the observed value and the actual value. It is used to measure the deviation between the observed value and the actual value. According to the CLSMM model equation, the sample is calculated using the endmember spectral data and the percentage of the pixel area. Regarding the reflectivity and RMSE at each band, RMSE is used to represent the decomposition error of each mixed pixel sample; the equations are as follows:

In equations (1) and (2), represents the spectral reflectance of the mixed pixel, n is the number of endmembers, ai is the area percentage of the ith endmembers in the pixel, fλi is the spectral reflectance of the ith endmember in the λ band, and ελi is the error value. In equation (3), RMSE is the root mean square error of the sample, m represents the number of bands, Ri is the calculated reflectivity of the ith sample, and represents the measured reflectivity of the ith sample.

2.3. Proportion of White Endmember and Surface Roughness

This article uses the proportion of white endmember (PWE) to characterize the mixing ratio of the mixed pixels; that is, the ratio of the white ends to the area of the mixed pixels, expressed in decimals:

In equation (4), PWE represents the mixing index of the sample, AB represents the area of the black endmembers in the sample, Aw represents the area of the white endmembers, PWE ∈ [0, 1], 0 represents pure black pixels, and 1 represents pure white pixels.

Surface roughness is represented by Ra, which refers to the small spacing and the unevenness of tiny peaks and valleys on the surface of an object. The smaller the surface roughness, the smoother the surface. The arithmetic mean deviation of the contour is the arithmetic mean of the absolute deviation of the contour within the sampled length. In other words, it is the arithmetic means of the absolute value of the distance from each point on the surface contour to the midline of the contour within the sample length. In the actual measurement, the more the number of measurement points, the more accurate Ra. The surface roughness test of this test sample uses a handheld high-precision roughness meter (Jitech TR200). The instrument gives the parameters which are not only in compliance with Chinese standards (China GB/T 6062-2002) but also compatible with ISO, DIN, ANSI, JIS, US, German, Japanese, and British national standards. The instrument is suitable for metal and nonmetal tested parts with rough surfaces. The degree of measurement is not affected by high temperature and high humidity working environment. The measurement range of the roughness meter is 0.005–16.00 μm, the measurement accuracy is 0.001 μm, and the measurement error is ≤±7%. Surface roughness test points at different photo papers. The numbering design is shown in Figure 3.

3. Results and Discussion

3.1. Influence of Mixing Ratio and Surface Roughness on the Reflectivity of Mixed Pixels

In order to reveal the influence of mixing ratio and surface roughness on the reflectivity of mixed pixels, the black and white endmembers and the reflectance data of three typical mixed pixel samples are PWE = 0.8, 0.5, and 0.2. These samples were selected for comparative analysis. It can be seen from Figures 4(a)4(d) that the reflectance of the white endmembers on the satin paper and matte paper is slightly higher than one at 440–450 nm. The phenomenon is caused by the Rayleigh scattering, but the small overflow amplitude will not affect the subsequent analysis results. Therefore, there is no need to modify it. The reflectance of the three typical mixed pixels gradually decreases with the decrease in the proportion of white endmembers. The three groups of spectral curves are evenly spaced, indicating that the effect of the experimental data better reflects the effect of the mixing ratio. Regarding reflectivity difference, further comparison shows that when the white area is large, the reflectance of the same mixed pixel on different photo papers varies greatly. This difference is particularly significant in the trough region of 550–600 nm. When PWE = 0.8, by mixing, the reflectivity of the pixels from high to low is satin paper (blue)> matte paper (green)> matte paper (red)> glossy paper (black). As the proportion of white in the mixed pixel decreases, the difference in reflectance of mixed pixels of different photo papers gradually decreases. When PWE = 0.1, when the peaks and troughs of the reflectance spectrum disappear, the spectral curve appears as a straight line close to 0 (380–700 nm band interval). When the proportion of white endmembers in the mixed pixel is close to 0, there is a reflection on the four kinds of photo paper. The rate curves are the same, indicating that the difference in reflectance caused by roughness at this time is negligible.

In order to quantitatively reveal the difference in reflectance of different photo paper mixed pixels caused by the change in the ratio of white endmembers, the average of the reflectance of all mixed pixels in the visible light band of 380–780 nm was statistically calculated (Table 1). The reflectivity of the mixed pixel on four kinds of paper is divided into four groups. Moreover, the absolute value of the range within the group is calculated, which can reflect the difference of reflectivity. As can be seen from the table, the average reflectance of the satin paper with large surface roughness is higher, and the roughness is small. The reflectivity of glossy paper is low. This trend is not apparent on some mixed pixels, which is mainly because the data of the 380–440 nm and 700–780 nm bands are included in the statistical average, and the two band ranges are affected by Rayleigh scattering and near-infrared effects. Although the difference will be lost, the integrity of the visible light data is retained to facilitate the analysis of the spectrum change law.

As can be seen from Figure 5, the reflectance range decreases with the decrease in the proportion of white endmember, and this trend is undeniable. To further analyze the change rule, we perform a Pearson correlation analysis of the white area ratio and the reflectance range and find that there are significant correlation characteristics and the correlation coefficient r is 0.796. The significance test Sig is less than 0.01, which is a significant positive correlation. We perform a function fit on the two and find that the linear equation y = ax + b has the best fitting effect and R2 reaches 0.634. In summary, we found that the reflectivity of mixed pixels in different photo papers will vary due to the surface roughness of the photo paper. This difference in reflectivity is closely related to the PWE of mixed pixels. There are some differences between pure black endmembers and pure white endmembers on different photo papers and whether the surface roughness will lead to the accuracy of linear spectral unmixing needs to be further verified by error contrast experiments.

3.2. Analysis of Test Data for Surface Roughness of Four Kinds of Photo Paper

Use a surface roughness tester to measure the roughness of four kinds of photo paper: glossy paper, matte paper, matte paper, and satin paper. Nine test points were selected for each photo paper. The roughness measurement data of 9 test points were counted, and the roughness of each paper was expressed by arithmetic mean value. The specific results of the roughness of the photo paper are shown in the following table. From Table 2, the test results of the three kinds of photo paper can be seen: glossy paper, suede paper, and matte paper are ideal. The difference is small. The satin paper has roughness data measured at only 3 of the 9 test points, and the difference is significant. The remaining 6 points show out of range. It is speculated that this may be due to glossy paper or suede paper. The surface texture particles of sandpaper and matte paper are relatively evenly distributed. In contrast, the satin paper is unevenly distributed, which causes the roughness to exceed the measuring range of the instrument in most cases.

3.3. Analysis of the Influence of Surface Roughness on Linear Spectral Unmixing of Mixed Pixels

The linear spectral mixing model is used to unmix all the mixed pixels, and the actual theoretical reflectance and the observed reflectance are calculated to obtain the root mean square error RMSE. The changes of the RMSE can intuitively reflect the effect of the linear unmixing of the mixed pixels. It can be seen from Figure 6 that the distribution of the unmixing error of the mixed pixels on each type of photo paper shows a parabolic function shape. With PWE = 0.5 being the symmetry of the axial direction, which gradually decreases on both sides, when PWE of the mixed pixel is 0.4, at 0.6, the unmixing error is more significant. When the PWE equals 0.05 or 0.95, the unmixing error is the smallest. The unmixing error of the mixed pixels on the satin paper is the largest, and the unmixing error of the mixed pixels on the suede paper is the smallest. However, on the whole, the mixed pixel unmixing error on the four kinds of roughness photo paper is less than 0.08, and in some cases, the unmixing error is 0.01, and the average unmixing error of 76 mixed pixels is 0.042. According to statistical principles, it can be known that for the prediction of continuous variables, we can measure it with the help of RMSE mean square error. When the RMSE is smaller, the result means that the model works better. In other words, the predicted value will be closer to the actual value. When RMSE approaches 0, the surface roughness has little effect on the resolution accuracy of mixed pixels.

The influence of surface roughness on the unmixing accuracy of the linear spectrum was further explored, and the average unmixing error, average reflectance, and unmixing error in the band range of 380–780 nm were calculated. The difference can be seen in Table 3. The error of the unmixing of mixed pixels on the same type of photographic paper is 0.92%, and the maximum error of unmixing of the same type of mixed pixels on different photographic papers is 0.69%. The average unmixing error of all mixed pixels is only spectral reflectance measurement error of 0.53%, less than 1%. Experimental results show that although the four kinds of photo papers with different surface roughness have different spectral reflectance, they will not affect the results of linear spectral unmixing. Combined with the roughness data, when the surface roughness Ra of the test sample is less than 3.339 μm, use a spectrometer to perform the spectral measurement and perform a linear spectral unmixing experiment. The effect of roughness on the unmixing accuracy can be ignored.

According to the basic rules of ground object radiation and reflection, the phenomenon of photon energy reflected by the object surface includes specular reflection and diffuse reflection, which is related to the roughness of the object surface. This phenomenon is related to the roughness of the surface of the measured object. When the surface roughness is less than the human wave, spectral reflection is formed when the wavelength is 5%, and the reflection angle is equal to the angle of incidence at this time. When the surface roughness is greater than the wavelength of the incident wave, diffuse reflection is formed. The reflection line is not unique and irregular. The structural complexity of the mixed pixel itself and the external changes in the environment lead to different situations in which light is scattered under different experimental scenarios, so it is particularly important to choose different models according to the actual situation. Studies show that, under stable light conditions and smooth surfaces, the waveband was received by the detector. The signal is mainly single-scattering. At this time, the linear spectrum mixing model has a better decomposition effect. When the surface of the pixel is rough or fluctuates, the multiple scattering effects between the internal components of the pixel are intensified. The model is more in line with the actual situation. Therefore, the influence factor of multiple scattering in the pixel is the key to model selection and unmixing accuracy. The experiment uses linear spectral mixing. The error between the results obtained by the linear spectral mixing model and the actual situation is less than 1%. The experimental results show that if the surface roughness of the measured sample is less than the maximum wavelength of the spectrometer, the photon energy on the surface of the measured object is mainly specular reflection, with almost no diffuse reflection. At this time, it is still an excellent choice to decompose the mixed pixel by linear spectral mixing model.

4. Conclusion

The papers with different surface roughness have a specific effect on the reflectivity of the mixed pixels. This phenomenon is more apparent when the proportion of white endmember (PWE) is significant, and as the white area ratio decreases, the reflectance difference gradually decreases. When the mixed pixels are close to the pure black endmembers, the reflectance of the mixed pixels on different photo papers is the same. When the surface roughness Ra of the sample is less than 3.339 μm, the mixed pixel decomposition experiment is carried out in the range of visible light band using the linear spectral mixing model. The average unmixing error of the mixed pixel with different roughness is no more than 0.53%, which is less than 1% of the standard spectral measurement error. When the surface roughness of the sample is controlled within a specific range, the effect of roughness on the unmixing accuracy of the mixed pixels is small, and this effect can be almost ignored. When the ASD FieldSpec3 spectrometer is used for spectral reflectance measurement and linear unmixing experiment, multiple scattering in the image element is the key to model selection and unmixing accuracy. If the surface roughness of the sample under test is less than the maximum wavelength of the spectrometer, the photon energy on the surface of the object is mainly specular reflection, and diffuse reflection hardly occurs. At this time, it is still an excellent choice to decompose the mixed pixel by linear spectral mixing model.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request after this project has been finished in 2021.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The study was supported by the remote sensing geochemistry disciplinary innovation team, Kunming University of Science and Technology, Kunming, China. Furthermore, the work was funded by the National Natural Science Fund of China (41101343) and 2018 Yunnan Province Ph.D. Academic Newcomer Award.