Multi-Sensor UAV Application for Thermal Analysis on a Dry-Stone Terraced Vineyard in Rural Tuscany Landscape

Italian dry-stone wall terracing represents one of the most iconic features of agricultural landscapes across Europe, with sites listed among UNESCO World Heritage Sites and FAO Globally Important Agricultural Heritage Systems (GIAHS). The analysis of microclimate modifications induced by alterations of hillslope and by dry-stone walls is of particular interest for the valuation of benefits and drawbacks of terraces cultivation, a global land management technique. The aim of this paper is to perform a thermal characterization of a dry-stone wall terraced vineyard in the Chianti area (Tuscany, Italy), to detect possible microclimate dynamics induced by dry-stone terracing. The aerial surveys were carried out by using two sensors, in the Visible (VIS) and Thermal InfraRed (TIR) spectral range, mounted on Unmanned Aerial Vehicles (UAVs), with two different flights. Our results reveal that, in the morning, vineyard rows close to dry-stone walls have statistically lower temperatures with respect to the external ones. In the afternoon, due to solar insulation, temperatures raised to the same value for each row. The results of this early study, jointly with the latest developments in UAV and sensor technologies, justify and encourage further analyses on local climatic modifications in terraced landscapes.


Introduction
Dry-stone wall terracing represents one of the most iconic features of agricultural landscapes across Europe, as well as in other rural regions of the world.Italian and other dry-stone walls were added in November 2018 to UNESCO's Representative List of the Intangible Cultural Heritage of Humanity (Decision 13.COM 10.b.10).Dry stone walling was recognized to be still in use in several countries, playing a fundamental role in maintaining the environment and landscape.
In addition, many terraced landscapes of the world are already inserted among UNESCO World Heritage Sites list, as well as in the newly established FAO Globally Important Agricultural Heritage Systems (GIAHS), including cases of Italian dry-stone wall terracing, such as "Portovenere, Cinque Terre, and the Islands terraced vineyards", (UNESCO Heritage Site), or the "Olive Groves of the Slopes between Assisi and Spoleto" (FAO GIAHS).
In Italy, and in particular in the Tuscany region, rural landscape is characterized by the presence of stone wall bench terraces [1].It is evident how terracing can alter the physical landscape, influencing factors such as soil moisture, erosion control and soil conservation [2].The different ecosystem services provided by agricultural terraces have been object of a large number of studies in recent years [2][3][4].These structures can induce significant modifications to local microclimate, such as a shading, as well as effects on local winds, inducing a potential thermal effect.
Due to the growing availability of short-(cameras) to long-range (satellite) sensors, remote sensing-based methodologies have been applied to the analysis of agricultural terraced fields, focusing on geo-morphological aspects, including analysis of preferential drainage patterns [5], and large scale [6] to medium scale terraces sites mapping [7].Furthermore, remote sensing methodologies can be used to detect relevant bio-geo physical phenomena in these unique agricultural landscapes, such as microclimatic effects induced by the geo-morphological shape of bench terraces, as well as the presence of anthropic artifacts, such as dry-stone walls.
Nowadays, Unmanned Aerial Vehicle (UAV) systems can compete with traditional remote-sensing platforms (i.e., aircrafts and satellites) for agricultural mapping, thanks to their high flexibility of use, low operational cost and very high spatial resolution (cm-level).In Matese et al. [8] a comparison between vegetation indices obtained by satellite, aircraft and UAV platforms for intra-vineyard variability monitoring, highlights pros and cons of the different technologies: (i) satellite survey allows to simultaneously cover a wider area but is affected by some limitations, such as low spatial resolution for precision agriculture purposes, meteorological-and orbit-dependent acquisition, high economic costs; (ii) aircraft surveys are characterized by higher planning flexibility although noticeable organizational efforts and costs needs to be considered; (iii) UAVs are suitable systems for small field applications, but are still limited on wider areas because of their short flight endurance and limited payload for sensor carriage.
UAV photogrammetry has been extensively used for mapping and 3D modeling in the geomatics field, as evidenced by the huge increase in this research field in the last few years [9][10][11][12].These platforms, equipped with digital RGB cameras, represent now a common alternative to traditional aerial manned photogrammetry.The application fields of image-based photogrammetric products from UAVs, such as 3D models or orthoimages, relate to forestry and agriculture [13], archaeology [14], architecture [15], cultural heritage [16][17][18], environmental surveying [19], traffic monitoring [20], and many others.Considering the above-mentioned features, and the possibility of carrying different type of payloads, ranging from RGB to Thermal InfraRed (TIR) sensors, UAVs appears to represent a key technology allowing a deeper understanding of small-to medium-scale bio-geophysical processes in agricultural landscapes.
The acquisition of the emitted radiation from natural surfaces in the range of TIR, thanks to airborne mounted sensors, allows us to obtain useful indication about soil and crops status by the analysis of temperature distribution across the field.For the specific case of vineyards, temperature patterns play a crucial role in regulating plant physiology, influencing phenology [21] and photosynthesis [22].Furthermore, temperature trends were found to significantly impact on the berry developmental processes, such as the accumulation of dry matter and sugar [23], or the accumulation of anthocyanins [24], thus potentially determining the quality of the final product, wine.Costa et al. [25] reported a research on the evaluation of canopy temperature to generate maps for irrigation scheduling; other thermal-related indices were proposed by Idso et al. [26] and Jones [27,28] for the evaluation of crop water stress and stomatal conductance, respectively.Applications of thermal imaging on vineyard has been tested by García-Tejero et al. [29], Grant et al. [30], Jones et al. [28,31] and Moller et al. [32], who researched the correlation between thermal indices and different water regimes of vineyards.Further tests with thermal camera mounted on UAV systems are reported in Baluja et al. [33], for the evaluation of water status variability; Berni et al. [34] estimated biophysical parameters and water stress from thermal imagery; Turner et al. [35] presented a UAV system equipped with a thermal camera for soil moisture assessment to improve irrigation efficiency; Bellvert et al. [36] reported the evaluation of crop water stress index (CWSI) for precision irrigation in a Spanish vineyard.
To our knowledge, few, early experiences of thermal monitoring of dry-stone agricultural terraces have been carried out.Warren et al. [37] tested the use of TIR imagery for the detection of dry-stone walls features.They showed how thermal imagery can be used to detect wall features such as the depth of backfill, areas of high moisture and the presence of particular structural elements not visible from conventional surveys.Barbera et al. [38] analyzed the microclimatic modifications induced by dry stone walls in particular structures called "jardinu" (circular dry-stone walls protecting single or groups of trees), located in the Pantelleria Island terraced fields.They found out that the thermal inertia of the stones may greatly enhance dewfall occurrence, both on stones and on leaves, suggesting a thermal gradient between the stones and the nearby environment.Moreover, they observed that the dissipation of the heat cumulated during the day by the stones reduced temperature excursion at night, positively affecting citrus adaptation.The importance of microclimatic monitoring through remotely sensed methods for a more reliable and advanced management of eco-, as well as agricultural systems, is advancing in parallel to technological innovations of sensors and UAV platforms, as pointed out by Zellweger et al. [39].For instance, the recent study of Romboli et al. [40] showed how microclimate spatial variability interacts with vine vigor in determining grape and wine characteristics.This latter one, however, was carried out in one non-terraced vineyard, with no anthropic modification to the surrounding landscape apart from cultivations.From the perspective of better understanding all the implications connected to dry-stone walls traditional structures over microclimate variations, further investigations are necessary.
Taking advantage by recent progresses in UAVs technologies, this paper presents a preliminary test carried out using a TIR sensor mounted on a UAV, coupled with a common RGB survey, to detect the thermal dynamics of a vineyard grown on a dry-stone wall terraced land, with the specific purposes of:

