**1. Introduction**

Remote sensing (RS) is the practice of obtaining information regarding an object, an area, or a phenomenon through the analysis of images acquired by means of a device that does not make physical contact with them [1]. RS investigations are mostly based on the development of a deterministic relationship between the amount of reflected emitted or backscattered electromagnetic energy, in specific bands or frequencies, and the chemical, biological, and physical characteristics of the studied phenomena. Conventionally, RS has been associated with satellites and manned aircrafts equipped with different sensors [2]. In recent years, considerable technological developments—mainly concerned with the use of unmanned aerial vehicle (UAV) platforms—have been registered [3]. Born for military applications, UAVs have become a common tool for use in geomatics for data acquisition in various research and operational fields: regional security, monitoring of structures and infrastructures, monitoring of archeological sites, environmental monitoring, application in agriculture, etc. [2,4,5]. Indeed, low-altitude RS by means of UAVs is one of the most powerful tools in precision agriculture (PA).

The International Society of Precision Agriculture (ISPAG) defines PA as "a management that gathers, processes and analyzes temporal, spatial and individual data and combines it with other information to support management decisions according to estimated variability for improved resource use efficiency, productivity, quality, profitability and sustainability of agricultural production" (www.ispag.org, last access 18 April 2020). The general stages of PA practice involve data collection, field

variability mapping through the use and development of algorithms, decision making, and management practice [6].

In this framework, thermal RS has been found to be a promising tool by measuring surface temperature [7]. In the last years, thermal sensors have also gained popularity because of the improvements in sensor technology and the reduction in costs. The surface temperature detected by thermal sensors has been found to be a rapid response variable to monitor plant growth and stress [7,8]. Indeed, the temperature is a fundamental environmental variable that plays an essential role in plant physiological processes, such as transpiration, leaf water potential, and photosynthesis [9].

The potentialities of the application of thermal UAV RS concern several mapping and monitoring issues, such as yield estimation, plant phenotyping [10,11], plant water stress detection [12,13], and plant disease detection [14], as explained in the further paragraphs of this review. When considering these last two applications, the capability to identify, using sensors mounted on UAVs, crop stress, and health before said crops are significantly damaged could be a crucial goal.

Concerning water stress, it has been found that thermal images show a correlation between minor changes in water stress that are undetectable by the normalized difference vegetation index (NDVI) [15]. In this regard, temperature-based indices represent a fast and practical way to evaluate and estimate crop water status, indicating plants' water content [16]. Among these indices, the crop water stress index (CWSI) [17], which is the most used, is often exploited to monitor plants' water status and, therefore, for the management of irrigation [18]. Such indices have been applied to different crops, both tree and herbaceous, such as olives [19], grapevines [20], sugar-beet [21], maize [22], rice [23], wheat and cotton [24].

However, the correct interpretation and use of a thermal image are based on basic knowledge of the nature of thermal radiation [1]. Thermal images can be influenced by various factors, including the characteristics of the thermal camera, meteorological conditions, and several sources of emitted and reflected thermal radiation [7]. For this reason, the aspects that are related to calibration, ground data collection, the step in which the use, and the measurement of the temperature of reference panels are recommended [25], and data processing must be carefully carried out for correct temperature retrieval.

The objective of this review is to cover the state of the art of thermal UAV RS in the framework of PA, outlining an overview of the latest applications. The structure of this paper is as follows. In Section 2, a brief description of the basic principles of thermography is provided. Section 3 deals with the characteristics of thermal cameras; the calibration, data collection, and data processing aspects are discussed. Section 4 is devoted to the main thermal RS in PA. Finally, in Section 5, we provide the present—as well as future—challenges in thermal RS.

#### **2. Basic Principles of Thermography**

Thermal RS uses the information at the emitted radiation in the thermal infrared (TIR) range (Figure 1) of the electromagnetic (EM) spectrum [26]. This information is converted into temperature [7]. Two categories can be distinguished within the IR region (0.7–100 μm), namely, the reflected-IR (0.7–3.0 μm) and TIR (3.0–100 μm). Generally, all of the elements of the landscape, such as vegetation, soil, water, and people, emit TIR radiation in the 3.0–14 μm portion of the EM spectrum [27].

In this range of the EM spectrum, part of the IR energy is transmitted to the Earth's surface through two so-called atmospheric windows that range from 3 to 5 μm and from 8 to 14 μm [28] (Figure 2). Almost all of the radiation between 5 and 8 μm is absorbed by atmospheric gasses (i.e., water, carbon dioxide (CO2), and ozone (O3) molecules), as indicated in Figure 2 [27] (Figure 2).

**Figure 1.** The electromagnetic (EM) spectrum. In evidence, the infrared (IR) region, in which the reflected-IR (0.7–3.0 μm) and the emitted-IR (3.0–100 μm) are further detailed.

**Figure 2.** Atmospheric transmittance in the thermal region with typical absorption bands induced by gasses and water (modified from [29]).

The physical laws of Planck, Wien, Stefan–Boltzmann, and Kirchhoff make it possible to better understand the behavior of EM radiation. According to Planck, each energy element (Q) is proportional to its frequency (ν), while the Planck's constant h is used to adjust this relationship (Equation (1)):

$$\mathbf{Q} = \mathbf{h}\mathbf{v} \tag{1}$$

