Next Article in Journal
The Urban Observatory: A Multi-Modal Imaging Platform for the Study of Dynamics in Complex Urban Systems
Next Article in Special Issue
Machine Learning Optimised Hyperspectral Remote Sensing Retrieves Cotton Nitrogen Status
Previous Article in Journal
Applications of Satellite Radar Imagery for Hazard Monitoring: Insights from Australia
Previous Article in Special Issue
3DeepM: An Ad Hoc Architecture Based on Deep Learning Methods for Multispectral Image Classification
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 13
Issue 8
10.3390/rs13081425
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting Water Stress in Wild Blueberry Fields Using Airborne Visible and Near Infrared Imaging Spectroscopy

by
Catherine Chan
1,*,
Peter R. Nelson
2,
Daniel J. Hayes
1,
Yong-Jiang Zhang
3 and
Bruce Hall
4
1
School of Forest Resources, University of Maine, Orono, ME 04469, USA
2
Schoodic Institute, 9 Atterbury Circle, P.O. Box 277, Winter Harbor, ME 04693, USA
3
School of Biology and Ecology, University of Maine, Orono, ME 04469, USA
4
Jasper Wyman & Son, P.O. Box 100, Milbridge, ME 04658, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(8), 1425; https://doi.org/10.3390/rs13081425
Submission received: 28 February 2021 / Revised: 30 March 2021 / Accepted: 3 April 2021 / Published: 7 April 2021
(This article belongs to the Special Issue Feature Extraction and Data Classification in Hyperspectral Imaging)
Browse Figures
Review Reports Versions Notes

Abstract

:
Water management and irrigation practices are persistent challenges for many agricultural systems, exacerbated by changing seasonal and weather patterns. The wild blueberry industry is at heightened susceptibility due to its unique growing conditions and uncultivated nature. Stress detection in agricultural fields can prompt management responses to mitigate detrimental conditions, including drought and disease. We assessed airborne spectral data accompanied by ground sampled water potential over three developmental stages of wild blueberries collected throughout the 2019 summer on two adjacent fields, one irrigated and one non-irrigated. Ground sampled leaves were collected in tandem to the hyperspectral image collection with an unoccupied aerial vehicle (UAV) and then measured for leaf water potential. Using methods in machine learning and statistical analysis, we developed models to determine irrigation status and water potential. Seven models were assessed in this study, with four used to process six hyperspectral cube images for analysis. These images were classified as irrigated or non-irrigated and estimated for water potential levels, resulting in an R2 of 0.62 and verified with a validation dataset. Further investigation relating imaging spectroscopy and water potential will be beneficial in understanding the dynamics between the two for future studies.

Graphical Abstract

1. Introduction

The application and use of hyperspectral imaging technology to plant (stress) health and vegetation traits has made great advances in the past decade [1]. Hyperspectral imaging, or imaging spectroscopy, is a method used in detecting and classifying objects based on the light reflectance (or spectral signature) across narrowly-resolved wavelengths within the optical portion of the electromagnetic spectrum [2]. Applications range from differentiating tree species [3] to aiding agricultural quality control [4]. One of imaging spectroscopy’s most valued benefits is its ability to collect information in a non-destructive manner as it does not require direct contact with the scanned object [5]. In this work, hyperspectral data analysis was used to classify irrigation status of blueberry crop areas, and determine water potential of the plant canopies.
Remote sensing has played an increasingly important role in monitoring and assessing plant health. Through its growing application (i.e., added satellite missions, sensor development, etc.) [6] and user-friendliness (i.e., programming, large data management, etc.) [7], remote sensing can provide timely and quality data over different cover types [8]. Targets of plant health include physiological responses detectable through remote sensing techniques such as hyperspectral imaging [9]. These responses include nutrient declines indicating disease [10], in addition to chlorophyll content estimated through the band ratio NDVI (normalized difference vegetation index), to show sugarcane health [11]. In relation to remote sensing, imaging spectroscopy can be used as a means to estimate water potential [12].
Water management and irrigation practices are a persistent challenge for many agricultural systems in the world [13]. The effects of water stress are complex and highly influential on plant growth and productivity [14]. With many factors involved in water stress, numerous methods in quantifying water status have been developed and implemented including evapotranspiration models, soil water balance measurements, and leaf water potential [15]. Water potential is the chemical potential reduction (force that causes water movement) relative to pure water at sea level, or is a measure of the driving force of water flow [16,17] and has been shown to be an effective indicator of water stress [18,19]. Measuring water stress can further our understanding of vegetation health within natural resources, including forested lands and agricultural crops [20].
Wild blueberries (Vaccinium angustifolium Ait.) (also referred to as lowbush) are grown commercially in only Quebec, Atlantic Canada, and the state of Maine [21]. Lowbush blueberries require unique growing conditions including acidic and infertile soils, unsuitable for many other crop types [22]. This distinctive environment increases vulnerability due to its uncommon land makeup, further amplified by climate change factors [23,24]. In recent years, wild blueberry production has faced challenges relating to warming, drought, freezing, and pathogens [25,26]. As a result, forecasting land conditions and taking prompt mitigative action, including water management, have become increasingly needed. We analyzed the reflectance properties of blueberry crops throughout different development stages to detect patterns temporally and spatially.
Our goal was to assess the utility of a high-resolution remote sensing system to classify and predict irrigation status and water stress in blueberry crops. To address this goal, we collected hyperspectral imagery over wild blueberry fields in Downeast Maine using an unoccupied aerial vehicle (UAV), a remotely controlled aircraft also referred to as a drone [27]. From this imagery we then developed, applied, and tested machine-learning models to associate reflectance measurements to (1) a categorical response (irrigated or non-irrigated) [28] and (2) water potential as a continuous variable [29]. We planned to achieve this through:
  • Collecting airborne data on an irrigated and non-irrigated field over three plant development stages.
  • Acquire water potentials of canopy leaf samples in each field.
  • Generate classified maps of irrigated vs. non-irrigated areas and estimated water potentials to determine locations of low or high plant canopy water stress.
The objectives were to configure methods and processes that will allow for more efficient detection via freely available code [30]. The process of configuration involves combining spatial measurements and ground data, in addition to approaches in UAV systems, image processing, and machine learning classifiers [31]. We sought to remotely assess irrigation status and water potentials in wild blueberry fields as a means to offer techniques in crop health monitoring.

2. Materials and Methods

2.1. Study Site