•
Investigate strengths and weaknesses of the use of TIR sensor mounted on UAV for thermal analysis of terraced crops; • conduct a preliminary test on the possible thermal effect that dry-stone walls can have on the vineyard microclimate, testing the hypothesis that stones have an influence on the temperature patterns of the field which can influence grape ripening and quality.

Study Area
The terraced area of Lamole, within the municipality of Greve in Chianti (Florence, Italy; 43 •  Despite not yet being present in FAO GIAHS and UNESCO Heritage lists, the site was included in the Rural Historic Landscape Catalogue [41] of the Italian Ministry of Agricultural, Food, Forestry Policies in February 2018.As reported in the related Ministerial Decree (n.6415 of 20 February 2018), the Lamole area was included in the Catalogue because it represents "an historical landscape where soil protection measures, represented by dry-stone wall terraces, together with policoltures represented by olive groves and vineyards, the woodland and the historical settlements preserve the Chianti cultural identity.Moreover, Lamole represents a remarkable example of dynamic recover and conservation of historical agricultural practices, which play nowadays a peerless added value for the development of the entire territory." The study has been carried out in the Grospoli terraced vineyard, a site that has been object of a significant number of studies related to hydrological processes and failures on terraces (among others, Preti et al. [42,43]).The area lies in the Chianti Classico production region.Sangiovese, Petit Verdot, Cabernet Sauvignon cultivars are grown in the vineyard (Figure 1).The vineyard has an area of 1.76 ha, with slopes ranging from from 12 to 22% [43] and an average altitude of 545 m a.s.l.The average terrace width is around 11 m, with wall ranging from 1 to 2.7 m height.For the remote sensing analysis, rows of vines were divided in two groups, the North and the South part, that are shown in Figure 2.  The vineyard has an area of 1.76 ha, with slopes ranging from from 12 to 22% [43] and an average altitude of 545 m a.s.l.The average terrace width is around 11 m, with wall ranging from 1 to 2.7 m height.For the remote sensing analysis, rows of vines were divided in two groups, the North and the South part, that are shown in Figure 2. The vineyard has an area of 1.76 ha, with slopes ranging from from 12 to 22% [43] and an average altitude of 545 m a.s.l.The average terrace width is around 11 m, with wall ranging from 1 to 2.7 m height.For the remote sensing analysis, rows of vines were divided in two groups, the North and the South part, that are shown in Figure 2.

Materials and Methods
An aerial survey campaign was performed on 8 September 2017 over the Grospoli terraced vineyard in Lamole, a few weeks before grape harvesting, by employing two different UAVs, equipped with different sensors.The operations were carried out by the research group of the Survey, Cultural Heritage, Engineering, Monitoring, Analysis (SCHEMA) Laboratory, a joint collaboration between the University of Florence and the Italian Military Geographic Institute (IGMI).The first photogrammetric test was performed with a classic aerial approach in the visible (VIS) range, to produce a Digital Elevation Model (DEM) of the area and the respective projected orthophoto.The second test was performed in the TIR range in order to study the influence of stone walls upon the thermal properties of the neighboring rows.The TIR campaign was originally designed by planning three flights, in the morning, half day and evening, to monitor the thermal range over the entire day.Two flights were finally performed, due to the adverse climate conditions, in the morning (08:50 CET) and at 15:00 CET in the afternoon.Each flight was performed in terrain-following mode, which ensures the correct value of overlap and Ground Sample Distance (GSD) on uneven ground is maintained.The surveys with UAVs respected the regulation of the Italian Civil Aviation Authority (Ente Nazionale Aviazione Civile) in terms of non-critical operations, thus in Visual Line of Sight (VLoS).

VIS Data Capture and Processing
A single flight was performed with a DJI Phantom 4 Pro over the entire area of the Grospoli vineyard, for the photogrammetric data capture in the VIS range.The used UAV platform is a multirotor quadcopter, which can be remotely controlled and programmed for automatic navigation, provided by Global Positioning System (GPS) waypoints.The aircraft is characterized by limited weight of about 1400 g (battery, propellers and sensor included) and size (350 mm diagonal size), with max flight time of approximately 30 minutes and max speed of 72 km/h.The system supports two operating frequencies, at 2.4 GHz and 5.8 GHz and the maximum transmission range is of 3.5 km (according to CE standard).The platform is also equipped with a sensor system, with dual rear front optical sensors and side infrared sensors, which allow a total of 5 directions obstacle sensing and 4 direction obstacle avoidance.Dual Inertial Measurement Units (IMU) and compasses design provides redundancy in case of errors.Other technical characteristics are reported in Table 1.

Materials and Methods
An aerial survey campaign was performed on 8 September 2017 over the Grospoli terraced vineyard in Lamole, a few weeks before grape harvesting, by employing two different UAVs, equipped with different sensors.The operations were carried out by the research group of the Survey, Cultural Heritage, Engineering, Monitoring, Analysis (SCHEMA) Laboratory, a joint collaboration between the University of Florence and the Italian Military Geographic Institute (IGMI).The first photogrammetric test was performed with a classic aerial approach in the visible (VIS) range, to produce a Digital Elevation Model (DEM) of the area and the respective projected orthophoto.The second test was performed in the TIR range in order to study the influence of stone walls upon the thermal properties of the neighboring rows.The TIR campaign was originally designed by planning three flights, in the morning, half day and evening, to monitor the thermal range over the entire day.Two flights were finally performed, due to the adverse climate conditions, in the morning (08:50 CET) and at 15:00 CET in the afternoon.Each flight was performed in terrain-following mode, which ensures the correct value of overlap and Ground Sample Distance (GSD) on uneven ground is maintained.The surveys with UAVs respected the regulation of the Italian Civil Aviation Authority (Ente Nazionale Aviazione Civile) in terms of non-critical operations, thus in Visual Line of Sight (VLoS).

VIS Data Capture and Processing
A single flight was performed with a DJI Phantom 4 Pro over the entire area of the Grospoli vineyard, for the photogrammetric data capture in the VIS range.The used UAV platform is a multirotor quadcopter, which can be remotely controlled and programmed for automatic navigation, provided by Global Positioning System (GPS) waypoints.The aircraft is characterized by limited weight of about 1400 g (battery, propellers and sensor included) and size (350 mm diagonal size), with max flight time of approximately 30 minutes and max speed of 72 km/h.The system supports two operating frequencies, at 2.4 GHz and 5.8 GHz and the maximum transmission range is of 3.5 km (according to CE standard).The platform is also equipped with a sensor system, with dual rear front optical sensors and side infrared sensors, which allow a total of 5 directions obstacle sensing and 4 direction obstacle avoidance.Dual Inertial Measurement Units (IMU) and compasses design provides redundancy in case of errors.Other technical characteristics are reported in Table 1.The DJI Phantom 4 Pro is equipped with DJI FC6310 integrated digital RGB camera, stabilized with a 3-axis (pitch, roll, yaw) gimbal.The onboard camera uses a 1-inch CMOS 20 MP sensor, which can take images of 5272 × 3648 pixels as maximum resolution.The fixed focal length is 8.8 mm (35 mm format equivalent of 24 mm) with 84 • FOV and f/2.8-f/11.The mechanical shutter, with 1/2000 s max speed, as well as the lens arrangement allows to obtain good image quality in terms of both sharpness and vividness.Table 2 lists the relevant technical properties of the RGB sensor (left column).The DJI Phantom 4 Pro is equipped with DJI FC6310 integrated digital RGB camera, stabilized with a 3-axis (pitch, roll, yaw) gimbal.The onboard camera uses a 1-inch CMOS 20 MP sensor, which can take images of 5272 x 3648 pixels as maximum resolution.The fixed focal length is 8.8 mm (35 mm format equivalent of 24 mm) with 84° FOV and f/2.8-f/11.The mechanical shutter, with 1/2000 s max speed, as well as the lens arrangement allows to obtain good image quality in terms of both sharpness and vividness.Table 2 lists the relevant technical properties of the RGB sensor (left column).Acronyms explication: 1 Field Of View (FOV), 2 Red, Green, Blue, 3 Joint Photographic Experts Group (JPEG), 4 Thermal InfraRed (TIR), 5 Focal Plane Array (FPA), 6 RAdiometric VIdeo sequence (RAVI).
The flight plan was realized with UgCS software (SPH Engineering UgCS), by setting the suitable parameters in order to have a GSD of approximately 2 cm.The flight was designed according to a classic photogrammetric criterion, with a flight pattern consisting of 9 swipes, with forward overlapping of 80% and sidelap of 70%.The images were acquired with nadiral camera at a constant speed of 5 m/s and altitude of 70 m Above Ground Level (AGL).More details of the flight plan are reported in Table 3 (left column).
The georeferencing of the photogrammetric survey was made by using 10 Ground Control Points (GCP), which guarantee also the metric accuracy of the survey.The GCPs were materialized as targets functional for both visible and thermal survey: 60 x 60 cm 2 panels covered with aluminium foil and marked with black sprayed triangles.These targets were homogeneously distributed around the surveyed vineyard, to cover all the interested surface.The positioning and measuring activities for each target, were carried out by the Geodetic Direction of the IGMI, with the Global Navigation Satellite System (GNSS) TOPCON GR3.Two multi-frequency receivers were used for the coordinates acquisition, one as Master (base receiver) and the other as Rover receiver, working in Real Time Kinematic (RTK) radio mode.The coordinate system used in all data processing is ETRS89/UTM32N (EPSG:25832).The GCPs accuracy is about 1 cm.
The 206 images acquired for the VIS flight were processed with the photogrammetric software Agisoft Photoscan (version 1.2.6 build 2834).After image alignment, thanks to common Structure from Motion (SfM) procedures, GCPs were collimated and associated to their GNSS coordinates, for georeferencing and scale the photogrammetric model.A set of 7 of the 10 total collimated targets were used as GCPs, while 3 of them were used as Check Points (CP).Multi-View Stereo (MVS) The DJI Phantom 4 Pro is equipped with DJI FC6310 integrated digital RGB camera, stabilized with a 3-axis (pitch, roll, yaw) gimbal.The onboard camera uses a 1-inch CMOS 20 MP sensor, which can take images of 5272 x 3648 pixels as maximum resolution.The fixed focal length is 8.8 mm (35 mm format equivalent of 24 mm) with 84° FOV and f/2.8-f/11.The mechanical shutter, with 1/2000 s max speed, as well as the lens arrangement allows to obtain good image quality in terms of both sharpness and vividness.Acronyms explication: 1 Field Of View (FOV), 2 Red, Green, Blue, 3 Joint Photographic Experts Group (JPEG), 4 Thermal InfraRed (TIR), 5 Focal Plane Array (FPA), 6 RAdiometric VIdeo sequence (RAVI).
The flight plan was realized with UgCS software (SPH Engineering UgCS), by setting the suitable parameters in order to have a GSD of approximately 2 cm.The flight was designed according to a classic photogrammetric criterion, with a flight pattern consisting of 9 swipes, with forward overlapping of 80% and sidelap of 70%.The images were acquired with nadiral camera at a constant speed of 5 m/s and altitude of 70 m Above Ground Level (AGL).More details of the flight plan are reported in Table 3 (left column).
The georeferencing of the photogrammetric survey was made by using 10 Ground Control Points (GCP), which guarantee also the metric accuracy of the survey.The GCPs were materialized as targets functional for both visible and thermal survey: 60 x 60 cm 2 panels covered with aluminium foil and marked with black sprayed triangles.These targets were homogeneously distributed around the surveyed vineyard, to cover all the interested surface.The positioning and measuring activities for each target, were carried out by the Geodetic Direction of the IGMI, with the Global Navigation Satellite System (GNSS) TOPCON GR3.Two multi-frequency receivers were used for the coordinates acquisition, one as Master (base receiver) and the other as Rover receiver, working in Real Time Kinematic (RTK) radio mode.The coordinate system used in all data processing is ETRS89/UTM32N (EPSG:25832).The GCPs accuracy is about 1 cm.
The 206 images acquired for the VIS flight were processed with the photogrammetric software Agisoft Photoscan (version 1.2.6 build 2834).After image alignment, thanks to common Structure from Motion (SfM) procedures, GCPs were collimated and associated to their GNSS coordinates, for georeferencing and scale the photogrammetric model.A set of 7 of the 10 total collimated targets were used as GCPs, while 3 of them were used as Check Points (CP).Multi-View Stereo (MVS) Acronyms explication: 1 Field Of View (FOV), 2 Red, Green, Blue, 3 Joint Photographic Experts Group (JPEG), 4 Thermal InfraRed (TIR), 5 Focal Plane Array (FPA), 6 RAdiometric VIdeo sequence (RAVI).
The flight plan was realized with UgCS software (SPH Engineering UgCS), by setting the suitable parameters in order to have a GSD of approximately 2 cm.The flight was designed according to a classic photogrammetric criterion, with a flight pattern consisting of 9 swipes, with forward overlapping of 80% and sidelap of 70%.The images were acquired with nadiral camera at a constant speed of 5 m/s and altitude of 70 m Above Ground Level (AGL).More details of the flight plan are reported in Table 3 (left column).
The georeferencing of the photogrammetric survey was made by using 10 Ground Control Points (GCP), which guarantee also the metric accuracy of the survey.The GCPs were materialized as targets functional for both visible and thermal survey: 60 × 60 cm 2 panels covered with aluminium foil and marked with black sprayed triangles.These targets were homogeneously distributed around the surveyed vineyard, to cover all the interested surface.The positioning and measuring activities for each target, were carried out by the Geodetic Direction of the IGMI, with the Global Navigation Satellite System (GNSS) TOPCON GR3.Two multi-frequency receivers were used for the coordinates acquisition, one as Master (base receiver) and the other as Rover receiver, working in Real Time Kinematic (RTK) radio mode.The coordinate system used in all data processing is ETRS89/UTM32N (EPSG:25832).The GCPs accuracy is about 1 cm.
The 206 images acquired for the VIS flight were processed with the photogrammetric software Agisoft Photoscan (version 1.2.6 build 2834).After image alignment, thanks to common Structure from Motion (SfM) procedures, GCPs were collimated and associated to their GNSS coordinates, for georeferencing and scale the photogrammetric model.A set of 7 of the 10 total collimated targets were used as GCPs, while 3 of them were used as Check Points (CP).Multi-View Stereo (MVS) algorithms were then used to obtain the dense cloud, from where the DEM of the surveyed area was generated.Finally, an orthophoto was projected as texture on the 3D model.

Flight Plans
algorithms were then used to obtain the dense cloud, from where the DEM of the surveyed area was generated.Finally, an orthophoto was projected as texture on the 3D model.

TIR Data Capture and Processing
The flight campaign for the thermal data capture was carried out on the same day with two distinct flights, at 08.50 CET and at 15:00 CET, by using an hexacopter customized by MicroGeo s.r.l. as UAV.
The flight plan was realized with the open source Mission Planner software, by setting the suitable parameters in order to have a GSD of 11 cm.The flights were designed according to a classic photogrammetric criterion, with a flight pattern consisting of 12 swipes, with forward overlapping of 80% and sidelap of 63%.The radiometric video sequences were acquired with nadiral orientation at a constant speed of 3 m/s and altitude of 40 m AGL.More details of the flight plan are reported in Table 3 (right column).
The UAV platform mounted the thermal camera OPTRIS PI450, with maximum optical resolution of 382 x 288 pixels at 80 Hz frame rate.The fixed focal length is 8 mm, with an angular FOV of 62° x 49°.The focal plane array (FPA) image sensor is made of uncooled microbolometers (25 x 25 μm each), with a spectral TIR response in the range 7.5-13 μm and accuracy of ± 2 °C.Table 2 (right column) shows a summary of thermal camera specifications.The camera can be connected to a computer via a USB 2.0 protocol to download the acquired data in radiometric video sequences format (.RAVI).Single frames had consequently to be extracted from the video to be suitable for a classical photogrammetric project.
The same targets used for RGB photogrammetry were used for georeferencing the two thermal orthomosaics: since aluminium foil has a very low emissivity (0.09 ε) it can be easily seen in the TIR range.Some other panels were painted in black and white, onto the aluminium foil, and distributed all around the test field, as delimited test areas for pre-flight calibration.In fact, before the TIR flights, a calibration with the thermal sensor was made by detecting the thermal response over the GCP targets and the b/w panels and by comparing them with temperature measurements taken with contact thermometers.
The final values of ambient temperature and emissivity were kept as default, Tamb 23 °C and ε = 1 (vegetation usually has high values of emissivity, around 0.98 [34]) during the data processing step.

VIS Range
algorithms were then used to obtain the dense cloud, from where the DEM of the surveyed area was generated.Finally, an orthophoto was projected as texture on the 3D model.

TIR Data Capture and Processing
The flight campaign for the thermal data capture was carried out on the same day with two distinct flights, at 08.50 CET and at 15:00 CET, by using an hexacopter customized by MicroGeo s.r.l. as UAV.
The flight plan was realized with the open source Mission Planner software, by setting the suitable parameters in order to have a GSD of 11 cm.The flights were designed according to a classic photogrammetric criterion, with a flight pattern consisting of 12 swipes, with forward overlapping of 80% and sidelap of 63%.The radiometric video sequences were acquired with nadiral orientation at a constant speed of 3 m/s and altitude of 40 m AGL.More details of the flight plan are reported in Table 3 (right column).
The UAV platform mounted the thermal camera OPTRIS PI450, with maximum optical resolution of 382 x 288 pixels at 80 Hz frame rate.The fixed focal length is 8 mm, with an angular FOV of 62° x 49°.The focal plane array (FPA) image sensor is made of uncooled microbolometers (25 x 25 μm each), with a spectral TIR response in the range 7.5-13 μm and accuracy of ± 2 °C.Table 2 (right column) shows a summary of thermal camera specifications.The camera can be connected to a computer via a USB 2.0 protocol to download the acquired data in radiometric video sequences format (.RAVI).Single frames had consequently to be extracted from the video to be suitable for a classical photogrammetric project.
The same targets used for RGB photogrammetry were used for georeferencing the two thermal orthomosaics: since aluminium foil has a very low emissivity (0.09 ε) it can be easily seen in the TIR range.Some other panels were painted in black and white, onto the aluminium foil, and distributed all around the test field, as delimited test areas for pre-flight calibration.In fact, before the TIR flights, a calibration with the thermal sensor was made by detecting the thermal response over the GCP targets and the b/w panels and by comparing them with temperature measurements taken with contact thermometers.
The final values of ambient temperature and emissivity were kept as default, Tamb 23 °C and ε = 1 (vegetation usually has high values of emissivity, around 0.98 [34]) during the data processing step.

TIR Data Capture and Processing
The flight campaign for the thermal data capture was carried out on the same day with two distinct flights, at 08.50 CET and at 15:00 CET, by using an hexacopter customized by MicroGeo s.r.l. as UAV.
The flight plan was realized with the open source Mission Planner software, by setting the suitable parameters in order to have a GSD of 11 cm.The flights were designed according to a classic photogrammetric criterion, with a flight pattern consisting of 12 swipes, with forward overlapping of 80% and sidelap of 63%.The radiometric video sequences were acquired with nadiral orientation at a constant speed of 3 m/s and altitude of 40 m AGL.More details of the flight plan are reported in Table 3 (right column).
The UAV platform mounted the thermal camera OPTRIS PI450, with maximum optical resolution of 382 × 288 pixels at 80 Hz frame rate.The fixed focal length is 8 mm, with an angular FOV of 62 • × 49 • .The focal plane array (FPA) image sensor is made of uncooled microbolometers (25 × 25 µm each), with a spectral TIR response in the range 7.5-13 µm and accuracy of ±2 • C. Table 2 (right column) shows a summary of thermal camera specifications.The camera can be connected to a computer via a USB 2.0 protocol to download the acquired data in radiometric video sequences format (.RAVI).Single frames had consequently to be extracted from the video to be suitable for a classical photogrammetric project.
The same targets used for RGB photogrammetry were used for georeferencing the two thermal orthomosaics: since aluminium foil has a very low emissivity (0.09 ε) it can be easily seen in the TIR range.Some other panels were painted in black and white, onto the aluminium foil, and distributed all around the test field, as delimited test areas for pre-flight calibration.In fact, before the TIR flights, a calibration with the thermal sensor was made by detecting the thermal response over the GCP targets and the b/w panels and by comparing them with temperature measurements taken with contact thermometers.
The final values of ambient temperature and emissivity were kept as default, T amb 23 • C and ε = 1 (vegetation usually has high values of emissivity, around 0.98 [34]) during the data processing step.For this reason, all the resulting thermal values must be considered only as qualitative information, because for absolute values a previous radiometric calibration would have been made.
The output of the thermal survey were, as mentioned, .RAVI radiometric videos, managed with the camera related software Optris PI Connect.The procedure for the frame extraction, in order to obtain pictures combined with the thermal information to use in the photogrammetric project, consisted of the following steps: • automatic extraction of single frames from the thermal video, by setting the proper time of acquisition to maintain a suitable overlapping (Auto Key Presser was used as software).The result is a series of .csvdatasheet files for each screenshot; • conversion of the .csvfiles into 16-bit TIF images, through a point-to-point conversion (software ThermoVision_JoeC v. 1.0.6.0), to obtain three kinds of files: i) temperature 16-bit raster images (thermograms) with black and white values scaled according to the min/max temperature values of the overall set of images, ii) pictures with color palette set according to their own min/max temperature values, iii) an overview .txtfile which reports the overall temperature range; • evaluation of thermal outliers for each frame: a temperature range too large results in an insufficient contrast for the identification of homologous points by the photogrammetric software.In this case, the frames containing people were removed, because body temperature (~36 • C) increased the temperature maximum overall value.On the contrary, the aluminium targets, which reflect sky temperature, gave outliers of about −30 • C. In fact, the solar rays enter the thermal camera as reflected instead of emitted radiation, thus compromising the measures.For this reason, a default threshold of 0 This procedure of data processing has been developed and tested at the Polytechnic University of Milan, by Sona, as reported in Reference [44].

GIS Processing
The obtained thermal orthomosaics, were then processed in GIS software (open source QGIS).Every single vineyard row included in the test area has been acquired as vector feature, by manual photointerpretation of the RGB orthophoto, in order to classify only the vegetation component.Each row has been identified considering the position on the right and on the left in respect to the dry-stone wall of that terraced bench.Then a statistical zonal analysis was performed on each row, resulting from the overlay with the two thermal orthomosaic raster, to classify each row by their own thermal characteristics.

Statistical Analysis
Vineyard row thermal behaviour was finally analysed starting from the differences in average temperatures between the ones in external and the ones in internal position.For each bench terrace, we considered as external row the one at the edge of the terrace, and as internal the ones closest to the dry-stone wall supporting the overlying bench, as shown in Figure 3.All the rows placed in between the external and internal rows of each bench are identified as intermediate.For both the morning and the afternoon flights, the mean and the standard deviation of the temperatures of the internal rows and the one of external rows were computed, and the difference of the means was then calculated and analysed.In both cases, the level of significance of the differences in mean temperatures of the two samples (internal rows and external rows) was tested with Student's t-test for the comparison of two means.

Results of VIS Image Analysis
The processing of the 206 pictures from the VIS flight resulted in a 2 mm GSD orthophoto with a pixel size of 2.41 x 2.41 μm 2 (see Figure 4).The first features-matching procedure implemented in the SfM algorithm, produced a point cloud consisting of 1,308,933 features in a local coordinate system over an area of 114,000 m 2 .By inserting the GCPs' coordinates, a bundle adjustment was performed to register the model in the UTM reference system, as reported in Section 3.1.The statistical values of the external orientation are listed in Table 4, reported as total RMSE values of GCPs and CPs.After the matching process which generated a sparse cloud, a denser cloud of 110,536,245 points was generated.Successively, a DEM surface model was created with a resolution of 4 cm/pix and point density of 772 points/m².

Results of VIS Image Analysis
The processing of the 206 pictures from the VIS flight resulted in a 2 mm GSD orthophoto with a pixel size of 2.41 × 2.41 µm 2 (see Figure 4).The first features-matching procedure implemented in the SfM algorithm, produced a point cloud consisting of 1,308,933 features in a local coordinate system over an area of 114,000 m 2 .By inserting the GCPs' coordinates, a bundle adjustment was performed to register the model in the UTM reference system, as reported in Section 3.1.The statistical values of the external orientation are listed in Table 4, reported as total RMSE values of GCPs and CPs.After the matching process which generated a sparse cloud, a denser cloud of 110,536,245 points was generated.Successively, a DEM surface model was created with a resolution of 4 cm/pix and point density of 772 points/m 2 .

Results of TIR Image Analysis
The different spectral range and the lower spatial resolution of the TIR sensor compared to the VIS one resulted in orthomosaics with ground resolution of 12 cm/pix both.In fact, the lack of information in the grey-scaled thermograms for the feature-matching procedure of the software gave a lower number of tie points as a result, compared to the VIS orthophoto, respectively of 76,287 for the morning orthomosaic and of 105,949 for the afternoon one.
The same GCPs used to scale and georeferencing the VIS 3D model, gave different values of RMSE for the TIR projects, as reported in Table 4.The higher error values associated to the TIR CPs

Results of TIR Image Analysis
The different spectral range and the lower spatial resolution of the TIR sensor compared to the VIS one resulted in orthomosaics with ground resolution of 12 cm/pix both.In fact, the lack of information in the grey-scaled thermograms for the feature-matching procedure of the software gave a lower number of tie points as a result, compared to the VIS orthophoto, respectively of 76,287 for the morning orthomosaic and of 105,949 for the afternoon one.
The same GCPs used to scale and georeferencing the VIS 3D model, gave different values of RMSE for the TIR projects, as reported in Table 4.The higher error values associated to the TIR CPs could be attributed to the difficulties in exactly picking the centre of the targets, mainly because of the low image quality, thus resulting in wrong identifications.In fact, the calculated Instantaneous Field Of View (IFOV) is of 91 mm, by considering the lens characteristics (FOV) of the thermal camera and the distance from the ground, thus giving a spatial resolution of about 1 cm, i.e., the minimum detail that can be seen on the thermal map.
The two generated dense clouds had similar sizes, of 964,066 points for the morning and 880,995 points for the afternoon.The DEM surface model was finally created with a resolution of 23 cm/pix and 18 points/m 2 for both.From a visible comparison, it appears that the external rows are warmer compared to the internal ones in the morning, while the same have similar temperature in the afternoon.The subtraction of the morning orthomosaic from the afternoon thermal map is shown in Figure 6, where the different temperature range of internal rows in the morning compared to the afternoon is clearly visible.

Results of Thermal Behaviour Analysis
The analysis of morning and afternoon temperatures showed how, on average, internal rows were colder then external ones of more than 2 °C, with a statistical significance of >99%, while temperatures recorded at 15:00 CET were almost equal for the two groups, around 23 °C (Table 5).Furthermore, the scatter plot of the morning and afternoon temperatures (Figure 7) showed how, in the morning, the whole population of internal rows exhibited temperatures lower than the external ones.On the other hand, afternoon temperatures were at the same level.The two populations appeared then as clearly separated.
To further investigate spatial and temporal temperatures patterns in the terraced hillslope, a full analysis of all rows of Grospoli vineyard was developed, considering also the values of temperatures of intermediate rows (Figure 8).The analysis showed how, for both the North and the South areas (see Figure 2), temperatures progressively decreased from the external edge to the internal position in the morning, while they grew up to the same value after 5 hours of insulation.

Results of Thermal Behaviour Analysis
The analysis of morning and afternoon temperatures showed how, on average, internal rows were colder then external ones of more than 2 • C, with a statistical significance of >99%, while temperatures recorded at 15:00 CET were almost equal for the two groups, around 23 • C (Table 5).Furthermore, the scatter plot of the morning and afternoon temperatures (Figure 7) showed how, in the morning, the whole population of internal rows exhibited temperatures lower than the external ones.On the other hand, afternoon temperatures were at the same level.The two populations appeared then as clearly separated.
To further investigate spatial and temporal temperatures patterns in the terraced hillslope, a full analysis of all rows of Grospoli vineyard was developed, considering also the values of temperatures of intermediate rows (Figure 8).The analysis showed how, for both the North and the South areas (see Figure 2), temperatures progressively decreased from the external edge to the internal position in the morning, while they grew up to the same value after 5 hours of insulation.

Discussion
Agricultural terraces are a fundamental feature of Mediterranean landscapes, which are threatened (among other causes) by the risk of abandonment and consequent degradation due to the progressive depopulation of marginal areas and the industrialization of agriculture [2].In some cases, these ancient structures allowed the settlement and the survival of entire peoples, which through the practice of terracing have allowed cultivation of the steepest slopes and, thanks to the enhanced slope stability, the building of villages at the valley bottoms, containing the risk of landslides.In these situations, the abandonment of such anthropic structures can even imply the increase of hazards related to geo-hydrological processes that could be triggered by extreme rainfall events.In this context, a dramatic example is represented by the tremendous event occurred in the Cinque Terre area in 2011 [45,46].However, the aspect of civil protection is limited to the few specific cases where settlements lie under terraced slopes.Beyond this, agricultural terraces represent a cultural heritage which should be maintained as part of a culture strongly connected to every territory where it is, or it was, present [47].In this perspective, an economical sustainability of agriculture on terraced land should be encouraged and helped to keep on going despite the higher costs that it involves.For this reason, understanding and quantifying all the benefits produced by a terraced landscape becomes of paramount importance to address effective policies for its conservation.
In this context, the availability of a cost-affordable and efficient methodology for monitoring of terraces evolutions, in terms of both spatial morphology (RGB) or thermal behaviour induced by the presence of dry-stone walls (TIR), can be a useful tool for risk prevention and managing purposes.
Our preliminary tests carried out by TIR surveys showed some influence on the temperature distribution across the vineyard due to the presence of terraces.
The thermal behaviour of the rows can be well explained by considering the sun position in respect to the vineyard.Since the slope of the terraced vineyard is aligned to the E-W position, at the sunset the internal rows are shaded by the presence of overlying dry-stone walls.From 08:50 CET to 15:00 CET, the sun moves along its trajectory, irradiating progressively the intermediate rows, in the middle of each bench, and then the internal ones.At 15:00 CET, the sun's position allowed the irradiation of all the rows, as reported in Figure 9.

Discussion
Agricultural terraces are a fundamental feature of Mediterranean landscapes, which are threatened (among other causes) by the risk of abandonment and consequent degradation due to the progressive depopulation of marginal areas and the industrialization of agriculture [2].In some cases, these ancient structures allowed the settlement and the survival of entire peoples, which through the practice of terracing have allowed cultivation of the steepest slopes and, thanks to the enhanced slope stability, the building of villages at the valley bottoms, containing the risk of landslides.In these situations, the abandonment of such anthropic structures can even imply the increase of hazards related to geo-hydrological processes that could be triggered by extreme rainfall events.In this context, a dramatic example is represented by the tremendous event occurred in the Cinque Terre area in 2011 [45,46].However, the aspect of civil protection is limited to the few specific cases where settlements lie under terraced slopes.Beyond this, agricultural terraces represent a cultural heritage which should be maintained as part of a culture strongly connected to every territory where it is, or it was, present [47].In this perspective, an economical sustainability of agriculture on terraced land should be encouraged and helped to keep on going despite the higher costs that it involves.For this reason, understanding and quantifying all the benefits produced by a terraced landscape becomes of paramount importance to address effective policies for its conservation.
In this context, the availability of a cost-affordable and efficient methodology for monitoring of terraces evolutions, in terms of both spatial morphology (RGB) or thermal behaviour induced by the presence of dry-stone walls (TIR), can be a useful tool for risk prevention and managing purposes.
Our preliminary tests carried out by TIR surveys showed some influence on the temperature distribution across the vineyard due to the presence of terraces.
The thermal behaviour of the rows can be well explained by considering the sun position in respect to the vineyard.Since the slope of the terraced vineyard is aligned to the E-W position, at the sunset the internal rows are shaded by the presence of overlying dry-stone walls.From 08:50 CET to 15:00 CET, the sun moves along its trajectory, irradiating progressively the intermediate rows, in the middle of each bench, and then the internal ones.At 15:00 CET, the sun's position allowed the irradiation of all the rows, as reported in Figure 9. Therefore, thermal differences appeared to be mainly influenced by the shading effect of the terrace risers, but we do not exclude the possibility that the walls themselves could have other effects related to the thermal behavior of stones.
A small fluctuation between the average temperature of internal and external rows can still be observed in the afternoon thermal data (see Figure 8).This variation might be explained by considering the high thermal inertia associated to stone wall.In fact, the medium value of specific heat for stone, is around 840 J/kg• • C (at 25 • C and atmospheric pressure) [48], thus making it a very common material in historical buildings of the Mediterranean area, typically characterized by thick and massive walls [49].In that context the main function was related to prevent summer overheating by keeping low and comfortable indoor temperature.The poor heat conductivity of rocks, which exact value depends on stone porosity, mineral composition and specific heat conductivity, usually entails low heat transfer from the irradiated surface to the inner part of the rock [48].These properties can be applied also on the dry-stone walls of the vineyard, which will warm up slowly during the day and they will passively release heat during the night, thus influencing also the facing internal rows temperature and microclimate.
Barbera et al. [38] analysed a wide range of environmental parameters to investigate the role of particular dry-stone wall structures in modifying the microclimate for crops.Among other interesting findings, they showed how the different thermal inertia of stones result in temperature differences (stones are colder in the morning and warmer in the evening, also showing very high temperatures when heated by the sun) which might influence the local microclimate.Working in a totally different context, Argyle and Stevens [50] investigated the distribution and growth patterns of a spontaneous shrub, the netleaf hackberry, growing around boulders in the cold desert biome of Wasatch Mountains (Utah).Unexpectedly, they found that hackberries associated with boulders were more likely to grow near the south side, suggesting that big stones may influence patterns of growth providing thermal radiation that melts snow thus providing a higher water availability for seedlings.Despite the fact that their research is far from the case of terraced vineyards, we hypothesize that the thermal inertia of dry-stone walls can have an influence, at least at the microclimatic scale, on temperature patterns of the crop and thus on its growth, phenology, ripening, and consequently on quality.However, neither works available in literature nor our findings have not yet shown it clearly.We believe that repeated measurements taken in different days at equal or comparable solar insulation conditions would provide a reliable time series to confirm these patterns of temperatures.
Nevertheless, the UAV-mounted TIR sensor, coupled with a RGB camera, proved to be a suitable monitoring method to reveal thermal variations induced by both natural sources (sun) as by anthropic artifacts (dry-stone walls).The main advantages of this approach consist in: (i) contained operative costs (ii) higher spatial resolution compared to traditional remote sensing platforms for TIR sensors (aircrafts and satellites); (iii) higher spatial covering in respect to ground-based TIR surveys; (iv) possibility of georeferencing spatial (RGB) and radiometric (TIR) information in a GIS software to maximize the data available for managing and prevention purposes; (v) reduced time required for data acquisition, thus allowing to repeat the survey multiple times during the day/month/year for monitoring of daily/seasonal variations; (vi) sufficient radiometric accuracy to discriminate qualitative daily thermal variations.Some limitations, however, need to be pointed out.As Conte et al. reported in their work [51], having both high spatial and radiometric resolution in aerial thermal imagery is currently not possible.Despite the above-mentioned advantages in using UAVs as platforms and the progressively miniaturizing process of multi-spectral sensors, the two resolutions are not comparable yet, even if an exponential technological growth is happening.
Another aspect to consider, is related to calibration.A radiometric calibration should be carried out in order to obtain absolute values of temperature.In this work, a relative thermal information was enough to obtain some preliminary information and to test the adopted methodology.For this reason, statistical analysis on exact temperature values (reported in Section 4.3) has to be considered as a relative information about the different thermal behaviour between morning and afternoon.Furthermore, the temperature resolution of the TIR camera does not allow us to make quantitative considerations.A proper radiometric calibration with a black body in laboratory, as well as setting adequately other factors as ambient temperature, emissivity, and spot size, might reduce the possible margin of error for more accurate temperature evaluations.
A geometric calibration, instead, concerns the metric accuracy of the orthomosaics, which is fundamental when RGB and TIR data are overlapping in GIS.In fact, TIR orthomosaics geometric resolution is significantly lower than the RGB, where accuracies of 0.1-0.2pix can be achieved.The main source of error is due to the radial distortion of the TIR sensor, deriving from refraction variations and other factors, as reported by Luhmann [52].For this reason, a pre-flight calibration for the thermal camera is suggested, with special targets to enhance contrast in thermal images.In fact, the reduced contrast and the blurring effects on objects edges make very difficult a precise collimation of a point in thermal frames.The authors suggest using algorithms for radiometric image enhancement and increasing the number of overlapping images (i.e., redundancy) to ease the photogrammetric workflow (SfM algorithm).Furthermore, proper accuracy evaluations should be made by considering fixed calibration certificates, obtained in laboratory, and the analysis of computed ones, in addition the sole evaluation of check point residuals.
Further analysis should be considered for evaluating the potential thermal inertia effect of dry-stone walls on the surrounding vegetation, as shown by Barbera et al. [38].Depending on its strength, this additional effect may have the potential of influencing local microclimate, as well as the characteristics of cultivations grown on bench terraces, such as grapevine.Considering the results of the present experiment, the use of TIR remote sensing based on UAV appears to be a reliable technological solution to assess this additional effect, for instance by mapping temperatures during or after the sunset.Further tests should be also performed to improve the photogrammetric process, as reported by Arandjelović et al. [53] for the registration of time-separated aerial images, to improve the overall geometric accuracy of the produced outputs (orthomosaics).Moreover, applications of TIR remote sensing that can be useful for analysis and monitoring of terraced landscape are not limited to thermal analysis.In fact, as shown by Glaser et al. [54], TIR analysis can highlight surface and subsurface water saturation dynamics.This application can provide additional insights to the studies carried out on the generation of preferential drainage patterns in terraced hillslopes, fundamental to understanding terraces degradation dynamics, that so far were based only on geo-morphological analyses and modelling [42,43,55].
Another future extension of this contribution might be the integration of different sensors on the same UAV platform, operating in different ranges of the electromagnetic spectrum, such as multi-or iper-spectral sensors [33,56,57].In this way, additional information coming from selected or wider band of the spectral range would complete and extend the data collected, with biophysical vegetation indexes useful also for precision agriculture purposes.

Conclusions
The present work describes an early application of the use of combined multi-sensor surveys, in the VIS and TIR range mounted of UAV platforms, to detect thermal dynamics of dry-stone terraced areas, which represent one of the most iconic features of agricultural landscapes in Italy and Europe.The use of UAVs platforms for remote sensing presents advantages in terms of operative cost, high planning flexibility and higher spatial resolution if compared to common remote-sensing platforms (such as aircrafts and satellites) [8].The methodology tested in this work, with an application in the historic Chianti area of Tuscany, proved to be an efficient tool for monitoring and managing these traditional rural landscapes.In particular, the UAV-mounted TIR sensor, coupled with an RGB camera approach allowed us to: • obtain sufficient radiometric accuracy to reveal daily thermal variations induced by both natural sources (sun) as by anthropic artifacts (dry-stone walls) • contain the operative costs and have higher spatial resolution compared to traditional remote sensing platforms for TIR sensors (aircrafts and satellites) • cover a wider area in respect to ground-based TIR surveys • georeferencing spatial (RGB) and radiometric (TIR) information in a GIS software • repeat the survey multiple times during the day, thanks to reduced time needed for data acquisition The preliminary results proposed in this contribution highlight a different thermal behavior of rows in the morning, when the internal rows show lower temperature than the external rows, due to solar position and shading effect.After some hours of direct irradiation, in the afternoon, temperature of both internal and external rows increased, reaching almost the same values.
These outcomes, as well as the rapid evolution in sensors and UAV technologies, encourage further tests with the above-mentioned methodology to detect microclimatic effects of dry-stone walls on the bench-growing vegetation.For instance, repeating thermal flights more frequently within the entire 24 hours could help to better understand temperature patterns across the entire terraced field.Other future activities might concern the integration with other sensors, operating in different bands of the electromagnetic spectrum to obtain biophysical indexes of vegetation, as well as monitoring at different time of the day/month/year.
32 34.73 N, 11 • 21 29.14 E), represents one of the most relevant examples of agricultural terracing sites of the country.

Figure 2 .
Figure 2. Location of the study area.Red border -North area; blue border -South area of the Grospoli vineyard in Tuscany region (Italy).

Figure 2 .
Figure 2. Location of the study area.Red border -North area; blue border -South area of the Grospoli vineyard in Tuscany region (Italy).

Figure 2 .
Figure 2. Location of the study area.Red border-North area; blue border-South area of the Grospoli vineyard in Tuscany region (Italy).

Table 3 .
The three flight plans in the visible (VIS) and Thermal InfraRed (TIR) ranges.The pictures represent camera locations (black dots) of each swipe and image overlap (colors).
ISPRS Int.J. Geo-Inf.2019, 8, x FOR PEER REVIEW 9 of 21 temperatures of the two samples (internal rows and external rows) was tested with Student's t-test for the comparison of two means.

Figure 3 .
Figure 3.Each vineyard row has been labelled as: internal (blue) for the closest rows to the dry-stone wall (dotted line), external rows (red) for the ones at the edge of the terrace (black lines).All the rows comprised between internal and external have been identified as intermediate (green).

Figure 3 .
Figure 3.Each vineyard row has been labelled as: internal (blue) for the closest rows to the dry-stone wall (dotted line), external rows (red) for the ones at the edge of the terrace (black lines).All the rows comprised between internal and external have been identified as intermediate (green).

Figure 4 .
Figure 4.The orthophoto of the study area produced by the photogrammetric processing of the images acquired from the RGB flight over the Grospoli vineyard.

Figure 4 .
Figure 4.The orthophoto of the study area produced by the photogrammetric processing of the images acquired from the RGB flight over the Grospoli vineyard.

Figure 5
Figure 5 reports the two orthomosaics scaled differently, in function of their maximum and minimum temperature, to enhance visualization.In particular the morning orthomosaic shows a temperature range of 14.8-18.3• C while in the afternoon the T range is between 21 • C and 25 • C.From a visible comparison, it appears that the external rows are warmer compared to the internal ones in the morning, while the same have similar temperature in the afternoon.The subtraction of the morning orthomosaic from the afternoon thermal map is shown in Figure6, where the different temperature range of internal rows in the morning compared to the afternoon is clearly visible.

Figure 5 .
Figure 5. Thermal maps of the Grospoli vineyard in (a) morning (08.50CET) and (b) afternoon (15:00 CET).The enlarged windows highlight internal and external rows, showing their respective thermal behaviours.

Figure 5 .
Figure 5. Thermal maps of the Grospoli vineyard in (a) morning (08.50CET) and (b) afternoon (15:00 CET).The enlarged windows highlight internal and external rows, showing their respective thermal behaviours.

Figure 6 .
Figure 6.Difference between the afternoon thermal map and the morning orthomosaics highlighting, in the enlarged windows, the behaviour of the internal and external rows.

Figure 6 .
Figure 6.Difference between the afternoon thermal map and the morning orthomosaics highlighting, in the enlarged windows, the behaviour of the internal and external rows.

Figure 7 .Figure 7 .
Figure 7. Scatter plot of morning and afternoon temperatures for external (circles) and internal (black squares) rows.

Figure 7 .
Figure 7. Scatter plot of morning and afternoon temperatures for external (circles) and internal (black squares) rows.

Figure 8 .
Figure 8. Temperature average values along the hillside terraced Grospoli vineyard.(a) North area; (b) South area.White and grey timeseries represent morning and afternoon temperatures respectively.Black dots represent internal rows temperature mean values.

Figure 8 .
Figure 8. Temperature average values along the hillside terraced Grospoli vineyard.(a) North area; (b) South area.White and grey timeseries represent morning and afternoon temperatures respectively.Black dots represent internal rows temperature mean values.

Table 1 .
Technical specifications of DJI Phantom 4 Pro.

Table 3 .
The three flight plans in the visible (VIS) and Thermal InfraRed (TIR) ranges.The pictures represent camera locations (black dots) of each swipe and image overlap (colors).

Table 3 .
The three flight plans in the visible (VIS) and Thermal InfraRed (TIR) ranges.The pictures represent camera locations (black dots) of each swipe and image overlap (colors).

Table 4 .
Statistics of Ground Control Points (GCPs) and Check Points (CPs) for RGB and TIR photogrammetric projects provided by Agisoft PhotoScan.

Table 5 .
Statistical analysis of averages of morning and afternoon temperatures, for external and internal rows (°C).

Table 5 .
Statistical analysis of averages of morning and afternoon temperatures, for external and internal rows ( • C).