When considering that the frequency of the wave (ν) is directly proportional to the speed of light (c) and inversely proportional to its length *(*λ), Equation (1) can be rewritten as follows Equation (2):

$$\mathbf{Q} = \frac{hc}{\lambda} \tag{2}$$

In other words, the energy of a quantum is inversely proportional to its wavelength. Thus, the longer the wave (larger wavelength), the lower its energy. Wien's and Boltzmann's laws describe the relationship of a black body's radiations (i.e., an ideal object that absorbs and reemits all of the incident energy) and the wavelength of the maximum emission with a black body's temperature [30].

Wien's displacement law explains the relationship between the true temperature of a black body, being expressed in degrees Kelvin, and its peak spectral exitance or dominant wavelength [27]. As the temperature increases, its maximum exitance shifts towards shorter wavelengths [31]:

$$
\lambda\_{\text{max}} = \frac{b}{T} \tag{3}
$$

Given that the Wien's constant *b* is equal to 2898 μm K, this formula indicates the wavelength at which the maximum radiant spectral exitance can be obtained. It is possible to observe such effects in nature. For example, a body with a very high absolute temperature like the sun (about 6000 K) has a λ*max* and, thus, a peak of emission in the visible part of the spectrum [32]. The result of the above formula is useful to indicate what should be the measurement range of the sensor used to measure the radiation emitted by a given body [27].

The Stefan–Boltzmann law states that the emittance of a black body is proportional to the fourth power of its absolute temperature (Equation (4)):

$$\mathbf{E} = \ \sigma \mathbf{T}^4\tag{4}$$

where E represents the spectral radiant exitance expressed in W·m<sup>−</sup>2, σ is the Stefan–Boltzmann constant, and T is the absolute temperature [K]. The formula clearly shows that the total EM radiation that is emitted by a black body is a function of its absolute temperature. Therefore, the radiation emitted by a body increases as its temperature increases, as mentioned above.

Kirchhoff's law (1860) (Equation (5)) states that, at a given wavelength, the emittance of a body is equal to its absorption capacity, which is:

$$
\varepsilon = \mathfrak{a} \tag{5}
$$

where ε represents the emittance and α the absorbance.

This is frequently formulated as "good absorbers are good emitters and good reflectors are poor emitters" [27]. The principle of energy conservation is defined by the following equation (Equation (6)):

$$
\varepsilon + \mathfrak{p} + \mathfrak{r} = 1 \tag{6}
$$

where ρ is the reflection and τ is the transmission. When considering that most objects are opaque to TIR radiation, the above equation becomes (Equation (7)):

$$
\varepsilon + \mathfrak{p} = 1 \tag{7}
$$

Materials with a high ε absorb a large quantity of incident energy and radiate large quantities of energy, while materials with low ε absorb and radiate less energy (Kirchhoff, 1860) [33]. All bodies with a temperature above absolute zero are characterized by random movement—i.e., the kinetic heat, whose measure is the kinetic temperature Tkin [27]. Besides, an object emits energy as a function of temperature, and the emitted energy is used in order to determine its radiant temperature Trad [28]. Although there is a strong positive linear correlation between Tkin and Trad, Trad is lower than Tkin

due to emissivity (ε) [34]. For this reason, the temperature that is measured by a sensor (Trad) will always be lower than the real temperature (Tkin) (Equation (8)) [31]. It follows from this, as explained by Kirchhoff's law (Equation (8)), that:

$$\mathbf{T\_{rad}} = \ \varepsilon^{\frac{1}{4}} \ \mathbf{T}^{\text{kin}} \tag{8}$$

The emissivity is the ratio between the radiation that is emitted by the surface and the radiation emitted by a black body at the same temperature [35]. Because the radiance of any real body, at the same temperature, is always lower than that of a black body (equal to 1), its emissivity has a value between 0 and 1 [27]. Anybody absorbs and emits radiation less effectively at a given temperature when compared to a black body, as explained by Planck's law. Practically, the sensed temperatures of materials with low emissivity appear to be much lower than those of nearby objects with the same temperatures, making the Tkin assessment less precise [31]. Many factors influence emissivity: color, chemical composition, surface roughness, moisture content, field of view, viewing angle, spectral wavelength, etc. [27,35,36]. The emissivity of materials is difficult to measure, although it is more or less constant in the region of the EM that ranges from 8 to 14 μm. The emissivity of vegetation ranges from 0.96 to 0.99, while that of soil is around 0.89, and the emissivity of the water is 0.99 (Table 1). In particular, focusing on vegetation, emissivity values of the leaves are available in literature for many plant species [37,38] and some studies were dedicated to the determination of emissivity in herbaceous and tree species [39,40].

**Table 1.** The emissivity of different surfaces over the range of 8–14 μm [1,36].


Regarding the intrinsic characteristics of a body, which affect its emissivity, thermal capacity (or heat capacity, measured in J kg<sup>−</sup>1K<sup>−</sup>1) measures the quantity of heat energy that is necessary for a body to increase its temperature by one degree—the lower the heat capacity, the less energy required. The thermal conductivity (measured in W m−<sup>1</sup> K<sup>−</sup>1), which measures the rate at which heat passes through a material, is another important parameter. This capacity is greater the higher the value of such a parameter in a material. These parameters can be integrated into the thermal inertia (expressed in J m−<sup>2</sup> K−<sup>1</sup> s−<sup>1</sup> <sup>2</sup> ), which measures a body's tendency to change in temperature or the rate of heat transfer between two substances put in contact. This concept is of great importance in the field of TIR, because the capacity of a body to quickly change its temperature (if the thermal inertia is low) depends on this parameter.