Our research was conducted on commercial blueberry fields in Deblois, Maine, owned by Jasper Wyman and Son (Wyman’s), a corporation specializing in frozen fruit. Maine’s coastal region has a four-season climate with an average low of −10.6 °C and high of 24.2 °C (lowest and highest month’s average over 2006–2010), and monthly average precipitation low of 85.1 mm and high of 136.4 mm (lowest and highest month’s average over 1981–2010) [32]. The site includes an irrigated and a non-irrigated [33] field, approximately 23 and 16 hectares, respectively (Figure 1 displays these two fields). The irrigated field is uniformly irrigated with Nelson Full-Circle Impact sprinklers (Walla Walla, WA, USA). Each irrigation event compensates for natural precipitation with a maximum irrigation regime of 2.5 cm of water per acre every 3 to 4 days as needed. The location is defined by Spodosols soils of the series Masardis consisting of a fine sandy loam with 0–3% slope. These crop fields tend to contain a number of different genotypes, growing in sections within a particular field [34,35]. Due to this pattern, a number of plots were selected from the area in order to accurately assess the locations and tendencies of water potential. The irrigated and non-irrigated experimental fields are co-located adjacent to one another to reduce climate and geologic site variability and to act as a comparison between a treated and untreated site.

2.2. Workflow Overview

The project workflow is outlined in Figure 2. Input data included the ground data (irrigated/non-irrigated classes and water potential) and the manually digitized polygons from the imagery that trained the predictor models. These points and polygons are displayed on the study site fields in Figure 1. The remainder of the nodes relate to model processing steps which are detailed further below. The models were applied to images, generating maps which were compared to validation ground measurements [36,37].

2.3. Ground Sampling

Measurements were taken over the two fields three times in the 2019 spring and summer including peak bloom (7 June), green fruit (3 July), and color break (25 July). These represent some of the development stages of wild blueberries throughout the summer growing season [38,39], where peak bloom is determined when 50% of flower buds have opened, green fruit by green color and initial fruit enlargement caused by active cell division, and color break as the period of fruits changing from green to pink to blue through cell expansion. We chose to collect data over the different stages to determine temporal variation in water conditions and spectral response over the seasonal cycle of blueberry growth [40].
We used 20 ground sample measurements on each field that were selected by a random sampling design, guided by genotypic distinction to cover spatial variations. The samples entailed a small branch of leaves attached to the stem. These were gathered as the drone captured image data. The samples were immediately stored in moist plastic Ziploc bags in coolers to minimize water loss, and measured approximately two hours later for water potentials. The water potential measurements were taken using a leaf pressure chamber (Model 1505D; PMS, Albany, OR, USA) in a climate-controlled laboratory space in Wyman’s facilities.

2.4. Image Data Collection and Sampling

We used a Micro A-Series Sensor hyperspectral imaging spectrometer manufactured by Headwall Photonics (Bolton, MA, USA). The device is attached to a DJI Matrice 600 Pro UAV, and operates as a line-scanning (pushbroom) instrument. The sensor captures the visible and near-infrared portion of the electromagnetic spectrum from 400 to 1000 nm and collects 324 spectral bands. The number of flight lines, heights, and speeds were determined by site area, takeoff distance from scanning area, land topography, and specific daylight conditions. Data collection with the UAV was conducted between 10:00 AM and 2:00 PM local time to avoid shadows in the data. Image collection was limited to these timeframes, battery capacity, and the amount of time access was allowed to the field sites.
Imagery was processed using the Headwall application Spectral View. A white tarp taken in the imagery and a spectralon panel were both considered as the white reference for image processing, however the spectralon provided more accurate spectra reflectances and greater consistency. Processing entailed transforming imagery from its raw form, to radiance, reflectance, and then orthorectification. The ground sampling distance was approximately 10 × 10 cm per pixel.
Using the collected imagery, we delineated pixels of ground sampled, or known irrigation status, blueberry leaves through the ENVI software program (version 5.5 64-bit). These delineations were extracted and used as samples, which were partitions in both the training and validation data. All downstream analyses from this point were performed in the R programming environment [30].
Training pixels in the irrigated/non-irrigated categorical model were digitized in large samples. With a binary response of irrigated and non-irrigated uninfluenced by the water potential ground reference data, we digitized large polygon areas within each field to gain a larger training size. Four samples from two images of each field were digitized for each stage.
The water potential digitizations were guided by ground reference samples. Four to eight pixels were digitized around the coordinates of where samples were collected. Pixels that were in clouded or shadowed images were discarded. Two samples were removed from the first scan date, and eight from both the second and third.

2.5. Model Development

A series of spectral indices were calculated from the reflectance values to use as predictors in our models [41,42]. One of these methods in deriving the variables included resampling at 5, 10, 50, and 100 nm [43,44]. These values were selected to roughly simulate wider bandpass, common in multi-spectral sensors that are commercially available. The other method of calculating vegetation indices capitalized on the ‘vegindex’ function from package ‘hsdar’ [45], which calculates almost 100 different indices depending on the bandpass of the data. The resampling and vegetation indices totaled 260 variables as predictors.
A random forest model (classifier which uses decision trees) was trained to classify images [46] with the calibration data (either water potential or irrigated/non-irrigated) using the ‘ranger’ package to determine the relationship between the water potentials of ground measurements and the spectral index predictors. Ranger was selected rather than the ‘Random Forest’ package due to its speed and improvements in variable importance measures and bias [47].
The process of developing the models (to classify and predict images) first entailed removing intercorrelated variables among all predictors using the ‘caret’ package in R through Pearson correlation coefficient [48]. After removing these variables, we applied ranger to determine which predictors are most important. Then, we again used ranger on the top 50 important variables. A list of models was produced, each utilizing a different number of the most important variables with accuracy rates. We selected the model that utilized the fewest number of variables, but produced a low error rate.
Prior to final model selection and image classification, we performed a validation test on the water potential models. We redeveloped the models with a split of 70% calibration and 30% validation from the total sample set, randomly chosen from each field. The split was chosen due to the small sample sizes with peak bloom having 139 pixel samples, green fruit 115, and color break 99. The global model prediction, however, used an 80% and 20% split of the total sample data due to the combined sample sizes.

3. Results

We developed models that used the extracted pixel samples and water potential measurements as inputs to predict characteristics of unknown areas of imagery [49]. Using functionality developed by related projects, the products of the models were generated maps of irrigated or non-irrigated areas, and water potential classifications [30].
The imagery dataset of irrigated and non-irrigated fields consists of 48 images collected in the three phenological stages, with a spatial resolution of approximately 10 cm. Ground referenced water potentials were measured as described above. One section on the irrigated field experienced winter damage, leading to a duplicated sample to detect potential differences.

3.1. Model

Different models were developed from the digitized samples of known conditions [50]. Two models were created for each field stage, with one as a binary classification (irrigated or non-irrigated), and the other a continuous estimation (water potential). Additionally, a global model for water potential was trained with combined samples from all three stages. Table 1 outlines the model information and error rates.
Model fit improved when highly intercorrelated predictors were removed. Errors in removing intercorrelated predictor variables were calculated at 0.9, 0.93, 0.96, and 0.99 pair-wise absolute correlation cutoff levels. Removing these at a 0.99 cutoff produced the highest accuracy rates.
Each model’s sample size was the number of digitized pixels. The categorical models had out-of-bag (OOB) errors from 0.003% to 1.402%. The stage-specific (or local) models for water potential had an R2 ranging from 0.437 to 0.487. The global model utilized all pixel samples as training data and resulted in an R2 of 0.554. For further analysis, predictions were also conducted for all four of these water potential models.

