Next Article in Journal
Heat vs. Health: Home Office under a Changing Climate
Previous Article in Journal
Identifying and Categorizing Sustainable Supply Chain Practices Based on Triple Bottom Line Dimensions: Evaluation of Practice Implementation in the Cement Industry
Previous Article in Special Issue
A Machine Learning-Enhanced 3D Reverse Design Approach to Personalized Garments in Pursuit of Sustainability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Grape Maturity Estimation for Personalized Agrobot Harvest by Fuzzy Lattice Reasoning (FLR) on an Ontology of Constraints

by
Chris Lytridis
1,
George Siavalas
1,
Theodore Pachidis
1,
Serafeim Theocharis
2,
Eirini Moschou
1 and
Vassilis G. Kaburlasos
1,*
1
HUMAIN-Lab, International Hellenic University (IHU), 65404 Kavala, Greece
2
Viticulture Laboratory, Department of Horticulture, Viticulture School of Agriculture, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(9), 7331; https://doi.org/10.3390/su15097331
Submission received: 31 March 2023 / Revised: 14 April 2023 / Accepted: 25 April 2023 / Published: 28 April 2023
(This article belongs to the Special Issue Computational Intelligence for Sustainability)

Abstract

:
Sustainable agricultural production, under the current world population explosion, calls for agricultural robot operations that are personalized, i.e., locally adjusted, rather than en masse. This work proposes implementing such operations based on logic in order to ensure that a reasonable operation is applied locally. In particular, the interest here is in grape harvesting, where a binary decision has to be taken regarding the maturity of a grape in order to harvest it or not. A Boolean lattice ontology of inequalities is considered regarding three grape maturity indices. Then, the established fuzzy lattice reasoning (FLR) is applied by the FLRule method. Comparative experimental results on real-world data demonstrate a good maturity prediction. Other advantages of the proposed method include being parametrically tunable, as well as exhibiting explainable decision-making with either crisp or ambiguous input measurements. New mathematical results are also presented.

1. Introduction

The exponentially increasing world population [1]—which has doubled every 50 years during the last 120 years—calls, in the short term, for a proportional increase in the food supplied by conventional agriculture. However, such an increase in food is hindered by various factors worldwide including labor shortages. At the same time, farmers are required to supply high-quality agricultural products using environmentally friendly methods.
To satisfy the above requirements, an extension of “Industry 4.0” digital practices to agriculture, namely “Agriculture 4.0”, is promising [2]. In the context of Agriculture 4.0, precision agriculture (PA) has been defined as “a management strategy that uses electronic information and other technologies to gather, process and analyze spatial and temporal data for the purpose of guiding targeted actions that improve the efficiency, productivity and sustainability of agricultural operations” [3]. Moreover, Agriculture 4.0 considers the introduction of robots in agriculture, or agrobots for short [4]. Reviews of state-of-the-art agrobots can be found in [5,6]. In particular, agrobots are developed in viniculture regarding monitoring [7], harvesting [8], pruning [9] and cluster detection [10]. An approach regarding two cooperative agrobots toward grape harvesting has been reported [11]. The interest here is in “personalized” grape harvesting as explained next.
To date, grape harvest has been carried out as outlined in the following. Human experts periodically sample grapes from a vineyard, analyze them chemically and, finally, decide when the harvest should commence. In conclusion, during harvest, human workers, either manually or mechanically, harvest all the grapes en masse regardless of their degree of maturity. Aiming at a higher quality of final produce [12], the harvest should be “personalized” in the sense that it is carried out conditionally on a per-grape-bunch basis to ensure that only ripe grapes are harvested. Note that such a method of personalized grape harvesting is especially advantageous for wine-making.
This paper presents a logic-based tunable methodology for personalized grape harvesting by an autonomous agrobot. Preliminary results have been presented in [13]. The significant enhancements and novelties of this work are summarized in the following. Real-world data are used here for regression, whereas artificial data were used in [13] for classification. Two parametrically optimized FLRule methods are employed here comparatively; moreover, in this current work, structurally optimized positive valuation functions are used. The work in [13] has comparatively demonstrated results by a conventional 4-layer backpropagation neural network, as well as a kNN with k = 1; whereas, the work here comparatively demonstrates results by a multilayer CNN, as well as a kNN with k > 1. A new theorem proof is also presented here.
Agrobots often need to take a binary decision, such as harvest (or not harvest) or prune (or not prune), conditioned on a number of constraints that need to be satisfied simultaneously, as well as adequately. The aim of this work is to propose a tunable rule-based method for explainable decision-making, under constraints, based on logic. The aforementioned aim was pursued by the development of logic-based decision-making, applicable on a mathematical lattice ontology of constraints. A specific agrobot application was dealt with, namely grape harvesting, based on either two or three constraints. Two decision-making models are presented below. Comparative experimental results on real-world data have been encouraging toward further extensions.
The paper is structured as follows: Section 2 describes the materials and methods including the physical problem, the datasets employed, mathematical background and computational algorithms. Section 3 presents computational experiments and results, as well as a discussion. Finally, Section 4 concludes by summarizing the contribution and proposing future work extensions.