#### **3. Thermal Cameras and Unmanned Aerial Vehicles (UAVs)**

Thermal cameras typically carry a sensor that detects the infrared radiation emitted by a body, displaying its temperature in a digital radiometric image. Two types of thermal cameras are currently available: scanning devices that allow for capturing a point or a line and those with a two-dimensional infrared focal plane array. The latter, which allows for capturing all of the elements of an image at once, is faster if combined with a better image resolution [41] and it is the most commonly used [42]. A further distinction can be made between thermal and photon (or quantum) detectors. The latter convert directly absorbed EM into a change in the distribution of electric energy in a semiconductor by varying

the concentration of free charge carriers [42] and require a cooling system, usually made while using helium or liquid nitrogen at a temperature of −196 ◦C [27]. The cryocooler lowers the temperature of the sensor with the aim to reduce the thermally induced noise to a level that is lower than that of the signal coming from the image [43]. The more efficient the cooling system, the more accurate the measurements of the instrument, which makes cooled thermal cameras more precise and accurate [44], thus allowing the detection of the slightest temperature differences in the image. These cameras generally work in the mid-wavelength infrared (MWIR) region (3–8 μm), where the thermal contrast is high [42]. The thermal radiance of a target is easier to detect once it is distinguishable (higher or lower) from the background [43]. Unfortunately, cooled sensors, besides being large and expensive, have higher energy consumption and they are not suitable for UAVs because of their weight [45]. Thermal detectors are less sensitive (±0.1 ◦C) and they are slower than quantum detectors, but, instead, have the advantage of not requiring cooling elements [46]. Different types of uncooled detectors are available, all of which are made of different and unconventional materials; the three most common types are VOx (vanadium oxide), amorphous silicon (α-Si) microbolometers, and ferroelectrics [47]. The operating principle of thermal detectors is based on the conversion of absorbed EM radiation into thermal energy [48]. Ferroelectric detectors are based on the ferroelectric phase transition that can be detected in some dielectric materials. The microbolometer is a resistor organized in arrays—called focal plane arrays, which are made up of VOx and α -Si—which is composed of a thermometer, integrated on a micro-bridge, and an adsorber. Temperature increases, which are caused by absorption of IR radiation, determine large fluctuations of its electrical resistance, which can be converted into electrical signals and processed in order to generate an image [42,49–52], whose geometric resolution depends on the number of detectors. The low values of temperature differences equivalent to noise in uncooled thermal sensors, which reach 20 m ◦C, allows it to be used in applications where previously only cooled thermal sensors could be used [44]. Furthermore, especially due to the rapid development of micro- and nanotechnology, microbolometers have become cheaper and more efficient [53]. The lenses are made of germanium—a shiny semi-metal, chosen because of its transparency to infrared radiation—and they are reflective of visible radiation [42]. Several parameters characterize a thermal camera; among these, worthy of a mention, include the temperature range measured, usually between −20 and +120 ◦C, and the thermal sensitivity, which determines the minimum value of temperature difference (ΔT) detectable in an image and generally ranges from 40 to 20 m ◦C for uncooled and cooled devices, respectively [42]. Regarding the geometric resolution, to date, it is still very low when compared to RGB cameras (currently, higher image resolutions ranging from 320 × 240 to 640 × 512 pixels), while the spectral resolution normally ranges from 7 to 14 μm [12]. As for the price of cameras, their cost can vary from €1.000 to more than €10,000, depending on sensor resolution and radiometric calibration accuracy [54,55].

Cooled thermal infrared cameras have the widest use in satellite and aerial RS, thanks to their thermal sensitivity and precision [56]; on the other hand, these devices are larger and more expensive, also in terms of energy, than the uncooled ones [57]. In contrast, uncooled thermal cameras are normally mounted on UAVs (Figure 3a,b), because they are smaller, lighter, and have lower energy consumption [58].

The limited payload is perhaps the principal limit of UAVs, together with the limited battery life, which affects the duration of the flight. For example, it is not unlikely that 45 minutes or more pass between the first and the final image of the dataset when larger areas need to be covered and several flights are required [59].

Payload integration considerably changes between cooled and uncooled TIR cameras, where different ventilation modes are a key factor in image quality; thus, generally, an uncooled microbolometer is preferred for its weight benefits [60].

**Figure 3.** DJI Phantom 4 pro (**a**) and (**b**) DJI Inspire 1, both equipped with a FLIR uncooled thermal camera (photos taken by the authors).

#### *3.1. Camera Calibration and Data Collection*

Low-cost thermal cameras, which are generally not radiometrically calibrated, can only provide information regarding relative temperature differences [54]. The data provided by such instruments are represented in the form of raw digital numbers (DNs) that express radiance. Even when using radiometrically calibrated UAV cameras, it is not easy to derive accurate surface temperature measurements, because their low accuracy is due to the presence of uncooled microbolometers [54]. The sensitivity and, therefore, the accuracy of each microbolometer, is affected by both the temperature of the focal plane array [61] and that of the other components of the thermal camera, body, and lens, such as to create a weak signal-to-noise ratio [53].