3.2. Validation

To further analyze the performance of the continuous models, predictions using calibration and validation samples were conducted. The increased split ratio for calibration data in the global water potential model was increased because it improved the model’s R2 while still maintaining a validation size of about 70. Table 2 outlines each model’s OOB error, R2, and calculated RMSE.
The table shows that all three local models had very similar R2 values although calculated RMSEs were dissimilar and increased with each model. The global model had a comparable RMSE to the third local model but also had the highest R2. As a result of the prediction, we decided to use only the global water potential model for our image classifications.

3.3. Variable Importance

Each model was developed through particular predictor variables and a certain number of the top predictors [51]. These predictor variables, or spectral derivatives, consist of vegetation indices and resampled bands. Models of differing numbers of top predictors (in multiples of five) were generated in the process along with accuracy rates, however the one with a lower number of variables with a comparable lower error rate was selected to improve efficiency but maintain efficacy. The plots of the top predictors show the diminishing importance of each variable and how, despite there being a level of relevancy, inclusion of more variables does not improve model accuracy.
Figure 3 is a plot of the 35 most important predictor variables in the global model with the importance levels, and Table 3 lists the top 20 predictors that were used in the model. Bandpasses are expressed with an ‘X’ preceding the bandpass wavelength, followed by the resampling range. All of the top five predictors involve bands in the 700 to 800 nm region, a range associated with chlorophyll measurements. The top predictor was the Datt index, which was developed to estimate chlorophyll content, particularly in higher plants or tree canopies [52]. The Giltelson index is a measure of chlorophyll fluorescence, which is proportional to actual chlorophyll content [53]. The third most important variable was TOCARI2OSAVI2, which incorporates soil adjustment in chlorophyll measurements and has been recommended for agricultural applications [54].
In the supplementary, Figure S1 and Figure S2 show the top predictor variables for each model of local classification and regression. The top predictors of the three categorical models showed very distinct variable selection results with no model containing a shared top five variable. The first shared predictor was CRI2 as the fifth for color break and the eighth of green fruit. This predictor was the 24th for peak bloom. Green fruit and color break had predictors focused on chlorophyll whereas the top predictor of peak bloom was bandpass X897.593 resampled at 100 nm. In all three individual models for water potential, MTCI (MERIS Terrestrial Chlorophyll Index) was a high predictor. Medium Resolution Imaging Spectrometer (MERIS) terrestrial chlorophyll index measures chlorophyll in the red edge region distinguishing it from other red edge position indices by its sensitivity to higher levels of chlorophyll [55]. The local regression plots are included in the supplementary materials to show the similarities and differences among local predictors, and in comparison to the global model.
A component of this study assessed variable importance for each model to understand how our imaging spectroscopy methods provide information. Most top predictor variables in the models focus on assessing chlorophyll levels. Some, however, include elements of soil adjustment such as TCARI2OSAVI2 and REPLi. Most of the vegetation indices and relevant bands range between 650 to 800 nm on the electromagnetic spectrum, or in the visible to near infrared region, which is displayed in the reflectance spectra (Figure 4).

3.4. Spectral Signatures

Various spectra of blueberry samples and pixels were observed to assess radiometric calibration and provide visual representation. The critical analysis, however, was using these spectra in development of the spectral variables, being the calculations of the vegetation indices and band resampling values. Nonetheless, the signatures can showcase information of interest. Figure 4 plots the median value signatures of combined samples, one from each field over each development stage. As can be seen in the figure, the spectra of peak bloom in the shades of pink closely align in the 500 to 700 nm range, but diverge after where the reflectance of the non-irrigated field shows higher than the irrigated. The reflectance of green fruit (shades of green) and color break (shades of blue) differ from the first stage where the irrigated signatures have higher reflectance rates than the non-irrigated, although the spectra of color break have similar values in the higher section of the spectrum. The difference of the irrigated and non-irrigated spectra of green fruit are also much larger than that of color break. These median reflectance spectra display the collected data in a raw form, but the transformations of these is what the model is developed from, and used to classify imagery.

3.5. Classification

Cube images of the same area that were well-illuminated for all three stages were classified for analysis which amounted to a total of 12. Figure 5, Figure 6 and Figure 7 each represent a separate stage, or date the image was taken with an irrigated image on top, and non-irrigated image on bottom. They first display NDVI images of the two sections, followed by a classified image as irrigated or non-irrigated, and then a predicted image of water potentials. The cube images from the different stages are of the same area and extent on the irrigated and non-irrigated fields.
All categorical images were classified accurately, meaning a majority of each image was classified as it actually was. All irrigated images were classified correctly, and the same result occurred with the non-irrigated. The global water potential model estimated large portions of the irrigated fields to have higher values than the non-irrigated fields of the same stage. This was especially evident in the peak bloom stage, which is consistent with the conditions of that year as precipitation was high that spring. The collection date also took place later in the season than scheduled due to unanticipated rain dates, resulting in the collected data being more representative of the transition stage between peak bloom and green fruit. The images are also shown to decrease in values throughout the season, or become more stressed throughout the summer.

4. Discussion

