Spatial Analysis of Soybean Plant Height and Plant Canopy Temperature Measured with On-the-Go Tractor Mounted Sensors

There is a growing interest in the Open Ag community to use inexpensive sensors controlled by open-source software to measure plant height and plant canopy temperature of agricultural crops. Plant height and plant canopy temperature are key indicators of plant health. This research study reports on an ongoing research initiative to test a compact and inexpensive mobile sensor to measure plant height and plant canopy temperature. The system is controlled by open source software and hardware. The specific objectives for this study were to analyze the relationship between plant height and plant canopy temperature of soybeans (Glycine max L.) measured with the mobile system and to analyze the spatial correlation of the plant height and plant canopy temperature measurements. Data were collected in a soybean plot in 2018 and 2019. Descriptive statistics, Pearson correlation, and geostatistical techniques were used to evaluate the data. A negative statistically significant (p ≤ 0.05) relationship was observed between the plant height and the plant canopy temperature measurements (r = −0.54, 2018; r = −0.37, 2019). Also, both parameters were spatially correlated; however, plant height had a greater spatial continuity than plant canopy temperature. Furthermore, similar patterns were observed for the in-field variability of the plant height and plant temperature maps derived via kriging. Similarities in plant height and plant canopy temperatures were observed from one year to the next, suggesting that the sensor technologies could be used as a historical record for monitoring growth patterns in soybean fields. The sensors and techniques used in this study can be easily adapted to other crops, thus providing two important


Introduction
Plant height and plant canopy temperature are two important plant parameters in agricultural production systems. They are key indicators of plant growth, development, and yields [1] [2] [3]. Plant height and/or plant canopy temperature sensors have been developed by researchers to record data while mounted on a mobile platform [4] [5] [6]. Mobile plant height and/or plant canopy temperature sensors have been used to obtain measurements for various crops, including sugarcane (Saccharum officinarum L.), cotton (Gossypium hirtusuum L.), corn (Zea mays L.), soybean (Glycine max L.), wheat (Triticum aestivum L.), and wild blueberry (Vaccinium angustifolium Aiton) [4] [6] [7] [8].
Recently, the agricultural community has become interested in the opensource market for sensor development because the user can share and modify system components to meet their needs. Open-source technologies have shown promise for the following agriculture applications, monitoring water levels in evaporation pans used in evaporation studies and irrigation scheduling [9], checking soil water content and multiple soil air-vegetation parameters [10], monitoring air humidity, air temperature, soil moisture, and soil temperature [11], and tracking meteorological conditions associated with downy mildew (Hyaloperonospora brassicae) development in vineyards and determining the most suitable time to treat downy mildew disease [12].
Recently, [5] developed a plant-canopy monitoring system with inexpensive sensors and open-source software, demonstrated how the system components worked, and discussed data output and potential applications of the data as a survey tool. The sensor system can be easily mounted on an agricultural vehicle such as a tractor to record plant height, plant canopy temperature, and location data as the device is moved throughout the field. The components for their four-row system only cost $292.00, and fabrication of the parts took only four hours. Currently, information is lacking on the relationship between plant height and plant temperature data recorded with the sensors on that system and on the potential of using the system to develop maps showing in-field variability of plant height and plant canopy temperature. This research study builds on the work of [5] and focuses on the following objectives: 1) analyze the relationship between plant height and plant canopy temperature of soybeans measured with on-the-go ultrasonic and infrared thermometers and 2) analyze the spatial variability of the plant height and plant canopy temperature measurements. Positive findings from this study will further support using inexpensive systems and  Figure 1). The study plot is referred to as Field 14 ( Figure 1). The field consisted of Sharkey clay, (0.5% to 2% slope, very-fine, smectitic, thermic Chromic Epiaquerts) and Tunica clay (0% to 2% slope, Clayey over loamy, smectitic over mixed, superactive, nonacid, thermic Vertic Epiaquepts).The average rainfall, maximum temperature, and minimum temperature of Washington County, MS are 135 cm, 33.6 degrees Celsius in July, and 2.2 degrees Celsius in January, respectively [13]. Major crops grown in Washington County include soybeans, corn, cotton, rice (Oryza sativa L.), and wheat [14].

Data Collection
The sensor described by [5] was used in this study. It consisted of a batterypowered, microcontroller-based system containing a Bluetooth radio, global positioning receiver, four infrared temperature sensors (plant canopy temperature measurements), four ultrasonic distance sensors (plant height measurement), micro secure digital card, and voltage regulation components. The sensors were installed in plastic enclosures, with one temperature and one ultrasonic sensor in each enclosure, with enclosures mounted at 1-m spacings along a horizontal aluminum channel. The system was then mounted on the front of an agricultural spray vehicle to collect data simultaneously from four crop rows. Prior to data collection in the field, the sensor system was driven to a bare soil area adjacent to the study plot; the distance between the plant height sensors and bare ground area was measured to establish the maximum height of the sensor above ground level, the reference distance. This information was used to calculate the plant height measurements obtained with the ultrasonic sensor in the field. Additionally, the distance between the sensor and bare ground level was measured with a tape measure to evaluate the accuracy of plant height sensors. The plant height sensors measurements were ±2 cm of the manual measurements. In postprocessing of the data, the plant height was determined by subtracting the measured distance obtained from the plant canopy sensor during field measurements from the bare ground reference-distance measurement.
The study plot was divided into four-row transects for data collection. Figure  2 shows an example of the sensor system collecting data from one of the