Besides this, other causes make it necessary to calibrate a thermal camera, as shown in [53] and [57]. The atmosphere affects the quality of the thermal image as it absorbs and emits IR. The atmospheric effects on UAVs at low-altitude measurements are considered to be negligible when compared to airborne or satellite measurements [31].

The effect of relative humidity, air density, and altitude can only be avoided by making measurements within about 10 m or less of the target's surface [62]; under different conditions, their effects must be taken into account. Meteorological conditions can have an indirect effect on the temperature measurements of uncooled thermal cameras [59]. In the field, some recommendations are useful for reducing the effects; for instance, the critical factors for data acquisition are the time of day, weather, and the knowledge of the surrounding environment [60]. Regarding the applications in PA, midday has been identified as the best time for flying in terms of thermal accuracy [63] and the reduction of background effects [64].

UAV thermal surveys should always be done in the absence of clouds, rain, snow, smoke, dust, or any other darkening agents, because all of these reduce atmospheric transition, as well as change the temperature of the background [62]. Before the flight, after powering up, the temperature of the camera sensor is expected to stabilize from 20 minutes [25] up to one hour [65]. It is also essential to perform a test flight to allow the camera to acclimatize to local weather conditions, as well as to protect the thermal camera with a casing (to reduce air temperature effects on the sensor) when mounted on a quadcopter UAV [66].

Besides, temperature calibration in the field is an essential step. The temperature of the presumed hottest and coldest objects within the area of image acquisition should be measured for ground truth calibration [31]. Field calibration should be performed while using the thermal measurements of target surfaces during the flight [20,67]. These are placed in the area of the study where the temperature is measured using thermocouples [59] or infrared thermo-radiometers [25], combined with data-loggers capable of recording the temperature throughout the flight. Targets can be made using black and white polypropylene panels, which represent the thermal extremes of the study site, and whose size is such as to represent several homogeneous pixels in the infrared thermal image [25].

In our experience, as is the case in Messina et al.'s [68] work, wet and dry reference surfaces were used during the thermal UAV survey by taking images of them before and after the flights, as indicated by [59]. Reference surfaces were placed close to UAV take-off and landing points, and their temperature was measured while using a handheld infrared thermometer (FLIR E6) (Figure 6b).

The temperature was measured at three moments of the flight: At take-off, during, and end of every flight. Four reference surfaces were used: three dry panels (black, grey, and white) and one created delimiting a square piece of the ground while using circular targets and covered with aluminum (Figure 4).

**Figure 4.** The temperature reference targets used during thermal unmanned aerial vehicle (UAV) surveys in an onion crop field in Calabria (Italy) (photo taken by the authors).

Portions of dry and wet soil can also be considered as targets, as described in the following paragraph. In addition to these, as in the case of multispectral sensors, there are ground control points (GCP)s—points that are marked on the ground that have a known geographical position—which improve the positioning and accuracy of the mapping outputs.

To make GCPs more visible in thermal images, they should have a low emissivity when compared to the adjacent vegetation and other bodies [69].

In the case of thermal measurements, GCPs are made of aluminum, exploiting the low emissivity that makes it appear as a cold object in the images [70] (Figure 5b). This GCP was made using aluminum [71]. Our proposed GCPs were made using 50 cm × 50 cm white polypropylene panels and covering two quadrants while using aluminum sheets (Figure 5a). Black cardboard was used to

partially cover the other two quadrants to locate the point and also to make the GCP clearly visible and usable in multispectral surveys (Figure 5c). Figure 5 shows the detectability of the target in thermal and spectral images.

**Figure 5.** Example of a homemade target for ground control points (GCPs) designed to be easily detected in thermal as well as in multispectral UAV surveys. (**a**). The GCP target (white dashed circle) as it appears in (**b**) thermal and (**c**) near-infrared (NIR) multispectral images.

A further effect of the weather conditions to be taken into account is that of attenuating the thermal radiance by the atmosphere [59]. The sensor registers, for every pixel, an at-sensor radiance (Lat-sensor), being expressed in Wm<sup>−</sup>2, which is determined by the following formula (Equation (9)) [72]:

$$\mathbf{L}\_{\text{at-searson}} = \pi \, \mathbf{L}\_{\text{surf}} + \mathbf{L}\_{\text{atm}} \tag{9}$$

In Equation (9), τ represents the atmospheric transmittance, while Latm, measured in Wm<sup>−</sup>2, is the upwelling thermal radiation, created as a result of particles in the atmosphere, both depending on the distance of the sensor from the object and the water content in the atmosphere [59]. Parameters τ and Latm can be obtained by exploiting one of the theoretical atmospheric models, such as MODTRAN

(MODerate resolution atmospheric TRANsmission) [73], which are widely used for thermal data that are acquired from satellites [7]. This model allows for the estimation of atmospheric emission, thermal scattering, and solar scattering by incorporating the effects of molecular absorbers and scatterers, aerosols, and clouds, when considering the wavelengths from the ultraviolet region to the far-infrared [36].