2. Materials and Methods

This section presents a summary of the physical problem, followed by a description of the datasets used below. In addition, it summarizes the mathematical background followed by a presentation of the computational algorithms employed below.

2.1. The Physical Problem

The requirement of an autonomous agrobot to harvest only ripe grape bunches involves estimation of the ripeness of each grape in the vineyard. More specifically, based on criteria established by experts, a grape’s ripeness is decided by the concentration of sugars and acids in the berry, as well as factors related to aromatic maturity, soluble content, etc., [14]. The datasets employed in this work are described next.

2.2. The Datasets

A grape’s maturity is decided by the values of three chemical indices, namely total soluble solids (TSS), titratable acidity (TA) and pH, which specify the acidity and basicity of the grape. More specifically, the index TSS specifies the sugar content in grape juice in degrees (°) Brix, which is a measure of how many grams of sugars are present per 100 g of juice; levels between around 18 and 24 °Brix indicate maturity, depending on both the grape variety and the desired wine style. The index TA is a measure of the acid content in grape juice. Organic acids are responsible for the sour taste of the wine and influence wine stability, color and pH. Acid concentration increases in the berry during the early period of berry growth, whereas it decreases at veraison, reaching as low as 0.6–0.8 g of titratable acids/100 mL (%TA) at harvest. Finally, the index pH is a measure of active acidity in the juice and wine, and it affects a range of factors including microbial and physical stability, oxidation level, color and flavor. Desirable pH values at the time of harvest vary depending on the berry type. In general, white grapes are harvested with a pH in the range of 3.1 to 3.3 and red grapes with a pH in the range of 3.3 to 3.5 [15].
Berries were collected in the year 2019 from three commercial vineyards located in the Eastern Macedonia region of Greece, and more specifically in Kavala (40°48′43″ N, 23°59′25″ E) for the Xinomavro variety in an area of 1.0 hectare (ha) and in two locations in Drama, (41°5.8′ N, 23°56.7′ E) and (41°5.5′ N, 23°55.80′ E), for the Syrah and Sauvignon Blanc varieties in areas of 2.5 and 1.7 ha, respectively. The three cultivars trained on a similar bilateral cordon and planting at a density of 3333 vines/ha (2.5 m between rows × 1.2 m within rows). Berry samplings took place on a weekly basis following veraison and continued until harvest depending on the grape variety. For each variety, 5 berries per vine were collected from 40 vines randomly within the vineyard, resulting in an amount of 200 berries per sampling date. The berries were pressed, and the must was chemically analyzed for soluble solids (Brix) by refractometry (HI96841, HANNA), titratable acidity (g L-1 tartaric acid) and pH (HI2020-02, HANNA). The degree of maturity was estimated by experts, and it was represented by a real number in the closed interval [0,1], with 1 indicating a 100% mature grape.
The datasets acquired and used in this work are available via the link at the end of the paper. Table 1 displays the ranges of the three indices for each of the three cultivars together with the corresponding unit, as well as the corresponding number of data records acquired per grape cultivar. Note that different grape varieties were developed in different vineyards by different owners who applied different sampling practices; the latter explains any differences in the recorded dataset sizes.

2.3. Mathematical Background

The mathematical background has been presented in [13]. For the reader’s convenience, this section repeats Definition 1 as well as Theorem 1. Τhe proof of Theorem 1 is introduced here in the Appendix A.
Definition 1.
Let (L,⊑) be a lattice. A function σ: L × L→[0,1] is called an inclusion measure if and only if the following equivalences hold: σ(u,w) = 1 ⇔ u ⊑ w ⇔ σ(x,u) ≤ σ(x,w).
Any use of an inclusion measure is called fuzzy lattice reasoning, or FLR for short. An inclusion measure supports two different modes of reasoning, namely generalized modus ponens and reasoning by analogy. Two inclusion measures can be defined by σ(x,u) = v(u)/v(xu) and σ(x,u) = v(xu)/v(x), namely “sigma join” and “sigma meet”, respectively, where the real function v: L→ℝ is a positive valuation in lattice (L,⊑).
It is known that if function vi: Li→ℝ is a positive valuation in lattice (Li,⊑i), i∈{1,…,N}, then the function v = v1 + … + vN is a positive valuation in the Cartesian product lattice (L,⊑) = (L1,⊑1)×…×(LN,⊑N), where a lattice (Li,⊑i), i∈{1,…,N}, is called constituent (lattice). A constituent lattice can be interpreted as a data dimension.
Theorem 1.
Let function σi: Li × Li→[0,1] be an inclusion measure in a constituent lattice (Li,⊑i), i∈{1,…,N}. Consider the Cartesian product lattice (L,⊑) = (L1,⊑1)×…×(LN,⊑N), where L = L1 × L2 × … × LN and ⊑ = ⊑1 × ⊑2 × … × ⊑N. Then, given u = (u1,…,uN), w = (w1,…,wN) and all λi ≥ 0 such that i = 1 N λ i = 1 , an inclusion measure function σ: L × L→[0,1] can be defined as either (A1) σΣ(u,w) = i = 1 N λ i σ i u i , w i  or (A2) σΠ(u,w) = i = 1 N σ i u i , w i .
The proof of Theorem 1 is presented in Appendix A.
A constituent lattice (Li,⊑i), i ∈ {1,…,N}, here was the lattice (ℝ,≤) or real numbers. Hence, a positive valuation function vi: ℝ→ℝ was a strictly increasing real function. Moreover, the strictly decreasing function θ(x)= −x was used here exclusively.