This study utilized machine learning and remote imaging spectroscopy to predict water irrigation status and water potentials in wild blueberry fields. The spectral indices produced using the ranger package were found to be effective as predictor variables in classifying hyperspectral cube images, where the indices used had a logical premise and displayed spectral sensitivity through band transformations. We analyzed our developed models to understand how the results were predicted.
The predictors in the local categorical models showed more variation among the three than those predicting local water potentials. This result is likely due to the prediction method in which the binary model target was more generalized than the continuous. Assigning a broad binary response to the large input sample data may have been too encompassing for a diverse sample, leading the model to use varying predictor variables [56]. Conversely, a specific water potential value assigned to fewer input pixels could have resulted in more refined and distinct predictor variables. Predictors in the local water potential models were found to be less variable than those of the categorical models, however both of these sets also showed differences as a result of the development stages.
The models in each set of local predictors exhibited variation, potentially a result of developmental growing patterns and a reflection of the distinct characteristics of each stage used for estimation. Our observations provide a framework in which each date represented a stage, where peak bloom is characterized as healthy and lush vegetation but lower leaf area, green fruit as higher vegetation area with additional purples, pinks, and whites from flowerings, and color break expressing darker shades from fruit maturity and the start of senescing leaves. The predicted images of water potential display temporal variance between the blueberry development stages, and both categorical and continuous models showed distinction in statistics and outputs. This was shown particularly between the first stage and the following two of the categorical models. These observations are supported through peak bloom having more distinguishable phenology from the other two, where leaf area was much smaller and ground soils were relatively more exposed [57]. Capturing the bare ground may have been a factor in the soil adjusted indices that were relevant in the first blueberry scanning [58]. In contrast, green fruit and color break had more vegetation coverage which can be observed in the NDVI images of Figure 5, Figure 6 and Figure 7 (with the exception of the non-irrigated image of the color break stage, possibly a result of drought stress or other factors). Due to the respective training inputs, differences between local models are to be expected [59], however the differences between local models and one that is global can be analyzed to understand how targets are estimated.
The global model used input data from all three development stages as training data, differing from the local models which used only those data from its respective date. Its estimation method was the result of a consolidation of characteristics from peak bloom, green fruit, and color break, resulting in important indices that will account for the variability [59]. It can be observed that the top three predictors of the global model had a significant level of importance over the following ones, or in other words, the importance level of the fourth predictor is much lower than the third. This differs from the local models where the levels had a more gradual decline in importance. The diverse training samples may have resulted in generalized predictors that emphasize traits related to chlorophyll [60,61], rather than a specific index, such as the first predictor of the local water potential model of color break, which was plant senescence reflectance index (PSRI) [62]. Understanding the manner in which the global water potential model makes estimates provides guidance in interpreting the predicted images.
The predicted water potential images for each field in each development stage show variations that correspond to conditions of that field or season. As the season progressed, estimated images exhibit increased water stress in all scenes. The irrigated fields show to be less stressed than the non-irrigated, supporting the irrigation management practices. This does however raise the question of what conditions are most conducive to blueberry yield, and what stresses might be permissible [63]. Understanding the baseline levels of water potential for each stage is critical in determining favorable or damaging conditions [64]. Typically, we have observed that the irrigated field generates a much higher yield than the non-irrigated, however in 2019, the irrigated field only produced a small amount more. The cause of this decreased yield however, was determined to be winter damage from a consistent freeze in the irrigated field [65]. Overall, the water potential model was found to estimate well across the development stages. The categorical models of predicting irrigated and non-irrigated also classified accurately, however these results raise the question of how useful the information actually is.
An objective of this study was to associate reflectance measurements with a categorical treatment (irrigated or non-irrigated) to make predictions. This was achieved, however there were areas that were misclassified. These misclassifications may show the regions of the field that are particularly water stressed in the irrigated field, or water sufficient in the non-irrigated field. Regardless, the method of classifying areas as irrigated or non-irrigated prompts questioning about how constructive the answer actually is. The local models were separated using only data taken from the respective date. This was due to the relativity in deeming a field irrigated or not irrigated throughout a changing season [66]. An area classified as irrigated at the beginning of the season may not be classified the same towards the end. The irrigated/non-irrigated distinction is a convenient method in classification and may be a useful precursor to mapping water potential, but finding utility in it alone may be a challenge.
A more sensitive approach is to predict and calculate water potential which provides relevant values [67]. Measuring water potential and how it would impact irrigation practices is a more appropriate pursuit in providing farmers and landowners greater use. Although the global water potential model predicted well against the validation dataset (R2 of 0.62), adding more samples in water potential and from more than a few select dates would assist in more applicable predictions. Another addition in strengthening the model and prediction process would be to include other types of variables from other functions or calculations. Spectral band derivatives and extensions of these, another common method in determining characteristics of vegetation [68], were not heavily incorporated in this process, with the exception of a few vegetation indices that include derivatives in the index formula. It may also be helpful in reducing the resampling predictors as it may not have had a significant effect on the predictors, especially as the number of band resampling increased.
In comparison to other studies that have used hyperspectral data and vegetation indices to predict crop water potential, this study shows similar promise in predictive ability [69,70,71]. In Zhang et al. (2019), it was found that select vegetation indices were sensitive to water status variables including canopy water content, but would benefit with additional data on biomass or vegetation structure [72]. With the more rudimentary methods in our statistical analysis and classifier model (as opposed to Pôças, 2017) [73], our study further supports spectral vegetation indices through utilizing a high number of predictors. It also contributes to predictive efforts in producing maps which show spatial variation on smaller scales. Given the finer resolution, this study’s application is mainly practical for limited areas, such as agricultural fields. Our methods however, were developed for efficient and applicable use.

5. Conclusions

The goal of this project was to use hyperspectral imaging processes in detecting water stress level practically, for the benefit of the agricultural industry. This was mainly achieved through the development of our models, although there are a number of actions that could strengthen the outcomes and better show the relationship between imaging spectroscopy and leaf water potential measurements. Models will be improved if larger training datasets are used. This is particularly relevant if other blueberry fields or varying crops are to be measured, such as fields in other areas. The process of streamlining by means of data collection, processing, and efficient programming, is also critical to provide timely and accurate data that can be efficiently used by the industry, which would require greater computing capacity and data management. With continued input additions and modifications, our methods can assist in improved agricultural practices.
Many societal structures are heavily reliant on natural and agricultural resources [74]. With the limited commercial production of wild blueberries, efficient practices and maintenance in the face of climate change are critical [24,26,75]. New technologies such as those pertaining to precision agriculture are becoming more widely used [76], necessitating the adoption of new methods to maintain techniques in a changing environment. With greater sampling and in-depth studies, hyperspectral imaging methods could have many suitable applications.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs13081425/s1, Figure S1. Plot of the 35 most important predictor variables in the local models and corresponding importance levels, Figure S2. Plot of the 35 most important predictor variables in the local water potential models and corresponding importance levels.

Author Contributions

Conceptualization, C.C. and P.R.N.; methodology, C.C., P.R.N. and D.J.H.; software, P.R.N.; validation, C.C., P.R.N. and D.J.H.; formal analysis, C.C., P.R.N.; investigation, C.C., P.R.N.; resources, P.R.N., Y.-J.Z. and B.H.; data curation, C.C., Y.-J.Z.; writing—original draft preparation, C.C.; writing—review and editing, C.C., P.R.N., D.J.H., Y.-J.Z., B.H.; visualization, C.C., P.R.N.; supervision, P.R.N. and D.J.H.; project administration, P.R.N.; funding acquisition, P.R.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Maine Economic Improvement Fund (MEIF) funded by the state of Maine as well as the Research Reinvestment Fund (RRF) through the University of Maine. This project was supported by the USDA National Institute of Food and Agriculture, Hatch (or McIntire-Stennis, Animal Health, etc.) Project number ME0-41907 through the Maine Agricultural & Forest Experiment Station. Maine Agricultural and Forest Experiment Publication Number 3806.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All code is available at the following two Github pages: https://github.com/nelsopet/lecospec (accessed on 28 February 2021) and https://github.com/catherinechan70/Blueberries (accessed on 28 February 2021), please use correspondence for any additional inquiries.

Acknowledgments