An example of MODTRAN application in UAV is shown by [65], which shows that it allows for obtaining the surface temperature by applying atmospheric correction methods that are based on entering the model as input data, i.e., data that are related to local atmospheric conditions. In the application on UAV data using MODTRAN, it is essential to use local measurements of temperature, relative humidity, and atmospheric pressure acquired from nearby meteorological stations (Figure 6a) placed in the field [39–43].

**Figure 6.** (**a**) The weather station used during the thermal surveys on an onion field. (**b**) Handheld infrared thermometer used to measure the temperature of reference surfaces (photos taken by the authors).

#### *3.2. Data Processing*

Several steps are needed to acquire georeferenced UAV imagery, starting from the flight planning that generates a navigation file guiding the UAV to capture images with the required overlaps and geometric resolution automatically. Thermal imagery can be processed with structure-from-motion (SfM) algorithm, which does not always work properly [59], as reported in [74,75]. Indeed, SfM can be unable to align thermal imagery, making it necessary to mosaic separate images and manually geo-reference them using GCPs. These issues are due to the reduced information contained in the thermal image, which complicates the identification of the common features that are needed for bundle adjustment [59]. Indeed, if compared to other types of imagery, such as RGB, thermal imagery is characterized by a lower geometric resolution and contrast and by a stronger optical distortion [69]. In addition to increasing vertical and horizontal overlaps, the presence of an incorporated multispectral (or RGB) camera, which is characterized by a higher geometric resolution, could help the alignment step of the photogrammetric process. A framework to process thermal imagery is shown in Turner et al.'s [76]

study. Firstly, image pre-processing provided the removal of blurry imagery and the subsequent conversion of images to a 16-bit file format. In this way, all of the images have the same dynamic scale, and a specific temperature value corresponds to the same digital number (DN) in all images. Secondly, image alignment was executed using GPS log files and the imagery's time stamps. Finally, the spatial image was co-registered to RGB images. Some improvements of this framework were provided in [59] in order to optimize the alignment and processing of thermal images by exploiting the increased information contained in the RGB data. Although recent UAV models are provided with real-time kinematic global navigation satellite systems (RTK-GNSS) that are able to reach a centimeter planimetric accuracy, to ensure high positional precision, the placement of GCP within the surveyed site is generally expected. Imagery must be geometrically aligned (orthorectification), calibrated, and corrected, while taking atmospheric effects into account before final orthomosaics are obtained [7] (Figure 7).

**Figure 7.** RGB (top) and thermal orthomosaic (bottom) images of an onion field (data collected by the authors).

#### **4. Thermal UAV Imagery Analysis and Applications in Precision Agriculture (PA)**

#### *4.1. Crop water Stress Monitoring*

Several detection systems have been developed for the production of specialty crops worldwide thanks to technological advances [77]. The collection of accurate information on the spatial variability of fields is essential in this context. A field's variability is affected by several factors, including crop yield and water content, and sensors that can be used to detect these factors include thermal cameras.

Plant water stress is one of the main critical factors of abiotic stress, as it limits the development of crops [78,79]. Therefore, among the applications of thermal UAV RS in agriculture, the identification of water stress from plant temperature data is of great importance, with irrigation resource management being a key issue for PA. The use of UAVs in the study of plant water status requires measurements of stem water potential (ψs) and stomatal conductance (gs), which are the most common physiological water stress indicators that are used to determine crop water status [15,76,80,81]. Both indicators can be measured in the field while using a pressure chamber and porometers, respectively, although the interpolation of such local observations is not straightforward.

The use of thermal images allows for detecting the water stress conditions of a plant, because, at the foliar level, stomata closure reduces transpiration and evaporative cooling, which results in an increased temperature of the leaf [82,83]. This increase in temperature can be detected by thermal cameras; therefore, thermal images can provide spatially continuous information concerning the water status of plants in a wider area than that obtained by local measurements [84,85]. The temperature of the plant is not only regulated by the water supply, but also by the micrometeorological conditions [86]. Among the climatic factors, atmospheric humidity plays a crucial role and, in environments with humid climates, cloudiness also becomes a critical factor [87]. Several indices have been developed in the past decades with the aim of compensating for the variation of these conditions [86].

The CWSI was developed precisely for that purpose. Jackson and colleagues formulated the CWSI while using portable IR thermometers on herbaceous crops [84,88], developing a normalized index in order to overcome environmental effects that may affect the relationship between plant temperature and water stress. This index, which can assume values that range between 0 and 1, results in being directly proportional to the water stress level of many species of interest. The CWSI is based on the normalization of the canopy–air temperature difference with evaporative demand by means of the vapor pressure deficit (VPD) of the air [14]. The formula to calculate the CWSI is as follows, according to the methodology that was proposed by Idso et al. [17] (Equation (10)):

$$\text{CWSI} = \frac{(\text{T}\_{\text{C}} - \text{T}\_{a}) - (\text{T}\_{\text{C}} - \text{T}\_{a})\,\text{LL}}{(\text{T}\_{\text{C}} - \text{T}\_{a})\,\text{UL} - (\text{T}\_{\text{C}} - \text{T}\_{a})\,\text{LL}} \tag{10}$$