2.4. Computational Algorithms

Consider the following set of inequalities (constraints) of chemical indices related to grape maturity: a: TSSmin ≤ TSS ≤ TSSMAX, b: TAmin ≤ TA ≤ TAMAX, c: pHmin ≤ pH ≤ pHMAX. The binary degree of satisfaction of those inequalities gives rise to a crisp Boolean lattice represented in Figure 1 by a Hasse diagram.
The bounds of the inequalities in Figure 1 depend on the grape cultivar such that, within them, grape maturity is maximized. More specifically, for Xinomavro, those bounds were [TSSmin, TSSMAX,] = [17.9, 24.7], [TAmin, TAMAX] = [4.1, 9.75] and [pHmin, pHmin] = [2.98, 3.42]; for Syrah, they were [TSSmin, TSSMAX,] = [23.6, 26.4] and [TAmin, TAMAX] = [3.6, 5.6]; whereas, for Sauvignon Blanc, they were [TSSmin, TSSMAX,] = [16.6, 21.6] and [TAmin, TAMAX] = [4.8, 8.1].
We remark that the Hasse diagram in Figure 1 is interpreted here as an ontology. Note that, from a data-processing point of view, ontologies have already been engaged in Industry 4.0, as well as in agricultural applications [16]. In this work, the crisp lattice in Figure 1 is fuzzified by a parametrically optimizable inclusion measure function σ(.,.). In conclusion, fuzzy lattice reasoning (FLR) is applied in the context of the lattice-computing (LC) information-processing paradigm [17,18].
Algorithms 1 and 2 from [13] were used here for regression. In particular, the mean (absolute error) in maturity estimation was computed as Q/ntst, where Q is the sum of absolute errors; moreover, ntst is the number of testing data x ¯ i , r i , i∈{1,…,ntst}. More specifically, Figure 2 explains in a flow chart the basic decision-making by the FLRule. In particular, partial decision-making is carried out by calculating the degree of satisfaction of a constraint by an inclusion measure function per grape maturity index. Then, the partial decisions are fused toward estimating the grape maturity. In conclusion, a binary decision, i.e., “harvest” or “not harvest”, is taken conditioned on a user-defined threshold toward a discriminative grape harvest. Note that the basic decision-making explained in Figure 2 is part of both the training phase toward optimal parameter estimation and the testing phase toward action.

3. Computational Experiments and Results