We acknowledge the land used and studied in this project is occupied territory of the native Penobscot and Passamaquoddy peoples. Additionally, we would like to thank field crew members Pratima Pahadi, Arin Chen, Aldous Hofmann, Sarah Marcotte, and Valeria Briones for help with field sample collection.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Stuart, M.B.; Mcgonigle, A.J.S.; Willmott, J.R. Hyperspectral Imaging in Environmental Monitoring: A Review of Recent Developments and Technological Advances in Compact Field Deployable Systems. Sensors 2019, 19, 3071. [Google Scholar] [CrossRef] [Green Version]
  2. Green, R.O.; Eastwood, M.L.; Sarture, C.M.; Chrien, T.G.; Aronsson, M.; Chippendale, B.J.; Faust, J.A.; Pavri, B.E.; Chovit, C.J.; Solis, M.; et al. Imaging Spectroscopy and the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Remote Sens. Environ. 1998, 65, 227–248. [Google Scholar] [CrossRef]
  3. Marrs, J.; Ni-Meister, W. Machine Learning Techniques for Tree Species Classification Using Co-Registered LiDAR and Hyperspectral Data. Remote Sens. 2019, 11, 819. [Google Scholar] [CrossRef] [Green Version]
  4. Nguyen-Do-Trong, N.; Dusabumuremyi, J.C.; Saeys, W. Cross-polarized VNIR hyperspectral reflectance imaging for non-destructive quality evaluation of dried banana slices, drying process monitoring and control. J. Food Eng. 2018, 238, 85–94. [Google Scholar] [CrossRef]
  5. Thomas, S.; Kuska, M.T.; Bohnenkamp, D.; Brugger, A.; Alisaac, E.; Wahabzada, M.; Behmann, J.; Mahlein, A.-K. Benefits of hyperspectral imaging for plant disease detection and plant protection: A technical perspective. J. Plant Dis. Prot. 2017, 125, 5–20. [Google Scholar] [CrossRef]
  6. Houborg, R.; Fisher, J.B.; Skidmore, A.K. Advances in remote sensing of vegetation function and traits. Int. J. Appl. Earth Obs. Geoinf. 2015, 43, 1–6. [Google Scholar] [CrossRef] [Green Version]
  7. Ma, Y.; Wu, H.; Wang, L.; Huang, B.; Ranjan, R.; Zomaya, A.; Jie, W. Remote sensing big data computing: Challenges and opportunities. Futur. Gener. Comput. Syst. 2015, 51, 47–60. [Google Scholar] [CrossRef] [Green Version]
  8. Kerr, J.T.; Ostrovsky, M. From space to species: Ecological applications for remote sensing. Trends Ecol. Evol. 2003, 18, 299–305. [Google Scholar] [CrossRef]
  9. Koch, B. Status and future of laser scanning, synthetic aperture radar and hyperspectral remote sensing data for forest biomass assessment. ISPRS J. Photogramm. Remote Sens. 2010, 65, 581–590. [Google Scholar] [CrossRef]
  10. Abdulridha, J.; Ehsani, R.; Abd-Elrahman, A.; Ampatzidis, Y. A remote sensing technique for detecting laurel wilt disease in avocado in presence of other biotic and abiotic stresses. Comput. Electron. Agric. 2019, 156, 549–557. [Google Scholar] [CrossRef]
  11. Rahman, M.R.; Islam, A.H.M.H.; Rahman, M.A. NDVI derived sugarcane area identification and crop condition assessment. Plan Plus 2004, 1, 1–12. [Google Scholar]
  12. Asner, G.P.; Nepstad, D.C.; Cardinot, G.A.; Ray, D. From The Cover: Drought stress and carbon uptake in an Amazon forest measured with spaceborne imaging spectroscopy. Proc. Natl. Acad. Sci. USA 2004, 101, 6039–6044. [Google Scholar] [CrossRef] [Green Version]
  13. Shao, H.-B.; Chu, L.-Y.; Jaleel, C.A.; Manivannan, P.; Panneerselvam, R.; Shao, M.-A. Understanding water deficit stress-induced changes in the basic metabolism of higher plants–biotechnologically and sustainably improving agriculture and the ecoenvironment in arid regions of the globe. Crit. Rev. Biotechnol. 2009, 29, 131–151. [Google Scholar] [CrossRef]
  14. Reddy, A.R.; Chaitanya, K.V.; Vivekanandan, M. Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants. J. Plant Physiol. 2004, 161, 1189–1202. [Google Scholar] [CrossRef]
  15. Ihuoma, S.O.; Madramootoo, C.A. Recent advances in crop water stress detection. Comput. Electron. Agric. 2017, 141, 267–275. [Google Scholar] [CrossRef]
  16. Kramer, P.J.; Boyer, J.S. Water Relations of Plants and Soils; Academic Press: New York, NY, USA, 1995. [Google Scholar]
  17. Robbins, N.E.; Dinneny, J.R. The divining root: Moisture-driven responses of roots at the micro- and macro-scale. J. Exp. Bot. 2015, 66, 2145–2154. [Google Scholar] [CrossRef] [Green Version]
  18. Améglio, T.; Archer, P.; Cohen, M.; Valancogne, C.; Daudet, F.-A.; Dayau, S.; Cruiziat, P. Significance and limits in the use of predawn leaf water potential for tree irrigation. Plant Soil 1998, 207, 155–167. [Google Scholar] [CrossRef]
  19. Jones, H.G. Physiological Aspects of the Control of Water Status in Horticultural Crops. HortScience 1990, 25, 19–25. [Google Scholar] [CrossRef] [Green Version]
  20. Borengasser, M.; Hungate, W.S.; Watkins, R. Hyperspectral Remote Sensing: Principles and Applications; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  21. Rose, A.; Drummond, F.A.; Yarborough, D.E.; Asare, E. MR445: Maine Wild Blueberry Growers: A 2010 Economic and Sociological Analysis of a Traditional Downeast Crop in Transition; DigitalCommons@UMaine: Orono, ME, USA, 2013. [Google Scholar]
  22. Glass, V.M.; Percival, D.C. Proctor tolerance of lowbush blueberries (Vaccinium angustifolium Ait.) to drought stress. II. Soil water and yield component analysis. Can. J. Plant Sci. 2005, 85, 911–927. [Google Scholar] [CrossRef]
  23. Obsie, E.Y.; Qu, H.; Drummond, F. Wild blueberry yield prediction using a combination of computer simulation and machine learning algorithms. Comput. Electron. Agric. 2020, 178, 105778. [Google Scholar] [CrossRef]
  24. Tasnim, R.; Drummond, F.; Zhang, Y.-J. Climate Change Patterns of Wild Blueberry Fields in Downeast, Maine over the Past 40 Years. Water 2021, 13, 594. [Google Scholar] [CrossRef]
  25. Whittle, P. Hard Times as Maine Wild Blueberry Industry in Decline. 2018. Available online: https://www.timesrecord.com/articles/maine-1/hard-times-as-maine-wild-blueberry-industry-in-decline/ (accessed on 24 November 2018).
  26. Tasnim, R.; Calderwood, L.; Annis, S.; Drummond, F.A.; Zhang, Y.J. The Future of Wild Blueberries: Testing Warming Impacts Using Open-Top Chambers; Spire: Orono, ME, USA, 2020. [Google Scholar]
  27. Johnston, D.W. Unoccupied aircraft systems in marine science and conservation. Annu. Rev. Mar. Sci. 2019, 11, 439–463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Zhou, X.; Sun, J.; Mao, H.; Wu, X.; Zhang, X.; Yang, N. Visualization research of moisture content in leaf lettuce leaves based on WT-PLSR and hyperspectral imaging technology. J. Food Process Eng. 2018, 41, e12647. [Google Scholar] [CrossRef]
  29. Rapaport, T.; Hochberg, U.; Shoshany, M.; Karnieli, A.; Rachmilevitch, S. Combining leaf physiology, hyperspectral imaging and partial least squares-regression (PLS-R) for grapevine water status assessment. ISPRS J. Photogramm. Remote Sens. 2015, 109, 88–97. [Google Scholar] [CrossRef]
  30. Nelson, P.R. Github Repository; LECOSPEC: San Francisco, CA, USA, 2020; Available online: https://github.com/nelsopet/lecospec (accessed on 28 February 2021).
  31. Nelson, P.; Thompson, N. Visible and Infrared Imaging Spectroscopy for High Resolution Mapping and Health Assessment of Maine’s Forest and Agricultural Resources; Maine Economic Improvement Funds (MEIF) Small Campus Initiative (SCI) Research Grant: Orono, ME, USA, 2018. [Google Scholar]
  32. Anthony, A.; Imke, D.; Scott, A.; Mike, S.; Russell, V.; Xungang, Y.; Rocky, B. NOAA’s U.S. Climate Normals (1981–2010); NOAA National Centers for Environmental Information: Asheville, NC, USA. [CrossRef]
  33. Farooque, A.A.; Zaman, Q.U.; Schumann, A.W.; Madani, A.; Percival, D.C. Response of wild blueberry yield to spatial variability of soil properties. Soil Sci. 2012, 177, 56–68. [Google Scholar] [CrossRef]
  34. Annis, S.L.; Stubbs, C.S. Stem and Leaf Diseases and Their Effects on Yield in Maine Lowbush Blueberry Fields. Small Fruits Rev. 2004, 3, 159–167. [Google Scholar] [CrossRef]
  35. Alghory, A.; Yazar, A. Evaluation of crop water stress index and leaf water potential for deficit irrigation management of sprinkler-irrigated wheat. Irrig. Sci. 2019, 37, 61–77. [Google Scholar] [CrossRef]
  36. Forkuor, G.; Hounkpatin, O.K.L.; Welp, G.; Thiel, M. High Resolution Mapping of Soil Properties Using Remote Sensing Variables in South-Western Burkina Faso: A Comparison of Machine Learning and Multiple Linear Regression Models. PLoS ONE 2017, 12, e0170478. [Google Scholar]
  37. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  38. Schilder, A. Growth Stages. Michigan State University Extension. April 2015. Available online: https://www.canr.msu.edu/blueberries/growing_blueberries/growth-stages (accessed on 20 November 2018).
  39. Hepler, P.R.; Yarborough, D.E. Natural Variability in Yield of Lowbush Blueberries. HortScience 1991, 26, 245–246. [Google Scholar] [CrossRef] [Green Version]
  40. Zarco-Tejada, P.J.; Ustin, S.L.; Whiting, M.L. Temporal and Spatial Relationships between Within-Field Yield Variability in Cotton and High-Spatial Hyperspectral Remote Sensing Imagery. Agron. J. 2005, 97, 641–653. [Google Scholar] [CrossRef] [Green Version]
  41. Xue, J.; Su, B. Significant remote sensing vegetation indices: A review of developments and applications. J. Sens. 2017, 2017, 1353691. [Google Scholar]
  42. Calderón, R.; Navas-Cortés, J.A.; Zarco-Tejada, P.J. Early detection and quantification of Verticillium wilt in olive using hyperspectral and thermal imagery over large areas. Remote Sens. 2015, 7, 5584–5610. [Google Scholar] [CrossRef] [Green Version]
  43. Poona, N.; Van Niekerk, A.; Ismail, R. Investigating the utility of oblique tree-based ensembles for the classification of hyperspectral data. Sensors 2016, 16, 1918. [Google Scholar] [CrossRef] [Green Version]
  44. Adjorlolo, C.; Mutanga, O.; Cho, M.A.; Ismail, R. Spectral resampling based on user-defined inter-band correlation filter: C3 and C4 grass species classification. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 535–544. [Google Scholar] [CrossRef]
  45. Lehnert, L.W.; Meyer, H.; Obermeier, W.A.; Silva, B.; Regeling, B.; Bendix, J. Hyperspectral Data Analysis in R: The hsdar Package. J. Stat. Softw. 2019, 89, 1–23. [Google Scholar] [CrossRef] [Green Version]
  46. Boonprong, S.; Cao, C.; Chen, W.; Bao, S. Random forest variable importance spectral indices scheme for burnt forest recovery monitoring—Multilevel RF-VIMP. Remote Sens. 2018, 10, 807. [Google Scholar] [CrossRef] [Green Version]
  47. Wright, M.N.; Ziegler, A. Ranger: A fast implementation of random forests for high dimensional data in C++ and R. J. Stat. Softw. 2017, 77, 1–17. [Google Scholar] [CrossRef] [Green Version]
  48. Rinnan, Å.; Bruun, S.; Lindedam, J.; Decker, S.R.; Turner, G.B.; Felby, C.; Engelsen, S.B. Predicting the ethanol potential of wheat straw using near-infrared spectroscopy and chemometrics: The challenge of inherently intercorrelated response functions. Anal. Chim. Acta 2017, 962, 15–23. [Google Scholar] [CrossRef]
  49. Ludovisi, R.; Tauro, F.; Salvati, R.; Khoury, S.; Mugnozza Scarascia, G.; Harfouche, A. UAV-based thermal imaging for high-throughput field phenotyping of black poplar response to drought. Front. Plant Sci. 2017, 8, 1681. [Google Scholar] [CrossRef] [PubMed]
  50. Zarco-Tejada, P.J.; Camino, C.; Beck, P.S.A.; Calderon, R.; Hornero, A.; Hernández-Clemente, R.; Kattenborn, T.; Montes-Borrego, M.; Susca, L.; Morelli, M. Previsual symptoms of Xylella fastidiosa infection revealed in spectral plant-trait alterations. Nat. Plants 2018, 4, 432. [Google Scholar] [CrossRef]
  51. Fan, L.; Zhao, J.; Xu, X.; Liang, D.; Yang, G.; Feng, H.; Yang, H.; Wang, Y.; Chen, G.; Wei, P. Hyperspectral-based estimation of leaf nitrogen content in corn using optimal selection of multiple spectral variables. Sensors 2019, 19, 2898. [Google Scholar] [CrossRef] [Green Version]
  52. Datt, B. A New Reflectance Index for Remote Sensing of Chlorophyll Content in Higher Plants: Tests using Eucalyptus Leaves. J. Plant Physiol. 1999, 154, 30–36. [Google Scholar] [CrossRef]
  53. Gitelson, A.; Buschmann, C.; Lichtenthaler, H.K. The Chlorophyll Fluorescence Ratio F735/F700 as an Accurate Measure of the Chlorophyll Content in Plants. Remot Sens. Environ. 1999, 69, 296–302. [Google Scholar] [CrossRef]
  54. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  55. Dash, J.; Curran, P.J. The MERIS terrestrial chlorophyll index. Int. J. Remote Sens. 2004, 25, 5403–5413. [Google Scholar] [CrossRef]
  56. Broge, N.H.; Leblanc, E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 2001, 76, 156–172. [Google Scholar] [CrossRef]
  57. Vidhya, R.; Vijayasekaran, D.; Farook, M.A.; Jai, S.; Rohini, M.; Sinduja, A. Improved Classification of Mangroves Health Status Using Hyperspectral Remote Sensing Data. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2014, XL-8, 667–670. [Google Scholar] [CrossRef] [Green Version]
  58. Thenkabail, P.S.; Smith, R.B.; De Pauw, E. Hyperspectral Vegetation Indices and Their Relationships with Agricultural Crop Characteristics. Remote Sens. Environ. 2000, 71, 158–182. [Google Scholar] [CrossRef]
  59. Ge, X.; Wang, J.; Ding, J.; Cao, X.; Zhang, Z.; Liu, J.; Li, X. Combining UAV-based hyperspectral imagery and machine learning algorithms for soil moisture content monitoring. PeerJ 2019, 7, e6926. [Google Scholar] [CrossRef]
  60. Liang, L.; Di, L.; Zhang, L.; Deng, M.; Qin, Z.; Zhao, S.; Lin, H. Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method. Remote Sens. Environ. 2015, 165, 123–134. [Google Scholar] [CrossRef]
  61. Rumpf, T.; Mahlein, A.-K.; Steiner, U.; Oerke, E.-C.; Dehne, H.-W.; Plümer, L. Early detection and classification of plant diseases with Support Vector Machines based on hyperspectral reflectance. Comput. Electron. Agric. 2010, 74, 91–99. [Google Scholar] [CrossRef]
  62. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Rakitin, V.Y. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening. Physiol. Plant. 1999, 106, 135–141. [Google Scholar] [CrossRef] [Green Version]
  63. Jones, H.G. Monitoring plant and soil water status: Established and novel methods revisited and their relevance to studies of drought tolerance. J. Exp. Bot. 2006, 58, 119–130. [Google Scholar]
  64. Chaerle, L.; Lenk, S.; Leinonen, I.; Jones, H.G.; Van Der Straeten, D.; Buschmann, C. Multi-sensor plant imaging: Towards the development of a stress-catalogue. Biotechnol. J. 2009, 4, 1152–1167. [Google Scholar] [CrossRef] [PubMed]
  65. Panta, G.; Rieger, M.; Rowland, L. Effect of cold and drought stress on blueberry dehydrin accumulation. J. Hortic. Sci. Biotechnol. 2001, 76, 549–556. [Google Scholar] [CrossRef]
  66. Zovko, M.; Žibrat, U.; Knapič, M.; Kovačić, M.B.; Romić, D. Hyperspectral remote sensing of grapevine drought stress. Precis. Agric. 2019, 20, 335–347. [Google Scholar] [CrossRef]
  67. Cohen, Y.; Alchanatis, V.; Meron, M.; Saranga, Y.; Tsipris, J. Estimation of leaf water potential by thermal imagery and spatial analysis. J. Exp. Bot. 2005, 56, 1843–1852. [Google Scholar] [CrossRef] [Green Version]
  68. Thorp, K.; Wang, G.; Bronson, K.; Badaruddin, M.; Mon, J. Hyperspectral data mining to identify relevant canopy spectral features for estimating durum wheat growth, nitrogen status, and grain yield. Comput. Electron. Agric. 2017, 136, 1–12. [Google Scholar] [CrossRef] [Green Version]
  69. González-Fernández, A.B.; Sanz-Ablanedo, E.; Gabella, V.M.; García-Fernández, M.; Rodríguez-Pérez, J.R. Field Spectroscopy: A Non-Destructive Technique for Estimating Water Status in Vineyards. Agronomy 2019, 9, 427. [Google Scholar] [CrossRef] [Green Version]
  70. Liu, S.; Peng, Y.; Du, W.; Le, Y.; Li, L. Remote Estimation of Leaf and Canopy Water Content in Winter Wheat with Different Vertical Distribution of Water-Related Properties. Remote Sens. 2015, 7, 4626–4650. [Google Scholar] [CrossRef] [Green Version]
  71. Osco, L.P.; Ramos, A.P.M.; Pinheiro, M.M.F.; Moriya, É.A.S.; Imai, N.N.; Estrabis, N.; Ianczyk, F.; De Araújo, F.F.; Liesenberg, V.; Jorge, L.A.D.C.; et al. A Machine Learning Framework to Predict Nutrient Content in Valencia-Orange Leaf Hyperspectral Measurements. Remote Sens. 2020, 12, 906. [Google Scholar] [CrossRef] [Green Version]
  72. Zhang, F.; Zhou, G. Estimation of vegetation water content using hyperspectral vegetation indices: A comparison of crop water indicators in response to water stress treatments for summer maize. BMC Ecol. 2019, 19, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Pôças, I.; Gonçalves, J.; Costa, P.M.; Gonçalves, I.; Pereira, L.S.; Cunha, M. Hyperspectral-based predictive modelling of grapevine water status in the Portuguese Douro wine region. Int. J. Appl. Earth Obs. Geoinform. 2017, 58, 177–190. [Google Scholar] [CrossRef]
  74. Scialabba, N. Factors influencing organic agriculture policies with a focus on developing countries. In Proceedings of the IFOAM 2000 Scientific Conference, Basel, Switzerland, 28–31 August 2000; Volume 28, p. 31. [Google Scholar]
  75. Ramoelo, A.; Dzikiti, S.; Van Deventer, H.; Maherry, A.; Cho, M.A.; Gush, M. Potential to monitor plant stress using remote sensing tools. J. Arid. Environ. 2015, 113, 134–144. [Google Scholar] [CrossRef]
  76. McBratney, A.; Whelan, B.; Ancev, T.; Bouma, J. Future Directions of Precision Agriculture. Precis. Agric. 2005, 6, 7–23. [Google Scholar] [CrossRef]