Where Tc − Ta is the canopy–air temperature difference, while LL refers to the Tc − Ta values for the lower limit and UL for the upper limit. The normalization related to the VPD considers the Tc − Ta difference of a canopy under two boundary conditions: (a) a lower limit when it transpires at its potential rate (i.e., under well-watered conditions), and (b) an upper limit under no transpiration. The lower limit is described by linear regression between Tc − Ta and the VPD, which is known as the non-water-stressed baseline (NWSB). The NWSB is empirically derived by measuring the leaf–air temperature difference for a well-watered crop in the experimental environment and provides the lowest temperature difference likely in that environment. Once the NWSB, (Tc − Ta)LL can be calculated by solving the baseline equation for the current VPD. The upper limit (Tc − Ta)UL, which is a constant, is obtained by solving the same NWSB equation for a hypothetical slightly negative VPD. It represents the vapor pressure difference that is generated by the temperature differential (Tc − Ta) when VPD is 0 [17]. The upper limit value is close to the NWSB interception *a* (depending on the temperature) and only corresponds to *a* when *a* is equal to 0 [89]. This method is site-dependent. A second approach for obtaining the upper and lower limits of (Tc − Ta) is theoretical and it foresees the combination of energy balance and diffusion equations. This method requires knowledge of challenging to obtain variables,

such as net radiation and aerodynamic resistance [19,89]. Another approach, defined statistical, was used with good results [15,63,90,91]. In brief, it foresees the use of the average temperature of the most cooling 5–10% of the canopy pixels as the wet reference to calculate the lower limit [92].

A further method that was proposed to obtain the two limits consists in the use of direct measurements over wet or dry reference surfaces, natural or artificial. The reference dry targets, as suggested by Jones [93,94], can be achieved, impeding leaf transpiration by covering the leaf surface (on one or both sides) with petroleum jelly. In this case, the temperature of the upper limit reference results from the measurement of the leaf temperature carried out about 30 minutes after the application of the petroleum jelly layer. The wet targets are obtained by spraying a thin layer of water on one or both leaf sides before taking the thermal image—this is valid on a single plant scale [94]. The advantage of using this method is that the stress levels are normalized to the actual response of the plants. However, the need to repeat the measurement for each test site after each flight of the UAV can be an obstacle to its applicability. In addition, another problem is related to the difficulty of identifying the leaves that are covered by petroleum jelly [95,96]. Maes proposed an example of a wet artificial reference target in [97] using a cloth and a steel wireframe. The target thus created, with an appearance, shape, and size similar to kiwi and grapevines leaves was kept wet for days by keeping the lower part of the cloth immersed inside a bottle filled with water. Artificial targets were also used in [98]. This study provides an approach in which the reference temperature of the upper limit is set at 5 ◦C above the air temperature, while that of the lower limit is derived from the temperature of the artificial target. Berni et al. [19] proposed an approach to monitor relatively large areas while using UAVs to detect water stress; the work, focused on the control of water stress in an olive grove, based on physical models for the estimation of input variables of energy balance equations, did not require the use of reference surfaces.

Two problems have prevented the widespread use of the CWSI: the first concerns its use under changing atmospheric conditions, and the second occurs when using lower resolution data from satellites or aircrafts platforms (as compared to UAVs), concerning, namely, the problem of "mixed-pixel value", which is where part of the temperature of the pixel comes from the background soil and not from the pure canopy, thus reducing the quality of the data [99].

As far as the first aspect is concerned, CWSI was found to work better in dry climates, while it showed considerable limitations in wet climates and in environments with substantial climatic variability [100]. In fact, it must be taken into account that the absolute in leaf-air temperature difference decreases as the atmospheric humidity deficit decreases, and the same goes for sensitivity to any measurements made. Furthermore, while taking the latter aspect into account in the calculation of the CWSI, as the humidity (and temperature) deficit decreases, the signal-to-noise ratio is reduced [101]. Finally, it must also be considered that the canopy temperature depends, in part, on variations in the roughness of the canopy, wind speed, and net radiation, all of which are more variable in humid climates [87].

As far as the problem of "mixed-pixel value" is concerned, while the resolution of satellite images is limited, the sensors mounted on UAVs, because of their sub-meter spatial resolution, make it easier to recover the pure canopy temperature while minimizing the thermal effects of the background soil [65,102,103]. Different approaches can be adopted, as shown in [99]. An approach to extracting the temperatures of canopies' sunlit leaves foresees distinguishing canopies' pixels from background pixels that are based on an analysis of the temperature distribution among the pixel population and then segment the image into two distinct classes. This solution is possible where there is a clear difference between the average temperature of the canopy and that of the soil/background. A valid solution also includes the use of additional information that is derived from RGB, multispectral and hyperspectral images (in cases where several sensors are used simultaneously), in order to identify plant's pixels [99]. In this case, the use of vegetation indices, like NDVI, to separate soil pixels and plant pixels can be useful.

In addition to the CWSI, several normalized thermal indices were developed between the 1970s and 1980s. Among these, there are the conductance index (IG) and the stomatal conductance index (I3), whose formulas (Equations (11) and (12)) are as follows:

$$IG = \frac{\mathcal{T}\_{dry} - \mathcal{T}\_{\mathcal{C}}}{\mathcal{T}\_{\mathcal{C}} - \mathcal{T}\_{\text{wet}}} \tag{11}$$