The dataset for each grape cultivar was partitioned randomly into 80% training data and 20% testing data. The mean error and standard deviation were calculated for 10 random permutations of the dataset for each cultivar. The mean error is defined as the average value difference between the actual maturity prediction (as defined in the dataset by the experts) and the calculated prediction. The inclusion measure “sigma join”, i.e., σ(x,u) = v(u)/v(xu), was used here exclusively because it calculates estimates beyond intervals.
Initial experiments were aimed at computing an optimal, parametric positive valuation function v(x). Alternative function types were tried under the constraint that v(x) must be a monotonically increasing function. More specifically, linear functions v(x) = ax + b, where a,b ∈ℝ with a > 0, were tried, as well as sigmoid functions A / 1 + e k x + x 0 , where A,k,x0 ∈ℝ with A,k > 0. Furthermore, convex sums v(x) including two, three and four sigmoid functions were tried denoted as sigmoid2, sigmoid3 and sigmoid4, respectively.
To estimate the values of the parameters that yield the most accurate maturity prediction, given values for TSS, TA and pH, a genetic algorithm (GA) was engaged. Figure 3 shows the structure of an individual chromosome in the genetic algorithm. The number of genes in each chromosome depends on both the cultivar and the chosen parametric function, the latter in turn dictates the number of parameters needed.
The algorithm included certain constraints in order to guarantee that the parameter values (individuals) yielded in each generation result in a ν(x) function that is strictly monotonically increasing. In addition, further constraints were applied to each individual so that the λ parameters were always in the range [0, 1]. For each optimization experiment, 1000 individuals were initially generated to serve as the initial population. The genetic algorithm was executed for 1000 generations. At the end of the optimization process, the parameters of the best individual were extracted and applied to the corresponding equations in order to determine the maturity percentage of the grapes, given three measurements for the TSS, TA and pH indices.
Figure 4 displays the prediction accuracy of an FLRule method for different v(x) functions for the grape cultivar Xinomavro.
The results of Figure 4 indicate that the best prediction accuracy was achieved using a v(x) parametric function involving the sum of two sigmoid functions. It can be seen that using more sigmoid functions for the definition of v(x) (and therefore increasing the number of tunable parameters) does not improve the prediction accuracy. Therefore, for the remainder of this section, where prediction accuracies for the FLRule (σΣ) and FLRule (σ) methods are presented, they are derived from calculations using a v(x) parametric function involving the sum of two sigmoid functions.
Having established the optimal parametric function v(x) for this problem, computational experiments for estimating maturity by the FLRule (σΣ) and FLRule (σ) methods for the three cultivars, Xinomavro, Syrah and Sauvignon Blanc, were carried out. Comparative results by both the convolutional neural network (CNN) and a kNN regressor are also demonstrated. On one hand, the CNN (sequential model) consisted of either two inputs (Xinomavro) or three inputs (Syrah, Sauvignon Blanc) representing the maturity indices. Furthermore, the CNN included one input layer, six hidden layers (500 neurons per layer with relu activation functions) and a single output layer (with a sigmoid activation function). On the other hand, the kNN regressor used k = 3 (for Xinomavro) and k = 2 (for Syrah and Sauvignon Blanc).
Two types of computational experiments were carried out including, first, trivial interval inputs and, second, non-trivial interval inputs. More specifically, first, a measured value “x” of a maturity index was represented by the trivial interval [x, x] for an FLRule method; second, a measured value “x” of a maturity index was augmented to the non-trivial interval [x − δ, x + δ] for an FLRule method in order to accommodate uncertainty/ambiguity. An arbitrary value of δ = 1 was selected in this work for demonstration reasons. Only an FLRule method can deal with non-trivial interval inputs with consistency.
Table 2 displays the results regarding the grape cultivar Xinomavro. For trivial interval (point) inputs, Table 2 shows that for the Xinomavro cultivar, the prediction accuracy of the FLRule (σΣ) and FLRule (σ) methods was comparable to that achieved with the CNN and the kNN regressor, when the testing set is considered. Additional experiments were carried out for non-trivial interval inputs, as detailed above, only for the methods FLRule (σΣ) and (σ).
Similar experiments were carried out for the grape cultivar Syrah and the results are presented in Table 3. In this case, calculations were performed involving only the TSS and TA maturity indices, as there were no data available for the index pH. The results show that while the FLRule (σΣ) and FLRule (σ) methods achieve a similar prediction accuracy to the CNN, the kNN regressor is noticeably more accurate. This can be explained by the significantly smaller size of the dataset for this cultivar (N = 85), as well as by the absence of one of the dimensions (pH), factors which make the kNN regressor more suitable. Additional experiments were carried out for non-trivial interval inputs, as detailed above, only for the methods FLRule (σΣ) and (σ).
Finally, the comparative maturity prediction for the grape cultivar Sauvignon Blanc is shown in Table 4. As was the case with the cultivar Syrah, there were no values available for the pH index, and so only TSS and TA were considered as inputs. Here, the FLRule (σΣ) method appears to be significantly better when compared to either the FLRule (σ) method or the CNN, and it is comparable with the kNN regressor. Additional experiments were carried out for non-trivial interval inputs, as detailed above, only for the methods FLRule (σΣ) and (σ).
Computational experiments were also carried out using the inclusion measure σΠ with the FLRule. Nevertheless, the results were considerably worse due to the calculation of inclusion measure σΠ as a product of inclusion measures per constituent lattice.

Discussion

The results presented in the previous section have demonstrated the effectiveness of both the FLRule (σΣ) and FLRule (σ) methods in estimating grape maturity on real-world data compared with traditional machine-learning algorithms including the CNN and kNN regressors.
Maturity estimation was pursued here based on logic, namely fuzzy lattice reasoning (FLR), techniques. In particular, an optimizable inclusion measure function σ(.,.) was used based on a parametric strictly increasing real function v: ℝ→ℝ per data dimension. In the first place, extensive computational experiments were carried out using the Xinomavro grape cultivar dataset, because it was the largest one, toward selecting an optimal positive valuation function v(x) per data dimension. Two different types (structures) of strictly increasing functions v(.) were considered, namely linear and sigmoid. The clearly superior performance of the sigmoid function was attributed to the fact that a sigmoid function has more tunable parameters, namely three, compared to a linear function which has only two tunable parameters. Next, experiments were carried out by considering convex sums of sigmoid functions in order to further increase performance by increasing the number of tunable parameters. The best results were obtained for the sum of two sigmoid functions; whereas, for more than two, the performance slightly deteriorated as demonstrated in Figure 4. The latter deterioration was attributed to the exponential increase in the search space due to an increase in parameters per data dimension.
Using the optimal sigmoid2 positive valuation per data dimension, selected as explained above, computational experiments were carried out to estimate the grape maturity prediction accuracy by the proposed methods for the three cultivars Xinomavro, Syrah and Sauvignon Blanc.
For the Xinomavro cultivar, there were measurements for all three indices TSS, TA and pH available in the dataset. In addition, more data records were available for this cultivar (N = 825), compared to the records available for the other two cultivars. Training the FLRule (σΣ) and FLRule (σ) methods yielded maturity predictions comparable with the ones produced by the CNN and the kNN regressor.
For the Syrah and the Sauvignon Blanc grape cultivars, values for only the TSS, TA indices were available. Moreover, the training datasets for those cultivars are much smaller than that of the Xinomavro cultivar. In conclusion, for the Syrah cultivar (N = 85 samples), the experimental results have demonstrated that the FLRule (σΣ), FLRule (σ) and CNN methods perform well in estimating maturity, with errors similar to those observed for the Xinomavro cultivar. However, in this case, the kNN regressor appears to be more accurate than the other methods. Note also that the standard deviation for the CNN is significantly larger than that of other methods, which implies that the CNN’s performance in this problem is not stable, possibly because more data are required for training. Finally, for the Sauvignon Blanc (N = 96 samples), the FLRule (σΣ) method yielded a comparable accuracy to the kNN regressor. The accuracy of the FLRule (σ) method deteriorated significantly between the training and testing sets. The CNN performed slightly worse than either the FLRule (σΣ) or the kNN regressor.
It is understood that in certain applications, especially real-world applications, few data are available. The proposed FLRule method here has demonstrated effectiveness in dealing with few data.
Regarding the accuracy errors observed during the computational experiments, it has to be pointed out that, in practice, after grape maturity estimation, the agrobot needs to take a binary decision, that is to “harvest” or to “not harvest” a grape conditioned on a user-defined threshold. Therefore, an approximation of the grape’s maturity with such a small error, as demonstrated by the above experiments, was acceptable to experts.
Two additional advantages of an FLRule method include, first, its capacity to data process intervals that accommodate uncertainty/ambiguity and, second, its capacity to explain its answers by a rule; the latter is the set of inequalities a, b and c in Figure 1.

4. Conclusions

During harvest, a personalized binary decision is taken repeatedly regarding every fruit, namely “harvest” or “not harvest”. Such a decision is conditioned on a number of constraints that need to be satisfied simultaneously, as well as adequately. The aforementioned constraints were represented here by a Boolean lattice ontology that defines a rule; the latter is tunable by a parametric positive valuation function.
Two different FLRule methods for regression by reasoning have been presented based on the functions σΣ and σ, respectively. In particular, the parametric FLRule (σΣ) and (σ) methods were engaged for estimating the maturity of grapes for the purpose of harvesting. The proposed FLRule (σΣ), as well as FLRule (σ), methods have demonstrated a good performance in predicting maturity in comparison to alternative machine-learning methods, namely the CNN and the kNN classifier. Several comparative advantages of an FLRule method have been discussed.
An FLRule method here was based on three grape maturity indices, namely TSS, TA and pH. Future work will consider including more maturity indices, as well as applications to a variety of agricultural tasks such as spraying, pruning and others. An estimation of maturity index values by machine-vision techniques will also be pursued by extending preliminary works [19,20]. Future work will also seek further improvements by including all-order data statistics represented by intervals’ numbers [21].

Author Contributions

Conceptualization, V.G.K. and C.L.; methodology, V.G.K.; software, C.L.; validation, C.L. and T.P.; formal analysis, V.G.K. and S.T.; resources, S.T. and G.S.; data curation, S.T. and E.M.; writing—original draft preparation, C.L., V.G.K. and T.P.; writing—review and editing, C.L., V.G.K., T.P., G.S., S.T. and E.M.; visualization, C.L., and G.S.; supervision, V.G.K. and T.P.; project administration, V.G.K. and E.M.; funding acquisition, V.G.K. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge support of this work by the project “Technology for Skillful Viniculture (SVtech)” (MIS 5046047) which is implemented under the action “Reinforcement of the Research and Innovation Infrastructure” funded by the Operational Program “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014–2020) and co-financed by Greece and the European Union (European Regional Development Fund).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available at https://github.com/lytridic/Grape-Maturity-Dataset (accessed on 26 March 2023) under a Creative Commons Attribution 4.0 international license.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The proof of Theorem 1 is presented next.
Proof of Theorem 1.
Let function σi: Li × Li→[0,1] be an inclusion measure in a constituent lattice (Li,⊑i), i ∈ {1,…,N}. Consider the Cartesian product lattice (L,⊑) = (L1,⊑1)×…×(LN,⊑N), where L = L1 × L2 × … × LN and ⊑ = ⊑1 × ⊑2 × … × ⊑N. Then, given u = (u1,…,uN), w = (w1,…,wN) and all λi ≥ 0 such that i = 1 N λ i = 1 , an inclusion measure function σ: L × L→[0,1] can be defined as either (A1) σΣ(u,w) = i = 1 N λ i σ i u i , w i or (A2) σΠ(u,w) = i = 1 N σ i u i , w i .
Proof. 
(A1) Let u,w,x ∈ L with u = (u1,…,uN), w = (w1,…,wN) and x = (x1,…,xN).
  • In the one direction, we prove that uw ⇒ σΣ(u,w) = 1. More specifically, uwuiwi, i ∈ {1,…,N} ⇒ σi(ui,wi) = 1, i ∈ {1,…,N} ⇒ i = 1 N λ i σ i u i , w i = σΣ(u,w) = 1.
    In the other direction, we prove that uw ⇐ σΣ(u,w) = 1. More specifically, σΣ(u,w) = 1 ⇒ i = 1 N λ i σ i u i , w i = 1 ⇒ σ(ui,wi) = 1, i ∈ {1,…,N} ⇒ uiwi, i∈{1,…,N} ⇒ uw.
  • uwuiwi, i∈{1,…,N} ⇔ σi(xi,ui) ≤ σi(xi,wi), i∈{1,…,N} ⇔ i = 1 N λ i σ i x i , u i i = 1 N λ i σ i x i , w i ⇔ σΣ(x,u) ≤ σΣ(x,w).
(A2) Let u,w,x ∈ L with u = (u1,…,uN), w = (w1,…,wN) and x = (x1,…,xN).
  • In the one direction, we prove that uw ⇒ σΠ(u,w) = 1. More specifically, uwuiwi, i∈{1,…,N} ⇒ σi(ui,wi) = 1, i ∈ {1,…,N} ⇒ i = 1 N σ i u i , w i = σΠ(u,w) = 1.
    In the other direction, we prove that uw ⇐ σΠ(u,w) = 1. More specifically, σΠ(u,w) = 1 ⇒ i = 1 N σ i u i , w i = 1 ⇒ σ(ui,wi) = 1, i∈{1,…,N} ⇒ uiwi, i∈{1,…,N} ⇒ uw.
  • uwuiwi, i∈{1,…,N} ⇔ σi(xi,ui) ≤ σi(xi,wi), i∈{1,…,N} ⇔ i = 1 N σ i x i , u i i = 1 N σ i x i , w i ⇔ σΠ(x,u) ≤ σΠ(x,w).
                                                                                                                                                                   □

References

  1. United Nations World Population Projected to Reach 9.8 Billion in 2050. Available online: https://www.un.org/development/desa/en/news/population/world-population-prospects-2017.html (accessed on 12 December 2022).
  2. De Clercq, M.; Vatz, A.; Biel, A. Agriculture 4.0: The Future of Farming Technology. In Proceedings of the World Government Summit, Dubai, United Arab Emirates, 11–13 February 2018; pp. 4–25. [Google Scholar]
  3. Lowenberg-DeBoer, J.; Erickson, B. Setting the Record Straight on Precision Agriculture Adoption. Agron. J. 2019, 111, 1552–1569. [Google Scholar] [CrossRef]
  4. Bechar, A. Innovation in Agricultural Robotics for Precision Agriculture; Progress in Precision Agriculture; Bechar, A., Ed.; Springer International Publishing: Cham, Switzerland, 2021; ISBN 978-3-030-77035-8. [Google Scholar]
  5. Fountas, S.; Mylonas, N.; Malounas, I.; Rodias, E.; Hellmann Santos, C.; Pekkeriet, E. Agricultural Robotics for Field Operations. Sensors 2020, 20, 2672. [Google Scholar] [CrossRef] [PubMed]
  6. Oliveira, L.F.P.; Moreira, A.P.; Silva, M.F. Advances in Agriculture Robotics: A State-of-the-Art Review and Challenges Ahead. Robotics 2021, 10, 52. [Google Scholar] [CrossRef]
  7. Roure, F.; Moreno, G.; Soler, M.; Faconti, D.; Serrano, D.; Astolfi, P.; Bardaro, G.; Gabrielli, A.; Bascetta, L.; Matteucci, M. GRAPE: Ground Robot for VineyArd Monitoring and ProtEction. In Advances in Intelligent Systems and Computing; Springer: Berlin/Heidelberg, Germany, 2018; Volume 693, pp. 249–260. ISBN 9783319708331. [Google Scholar]
  8. Vrochidou, E.; Tziridis, K.; Nikolaou, A.; Kalampokas, T.; Papakostas, G.A.; Pachidis, T.P.; Mamalis, S.; Koundouras, S.; Kaburlasos, V.G. An Autonomous Grape-Harvester Robot: Integrated System Architecture. Electronics 2021, 10, 1056. [Google Scholar] [CrossRef]
  9. Botterill, T.; Paulin, S.; Green, R.; Williams, S.; Lin, J.; Saxton, V.; Mills, S.; Chen, X.; Corbett-Davies, S. A Robot System for Pruning Grape Vines. J. Field Robot. 2017, 34, 1100–1122. [Google Scholar] [CrossRef]
  10. Luo, L.; Tang, Y.; Zou, X.; Wang, C.; Zhang, P.; Feng, W. Robust Grape Cluster Detection in a Vineyard by Combining the AdaBoost Framework and Multiple Color Components. Sensors 2016, 16, 2098. [Google Scholar] [CrossRef] [PubMed]
  11. Lytridis, C.; Bazinas, C.; Pachidis, T.; Chatzis, V.; Kaburlasos, V.G. Coordinated Navigation of Two Agricultural Robots in a Vineyard: A Simulation Study. Sensors 2022, 22, 9095. [Google Scholar] [CrossRef]
  12. Bac, C.W.; van Henten, E.J.; Hemming, J.; Edan, Y. Harvesting Robots for High-Value Crops: State-of-the-Art Review and Challenges Ahead. J. Field Robot. 2014, 31, 888–911. [Google Scholar] [CrossRef]
  13. Kaburlasos, V.; Lytridis, C.; Siavalas, G.; Pachidis, T.; Theocharis, S. Fuzzy Lattice Reasoning (FLR) for Decision-Making on an Ontology of Constraints toward Agricultural Robot Harvest. In Machine Learning, Multi Agent and Cyber Physical Systems; World Scientific: Tianjin, China, 2023; pp. 80–87. [Google Scholar]
  14. Rabot, A. Using of a Combined Approach by Biochemical and Image Analysis to Develop a New Method to Estimate Seed Maturity Stage for Bordeaux Area Grapevine. OENO One 2017, 51, 29–35. [Google Scholar] [CrossRef]
  15. Dami, I. Determining Grape Maturity and Fruit Sampling. OHIO STATE Univ. Hortic. Crop Sci. OARDC 2009, 1–4. [Google Scholar] [CrossRef]
  16. Sampath Kumar, V.R.; Khamis, A.; Fiorini, S.; Carbonera, J.L.; Olivares Alarcos, A.; Habib, M.; Goncalves, P.; Li, H.; Olszewska, J.I. Ontologies for Industry 4.0. Knowl. Eng. Rev. 2019, 34, e17. [Google Scholar] [CrossRef]
  17. Kaburlasos, V.G. The Lattice Computing (LC) Paradigm. In Proceedings of the 15th International Conference on Concept Lattices and Their Applications (CLA 2020), Tallinn, Estonia, 29 June–1 July 2020; Volume 2668, pp. 1–7. [Google Scholar]
  18. Sussner, P.; Caro Contreras, D.E. Generalized Morphological Components Based on Interval Descriptors and N-Ary Aggregation Functions. Inf. Sci. 2022, 583, 14–32. [Google Scholar] [CrossRef]
  19. Bazinas, C.; Vrochidou, E.; Kalampokas, T.; Karampatea, A.; Kaburlasos, V.G. A Non-Destructive Method for Grape Ripeness Estimation Using Intervals’ Numbers (INs) Techniques. Agronomy 2022, 12, 1564. [Google Scholar] [CrossRef]
  20. Chariskou, C.; Bazinas, C.; Daniels, A.J.; Opara, U.L.; Nieuwoudt, H.H.; Kaburlasos, V.G. Variable Selection for the Prediction of TSS, PH and TA of Intact Berries of Thompson Seedless Grapes from Their NIS Reflection. In Proceedings of the 2022 International Conference on Software, Telecommunications and Computer Networks (SoftCOM), Split, Croatia, 22–24 September 2022; IEEE: Piscataway, NJ, USA, 2022; pp. 1–5. [Google Scholar]
  21. Lytridis, C.; Lekova, A.; Bazinas, C.; Manios, M.; Kaburlasos, V.G. WINkNN: Windowed Intervals’ Number KNN Classifier for Efficient Time-Series Applications. Mathematics 2020, 8, 413. [Google Scholar] [CrossRef]
