Digital Commons @ Michigan Tech Digital Commons @ Michigan Tech Utilizing hyperspectral remote sensing for soil gradation Utilizing hyperspectral remote sensing for soil gradation

: Soil gradation is an important characteristic for soil mechanics. Traditionally soil gradation is performed by sieve analysis using a sample from the ﬁeld. In this research, we are interested in the application of hyperspectral remote sensing to characterize soil gradation. The speciﬁc objective of this work is to explore the application of hyperspectral remote sensing to be used as an alternative to traditional soil gradation estimation. The advantage of such an approach is that it would provide the soil gradation without having to obtain a ﬁeld sample. This work will examine ﬁve di ﬀ erent soil types from the Keweenaw Research Center within a laboratory-controlled environment for testing. Our study demonstrates a correlation between hyperspectral data, the percent gravel and sand composition of the soil. Using this correlation, one can predict the percent gravel and sand within a soil and, in turn, calculate the remaining percent of ﬁne particles. This information can be vital to help identify the soil type, soil strength, permeability hydraulic conductivity, and other properties that are correlated to the gradation of the soil.


Introduction
The gradation (grain size distribution) is an important characteristic that influences the physical and mechanical properties of the soil. Studies show that the gradation of soils carries important details of the soil's mechanical properties and that one can create an index to find correlations to these values [1]. Kuenza [2] has looked into gravel and sand composition of soils in mountainous areas with slope instability; based on the amount of sand versus the amount of gravel present it was determined which particle size has a more dominating influence over the soil's stress-strain behavior for an undrained torsional shear test. The percentage of gravel and sand content in soil influences the shear strength and dilatancy of sand-gravel mixtures. Therefore, knowing the gradation of the soil can help estimate the soil's properties [3]. While gradation is a useful metric for geotechnical engineering, its estimation requires an in situ sample collection. Presently, soil gradation is primarily done via sieve analysis (also called a gradation test). The traditional approach of in situ sample collection and sieve analysis can be laborious, costly, and challenging to attain in remote or unfriendly territories. Therefore, ways to remotely characterize soil gradation can be of significant benefit for many applications.
One such application is in the mobility of military vehicles in unknown territories. There is present work being done where terrain strength characteristics and land mapping for autonomous mobility and simulation modeling are quintessential. The methods for the past models have relied heavily on in situ measurements, which requires one to gather the soil from an area of interest, which can be hard to obtain or impossible in war-zones and remote territories, and they are looking to be improved [4]. Therefore, we present an alternative method for predicting the gradation of the soil using remote sensing that would be particularly beneficial for remote territories.
With the rapid increase of computational capabilities and sensor improvements, remote sensing has developed into a more time-efficient and accurate tool, which will continue to grow with new and emerging ideas [5][6][7][8][9]. Remote imagery provides a means of quick and efficient data collection for rapid assessment of an area. Previous work has been done on analyzing the use of different means via cameras, unmanned aerial vehicles (UAVs), and satellites to provide a way to measure sand lower bound friction angles and bearing strength [10]. The agricultural industry has recently utilized techniques such as laser-induced breakdown spectroscopy (LIBS) for quantifying nutrients within soil [11]. Other techniques, such as X-ray diffraction [12][13][14] and diffuse reflectance spectroscopy [15] have been used for soil mineral classification. These methods are limited as they are either destructive (LIBS) or costly/time consuming (X-ray diffraction). Hyperspectral remote sensing is both a nondestructive and time efficient testing method. Hyperspectral imaging has already worked towards predicting heavy metal soil contamination [16], improving our ability to do anomaly detection, target recognition, and mineral characterization [17]. Studies using both multispectral and hyperspectral sensors show promise in the use of optical remote sensing to provide a novel method for improving the mapping of soils for discriminating between variational levels of sand and clay content [18]. Using hyperspectral data, one can predict different features/characteristics of bare soil, and certain wavelengths may provide a means of identifying/distinguishing certain features [19,20]. The agriculture field has utilized soil spectral reflectance for many purposes, such as identifying reflectance effects based on moisture content, soil tilling, organic material, iron content, and minerals within the soil [21][22][23][24][25][26][27][28] and how the open soil surface patterns impact that reflectance [29]. Visible and near-infrared spectroscopy has already been used for identifying clay mineral contents within the soil [30], and shown ease of use and good agreement with identifying soil mineral composition to that of X-ray diffraction analysis [31]. In contrast, this study focuses on characterizing the soil from a soil mechanics perspective.
Hyperspectral remote sensing has the capability of determining attributes such as clay and sand content of soil within the laboratory and field. This could be especially useful when a large number of soil analysis needs to be done and may provide a real-time method for field data collection [32]. Expanding beyond this, studies have already shown promise in developing a soil index, to help improve the classification of different colored sands based on both color and curvature of reflectance spectrums. However, work still needs to be done for differing water conditions, grain size distribution, and other factors [33]. Currently, there is research being done on the best algorithms to implement for processing hyperspectral data to classify different land cover types better. Considerable success has already been found using support vector machine (SVM) methods for classifying hyperspectral data over multiple land cover types [34]. Hyperspectral imaging of topsoil with the use of regression modeling has proved to be able to predict sand and clay percentages in a bare field [35]. Hyperspectral remote sensing shows promise as a rapid and non-destructive means to measure the soil gradation. In this study, we propose the use of hyperspectral remote sensing utilizing the visible and near-infrared (NIR) range (400-1000 nm) as this range of wavelengths is more readily available and will deliver a cost-effective solution for field implementation. In our paper, we verify the applicability of hyperspectral remote sensing to characterize the gradation of five different soil types.

Materials and Methods
In this study, we analyzed five different soil types obtained from the Keweenaw Research Center (KRC) that are used for testing different terramechanics properties. These soils are labeled as Fine, Coarse, Rink, Stability, and 2NS, which have varying amounts of sand, gravel, and fine particles. The soils we examined are coarse-grained or fine-grained and have no clay or organic content. Lab analysis of the soils was done prior to testing to calculate the soil parameters such as D10, D30, and D60 (the grain size at 10%, 30%, and 60% of the sample passing through a specific sieve size), coefficient of uniformity (Cu) and coefficient of curvature (Cc) shown in Equations (1) and (2) respectively, Unified Soil Classification for Soil (USCS) [36], and confirm there were no clays or organics present. The specifics for these soils are listed in Table 1. The locations of these soils at KRC are shown in Figure 1.   The study was conducted in the lab, with each of the five soils placed into their bin with dimensions roughly 61 cm × 61 cm × 30.5 cm (Length × Width × Height). These were then set in a row under a series of work lights (1 m above soil surface) and a track to have the hyperspectral camera move over the soils in a straight and level manner, as shown in Figure 2. These work lights provide, similar to all incandescent light bulbs, a continuous spectrum of wavelengths from the near ultraviolet to the far infrared region.   2NS, Stability, Rink, Coarse, and Fine. Above each bin are two 500-Watt work lights whose field of view cover the entire soil bin, and the middle bar is the track; the hyperspectral camera is moved along.

Remote Sensing Background
Most everyday cameras operate within three bands, red, green, and blue. These cameras provide us the color images we see in a regular photograph. However, there are many bands in the electromagnetic (EM) spectrum to which certain cameras can distinguish. If the light were to reflect from a smooth white surface, one would receive a 100% reflectance since nothing is absorbed into the material. Most objects have color, roughness, texture, which all change the amount of light reflected to the camera. The reflectance values for each different wavelength is then recorded by the camera, which is then processed and analyzed to produce a multitude of information. Hyperspectral and multispectral cameras work by taking this reflected EM energy from the surface and record the response.

Sensors and Hardware
In this study, different spectral resolution hyperspectral sensors were used: a BaySpec OCI-F hyperspectral camera and a very high-resolution Analytical Spectral Devices, Inc. (ASD) Spectroradiometer. This was done as a way of checking if more economically affordable hyperspectral sensors with a slightly lower spectral resolution could also be used to detect the similar features that the ASD Spectroradiometer identified. This BaySpec camera (shown in Figure 3) has a range of 400-1000 nm covering the entire visible (VIS) spectrum and some of the near infrared (NIR) region with roughly 5-7 nm full width at half max (FWHM) spectral resolution and 120 different bands. Only the 450-1000 nm range was used, as the initial 400-449 nm had excessive noise rendering it useless. This is also typically a very low spectral reflectance range with little variations in initial values for soils, making this range very inefficient for differentiation between soils. It is a push-broom camera, meaning it collects images as it moves in a straight line. The width of the scan line is known as the swath width. Taking images at rates of up to 60 frames per second (fps), it also has a spatial resolution of 800 pixels (swath width) by the scan-length. The camera is factory calibrated (calibration fixed permanently).
It is equipped with a 16 mm (21 • FOV) lens. The full technical aspects of this camera are listed in the BaySpec manual [37]. region with roughly 5-7 nm full width at half max (FWHM) spectral resolution and 120 different bands. Only the 450-1000 nm range was used, as the initial 400-449 nm had excessive noise rendering it useless. This is also typically a very low spectral reflectance range with little variations in initial values for soils, making this range very inefficient for differentiation between soils. It is a push-broom camera, meaning it collects images as it moves in a straight line. The width of the scan line is known as the swath width. Taking images at rates of up to 60 frames per second (fps), it also has a spatial resolution of 800 pixels (swath width) by the scan-length. The camera is factory calibrated (calibration fixed permanently). It is equipped with a 16 mm (21° FOV) lens. The full technical aspects of this camera are listed in the BaySpec manual [37]. A high-resolution ASD Spectroradiometer was also used, which has a useable band range from 400 to 2500 nm. The range from 350 to 399 nm has a lot of noise, which makes the data unusable [38]. The device has a total of 512 bands in the VIS-NIR range (350-1000 nm), roughly 600 bands in the short wave infrared (SWIR) 1 range (1000-1830 nm), and similarly in the SWIR 2 range (1830-2500 nm). The spectral resolution is 3 nm (FWHM) at 700 nm, 10 nm (FWHM) at 1400 nm, and 10 nm (FWHM) at 2100 nm. Square spatial pixels are 1.4 nm for 350-1000 nm and 2 nm for 1000-25000 nm. The fiber optic cable is run from the device into a pistol grip, and then a (1 degree) scope was attached A high-resolution ASD Spectroradiometer was also used, which has a useable band range from 400 to 2500 nm. The range from 350 to 399 nm has a lot of noise, which makes the data unusable [38]. The device has a total of 512 bands in the VIS-NIR range (350-1000 nm), roughly 600 bands in the short wave infrared (SWIR) 1 range (1000-1830 nm), and similarly in the SWIR 2 range (1830-2500 nm). The spectral resolution is 3 nm (FWHM) at 700 nm, 10 nm (FWHM) at 1400 nm, and 10 nm (FWHM) at 2100 nm. Square spatial pixels are 1.4 nm for 350-1000 nm and 2 nm for 1000-25000 nm. The fiber optic cable is run from the device into a pistol grip, and then a (1 degree) scope was attached to the grip. The device was facing at nadir to the soil, during data collection, at the same location as the BaySpec camera at the center of each bin. The reflectance values of the region being scanned are collected in CSV format. Full technical aspects of the ASD Spectroradiometer are listed in the manual [39].
To collect the volumetric water content, a CR1000X is used with the CS615 probes. The CS615 water content reflectometer is connected to the CR1000X datalogger which is used to collect the measured volumetric water content. This is measured through the effect of a changing dielectric constant on the EM waves propagating along a waveguide. This method is called time domain reflectometry. The technical details of the CR1000X and CS615 are available from the manuals [40,41]. The CS615 water content probes are set to record the water content average over a 15-min interval. The probes being 30.5 cm in length are placed into the soil at around an angle to be fully submerged into the soil and reach the bottom of the bin.

Laboratory Testing Methodology
First, the volumetric water content probes were inserted into the bins, and the water content was collected for 15 min. The average moisture content was then recorded for each of the soil bins. After the time allotted, the probes were removed, and the work lights were turned on. The BaySpec camera was loaded onto a trolley, which was then moved at nadir across all five bins. Once images were collected, the BaySpec CubeCreator automatically combines all the individual slices and stitches them into one wide image. This is then used to highlight each individual bin and see the average reflectance curve over that area of interest. Each wavelength recorded will have its own specific/unique value. For example, red reflectance at 700 nm wavelength could be 20%, 0.2, while blue reflectance at 500 nm wavelength is 30%, 0.3. Figure 4 shows the soil bins (color image), the area captured by the camera (greyscale), and the hyperspectral plot after the software averages the pixel values for one of the trials done with the Remote Sens. 2020, 12, 3312 6 of 13 BaySpec camera. The BaySpec camera tested the soils at varying water contents ranging from 2 to 5% moisture content, with the Stability soil having one trial with 8%.
was loaded onto a trolley, which was then moved at nadir across all five bins. Once images were collected, the BaySpec CubeCreator automatically combines all the individual slices and stitches them into one wide image. This is then used to highlight each individual bin and see the average reflectance curve over that area of interest. Each wavelength recorded will have its own specific/unique value. For example, red reflectance at 700 nm wavelength could be 20%, 0.2, while blue reflectance at 500 nm wavelength is 30%, 0.3. Figure 4 shows the soil bins (color image), the area captured by the camera (greyscale), and the hyperspectral plot after the software averages the pixel values for one of the trials done with the BaySpec camera. The BaySpec camera tested the soils at varying water contents ranging from 2 to 5% moisture content, with the Stability soil having one trial with 8%. The ASD Spectroradiometer is then placed with the optical gun facing directly down at nadir to the soil in the center. It was placed in such a manner that there were no shadows being cast onto the The ASD Spectroradiometer is then placed with the optical gun facing directly down at nadir to the soil in the center. It was placed in such a manner that there were no shadows being cast onto the soil. A series of 10 consecutive tests were done over this area and then averaged to get the final reflectance for each soil type.
Once all the soils had been recorded on both devices, the lights were shut off, and water was poured onto the soils via a watering can to simulate rainfall. The BaySpec camera was not used after the first test due to the processing time of the images collected, so it will only be examining drier conditions. Different moisture content levels were tested, ranging from 2 to <15% by adding various amounts of water in between trials and allowing around 10-15 min for the water to be absorbed by the soil to see if the moisture content influenced the results. Then the process was repeated using only the ASD Spectroradiometer to collect this data. The soil images and values from both sensors are then processed, and reflectance curves are gathered for each soil bin.

Results
Based on these hyperspectral curves, a novel soil classification index was developed for each soil, Equation (3). The developed index provides a means to identify the percent sand and gravel within a soil using hyperspectral remote sensing and leaving the remaining amount of soil to be comprised of fine particles. Needing only the reflectance values at 550, 600, and 650 nm range of a hyperspectral image, the range is within the visible spectrum ( Figure 5). These values we feel are representative of roughly the first 10-15 cm of the soil as the exposed surface is the topsoil which should have approximately the same gradation in this range. The index did not vary much with water being added, in other words increasing the moisture content of the soil, as this is a normalized value and is concerned with the slope of the reflectance curve, not the exact reflectance values. After completing multiple trials, these bands provided max spectral variance among the classes (gravel, sand, and fine). A more extensive search Remote Sens. 2020, 12, 3312 7 of 13 using machine learning will be carried out in the future. The index was checked at different moisture contents by adding water to the soil ranging from 0 to 4000 mL, see Figure 6.
representative of roughly the first 10-15 cm of the soil as the exposed surface is the topsoil which should have approximately the same gradation in this range. The index did not vary much with water being added, in other words increasing the moisture content of the soil, as this is a normalized value and is concerned with the slope of the reflectance curve, not the exact reflectance values. After completing multiple trials, these bands provided max spectral variance among the classes (gravel, sand, and fine). A more extensive search using machine learning will be carried out in the future. The index was checked at different moisture contents by adding water to the soil ranging from 0 to 4000 mL, see Figure 6.  Apart from Stability, due to a few outliers in the dry conditions, the soil classification index remained fairly stable. Figure 7 shows the relationship for both % gravel and % sand versus the soil classification index when the soil moisture content ranges from 2% to less than 15% moisture content when using the ASD FieldSpec Spectroradiometer. Results show that we only had difficulty with one type of soil for predicting the % gravel, but the trend did not fail even when the soil moisture was increased. The prediction for % sand, however, worked for all the five soil types at varying moisture contents. No direct relation was found between the soil classification index and the % fine content, so instead knowing the remainder after excluding % gravel and % sand would be comprised of solely the % fine content. Next, the relationships for when conditions are known about the moisture content: 2 to <3%, 3 to <5%, 5 to <10%, and 10 to <15% moisture content are shown in Figure 8, with % gravel in orange and % sand in blue. The BaySpec hyperspectral camera was also used to find similar trends. The BaySpec camera was tested for values ranging from 2 to 5% moisture content, with one soil (Stability) having a single trial at 8%, and the plots are displayed in Figure 9. Apart from Stability, due to a few outliers in the dry conditions, the soil classification index remained fairly stable. Figure 7 shows the relationship for both % gravel and % sand versus the soil classification index when the soil moisture content ranges from 2% to less than 15% moisture content when using the ASD FieldSpec Spectroradiometer. Results show that we only had difficulty with one type of soil for predicting the % gravel, but the trend did not fail even when the soil moisture was increased. The prediction for % sand, however, worked for all the five soil types at varying moisture contents. No direct relation was found between the soil classification index and the % fine content, so instead knowing the remainder after excluding % gravel and % sand would be comprised of solely the % fine content. Next, the relationships for when conditions are known about the moisture content: 2 to <3%, 3 to <5%, 5 to <10%, and 10 to <15% moisture content are shown in Figure 8, with % gravel in orange and % sand in blue. The BaySpec hyperspectral camera was also used to find similar trends. The BaySpec camera was tested for values ranging from 2 to 5% moisture content, with one soil (Stability) having a single trial at 8%, and the plots are displayed in Figure 9.
contents. No direct relation was found between the soil classification index and the % fine content, so instead knowing the remainder after excluding % gravel and % sand would be comprised of solely the % fine content. Next, the relationships for when conditions are known about the moisture content: 2 to <3%, 3 to <5%, 5 to <10%, and 10 to <15% moisture content are shown in Figure 8, with % gravel in orange and % sand in blue. The BaySpec hyperspectral camera was also used to find similar trends. The BaySpec camera was tested for values ranging from 2 to 5% moisture content, with one soil (Stability) having a single trial at 8%, and the plots are displayed in Figure 9.  contents. No direct relation was found between the soil classification index and the % fine content, so instead knowing the remainder after excluding % gravel and % sand would be comprised of solely the % fine content. Next, the relationships for when conditions are known about the moisture content: 2 to <3%, 3 to <5%, 5 to <10%, and 10 to <15% moisture content are shown in Figure 8, with % gravel in orange and % sand in blue. The BaySpec hyperspectral camera was also used to find similar trends. The BaySpec camera was tested for values ranging from 2 to 5% moisture content, with one soil (Stability) having a single trial at 8%, and the plots are displayed in Figure 9.
(a) (b) Figure 9. Trend lines for % gravel (a) and % sand (b) content among the five soils as recorded by the BaySpec hyperspectral camera.
When compared to the ASD, the trendlines for both % gravel and % sand versus the soil classification index, have similar trends as shown in Figure 10. The BaySpec and ASD trends slightly differ, while maintaining an accuracy of 86.2% (BaySpec) and 80% (ASD) for percent gravel When compared to the ASD, the trendlines for both % gravel and % sand versus the soil classification index, have similar trends as shown in Figure 10. The BaySpec and ASD trends slightly differ, while maintaining an accuracy of 86.2% (BaySpec) and 80% (ASD) for percent gravel estimation when under 5% moisture content. The % sand trends between the ASD and BaySpec devices are nearly parallel, while the prediction accuracy is either 73.3% (ASD) or 77.9% (BaySpec), also under 5% moisture content. Note that each trendline and its corresponding accuracy were made in accordance with each device's collected set of data, as seen in Figure 10.
Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 14 estimation when under 5% moisture content. The % sand trends between the ASD and BaySpec devices are nearly parallel, while the prediction accuracy is either 73.3% (ASD) or 77.9% (BaySpec), also under 5% moisture content. Note that each trendline and its corresponding accuracy were made in accordance with each device's collected set of data, as seen in Figure 10.

Discussion
Trends for both % gravel and % sand versus the soil classification index have been shown in Figures 7-10. Similar to Lucas' work [33] who found a brightness index to distinguish between 22 different colored natural sands, a similar range of spectral values for our index to identify the gradation was used. Moreover, in agreement with Sousa [18] we found using hyperspectral imagery

Discussion
Trends for both % gravel and % sand versus the soil classification index have been shown in Figures 7-10. Similar to Lucas' work [33] who found a brightness index to distinguish between 22 different colored natural sands, a similar range of spectral values for our index to identify the gradation was used. Moreover, in agreement with Sousa [18] we found using hyperspectral imagery to identify our needed bands for the soil classification index very useful, this does turn the final model into a more "multispectral approach" yet still has great promise for soil mapping. Both the BaySpec and ASD Spectroradiometer detected these trends with Rˆ2 values of 0.86 and 0.80 for % gravel and 0.78 and 0.73 for % sand. Other researchers [32] have been able to identify sand content using seven bands (in the range of 950-2400 nm) and getting an Rˆ2 value of 0.85 within their laboratory study. This group had also found correlations to clay content in contrast to us identifying the % gravel and % fine content of the soil. Selige's work [35], however, was able to utilize HyMap imagery with partial least squares regression (PLSR) and multilinear regression (MLR) to get Rˆ2 values of 0.94 and 0.87 for sand content, respectively. That process used more predictors, nine factors for PLSR and four bands for MLR and the wavelengths used were in the range of 2200-2400 nm. While work is being done in the wavelength range of 2000-2500 nm for clay and mineral identification research, this study is focused on simple gradation of the soil (% gravel, % sand, and % fine content) as the BaySpec camera does not reach that wavelength range. Both these groups use wavelength ranges past the visible spectrum, use more bands thereby requiring a larger computational time than using our index, and would require a more expensive sensor to collect the bandwidths. Figure 6 shows that various amounts of water being added 0-4000 mL to the soil bins had little effect on the soil classification index. If the water content of the soil is known, then Figure 8 allows for better prediction accuracy in all but the 3 to <5% water content for sand, and 0 to <3% and 5 to <10% water content for gravel. In other cases, the better prediction model comes from the 2% to <15% water content model in Figure 7.
Furthermore, for Figure 6, with the soil classification index remaining constant throughout adding varying milliliters of water to the soil, we assumed that 5-10 min was long enough for the water to soak into the soil. This may have led to skewing the moisture content: high water content in the top layer and drier in lower depths of the soils. This may be especially true within the Fine pit, with the highest fine particle content among the soil types making it the least permeable and potentially influencing its reflectance values. The moisture content probes, however, do an average over the entire soil depth, and the water had clearly soaked in beneath the visible surface of the soil.
Another potential reason for the fine soil not working for the soil classification index estimating % gravel in both devices could be due to the gravel content being zero, rather than having at least some gravel present as with the other four soils. Once the values of % gravel and % sand are known, the remaining makeup of the soil will consist of fine particles. Thus, giving one the entire particle size composition of the soil under examination. Specifics for the five soils tested are listed in Table 2. All the testing was done using pure soil samples and no organics present, which could skew the results when tested in the field. If prior knowledge of the site about whether or not organics are present would also allow for further classification to the soil type at hand to help get a better classification of the study site. A relation for remotely sensed data and the Cc and Cu could also be created to allow for a full USCS classification [36] of the soil down to specific soil type.

Conclusions
In conclusion, the present approach for soil gradation using sieve analysis of in situ soil samples can be costly or difficult to obtain. We have shown that using both a high resolution and a more cost-efficient hyperspectral sensor, one can pick up specific gradation details about each of the five different soils used in this study.
With the improvement of cameras and cost to develop hyperspectral and multispectral sensors reducing, these devices are going to become a much more readily available tool. These improvements in sensors have allowed us to create a soil classification index using values from the visible region of the reflectance spectrum. Utilizing the soil classification index, it is possible to remotely characterize the percent gravel, sand, and the remaining being fine soil. This method only requires the reflectance values at 550, 600, and 650 nm within a visible/near infrared hyperspectral or multispectral image, given the sensor collects these particular bands. This soil classification index represented approximately the top 10-15 cm of the soil and remained relatively constant for each soil with increasing the moisture content. Knowing the approximate or actual moisture content of the soil helps improve the accuracy of the model.
Keeping in mind the limitations, being that of different sensor resolutions and what bands the hyperspectral sensors can pick up may not precisely be 550, 600, and 650 nm wavelengths. Different platforms and the presence of organics in the soil may also have different results than that which this study has found. Different minerals will also have an influence on the soil classification index, such minerals like Fe-Oxide minerals where goethite and hematite both have an absorption band around 500 nm [30]. In future work, we hope to expand this method into testing many different soil types outside the five used in this study. We also hope to investigate variations of the predictions of percent gravel and sand, based on organics being present in the field. Finally, we also would like to take this study to the UAV platform in the field and see how this study holds up outside of a constrained environment at differing spatial resolutions, potentially allowing for another index to identify the Cu and Cc of the soil to classify using the USCS system.