$$I3 = \frac{\mathcal{T}\_{\mathcal{C}} - \mathcal{T}\_{\text{wet}}}{\mathcal{T}\_{dry} - \mathcal{T}\_{\mathcal{C}}} \tag{12}$$

where Tc represents the surface temperature of the canopy, while Tdry and Twet are entirely dry or wet reference surfaces to simulate leaf temperature under conditions of minimum and maximum transpiration.

The main aspects of these two indices are that IG increases with stomatal conductance, while I3 is positively correlated with stomatal resistance [94,96]. The three indices (CWSI, IG, and I3) need the knowledge of canopy temperature under both wet and dry conditions; however, the concept of the CWSI remains the most widely used, as it is best known [94].

Research carried out thus far has concerned the monitoring of water stress in both herbaceous and tree crops. In Sullivan et al.'s [24] and Bian et al.'s [104] research, thermal UAVs were used, respectively, to monitor *Gossypium hirsutum* L. crop residue management and the response to different irrigation treatments, calculating the CWSI. Martinez et al. [105] worked on monitoring sugar beet water stress by comparing thermal data, which were obtained using a conventional thermal camera, with those obtained using a low-cost infrared sensor.

Zhang et al. [106] monitored, at the farm scale, maize water stress using RGB and thermal images. The study of Crusiol et al. [107] dealt with the evaluation of the water status of soybean plants under different water conditions. Other works [108,109] focused on the monitoring of herbaceous crops while using different sensors (RGB, multispectral, and thermal) with the aim of producing orthomosaics and vigor maps.

The CWSI is also widely used for perennial crops. Bellvert et al. [110] calculated the index in peach trees (*Prunus Persica*), mapping the internal spatial variability of the orchard by thermal UAV and verifying the relationships between the index and the leaf water potential in different growing seasons. Gonzalez-Dugo et al. [91] studied the spatial variations in the water status of five different fruit tree species: almond (*Prunus dulcis*), apricot (*Prunus armeniaca*), peach (*Prunus persica*), lemon (*Citrus x limon*), and orange (*Citrus sinensis*). In this respect, it is important to remember that the stomatal response to environmental conditions can vary from species to species [111]. Olive and citrus show relevant stomatal closure at midday [89,91], while, in other species, such as almond trees, the stomatal behavior changes between the different cultivars [112].

Some studies have applied thermal UAV RS to citrus orchards (orange, *Citrus sinensis*, and mandarin, *Citrus reticulata*) in order to extract the temperature of the crowns for water stress detection [113,114]. Research has also been carried out on olive crops, a species of great importance in the countries of the Mediterranean basin [115], in order to verify the plant's behavior in response to various irrigation treatments. In particular, Poblete-Echeverría et al. [116] showed that the temperature difference between the canopy (Tc) and air (Ta) is related to the difference in water potential under different irrigation treatments when the plants are under water stress conditions. Berni et al. [19] used very high-resolution UAV thermal imagery to calculate and map the tree canopy conductance and the CWSI in a heterogeneous olive orchard.

Egea et al. [117] demonstrated the usefulness of the CWSI for monitoring water stress in a dense olive orchard by verifying sound relationships between the index and water stress indicators, such as stomatal conductance, stem water potential, and leaf transpiration rate. Ortega-Farías et al. [118], while using thermal and multispectral cameras that were mounted on UAVs, estimated the energy balance components on a drip irrigation olive grove, acquiring high-resolution images to evaluate intra-field spatial variability.

Several works have used thermal UAVs for the calculation of the CWSI in vineyards, such as that of Zarco-Tejada et al. [81]. In Baluja et al.'s [15] study, the relationships between the temperatures or indices that were derived from thermal and multispectral imagery and stomatal conductance and water potential were determined. In particular, air and leaf temperatures were recorded with a handheld thermometer, and stomatal conductance was measured with a leaf gas exchange system, while the stem water potential was measured using a Scholander pressure bomb; the CWSI, IG, and I3 indices were also calculated.

Similarly, the work by Bellvert et al. [20] related the CWSI with leaf water potential and that by Santesteban et al. [119] related the CWSI with stem water potential and seasonal leaf stomatal conductance to estimate the variability of plant water status in a vineyard. Matese et al. [120] and Pàdua et al. [16,121] exploited the use of different sensors (RGB, multispectral, and thermal) for several applications in precision viticulture: the production of vigor maps, multi-temporal analysis of vigor maps, and water stress detection. Gómez-Candón et al. [25], by testing in an apple orchard, proposed a methodology to derive thermal orthomosaics, including a method for the radiometric correction of thermal UAV imagery.

#### *4.2. Plant Disease Detection, Phenotyping, Yield Estimation, and Vegetation Status Monitoring*

Calderón et al. provides an interesting example of application in the framework of plant pathology of thermal UAV RS [14]. In particular, this research aimed to evaluate the use of thermal imagery and physiological indices derived from other types of sensors to detect the presence of infection that is caused by the soil-borne fungus *Verticillium dahliae*. The role of thermal RS in the diagnosis of the pathology is due to the water stress of plants, caused by the fungus, or by the stomatal closure, which determines the reduction of the transpiration rate. Therefore, decreasing evaporative cooling increases leaf temperature. By conducting field measurements, it was demonstrated that crown variations of temperature are higher and stomatal conductance lower as the severity level of the disease increased. Besides, crown temperature and the CWSI index were shown to be among the best indicators to detect *Verticillium dahliae* in the early stages of disease development.