Post Processing
For further analysis, we only used data collected from the two center sensors.
Measurements from these sensors coincide with other field sensors that are being tested on this plot. Data were cleaned prior to further analysis including removal of inaccurate values and points that were collected outside of the field boundary. Then, the plant height and plant canopy temperature data were averaged for the two center sensors.

Statistical Analysis
Descriptive statistics (i.e., minimum, maximum, mean, standard deviation, skewness), histograms, and frequency distribution plots were used to evaluate the normality and variability of the plant height and the plant canopy temperature datasets. Pearson correlation (p ≤ 0.05) [15] was employed to determine the relationship between the plant height and plant canopy temperature measurements. The R software ["Action of the Toes" [16]] was used to calculate the descriptive statistics and the Pearson correlation coefficient.

Summary Statistics and Correlation Measurements
Descriptive statistics are summarized in Table 1  Moderate to weak statistically significant (p ≤ 0.05) negative correlations were observed between the plant height and the plant canopy temperature ( Table 2).
The best correlation occurred for the 2018 data measurements, indicating the strength of the relationship can change from one year to the next.

Variogram Models
The plant height and the plant canopy temperature data were spatially correlated (Table 3). Isotropic exponential models best described the plant height spatial correlation for the 2018 and 2019 field measurements and the plant canopy temperature measurements for the 2019 dataset. An isotropic spherical model best described the spatial variability of the plant canopy temperature in 2018.

Plant Height and Plant Canopy Temperature Maps
Plant height and plant canopy temperature maps derived from block kriging are shown in Figure 3. In 2018, the tallest soybean plants were observed in the southwest region of the plot; shorter plants were detected in the mid central and southeast sections of the field (Figure 3(a)). The plant canopy temperature in-field variability pattern was similar to the plant height measurements' in-field variability pattern (Figure 3(b)). The lowest canopy temperature occurred in the southwest region of the plot where the tallest soybean plants were mapped. For the 2019 results, the tallest soybean plants occurred in the eastern and the southwestern sections of the plot (Figure 3(c)). There was a general agreement in the variability patterns of plant height and plant canopy temperature in that increases in plant height were associated with decrease in plant canopy temperature (Figure 3(c) and Figure 3(d)). The plant height differences observed in the southwest section of the field were consistent with the plant height differences observed in 2018. Furthermore, in 2019, the shortest plants were observed in the northern section of the field plot.

Discussion
The location and number of samples are important for understanding and mapping the spatial variability of plant and soil parameters. Within approximately 1.5 hrs, the mobile system acquired approximately 1000 plant height and plant canopy temperature measurements to use for studying the variability of those components. Additionally, the sensors can acquire measurements at any location within a field that is accessible by tractor.
Overall a negative relationship was observed between plant height and plant canopy temperature (Table 2). Reference [20] also observed a negative relationship between plant height and plant canopy temperature for similar sensors used for a soybean and wheat phenotyping study focusing on plant breeding. The difference between our study and their study was that our goal was to simulate using the plant height and plant canopy temperature system in normal farming operations. With that aspect in mind, the average speed of the tractor recorded by the global positioning system sensor was 5 km•hr −1 during data collection.
Future research initiatives need to focus on speeds commonly used by agricultural producers because the final goal is to be able to use this system simultaneously with other farming operations implements.
For geostatistical applications, [17] [18] indicated that 100 to 150 points are needed to develop good variogram models for creating maps via kriging. The sensing system collected more than enough points to meet that criteria. Geostatistical analysis showed that plant height and plant canopy temperature variation was spatially dependent for both years of data collection. The plant height was spatially correlated for longer distances than the plant canopy temperature for both years (see range parameters in Table 3). Maps derived from kriging indicated that the plot could be divided into zones, which was consistent with other researchers using on-the-go sensors to measure plant canopy height and volume [21]. The differences in plant height observed on the maps were verified during ground-truthing of the field, thus justifying accuracy and potential use of the system as a survey tool. Furthermore, we observed consistency in the plant height and the plant temperature variability from one year to the next. Thus, measurements obtained with this system can serve as a historical record of temperature and plant growth patterns observed in soybean fields.

Conclusion
The tractor is the primary tool used to get things done on a farm.