Figure 1. A Hasse diagram representing a Boolean lattice ontology of the three inequality constraints a: TSSmin ≤ TSS ≤ TSSMAX, b: TAmin ≤ TA ≤ TAMAX, c: pHmin ≤ pH ≤ pHMAX, regarding the grape maturity indices TSS, TA and pH.
Figure 1. A Hasse diagram representing a Boolean lattice ontology of the three inequality constraints a: TSSmin ≤ TSS ≤ TSSMAX, b: TAmin ≤ TA ≤ TAMAX, c: pHmin ≤ pH ≤ pHMAX, regarding the grape maturity indices TSS, TA and pH.
Sustainability 15 07331 g001
Figure 2. Basic decision-making by the FLRule for the maturity index values TSS, TA and pH. Partial decision-making is carried out by calculating the inclusion measure function σ(.,.) values σ([TSS, TSS]⊑[TSSmin, TSSMAX]), σ([TA, TA]⊑[TAmin, TAMAX]) and σ([pH, pH]⊑[pHmin, pHmin]). Then, the partial decisions are fused, by either σΣ(,.) or σ(.,.) toward grape maturity estimation. In conclusion, a binary decision, i.e., “harvest” or “not harvest”, is taken conditioned on a user-defined threshold.
Figure 2. Basic decision-making by the FLRule for the maturity index values TSS, TA and pH. Partial decision-making is carried out by calculating the inclusion measure function σ(.,.) values σ([TSS, TSS]⊑[TSSmin, TSSMAX]), σ([TA, TA]⊑[TAmin, TAMAX]) and σ([pH, pH]⊑[pHmin, pHmin]). Then, the partial decisions are fused, by either σΣ(,.) or σ(.,.) toward grape maturity estimation. In conclusion, a binary decision, i.e., “harvest” or “not harvest”, is taken conditioned on a user-defined threshold.
Sustainability 15 07331 g002
Figure 3. Structure of a chromosome containing the parameters represented as genes.
Figure 3. Structure of a chromosome containing the parameters represented as genes.
Sustainability 15 07331 g003
Figure 4. Maturity prediction accuracy for the grape cultivar Xinomavro depending on v(x) selection for the FLRule (σΣ) and FLRule (σ) methods.
Figure 4. Maturity prediction accuracy for the grape cultivar Xinomavro depending on v(x) selection for the FLRule (σΣ) and FLRule (σ) methods.
Sustainability 15 07331 g004
Table 1. Ranges of maturity indices for the three grape cultivars.
Table 1. Ranges of maturity indices for the three grape cultivars.
Grape CultivarMaturity
Index
UnitRangeNo. of Data Records
XinomavroTSSDegrees Brix (°Bx)9.2–24.7
TAg/L4–24.6825
pH-2.35–3.42
SyrahTSSDegrees Brix (°Bx)10–26.4
TAg/L3.6–15.685
pH--
Sauvignon BlancTSSDegrees Brix (°Bx)9.1–21.6
TAg/L4.8–21.0796
pH--
Table 2. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding grape cultivar Xinomavro for both trivial and non-trivial (interval) inputs.
Table 2. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding grape cultivar Xinomavro for both trivial and non-trivial (interval) inputs.
Trivial Interval InputsNon-Trivial Interval Inputs
MethodMeanStdMeanStd
FLRule (σΣ)Training0.0540.0020.0560.001
Testing0.0570.0020.0570.001
FLRule (σ)Training0.0540.0020.0800.001
Testing0.0560.0020.0860.001
CNNTraining0.0400.012--
Testing0.0500.010--
kNNTraining0.0340.001--
Testing0.0530.004--
Table 3. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding the grape cultivar Syrah for both trivial and non-trivial (interval) inputs.
Table 3. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding the grape cultivar Syrah for both trivial and non-trivial (interval) inputs.
Trivial Interval InputsNon-Trivial Interval Inputs
MethodMeanStdMeanStd
FLRule (σΣ)Training0.0460.0040.0580.004
Testing0.0570.0020.0690.009
FLRule (σ)Training0.0440.0040.0560.006
Testing0.0560.0030.0640.010
CNNTraining0.0450.010--
Testing0.0520.020--
kNNTraining0.0160.004--
Testing0.0380.014--
Table 4. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding the grape cultivar Sauvignon Blanc for both trivial and non-trivial (interval) inputs.
Table 4. Maturity prediction statistics, i.e., mean (absolute error) and std (standard deviation), regarding the grape cultivar Sauvignon Blanc for both trivial and non-trivial (interval) inputs.
Trivial Interval InputsNon-Trivial Interval Inputs
MethodMeanStdMeanStd
FLRule (σΣ)Training0.0460.0060.0620.013
Testing0.0540.0030.1030.013
FLRule (σ)Training0.0410.0040.0490.006
Testing0.1000.0180.1010.012
CNNTraining0.0550.023--
Testing0.0720.029--
kNNTraining0.0220.003--
Testing0.0570.009--
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lytridis, C.; Siavalas, G.; Pachidis, T.; Theocharis, S.; Moschou, E.; Kaburlasos, V.G. Grape Maturity Estimation for Personalized Agrobot Harvest by Fuzzy Lattice Reasoning (FLR) on an Ontology of Constraints. Sustainability 2023, 15, 7331. https://doi.org/10.3390/su15097331

AMA Style

Lytridis C, Siavalas G, Pachidis T, Theocharis S, Moschou E, Kaburlasos VG. Grape Maturity Estimation for Personalized Agrobot Harvest by Fuzzy Lattice Reasoning (FLR) on an Ontology of Constraints. Sustainability. 2023; 15(9):7331. https://doi.org/10.3390/su15097331

Chicago/Turabian Style

Lytridis, Chris, George Siavalas, Theodore Pachidis, Serafeim Theocharis, Eirini Moschou, and Vassilis G. Kaburlasos. 2023. "Grape Maturity Estimation for Personalized Agrobot Harvest by Fuzzy Lattice Reasoning (FLR) on an Ontology of Constraints" Sustainability 15, no. 9: 7331. https://doi.org/10.3390/su15097331

APA Style

Lytridis, C., Siavalas, G., Pachidis, T., Theocharis, S., Moschou, E., & Kaburlasos, V. G. (2023). Grape Maturity Estimation for Personalized Agrobot Harvest by Fuzzy Lattice Reasoning (FLR) on an Ontology of Constraints. Sustainability, 15(9), 7331. https://doi.org/10.3390/su15097331

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

Article Metrics

Back to TopTop