Remotesensing 13 01425 g001 550
Figure 1. Study site of the irrigated field (upper light blue polygon) and non-irrigated field (lower light green polygon) with field sampling points and polygons.
Figure 1. Study site of the irrigated field (upper light blue polygon) and non-irrigated field (lower light green polygon) with field sampling points and polygons.
Remotesensing 13 01425 g001
Remotesensing 13 01425 g002 550
Figure 2. Process flow from data collection to classified and predicted images.
Figure 2. Process flow from data collection to classified and predicted images.
Remotesensing 13 01425 g002
Remotesensing 13 01425 g003 550
Figure 3. Chart of the 35 most important predictor variables in the global water potential model, the first 20 of which were included in the selected classifier, and corresponding importance levels.
Figure 3. Chart of the 35 most important predictor variables in the global water potential model, the first 20 of which were included in the selected classifier, and corresponding importance levels.
Remotesensing 13 01425 g003
Remotesensing 13 01425 g004 550
Figure 4. Spectral plots of median reflectance for all sampled blueberry pixels from each field and stage.
Figure 4. Spectral plots of median reflectance for all sampled blueberry pixels from each field and stage.
Remotesensing 13 01425 g004
Remotesensing 13 01425 g005 550
Figure 5. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over peak bloom. (a) Normalized difference vegetation index (NDVI) images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Figure 5. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over peak bloom. (a) Normalized difference vegetation index (NDVI) images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Remotesensing 13 01425 g005
Remotesensing 13 01425 g006 550
Figure 6. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over green fruit. (a) NDVI images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Figure 6. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over green fruit. (a) NDVI images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Remotesensing 13 01425 g006
Remotesensing 13 01425 g007 550
Figure 7. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over color break. (a) NDVI images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Figure 7. Processed scene images of irrigated (top images) and non-irrigated (bottom images) field sections over color break. (a) NDVI images of scenes. (b) Classified images of irrigated and non-irrigated scenes. (c) Predicted images of water potential.
Remotesensing 13 01425 g007
Table
Table 1. Model information used in image classification with each collection stage having a separate model, and combined stages for global model regression. The table also includes out-of-bag prediction errors as percentages for classification, and mean squared error and R2 values for regression.
Table 1. Model information used in image classification with each collection stage having a separate model, and combined stages for global model regression. The table also includes out-of-bag prediction errors as percentages for classification, and mean squared error and R2 values for regression.
Irrigated/Non-Irrigated
Local Classification		Water Potential
Local Regression		Water Potential Global Regression
Peak
Bloom		Green
Fruit		Color
Break		Peak
Bloom		Green
Fruit		Color
Break
Sample Size (pixels)47,75832,018103,13513911599353
Independent variables25252515252520
Out-of-bag (OOB) prediction error/MSE1.346%0.003%1.402%320679754709
R2 (OOB)NA0.4870.4580.4370.554
Table
Table 2. Outline of global and local regression model statistics for predicting water potential.
Table 2. Outline of global and local regression model statistics for predicting water potential.
Validation MetricLocal Regression PredictionGlobal Regression Prediction
Peak
Bloom		Green
Fruit		Color
Break
OOB prediction
error (MSE)		237533556616
R2 (OOB)0.5540.5630.5640.617
Calculated RSME29.836.945.246.4
Table
Table 3. Top 20 predictor variables used in the global water potential model with description and band formula.
Table 3. Top 20 predictor variables used in the global water potential model with description and band formula.
AbbreviationNameFormula
1Datt‘Chlorophyll & height’(R749 − R720) − (R701 − R672)
2Gitelson‘Chlorophyll’1/R700
3TCARI2OSAVI2Transformed Chlorophyll Absorption Ratio 2/Optimized Soil Adjusted Vegetation Index 2(3 * ((R750 − R705) − 0.2 * (R750 −R550) * (R750/R705)))/
((1 + 0.16) * (R750−R705)/(R750 + R705 + 0.16))
4X692.593_wvl_005nm ‘Bandpass 692.593 resampled at 5 nm’
5Maccioni‘Chlorophyll’(R780 − R710)/(R780 − R680)
6Datt2‘Chlorophyll & height’R850/R710
7MTCIMERIS Terrestrial Chlorophyll Index(R754 − R709)/(R709 − R681)
8DDDouble Difference Index(R749 − R720) − (R701 − R672)
9DWSI4Disease Water Stress Index 4R550/R680
10GDVI4Green Difference Vegetation Index 4(R4800 − R4680)/(R4800 + R4680)
11GIGreenness IndexR554/R677
12CARIChlorophyll Absorption Ration IndexR700 * abs(a * 670 + R670 + b)/R670 * (α2 + 1)0.5
α = (R700 − R550)/150
b = R550 − (550 *α)
13TCARI2Transformed Chlorophyll Absorption Ratio 2(3 * ((R750 − R705) − 0.2 * (R750 − R550) * (R750/R705)))
14Boochs2Single Band 703 BoochsD703
15mSRmodified Simple Ratio(R800 − R445)/(R680 − R445)
16Ctr5Carter 5R695/R670
17EVIEnhanced Vegetation Index2.5 * ((R800 − R670)/(R800 − (6 * R670) − (7.5 * R475) + 1))
18SumDr1‘LAI & % green cover’ ∑ i = 626 795 D 1 i
19TCARIOSAVITransformed Chlorophyll Absorption Ratio/Optimized Soil Adjusted Vegetation Index(3 * ((R700 − R670) − 0.2 * (R700 − R550) * (R700/R670)))/
((1 + 0.16) * (R800 − R670)/(R800 + R670 + 0.16))
20X397.593_wvl_100_nm ‘Bandpass 397.593 resampled at 100 nm’
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chan, C.; Nelson, P.R.; Hayes, D.J.; Zhang, Y.-J.; Hall, B. Predicting Water Stress in Wild Blueberry Fields Using Airborne Visible and Near Infrared Imaging Spectroscopy. Remote Sens. 2021, 13, 1425. https://doi.org/10.3390/rs13081425

AMA Style

Chan C, Nelson PR, Hayes DJ, Zhang Y-J, Hall B. Predicting Water Stress in Wild Blueberry Fields Using Airborne Visible and Near Infrared Imaging Spectroscopy. Remote Sensing. 2021; 13(8):1425. https://doi.org/10.3390/rs13081425

Chicago/Turabian Style

Chan, Catherine, Peter R. Nelson, Daniel J. Hayes, Yong-Jiang Zhang, and Bruce Hall. 2021. "Predicting Water Stress in Wild Blueberry Fields Using Airborne Visible and Near Infrared Imaging Spectroscopy" Remote Sensing 13, no. 8: 1425. https://doi.org/10.3390/rs13081425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop