**1. Introduction**

Climate change and related natural disasters affect several sectors [1]. Agriculture is one of the most vulnerable [2,3]. Between 2005 and 2015, the impact of natural disasters on the agricultural sector was estimated to be 96 billion dollars in damaged, or completely lost, crops [4]. Climate change-related effects (e.g., temperature and precipitation increasing in terms of level, time, and variability) are expected to reduce the yield and quality of many crops, especially cereals and fodder cereals [5].

Storms and hail also can cause serious damage to crops [6]. Hurricanes can cause much damage, with grass lodging, uprooting of orchards, and falling trees [7,8]. These critical events, potentially highly impacting farmers' income, must be carefully accounted for in the context of risk management in agriculture.

Fruits and vegetables represent (year 2018) about 14% of the total value of European (EU) agricultural production [9,10]. These crops are very important for many EU member states, in particular for Mediterranean countries such as Spain, Italy, and France. Italy is one of the main European leaders in the apple sector [11]. Consequently, the yield loss risks concerning the fruit and vegetable sector must be minimized. Major threats concern diseases, insects, and natural disasters such as hail, drought, frost, and storms.

**Citation:** De Petris, S.; Sarvia, F.; Gullino, M.; Tarantino, E.; Borgogno-Mondino, E. Sentinel-1 Polarimetry to Map Apple Orchard Damage after a Storm. *Remote Sens.* **2021**, *13*, 1030. https://doi.org/ 10.3390/rs13051030

Academic Editor: Mario Cunha

Received: 9 February 2021 Accepted: 5 March 2021 Published: 9 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

Apple cultivation is very intensive today, with a plant density around 2000 plants per hectare [12]. Such density allows a very high yearly production (about 45 tons per hectare) [13], which is obtained by a row-based cultivation strategy where young plants begin to be productive after the third year. The adoption of low-vigor rootstocks enables an increase in planting density and rapid fruiting. Unfortunately, this kind of cultivation determines a very underdeveloped root system, not enough to guarantee plant stability under unfavorable conditions. The situation is more critical during extreme weather events, especially when there are many weighty fruits, i.e., before harvesting [14,15]. Steel cables anchored to concrete or wooden poles are used to improve row stability.

Within this context, when a stormy event occurs, it is important to assess the spatial level and extent of damage to start remedial actions and minimize crop loss. Farmers are interested in damage estimation especially when a refund is due by insurance companies [16,17]. In this case, damage is assessed through on-the-spot checks by an expert surveyor from the insurance company, who determines the extent, type, and quality of damage. Such an approach depends on a high level of subjectivity related to the expert's skill and experience. Moreover, these operations require a lot of time and are expensive, especially where large areas have been affected by the event.

In this operative context concerning crop damage analysis, a more objective monitoring could play a key role, providing more robust forecasts about potential yield or yield losses. Many agricultural stakeholders, such as farmers, consortia, agronomists, insurance companies, and local administrations, require a continuous monitoring of crops over large spatial extents.

A method based on free Earth Observation (EO) data can certainly represent an effective support [18] and the consequent technological transfer desirable [19–24].

In particular, optical remote sensing data have been successfully used in several operational frameworks, as proved by many works [25–32]; unfortunately, cloud conditions limit the nominal temporal resolution of this type of data, especially while monitoring natural disasters (e.g., floods or storms) that ordinarily occur when clouds are present. Data from synthetic aperture radar (SAR) systems can operate during all-weather conditions, and, while exploring agronomical issues, they can be used to analyze the moisture and geometrical conditions of crops [33–35]. In particular, dual-polarimetric SAR acquisitions from Copernicus Sentinel-1 mission (S1) provide unique opportunities to disseminate operational monitoring for several application communities [36,37]. Dual-pol acquisition mode has a larger swath and a lower data volume compared with full-pol acquisitions, thus improving data collection and processing for operational activities [38,39]. Polarimetric data can provide information about polarization amplitude and phase, allowing scattering mechanism definition (i.e., single-bounce, double-bounce, or volume scattering) induced by target properties. SAR polarimetry (PolSAR) is a technique that analyzes SAR polarization with respect to the vector of polarized electromagnetic waves. When a signal passes through a medium, the refraction index changes, or when it strikes an object, it is reflected; the so-called backscattering matrix [40] contains information about the reflectivity, shape, and orientation of the reflecting target. An important improvement in the extraction of physical information from the ordinary coherent backscattering matrix was achieved by Cloude and Pottier [41,42], who proposed the composition of system vectors. Most studies have assessed the sensitivity of polarimetric indicators derived from the C-band space-borne SAR to derive crop parameters [43]. The PolSAR technique was successfully applied to monitor crop growth and give estimates of yield. For example, Betberder [44] analyzed temporal trends of polarimetric indicators, proving their high potential to detect crop growth changes. Valcarce [45] used polarimetric data time series for land-cover classification, adopting a decision tree classification algorithm performing high crop class detection accuracies. Mercier and Qi [46,47] used PolSAR to support/integrate vegetation phenology monitoring based on optical data.

Only few works referring to PolSAR application in crop damage analysis are present in the literature [48,49], denoting a lack of scientific production about this issue. Nevertheless, hailstorms and storms are known to change vegetation structure, resulting in lodging or tree uprooting/breaking. Therefore, this peculiar effect changes polarimetric response and could be used to detect and characterize tree structure [50]. In general, it can be said that decomposition techniques offer a new insight into PolSAR data for describing vegetation structural proprieties [51].

The polarimetric decomposition technique decomposes the signal into its individual scattering components, permitting identification of the dominant scattering type [42,52]; this information is related to the target structural properties [18,53,54]. Various decomposition techniques have been proposed, and Lee and Cloude provided a comprehensive review about this topic [42,55]. Model-based [56] and eigenvector-based [41] algorithms have been preferred by many researchers [51]. According to Ji and his collaborators [57], the Cloude–Pottier H-α-A decomposition seems to be the most promising approach. It is based on second-order statistics extracted by a set of neighbor pixels that are used to calculate the local entropy H and the α angle (related to average scattering mechanisms). These are used to define a Cartesian space, H-α, that is linearly divided into nine zones describing the main scattering mechanisms. Recently, eigenvector decomposition has been widely applied in several applications [55,58–61]. The method was originally developed for quadpolarization data. Nevertheless, it was also adapted to work with dual-pol data [57,62,63], and consequently, it can be successfully used to retrieve polarimetric information also from S1 data that are unable to collect quad-pol data.

In this work, the applicability of the H-α-A polarimetric decomposition technique to the detection and mapping of damages from storms affecting fruit orchards was tested. In particular, the proposed case study refers to the stormy event that occurred in Northwest Italy on 12 August 2020. Consequently, a map of potentially damaged orchards was generated with the aim of supporting insurance companies and local administrations to address their funding restoration policies.

## **2. Materials and Methods**

#### *2.1. Study Area*

On 12 August 2020, an exceptional storm affected the Northwest of Italy. In particular, the storm uprooted many apple orchards in the province of Cuneo (Piemonte region, NW Italy). Moreover, it occurred in a critical period of the year, when the main fruits (apples, pears, and peaches) were still to be harvested (Figure 1). Because in this period the farmers are focused on harvesting, no early recovery efforts were performed in the damaged fields. Therefore, the majority of the uprooted trees were not removed until October.

**Figure 1.** An apple orchard (cultivar "Gala") with hail nets uprooted by the storm on 12 August 2020. At the bottom, many mature apples can be noted, suggesting the economic loss caused by the storm.

The study area includes four municipalities: Saluzzo, Verzuolo, Manta, and Lagnasco (Figure 2). The area of interest (AOI) is sized about 132.23 km2. It plays a crucial economic role in Piemonte fruit production. In fact, this zone is suitable for this cultivation: the loose soil without water stagnation, sunny and dry atmosphere, and strong temperature difference between day and night allow the correct ripening and coloring of fruits. Apples represent the primary crop in Manta. Since August is a droughty period in the AOI, no significative previous precipitations had occurred before the event; 1.2 mm had cumulated in the previous week, as reported by the regional environmental agency (www.arpa. piemonte.it). Therefore, the authors supposed that moisture-related conditions cannot significantly affect the SAR signal.

**Figure 2.** Italian regions (light gray) and the Piemonte region (dark gray). (Red) The AOI includes the Saluzzo, Manta, Lagnasco, and Verzuolo municipalities (reference frame: WGS84 UTM32N).

#### *2.2. Data and Data Collection*

#### 2.2.1. Sentinel-1 Data

Sentinel-1 is currently one of the largest space-borne missions providing free and openly accessible SAR data. The S1 mission relies on a constellation of two satellites (Sentinel-1A and Sentinel-1B) operating in the C-band (5.54 cm wavelength). The main acquisition mode over land is the Interferometric Wide (IW) swath, recording approximately 250 km in length at 5 × 20 m spatial resolution in a single look. Ordinarily, S1 records data in a dual pole mode (VV and VH), where electromagnetic waves are polarized vertically (V) for transmission and horizontally/vertically for reception. The data are recorded as complex values (I/Q components) and in SAR geometry (range and azimuth). A descending single-look complex (SLC) IW image (relative orbit no. 139), acquired after the storm (14 August 2020), was obtained from the Copernicus Open Access Hub (https://scihub.copernicus.eu/dhus/#/home, accessed: 20 December 2020).

#### 2.2.2. Cadastral Data

A cadastral map coupled with farmers' applications for EU Common Agricultural Policy (CAP) incentives was used in this work to classify the orchards in the AOI. The correspondent map (hereafter called orchard map (OM)) was consequently generated. The damaged orchards were analyzed at cadastral parcel level. The cadastral map was obtained for free from the regional geoportal in vector format georeferenced in the WGS84 UTM zone 32N reference frame and updated in 2018 (nominal scale was 1:2000). Databases containing farmers' applications for EU CAP incentives of 2019 were used to map orchard types in the AOI (2020 data are not yet available). Every year, farmers support their activities with CAP

incentives. These data were obtained for free from the regional public information system for agriculture. CAP applications contain the cadastral parcel code and the declaration of the most relevant crops as communicated by farmers. In this way, it is possible to couple the cadastral map with crop type information at parcel level by an ordinary join operation available in the Geographical Information System (GIS) software. In this work, 2040 (about 1136 ha) apple orchards were selected from the joined data to test the procedure.

#### 2.2.3. Ground Dataset

A ground survey was conducted to gather the field data needed to calibrate and validate the PolSAR-based mapping procedure. In total, 72 apple orchards were surveyed (about 3.5% of the apple orchards in the AOI) during a ground campaign aimed at labeling damaged (22) and undamaged (50) fields. Specifically, the surveyed fields have an average size of about 0.92 ha, fitting well with the S1 geometrical resolution. In fact, about 40 S1 pixels can characterize each field. In particular, a visual assessment aimed at recognizing the following conditions was performed: if the majority of the trees were uprooted, the field was labeled as damaged; otherwise, it was labeled undamaged, and the related cadastral parcel was selected from the OM layer.

The dataset was split in a training (60%) and a test set (40%) by random selection from the surveyed parcels. In total, 13 damaged fields (hereafter called DTFs) and 28 undamaged ones (hereafter called UTFs) were assigned to the training set. Conversely, 10 damaged fields (hereafter called DVFs) and 21 undamaged ones (hereafter called UVFs) were assigned to the test set. The training and test set parcels are shown in Figure 3.

**Figure 3.** Parcels belonging to the training and test sets. Colors (see legend) define the state of the surveyed parcel (damaged/undamaged). Reference frame is WGS84 UTM 32N.

This dataset was provided by local farmers. The authors found that the supplied sample includes 72 fields corresponding to about 3.5% of the apple orchards in the AOI. The authors had just the opportunity of comparing the sample size with the expected total number of apple orchards in the AOI (about 2050). The authors are aware that this sample size does not perfectly fit statical requirements. Nevertheless, it well represents ordinary availability of ground data from farmers when working with actual data not directly managed by scientists. This situation well represents a common operational condition

when working with technology transfer issues, especially in the agronomic sector. In fact, the most of data from farmers, generally, rely on their autonomous collections and decision of making them public. Moreover, the private property of parcels is an objective limiting factor for all the analyses, since free access is not guaranteed. With these premises, we proceed to process the data.

A preliminary economic assessment was also performed since the storm occurred close to the apple harvesting period, determining a significant problem for local apple yield in 2020. This was obtained considering, for damaged parcels, a potential yield equal to the average one in the Piemonte region (31 t·ha−1) and a reference unitary price of 380 €·<sup>t</sup> <sup>−</sup>1. These values were obtained from the Italian Statistics Institute (ISTAT) [64].

#### *2.3. Data Processing*

#### 2.3.1. Polarimetric Decomposition

The available S1 IW SLC image was processed to compute the polarimetric decomposition parameters. The adopted workflow is shown in Figure 4 and proposed by [65]. The target polarimetric analysis is ordinarily performed starting from the coherency matrix [66,67] or from the 2×2 covariance matrix (**C2**). Preprocessing steps were managed using the ESA SNAP v. 7.0.0 software [68].

**Figure 4.** The adopted workflow. All steps were managed in SNAP ESA v. 7.0.

First, the precise orbit state vector data were downloaded from the ESA archive (https://qc.sentinel1.eo.esa.int/, accessed: 20 December 2020) and applied to refine the satellite position. Precise orbit files are delivered within 20 days after data acquisition and provide accurate satellite position and velocity information. Using the TOPS split module, 1 sub-swath and 2 bursts were selected based on AOI coverage. A radiometric calibration was applied and the result saved in a complex-valued format needed to compute **C2**. TOPS deburst was applied by merging different bursts into a single SLC image. A spatial subset was then generated covering the AOI. The subset was multi-looked by 4 × 1 (range and azimuth direction, respectively) to generate squared pixels. The resulting multi-looked image, with a geometrical resolution equal to 15 m, was used to generate the local **C2** at pixel level. With respect to quad polarization, dual-polarimetric SAR sensors generate a matrix showing the half of the totally occurring scattering components involved in fully polarimetric imagery [69]. In particular, the covariance matrix for dual polarization (e.g., Sentinel-1) is often calculated with reference to a second-order scattering information [18] generated from the spatial averaging of the scattering vector k = [SVV, SVH] <sup>T</sup> as expressed in Equation (1):

$$\mathbf{C}\_{2} = \begin{bmatrix} \mathbf{C}\_{11} & \mathbf{C}\_{12} \\ \mathbf{C}\_{21} & \mathbf{C}\_{22} \end{bmatrix} = \begin{bmatrix} \langle |\mathbf{S}\_{VV}| \rangle & \langle \mathbf{S}\_{VH}\mathbf{S}\_{VV} \rangle^{\*} \\ \langle \mathbf{S}\_{VH}\mathbf{S}\_{VV} \rangle^{\*} & \langle |\mathbf{S}\_{VH}| \rangle \end{bmatrix} \tag{1}$$

where ∗ denotes the complex conjugate and - the local mean value in a 5 × 5 moving window. Each **C2** element (C11, C22, (C12), and (C12)) is stored individually and successively refined by Lee filtering (5 × 5 kernel size) to minimize speckle-related noise. H-α-A polarimetric decomposition was obtained by eigenvector computation as proposed by different authors [57,62,63]. The modified formula for dual-pol data, as proposed by [66], is reported in Equations (2) and (3).

$$
\langle \mathcal{C}\_{\mathbf{2}} \rangle = \left[ \mathcal{U} \right] \begin{bmatrix} \lambda\_1 & 0 \\ 0 & \lambda\_2 \end{bmatrix} \left[ \mathcal{U} \right]^{\*T} = \left. \lambda\_1 \mu\_1 \mathbf{u\_1}^{\*T} + \lambda\_2 \mu\_2 \mathbf{u\_2}^{\*T} \right. \tag{2}
$$

$$\begin{bmatrix} \mathbf{U} \end{bmatrix} = \begin{bmatrix} \begin{array}{cc} \mathcal{U}\_{11} & \mathcal{U}\_{12} \\ \mathcal{U}\_{21} & \mathcal{U}\_{22} \end{array} \end{bmatrix} = \begin{bmatrix} \mathbf{u}\_{1} \ \mathbf{u}\_{2} \end{bmatrix} = \begin{bmatrix} \cos a & -\sin a e^{-j\delta} \\ \sin a e^{j\delta} & \cos a \end{bmatrix} \tag{3}$$

where λ<sup>1</sup> ≥ λ<sup>2</sup> ≥ 0 are the local eigenvalues, [**U**] is the orthogonal unitary matrix, \* and T represents the complex conjugate and transpose matrices, respectively. The angles α and δ define the orientation and size of the polarization ellipse of the recorded signal [62]. The eigenvector dual-pol decomposition results in three roll-invariant parameters: polarimetric scattering entropy (H), mean scattering angle (α), and scattering anisotropy (A).

H was calculated from Equation (4):

$$H = \ -\sum\_{i=1}^{2} (-P\_i \ \log\_2 P\_i) \tag{4}$$

where

$$P\_i = \frac{\lambda\_i}{\lambda\_1 + \lambda\_2}$$

H defines scatter randomness; it can vary between 0 and 1 and is related to the number of dominant scattering mechanisms, being proportional to the degree of depolarization [70]. H = 0 means that the coherency matrix shows only one eigenvalue and, therefore, the relative orientation of the correspondent pixel elements is quite simplified (e.g., singlebounce reflection).

Anisotropy A (Equation (5)) provides additional information about H in terms of the difference between scattering mechanisms.

$$A = \frac{\lambda\_1 - \lambda\_2}{\lambda\_1 + \lambda\_2} \tag{5}$$

The anisotropy quantifies the relative strength between first and second dominant scattering mechanisms. It is strictly related to the degree of signal polarization [18,71,72]. According to Mandal [18], the state of polarization of an electromagnetic (EM) wave is characterized in terms of the degree of polarization (0 ≤ A ≤ 1). The latter is defined as the ratio between the average intensity of the polarized portion of the signal and its total intensity [73]. A = 1 and A = 0 for a completely polarized and completely unpolarized wave, respectively. The unpolarized part of the received wave, (1 − A), is assumed to represent the volume scattering component from the distributed targets [74].

Average scattering mechanisms (i.e., surface, double-bounce, and volume scattering) can be identified with respect to the α parameter, which is computed according to Equation (6):

$$\alpha = \sum\_{i=1}^{2} P\_i \cos^{-1} \left( \frac{\left| \lambda\_1 + \lambda\_2 \right|}{\sqrt{2} \sqrt{\left| \lambda\_1 \right|^2 + \left| \lambda\_2 \right|^2}} \right) \tag{6}$$

The α angles close to 0◦ denote a diffuse surface scattering, α close to 45◦ means dipole scattering (caused by volumes), and α close to 90◦ means double-bounce scattering mechanisms.

With these premises, the raster layer mapping local H, α, and A values was computed from the pre-processed SLC image. It was projected onto the WGS84 UTM 32N reference frame, applying the range–Doppler terrain correction. The adopted digital terrain model (DTM) needed for this step was the one freely obtainable from the Piemonte region geoportal [75]. It is supplied with a 5 m grid size and a height accuracy of ±0.30 m and was generated in 2011. The nearest-neighbor resampling method was adopted during the range–Doppler terrain correction.

#### 2.3.2. Testing H-α-A Values after the Storm

To assess how the storm changed the orchards' polarimetric behavior, a preliminary analysis was performed with reference to the training set. In particular, DTF and UTF pixels distributions were compared using the Mann–Whitney (MW) nonparametric test (onetailed) [76]. The MW null hypothesis is that DTFs and UTFs have an identical distribution. The one-sided alternative "greater" was set, assuming that the DTF cumulated frequency distribution was expected to have shifted to the right of the UTF one (i.e., DTFs were greater than UTFs) [77].

The authors preliminary explored the polarimetric indices' behavior using reference ground data. In particular, the frequency distributions were perceptively assessed using boxplots (see Section 3.1). The median value of distribution highlights a shift between damaged and undamaged fields. Therefore, to test these perceptive differences, the authors performed one tail test since the direction of changes is a priori known.

Three MW tests were performed to test if the DTF distributions of the H-α-A pixels within the parcels were statistically different from the UTF ones. All statistical analyses were performed using R software v. 3.6.3 [78]; conversely, spatial analysis was done using SAGA GIS 7.0 [79].

#### 2.3.3. Detection of Damaged Orchards

The main goal of this work was to test the capability of the PolSAR technique to recognize damaged orchards. For this task, UTFs were assumed as representatives of the state of undamaged orchards. Samples were sized about 23 ha and represented about 2% of OM. In spite of this small sample size, the UTFs preliminarily resulted in a good dataset, whose reliability was confirmed by ground surveys. With these premises, the H-α-A distributions within UTFs were used to represent the reference distributions of the undamaged orchards. All H-α-A distributions from the AOI mapped parcels were tested against undamaged ones by the MW test, checking the following conditions: (i) parcel H distribution was greater than that of the UTFs; (ii) parcel α distribution was greater than that of the UTFs; (iii) parcel A distribution was lower than that of the UTFs. The resulting MW U-statistic and related *p*-value were then mapped for each orchard parcel. Moreover, the compound probability (CP) [80] was also calculated according to Equation (7) using R software v. 3.6.3. CP represents the probability that the previously mentioned three conditions were simultaneously satisfied.

$$CP = (1 - p\_H)(1 - p\_a)(1 - p\_A) \tag{7}$$

where *pH* is the *p*-value resulting from the MW test under condition (i), *p<sup>α</sup>* is the *p*-value resulting from the MW test under condition (ii), and *pA* is the *p*-value resulting from the MW test under condition (iii). The resulting CP was then mapped for all OM parcels, representing its compound probability to have been damaged by the storm. A threshold value of CP able to separate damaged fields from undamaged ones has to be necessarily selected by final users, e.g., the insurance company or local public administration, according to their specific policies and strategies. Nevertheless, a possible solution is proposed here, relying on the standard error of the mean (SEM) of the CP distributions of the DTFs and UTFs. The estimated threshold value was used to generate the map of damaged orchards (DM): parcels showing a CP value lower than the threshold was classified as "undamaged," otherwise as "damaged." The DMs were then tested against the previously mentioned test set and the correspondent confusion matrix calculated to assess the accuracy of detection.

#### **3. Results**

*3.1. H-α-A Analysis*

The statistical distributions of H-α-A were computed with reference to DTFs and UTFs (Figure 5).

**Figure 5.** Boxplots of H-α-A distributions for UTFs and DTFs. The boxplot values are from bottom to top, respectively, 5th, 25th, 50th—cross is mean value—75th, and 95th percentiles. (**a**) Entropy pixel distribution; (**b**) alpha angle pixel distribution; (**c**) anisotropy pixel distribution.

The MW test results (Table 1) show that the H and α distributions of DTFs presented values significantly greater than UTFs; conversely, the A distribution of the DTFs was lower than that of the UTFs.

**Table 1.** MW test results obtained by comparing the H-α-A pixel distributions of DTFs and UTFs.


#### *3.2. Damaged Orchards' Mapping*

Based on the assumption that a storm can change the polarimetric behavior of orchards according to previously mentioned dynamics, a map of CP representing the parcel probability of being recognized as damaged was generated using the UTF dataset as reference (Figure 6).

With reference to CP, a threshold value was estimated to separate damaged fields from undamaged ones based on the SEM of CP statistic distributions of the DTFs and UTFs (Figure 7b). The DTFs showed a CP mean and SEM value of 0.715 and 0.125, respectively; consequently, one can assume that the CP mean value of all damaged orchards reasonably falls in the range 0.715 ± 0.125, about 0.6 being the lower boundary. A threshold equal to 0.6 was therefore selected to generate the DM binary classification (Figure 7a).

A total of 217 ha (430 orchards) of potentially damaged apple orchards were detected in the AOI. According to the OM layer, 19% of the apple orchards were damaged after the event.

**Figure 6.** A CP map of apple orchards in the AOI (Reference frame: WGS84 UTM32N).

**Figure 7.** (**a**) DM binary classification of the OM in the AOI (reference frame is WGS84 UTM 32N); (**b**) bar chart representing mean and 1 SEM of the CP for DTFs and UTFs.

The DM was validated with respect to the test set, and the correspondent confusion matrix computed (Table 2). Classification accuracy is defined here as the one for binary classification of imbalanced data [81–83] since, in the test set, the number of undamaged fields was significantly greater than that of damaged fields. The resulting precision and specificity were pretty high (0.80 and 0.71, respectively), while balanced accuracy was found to be 0.75. Overall accuracy was 0.74, while F1 score (harmonic mean of the precision and recall) and G-mean (geometric mean of sensitivity and precision) were both about 0.67.


**Table 2.** Metrics derived from the confusion matrix of the DM with respect to the test set. True positives (TPs): number of damaged elements predicted as damaged; false positives (FPs): number of undamaged elements predicted as damaged; false negatives (FNs): number of damaged ele-

Furthermore, it is worth stressing that the storm occurred close to the harvesting period, determining a significant problem for local apple yield in 2020. With reference to the AOI, a preliminary estimate of economic loss was computed to be about €2,500,000. Reported estimates could certainly vary according to the apple orchards' age, apple variety, plant density, agronomic management, and local soil properties. Nevertheless, these estimates constituted a preliminary assessment of storm damage that occurred on 12 August 2020. Future validation is expected to test these economic deductions.

Overall Accuracy 0.74 OA = (TP + TN) / (TP + TN + FP + FN) F1 Score 0.66 F1 = 2TP/ (2TP + FP + FN) Balanced Accuracy 0.75 BA = TPR + TNR/2 G-Mean 0.67 G-mean= sqrt (TPR \*PPV)