With regard to the issue of yield estimation and related aspects, thermal UAV RS was exploited to research cotton, soybean, and rice crops [23,122,123]. Feng et al. [122] used three types of sensors (RGB, multispectral, and thermal) in order to analyze several features to see which one had the best performance for the yield prediction and at which stages of the crop cycle: Four vegetation indices, canopy cover, plant height, temperature, and a cotton fiber index.

The study of Maimaitijiang et al. [123] on soybean concerned testing the potentialities of using different sensors (including the thermal camera) in the framework of multimodal data fusion and deep neural network to derive useful information for the yield prediction model. Liu et al. [23] investigated rice lodging while using RGB and thermal sensor UAV images. The objectives of the research were to measure the daily temperature differences between non-lodged and lodged rice crops, as well as to detect the optimal time window. Because the differences in temperature between 10 a.m. and 4 p.m. were significant, thermal camera devices allowed for identifying lodged rice plants.

Some other works have focused on the monitoring of vineyards and related landscape elements, such as agricultural terraces [124–127]. In particular, Tucci et al. [125,126] used RGB and thermal cameras to investigate the thermal dynamics of a vineyard that was grown on dry-stone wall terraced land. These dry-stone terraces are part of the UNESCO's Representative List of the Intangible Cultural Heritage of Humanity and are typical and iconic features of the agricultural landscapes across Europe, including Italy [128]. However, the risk of abandonment and degradation threaten agricultural terraces, which can involve the increase of hazards that are linked to geo-hydrological processes triggered by rainfall events. In this framework, the use of RGB and thermal sensors mounted on UAVs represents an efficient and cost-effective monitoring methodology, also given the high resolution of the images

and the reduced time for their acquisition. Tucci et al. [125], besides verifying the effectiveness of the used sensors in monitoring the terraced crops, investigated the influence that dry-stone walls can have on the microclimate of the vineyard and, consequently, on the quality of the obtained grapes. The results of the study showed differences in the temperature between the plants of the internal rows (lower temperatures) and the external rows (higher temperatures).

Phenotyping has an important role in crop science research. The acquisition of crop phenotypic information in different environments allows the association of genomic and phenotypic information useful for increasing yield [129]. In this framework, recently, the importance of the use of UAVs has increased, as they provide a rapid and non-destructive approach to phenotyping, also allowing the use of high spatial resolution images [129,130].

The use of UAVs has made it possible to overcome various limits that are linked to the use of traditional methods, such as the difficulty in making simultaneous measurements on different plots [131]. In this context, some works have been carried out involving the use of thermal UAV RS, thus proving its effectiveness, as in [9]. Natarajan et al. [132] employed different sensors (visible, multispectral, and thermal) that were mounted on UAVs for the phenotyping of indirect traits (including canopy temperature) for early-stage selection in sugarcane breeding.

Gracia-Romero et al.'s [133] work concerned the comparison of the performance of RGB, multispectral, and thermal data, which were derived from ground-based and UAV surveys, with the aim to assess genotypic differences in durum wheat's yield, under different growing conditions. In this case, the measurement of the canopies' temperature was an alternative valid to that of leaf stomatal conductance. Perich et al. [134] exploited thermal UAV RS to measure canopy temperature in wheat. Canopy temperature has a strong association with water status and stomatal conductance in wheat [135–137]; in particular, low canopy temperature can be associated with a 30% increase in yield, as well as an increase in water absorption by the deeper roots when measured during the grain filling phase [138]. Maimaitijiang et al.'s [139] study, which was performed on soybean, showed that the fusion of thermal and multispectral data could provide the best estimate of the biochemical traits of crops (chlorophyll content and N concentration) and biophysics (LAI, fresh and dry biomass).

The use of thermal UAV RS has proven to be very promising as a less-expensive way for mapping drainage systems. Subsurface tile drainage is of great importance in the Midwest of the United States [140]. Subsoil drainage allows for the level of groundwater to be lowered by removing excess water, as well as reducing soil erosion and increasing the aeration [141,142] and infiltration capacity of water derived from precipitation.

Knowing the precise layout of the drainage networks, as well as bringing benefits for their better functionality, favors soil and water conservation practices [143]. Additional benefits include increases in the soil productivity and yield of the crops, most of which do not tolerate excess water [143–145]. As the soil is directly above a drainage line, it is often drier than the soil between the drainage lines [146]. Differences in emissivity between dry and wet soil surfaces can be detected while using a thermal camera [147], which has been demonstrated in several studies using thermal cameras mounted on UAVs, even in combination with visible and multispectral cameras, showing the significant performance of thermal images and potentialities for their use in mapping agricultural drainage pipe systems [143,146,148–150].

Table 2 reports the dataset of all analyzed studies concerning the application of thermal UAV RS in precision agriculture (PA) and organized based on the aim of the study, the type of the used camera sensor, and the analyzed crops.

**Table 2.** Dataset compilation of studies concerning the application of thermal UAV remote sensing (RS) in precision agriculture (PA).

