Next Article in Journal
Plasma Amino Acids and Acylcarnitines Are Associated with the Female but Not Male Adolescent Swimmer’s Performance: An Integration between Mass Spectrometry and Complex Network Approaches
Previous Article in Journal
Spatiotemporal Variability of Trace Elements Fingerprints in Otoliths of Japanese Eel (Anguilla japonica) and Its Use in Tracing Geographic Origin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Artificial Intelligence-Based Robust Hybrid Algorithm Design and Implementation for Real-Time Detection of Plant Diseases in Agricultural Environments

Department of Electrical Electronics Engineering, Zonguldak Bülent Ecevit University, Zonguldak 67100, Turkey
*
Author to whom correspondence should be addressed.
Biology 2022, 11(12), 1732; https://doi.org/10.3390/biology11121732
Submission received: 29 October 2022 / Revised: 24 November 2022 / Accepted: 25 November 2022 / Published: 29 November 2022

Abstract

:

Simple Summary

Plant disease, defined as an abnormal condition that disrupts the normal growth of the plant, is one of the main causes of economic losses in the agricultural industry. Early diagnosis of plant disease is critical to increasing agricultural crop productivity. In this paper, a new robust hybrid classification model based on swarm optimization-supported feature selection, including machine learning and deep learning algorithms, that allows real-time classification of diseases in apple, grape, and tomato plants has been developed. In this way, it will be possible to diagnose the plant disease at an early phase and apply the appropriate treatment.

Abstract

The early detection and prevention of plant diseases that are an important cause of famine and food insecurity worldwide are very important for increasing agricultural product productivity. Not only the early detection of the plant disease but also the determination of its type play a critical role in determining the appropriate treatment. The fact that visual inspection, which is frequently used in determining plant disease and types, is tiring and prone to human error, necessitated the development of algorithms that can automatically classify plant disease with high accuracy and low computational cost. In this study, a new hybrid plant leaf disease classification model with high accuracy and low computational complexity, consisting of the wrapper approach, including the flower pollination algorithm (FPA) and support vector machine (SVM), and a convolutional neural network (CNN) classifier, is developed with a wrapper-based feature selection approach using metaheuristic optimization techniques. The features of the image dataset consisting of apple, grape, and tomato plants have been extracted by a two-dimensional discrete wavelet transform (2D-DWT) using wavelet families such as biorthogonal, Coiflets, Daubechies, Fejer–Korovkin, and symlets. Features that keep classifier performance high for each family are selected by the wrapper approach, consisting of the population-based metaheuristics FPA and SVM. The performance of the proposed optimization algorithm is compared with the particle swarm optimization (PSO) algorithm. Afterwards, the classification performance is obtained by using the lowest number of features that can keep the classification performance high for the CNN classifier. The CNN classifier with a single layer of classification without a feature extraction layer is used to minimize the complexity of the model and to deal with the model hyperparameter problem. The obtained model is embedded in the NVIDIA Jetson Nano developer kit on the unmanned aerial vehicle (UAV), and real-time classification tests are performed on apple, grape, and tomato plants. The experimental results obtained show that the proposed model classifies the specified plant leaf diseases in real time with high accuracy. Moreover, it is concluded that the robust hybrid classification model, which is created by selecting the lowest number of features with the optimization algorithm with low computational complexity, can classify plant leaf diseases in real time with precision.

1. Introduction

Due to the rapid increase in the world population and the current population size, efficient production of agricultural products is needed, as well as the efficient use of agricultural lands, which are limited. In the efficient production of agricultural products, it is critical to detect plant diseases at an early stage, as well as spraying, fertilizing, and weed detection processes. It is known that losses in agricultural product yield due to plant diseases range between 20% and 40%. This heavy yield loss also means a reduced market volume for buyers of that product [1]. A delay in the diagnosis of viral, pathogenic, or plague-borne plant diseases requires the application of more doses of pesticides to the diseased plant, thus resulting in a decrease in crop quality [2]. Detecting and correctly classifying plant diseases at an early stage not only contributes to improving the quality of agricultural products but also allows for the reduction of undesirable chemical spray applications such as fungicides and herbicides [3]. In this study, we propose a robust hybrid plant disease classification model based on metaheuristic optimization-assisted feature selection that includes machine learning and deep learning algorithms for real-time early detection of diseases in apple, grape, and tomato plants. The performance of the model embedded in the NVIDIA Jetson Nano developer kit on an unmanned aerial vehicle (UAV) has been experimentally tested on apple, grape, and tomato plants.
Plant disease, which is one of the main causes of economic losses in the agricultural industry in the world, is defined as an abnormal state that disrupts the normal growth of the plant [4]. Symptoms related to the abnormal state can usually be seen in the leaf, stem, and root parts of the plant. Leaf images are a good source of information in the classification of plant diseases. Therefore, researchers focus on the classification of plant diseases through leaf images [5]. In this study, plant diseases are classified using leaf images of apple, grape, and tomato plants.
Although traditional methods based on visual inspection are used in the detection of plant diseases in small-scale agricultural lands, it is very difficult to apply these methods in large-scale agricultural lands as they require tiring and continuous monitoring [6]. Especially in diagnosing both the species and the disease of plants with similar leaves, this process becomes more difficult and leads to visual errors. To cope with these challenges, various image processing and artificial intelligence techniques such as machine learning and deep learning are used in the real-time classification of plant diseases based on plant leaf images [7,8,9].
In the literature, there are image processing-based studies for real-time prediction of plant diseases and disease levels on many plant species. In [3], the severity level of white-tip disease (Aphelenchoides besseyi Christie), commonly seen in paddy crops, is estimated based on image processing techniques. The diseased area in the paddy crop is determined with the help of a color image segmentation algorithm based on chromatic aberration. The severity of the disease is estimated by proportioning the area of the determined region to the entire area of the leaf. Based on the estimation results obtained, a variable rate chemical spray system has been developed for the precise application of agrochemicals in real time. In [10], a ΔE (Delta E) segmentation-based feature extraction algorithm is proposed to classify diseases in images of citrus plants. The diseased region is identified by the ΔE method and hue saturation value (HSV), local binary patterns (LBP), and red green blue (RGB) histogram models are extracted from this region. The features obtained with the HSV, LBP, and RGB descriptors are combined, and a hybrid feature set is created. Principle component analysis (PCA) is used to reduce dimensionality. Fine k-nearest neighbor (kNN), cubic support vector machine (SVM), boosted tree, and bagged tree ensemble classifiers are tested. It is stated that the best performance is achieved with the bagged tree ensemble classifier. In [11], the temperature difference between the healthy and infected parts of the cucumber leaf is measured and downy mildew disease is detected using Fourier transform infrared (FTIR) spectroscopy. It is reported that the FTIR technique is effective in the pre-symptomatic detection of downy mildew disease on cucumber leaves. In [12], feature extraction and segmentation algorithms are used together for detection and classification of early scorch disease in banana leaves, fungal disease in bean leaves, sunburn disease in lemon leaves, and bacterial disease in roses. Segmentation of diseased regions has been carried out using a genetic algorithm (GA). It is also emphasized that the optimization algorithm used reduces the computational complexity.
In the image processing approach, the classification accuracy of plant diseases directly depends on the performance of feature extraction and segmentation. Especially, image segmentation plays a key role in detecting the severity of leaf disease, and the segmentation precision of leaf disease determines the accuracy of the classification of disease severity [13]. Although image processing techniques give good results in determining the severity of a disease in plant leaves, they cannot show the same success in classifying disease types in plant leaves where there is more than one disease type. Instead, artificial intelligence-based approaches, which use features with higher recognition accuracy through iterative learning without the need for specific features extracted by image processing techniques, have recently been frequently used in the classification of plant leaf diseases [14]. In [15], multiple types of tomato leaf diseases are classified using both deep learning and machine learning algorithms. To apply the machine learning models, a total of 52 texture features are extracted using a gray level co-occurrence matrix (GLCM) and local binary pattern (LBP) techniques, and 105 color features are extracted using color histograms and color moment approaches. The performances of kNN, SVM, and random forest (RF) from machine learning algorithms and AlexNet, VGG16, ResNet34, EfficientNet-b0, and MobileNetV2 architectures from deep learning algorithms are tested for 10 different types of tomato leaf disease. In [8], a set of images containing 12 plant species is separated into training and testing parts, and a convolutional neural network (CNN) is trained with the help of training data for the classification of plant diseases. It has been reported that the accuracy performance of the obtained model varies between 60% and 100% depending on the disease characteristics of each plant species. In [16], a deep learning model based on the pruned version of MobileNet architecture is proposed for efficient classification of multiple plant diseases based on leaf images of healthy and diseased plants. The proposed model is compared with VGG and classical MobileNet architectures. It is stated that the proposed architecture achieves 98.34% classification accuracy with six times fewer parameters compared to MobileNet and 29 times fewer than VGG. The performance of different optimizers is tested to reduce the computational overhead of the model. It is indicated that Adam and Nadam optimizers’ convergence rates are faster than stochastic gradient descent (SGD) and contribute to the reduction of computational cost. Besides these studies, there are many studies on the recognition and classification of leaf diseases of many plant species with artificial intelligence algorithms such as machine learning and deep learning techniques [17,18,19,20].
Although the classification accuracy of the techniques mentioned in the classification of plant leaf diseases belonging to more than one species is high, the high model complexity and high computational costs make their use difficult in real-time applications. In the real-time classification problem of plant leaf diseases, optimal values of the hyperparameters of artificial intelligence-based models built with specific architectures are determined using swarm-based optimization algorithms to overcome these challenges [21]. CNN models, which have been frequently used in the automatic diagnosis and classification of plant leaf diseases in recent years, try to solve the hyperparameter problem with various metaheuristic optimization algorithms. Although this contributes to reducing the parameter complexity of the model, it increases the computational load of the model.
This study is robust and can classify plant leaf diseases in real time with high accuracy, a low parameter complexity, and a low computational cost, using the feature-based robust hybrid classification technique, including the wrapper approach consisting of the flower pollination algorithm (FPA), SVM, and a CNN classifier. The study considers feature extraction with signal processing techniques as a preprocessing step, so that the computational load of the model is not increased. Until now, features have been extracted with many signal processing techniques in the classification of plant leaf diseases, but feature extraction has not been performed with the two-dimensional discrete wavelet transform (2D-DWT) signal processing techniques based on various wavelet families. The distinguishing features defining the characteristic features of plant leaf diseases are extracted for the first time in this study with the 2D-DWT technique using wavelet families, such as biorthogonal, Coiflets, Daubechies, Fejer–Korovkin, and symlets. The main contributions of this research can be summarized as follows:
  • A robust hybrid model based on 2D-DWT is proposed for the real-time classification of plant leaf diseases with high accuracy.
  • Feature groups are extracted for each family by applying 2D-DWT with the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families to the image dataset consisting of apple, grape, and tomato plants. The extracted feature groups for each wavelet family consist of distinctive features representing each plant leaf disease.
  • The features that keep classifier performance high for each wavelet family are selected by the wrapper approach, consisting of the population-based metaheuristic FPA and SVM algorithms. The fitness function is computed by considering both the number of features used in the model and the model’s performance in order to keep the model’s complexity and computation cost at a minimum level.
  • The efficiency of the proposed optimization algorithm is determined by comparing it with the particle swarm optimization (PSO) algorithm.
  • To overcome the model hyperparameter problem, the CNN classifier is used, which only has a classification layer without a feature extraction layer and uses the lowest number of features that can keep classification performance high.
  • For the real-time plant leaf disease classification problem, the model with the best performance is proposed, which includes the 2D-DWT signal processing method based on the “sym7” wavelet family, the wrapper approach consisting of FPA and SVM, and a CNN classifier.
  • The proposed model is embedded in the NVIDIA Jetson Nano developer kit on the UAV. Real-time classification tests have been performed on apple, grape, and tomato plants to demonstrate that the proposed model can classify plant leaf diseases in real time with high accuracy.
  • The experimental results obtained show that the model has low computational complexity and a minimum computational load; therefore, it can be used in real-time applications that require high classification accuracy.
The rest of this paper is structured as follows. The entire methodology that makes up the plant leaf disease classification model, including 2D-DWT, wavelet families, the wrapper approach consisting of the population-based metaheuristic FPA and SVM algorithms, CNN classifier, and classification performance metrics, is introduced in Section 2. The experimental results of the study, along with their detailed discussion, are presented in Section 3. Finally, Section 4 summarizes the concluding remarks and points to future work for research.

2. Framework of the Plant Diseases Detection Algorithm

This section describes our model structure that can classify leaf diseases of apple, grape, and tomato plants in real time with high accuracy. In our proposed model, features are extracted for each family by applying 2D-DWT signal processing techniques with the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families to the image dataset consisting of apple, grape, and tomato plants. The features that provide the highest model performance and the lowest model complexity are selected with the wrapper approach, consisting of the FPA and SVM algorithms. The CNN classifier, which is used to overcome the model hyperparameter problem, classifies plant leaf diseases with the help of selected features. All the methodology used in our model is described in detail in the following subsections.

2.1. Discrete Wavelet Transform

The Fourier transform (FT), which has been successfully applied for stationary signals, cannot be used for non-stationary signals. Since the spectrum changes with time, FT is insufficient to reveal the correct spectrum for non-stationary signals. Non-stationary signals are divided into sufficiently small pieces by short-term Fourier transform (STFT), and these small pieces are considered stationary. A time-frequency representation is needed to know what frequency components are present at different times and how they change as time passes. The sinusoidal frequency and phase content of a time-varying signal are found by STFT, and the relationship between frequency and time variation is defined with the help of a moving window. However, this technique is quite complicated to study non-stationary signals. Wavelet transform (WT) offers an alternative approach to overcome such challenges of FT and STFT with its variable-size windowing techniques [22,23]. The frequency and time information about a signal is obtained with this approach when short time intervals at high frequencies and long-time intervals at low frequencies can be selected. The wavelet transform is a transform technique that separates signals into different frequency components and examines each component with its resolution at that scale. Wavelets provide a good tool for time-frequency analysis [24].
Definition 1.
If  ψ ( t )  is a real-valued function whose Fourier spectrum,  ψ ( s )  satisfies the admissibility criterion.
C ψ = + | ψ ( s ) | 2 | s | d s <
then, ψ ( t ) is called a basic wavelet. Notice that, due to the s in the denominator of the integrand, it is necessary that:
ψ ( 0 ) = 0 + ψ ( t ) d t = 0
Furthermore, since ψ ( ) = 0 as well, we can see that the amplitude spectrum of an admissible wavelet is similar to the transfer function of a bandpass filter.
A set of wavelet basis functions, { ψ a , b ( t ) } , can be generated by translating and scaling the basic wavelet, ψ ( t ) , as:
ψ a , b ( t ) = 1 a   ψ ( t b a )
where  a > 0  and  b  are real numbers. The variable  a  reflects the scale of a particular basis function, while  b  specifies its translated position along the  t .
The continuous wavelet transform (CWT) of x ( t ) with respect to the wavelet ψ ( t ) is then:
W ( a , b ) = x , ψ a , b = + x ( t ) ψ a , b * ( t ) d t
The wavelet transform coefficients are given as inner products of the function being transformed with each of the basis functions.
The inverse continuous wavelet transform is mathematically defined as [25]:
f ( t ) = 1 C ψ = 0 + + W ( a , b ) ψ a , b ( t ) d b d a a 2
where C ψ is a constant of the wavelet used.
When we take the scaling and shifting wavelet parameters as dyadic variables m and n in discrete time, the discrete wavelet transform (DWT) of the signal is defined as:
W ( m , n ) = x , ψ m , n = a 0 m / 2 k f [ k ] ψ * ( k n a 0 m b 0 a 0 m )
where m and n denote frequency and time localizations, respectively. Here, a 0 > 1 , b 0 > 0 and m ,   n     Z , the scaling and shifting parameters are represented as a 0 m , n a 0 m b 0 , respectively.

2.2. Multiresolution Analysis

Not only the time domain features of the signals, but also the frequency domain features can be extracted and classified. DWT-based multi-resolution analysis is very useful in feature extraction applications from image signals. Undesirable components, such as noise and trend in the signal, can be separated by multiresolution analysis (MRA) [26].
Definition 2.
Analysis of signals in terms of both time and frequency is performed by computing scaling and wavelet functions. The scaling { φ m , n [ k ] } and wavelet { ψ m , n [ k ] } functions are given by:
φ m , n [ k ] = 2 m / 2 φ ( 2 m k n )
ψ m , n [ k ] = 2 m / 2 ψ ( 2 m k n )
where m ,   n     Z . The high-scale, low-frequency components correspond to approximation coefficient { φ [ k ] } that bridge the wavelets and filter banks and is denoted by:
φ [ k ] = k f [ k ] h [ k 2 n ]
where h [ n ] corresponds to the low pass filter that is associated with the scaling function. The low-scale high-frequency components correspond to detail coefficient { ψ [ k ] } is expressed by:
ψ [ k ] = k f [ k ] g [ k 2 n ]
where g [ n ] is the complementary high pass filter in this orthogonal filter bank. The scaling function, defined by the filter coefficients h [ n ] , provides approximation coefficients { φ [ k ] } , which are also referred to as low-pass output. The wavelet function, defined by the filter coefficients g [ n ] , provides the detailed coefficients { ψ [ k ] } , or alternatively the high-pass output.

2.3. Wavelet Families

The biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families used in the calculation of high and low pass filter coefficients of the DWT are introduced in the following subsections.

2.3.1. Biorthogonal Wavelet

Biorthogonal wavelets that are not based on vanishing moments have a compactly supported symmetrical structure. In the biorthogonal case, rather than having one scaling and wavelet function, there are two scaling functions { φ ,   φ ˜ } that may generate different MRA, and accordingly two different wavelet functions { ψ ,   ψ ˜ } . For orthogonal wavelets, recursively computations of the scaling function and mother wavelet are presented by [27]
φ ( t ) = 2 k h 0 ( k ) φ ( 2 t k )
φ ˜ ( t ) = 2 k h ˜ 0 ( k ) φ ˜ ( 2 t k )
ψ ( t ) = 2 k g 1 ( k ) φ ( 2 t k )
ψ ˜ ( t ) = 2 k g ˜ 1 ( k ) φ ˜ ( 2 t k )
where h 0 ( k ) , h ˜ 0 ( k ) , g 1 ( k ) , and g ˜ 1 ( k ) are dual filter coefficients.

2.3.2. Coiflet Wavelet

The Coiflet wavelet, designed by Ingrid Daubechies at the suggestion of Ronald Coifman, are discrete wavelets with vanishing moments and scaling functions. The idea behind Coiflet wavelets is to specify moment conditions approaching zero in the associated scaling function. This situation also causes the creation of a wavelet filter with descriptive properties [28]. The Coiflet wavelet function has 2 N moments equal to zero, and the scaling function has 2 N 1 moments equal to zero. The two functions have a support of the length of 6 N 1 . It is more symmetrical than compared to Daubechies wavelets. The general characteristic of the Coiflet wavelet is that for any given support width, it has the highest number of vanishing moments for both the scaling { φ } and the wavelet function { ψ } . The lengths of the scaling and wavelet functions for Coiflet wavelets are calculated as (L = 6, 12, 18, 24, 30). The approximation features also depend on the number of vanishing wavelet moments [29].
It is sufficient to define the Coiflets in terms of the H filter. Since for any integer n the filter z n H ( z ) generates the same MRA as H , we always assume the coefficients { h k } of H to be zero for k < 0 .
Definition 3.
Let { h j } j = 0 L 1 be the coefficients of a real quadrature mirror filter H . We say that H is a Coiflet of shift γ and moments M ,   N if the conditions in Equations (15)–(17) are satisfied:
j = 0 L 1 ( 1 ) j j k h j = 0           f o r     0 k < M
j = 0 L 1 j k h j = γ k                         f o r     0 k < N
3 M > L 1                           a n d       3 N L 1
If a Coiflet  H  also satisfies Cohen’s condition that  H  is non-zero in certain locations on the unit circle [30], its associated wavelet and scaling functions will have vanishing moments  M  and  N 1 , respectively [28]. The normalization  H ( 1 ) = 1  corresponds to Equation (16) with  k = 0 .

2.3.3. Daubechies Wavelet

Daubechies wavelets represent a wavelet summation that improves the frequency domain properties of Haar wavelets. These wavelets are a family of orthogonal wavelets indexed with N to represent the number of N vanishing wavelet moments. The Daubechies wavelets have two important properties. First, it has a finite number of non-zero p m scaling coefficients. This means that the scaling and wavelet functions are compact. The N -order Daubechies scaling function has 2 N nonzero scaling coefficients, and the support size width of the scaling and wavelet functions is in the range [ 0 ,   2 N 1   ] . The second property of the N -order Daubechies wavelet is that the initial N 1 moments of the wavelets are zero [31].
The moments of the scaling and wavelet functions are defined by:
+ φ ( k ) k j d ( k )
and:
+ ψ ( k ) k j d ( k )
respectively, where φ ( k ) is the scaling function and ψ ( k ) is the wavelet function. The moment of the Daubechies wavelets is expressed as:
+ ψ N ( k ) k j d ( k ) = 0           f o r   j = 0 , ,   N 1
The zeroing the first N moment of the ψ N wavelet is often abbreviated as ψ N , and has N zeroing moments. The zeroed moments mean that every N -order Daubechies wavelet basis function is orthogonal to all polynomials of degrees less than N . Daubechies wavelets are classified according to the number of zero moments they have. The smoothness of the scaling and wavelet functions increases in direct proportion to the number of zero moments [32].
φ ( k ) = 2 m = 0 N 1 p m φ ( 2 k m )
The wavelets ψ ( k ) are defined in terms of the scaling function. The expression relating the wavelet to the scaling function is:
ψ ( k ) = m = 0 N 1 ( 1 ) m r N 1 m φ ( 2 k m )
where N represents the order of the system, m denotes the localization parameter, p m and r N 1 m are the vector of scaling and wavelet filter coefficients, respectively.

2.3.4. Fejer–Korovkin Wavelet

The kernel function of the Fejer–Korovkin wavelet, which is more symmetrical but less smooth than the Daubechies filters, is calculated with the help of Equations (23)–(25). A special class of filters is the family of filters associated with MRA filters. An MRA filter m 0 generates the scaling function, associated with the MRA by [33]:
φ ^ ( ξ ) = j = 1 m 0 ( 2 j ξ )
A sufficient condition for a filter m 0 satisfying Equation (24), and taking on the value 1 at 0 , to be an MRA filter is that m 0 does not vanish on [ π / 2 ,   π / 2 ] .
| m 0 ( ξ ) | 2 + | m 0 ( ξ + π ) | 2 = 1
The kernel function { K ( ξ ) } is expressed by:
K ( ξ ) = 1 + π l = 0 N 1 ( 1 ) l ( 2 l + 1 ) a 1 c o s ( ( 2 l + 1 ) ξ )
where a 1 is sequence of coefficients and ξ is index of vector. The Fejer–Korovkin filters are defined by:
| m 0 n ( ξ ) | 2 = 1 2 π π / 2 π / 2 K ( ξ u ) d u
where m 0 n has length n if n is even and length n + 1 if n is odd.

2.3.5. Symlets Wavelet

Symlets wavelets, which are similar to the structure of Daubechies wavelets, are called orthogonal, biorthogonal, and least asymmetrical wavelets, and were introduced to the literature by Daubechies [29]. Daubechies wavelets have a maximum phase, while Symlets wavelets have a minimum phase. Unlike Daubechies wavelets, Symlets wavelets have a smoothed wavelet function with near-zero moments [34]. Symlets are more symmetrical than Daubechies, even though they have a 2 N 1 support size with N vanishing moments. Symlet wavelet coefficients for different filter lengths are calculated in [29].

2.4. Two-Dimensional Discrete Wavelet Transform (2D-DWT)

To simplify the vortex dynamics within a given implementation of the vorticity field, 2D-DWT is used, assuming that the vorticity field has finite energy; i.e., ζ 2 ( x ) d x < . Let φ be a scaling function and ψ the corresponding wavelet forming an orthonormal basis on L 2 ( R ) .
Three separable, “directionally sensitive” 2D-wavelet functions ψ H ( x , y ) , ψ V ( x , y ) , and ψ D ( x , y ) can be defined as
ψ H ( x , y ) = ψ ( x ) φ ( y )
ψ V ( x , y ) = φ ( x ) ψ ( y )
ψ D ( x , y ) = ψ ( x ) ψ ( y )
corresponding to horizontal, vertical, and diagonal directions, respectively. The directional sensitivity is a natural consequence of the separability in Equations (27)–(29); it does not increase the computational complexity of the 2D-DWT with the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families used in this study.
The separable scaling function φ ( x , y ) = φ ( x ) φ ( y ) is associated with the approximation space. The separable 2D-scaling and frequency dependent wavelet functions are defined as [35]:
φ j , m , n ( x , y ) = 2 j / 2 φ ( 2 j x m ,   2 j y n )
ψ j , m , n i ( x , y ) = 2 j / 2 ψ ( 2 j x m ,   2 j y n )
where i term denotes one of H , V , or D . The 2D-DWT calculation for the function f ( x , y ) in M and N dimensions is expressed by:
W φ ( j 0 , m , n ) = 1 M N x = 0 M 1 y = 0 N 1 f ( x , y ) φ j 0 , m , n ( x , y )
W ψ i ( j , m , n ) = 1 M N x = 0 M 1 y = 0 N 1 f ( x , y ) ψ j 0 , m , n i ( x , y )
for j = 0 , 1 ,   , J 1 and m , n = 0 , 1 , 2 ,   ,   2 j 1 . Note that j 0 is an arbitrary starting scale. The wavelet transform coefficients defined by Equations (32) and (33) are called approximation and detail coefficients, respectively. W φ ( j 0 , m , n ) coefficients describe an approximation of f ( x , y ) at this scale, and W ψ i ( j , m , n ) coefficients compose diagonal, vertical, and horizontal details for scales j j 0 .
In this study, 2D-DWT with the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families is applied to the image dataset consisting of apple, grape, and tomato plants. The feature groups extracted at various filter lengths for each wavelet family are used in the proposed plant leaf disease classification model structure. It should be noted that all of the wavelets used in the work avoid the computational complexity of complex numbers.
The scaling function, horizontal, vertical, and diagonal wavelet functions of wavelets that provide the best classification performance for each wavelet family are presented in Table 1. In the study, the wavelets that provide the best classification performance are obtained as “bior2.4” for the biorthogonal spline wavelet family, “coif1” for the Coiflet wavelet family, “db5” for the Daubechies wavelet family, “fk18” for the Fejer–Korovkin wavelet family, and “sym7” for the symlet wavelet family. Moreover, the best classification performance is obtained by the features extracted with “sym7” in the study, where there are twelve classes consisting of one healthy and three disease classes for each of the three plant species consisting of apple, grape, and tomato.

2.5. Flower Pollination Algorithm

The FPA introduced by Xin-She Yang in 2012 is a new metaheuristic inspired by the reproductive process of flowering plants to achieve the best result in the shortest time in solving global optimization problems [36]. Briefly, the algorithm mimics pollination in plants, including all the mechanisms for pollinators such as pollen transfer insects, bats, birds, wind, and water. The main purpose of flower pollination is to provide optimal vitality and the optimum biological reproduction phase. Pollination and other factors interact best to reproduce plants. Flowers need pollinators to carry out their pollination guidelines. There are two types of pollinators. The first of these are biotic pollinators, such as bats, flies, and bees. These pollinators can carry the pollen of flowers over great distances. In addition, the flight paths of these pollinators can be modeled with the Lévy distribution. The second type of pollinator is abiotic pollinators, such as water and wind. This form of pollination is considered to occur over short distances [37]. The pollination process is modeled by dividing it into two groups as cross-pollination and self-pollination. Cross-pollination is considered as the flowers of different plants pollinate each other. Self-pollination is defined as the pollination of the same flowers or different flowers of a plant. Self-pollination is used as the local search directive, while biotic cross pollination is defined as the global search directive. In the self-pollination process, there are some factors, such as the proximity of flowers to each other and the wind. Because of these factors, the local search process has more weight in the general search process. The pollinators and pollination types used in FPA are shown in Figure 1. The pollination process is mimicked using the following four basic rules:
(i)
Global pollination processes are carried out biotically, and pollinators carry pollen in the form of cross-pollination according to their Lévy flight.
(ii)
Abiotic pollination can occur in abiotic conditions, such as self-pollination and wind diffusion as local pollination.
(iii)
The coefficient called flower constancy is expressed as the probability of reproduction and varies in proportion to the similarity of flower species.
(iv)
Global pollination and local pollination are controlled by a switch probability p     [ 0 ,   1 ] . It is noted that p indicates the percentage balance between local and global search in the optimization search field.
In global pollination, pollen can be carried long distances because insects can fly for a long time. This guarantees the best reproduction. The mathematical expression of global pollination and flower constancy is defined by:
X i t + 1 = X i t + L ( X i t g * )
where X i t is the solution vector at iteration t , g * is the current best solution among all solutions in the current iteration, i is the pollen bundle or solution vector index, L is the step size representing the strength of pollination.
As insects and birds fly for a long distance, their motion can be shown according to the Lévy distribution. The long distances taken by pollinators are imitated by Lévy flight. Lévy flight distribution is defined as:
L ~ λ Γ ( λ ) sin ( π λ 2 ) π 1 s 1 + λ   ( s s 0 > 0 )
where Γ ( λ ) is the standard gamma function with an index λ and Lévy distribution is valid for large steps s > 0 . The abiotic pollination and flower constancy in the algorithm can be represented as:
X i t + 1 = X i t + ε ( X j t X k t )
where X j t and X k t represent pollen bundles of different flowers of the same plant; that is, different solutions of the solution set, ε denotes a random local pollination distance and has a normal distribution between [0, 1].
The remarkable feature of this algorithm is that many solution points are searched in the search space using the Lévy distribution. Determining the solution points in the solution space with the global search and searching the neighborhood of the solution points through the local search constitute the optimization logic of the algorithm.
The computational complexity of the wrapper approach, which consists of the FPA and SVM algorithms proposed from the study, is O ( N ( n × p + f n × p ) ) . Here, N denotes the number of iterations, n indicates the number of features, and f n defines fitness function value, p represents the number of selected features. Since FPA has fast convergence, low complexity, and good optimization performance, SVM has low computational complexity and is easily repeatable, the wrapper approach consisting of the FPA and SVM algorithms is used for feature selection in this study.

2.6. Convolutional Neural Network Classifier

Traditional neural networks are heavily inspired by the way biological neural systems work and consist of a high number of interconnected computational nodes working in a distributed manner to collectively learn from the input to optimize the final output. CNNs are similar to traditional neural networks in that they consist of neurons that self-optimize through learning. CNNs, which are used in many computer-vision fields, especially classification, include convolutional feature extraction and classification layers, unlike traditional neural networks [39]. In this study, the CNN classifier, which has only one classification layer without a feature extraction layer, is used to overcome the model hyperparameter problem and minimize the model complexity.
CNNs consist of three types of layers, the convolution layer, the pool layer, and the fully connected layer. The feature maps created in the final convolution and pooling layers are converted to a one-dimensional array of numbers and mapped to the final output of the network by connecting to one or more fully connected layers. In this study, the CNN classifier with two hidden layers consisting of flattened and fully connected layers is used. The features selected by the wrapper approach consisting of the FPA and SVM algorithms are given as input to the flatten layer of the CNN. Thus, the computational load of the model is prevented from increasing.

2.7. Performance Metrics for Classification

The whole dataset used in the study is randomly separated into two independent datasets, 80% and 20%, for the training and validation phases, respectively. The training performance of the models created in the study for plant leaf disease classification is measured by the 10-fold cross-validation method. The training performance of each classification model is calculated by averaging the accuracy values for each fold. The performances of all classification models used in the study are measured on the test data with the help of the metrics in Equations (37)–(40). The metrics given in Equations (37)–(40) are calculated based on the confusion matrix. Here, T P , F P , T N , and F N denote true positives, false positives, true negatives, and false negatives, respectively.
A c c u r a c y = T P + T N T P + T N + F P + F N   ,  
P r e c i s i o n = T P T P + F P   ,
R e c a l l = T P T P + F N   ,
F 1   s c o r e = 2 × R e c a l l × P r e c i s i o n R e c a l l + P r e c i s i o n   .
A c c u r a c y is generally a measure of how often the classifier estimates correctly. P r e c i s i o n is the ratio of how many positively predicted samples are predicted correctly. R e c a l l is an indication of how many of the samples that should be predicted positively were predicted correctly. F 1   s c o r e is the harmonic mean of recall and precision. It is a measure of how well the classifier is performing and is often used to compare classifiers [40].

2.8. Framework of the Proposed Methodology

In this study, a new robust hybrid model whose features are extracted with 2D-DWT and selected with the FPA-SVM wrapper approach and include a CNN classifier is proposed to classify plant leaf diseases in real time with high accuracy, low computational cost, and low parameter complexity. The framework of the proposed approach for real-time classification of plant leaf diseases is presented in Figure 2. The phases of the proposed approach are briefly summarized below:
In the data preparation phase, the image dataset consisting of apple, grape, and tomato plant diseases is randomly divided into two independent datasets, 80% and 20%, respectively, for the training and validation phases.
At the phase of applying 2D-DWT with wavelets, the distinguishing features defining the characteristic features of plant leaf diseases are extracted with the 2D-DWT using wavelet families such as biorthogonal, Coiflets, Daubechies, Fejer–Korovkin, and symlets.
In the feature extraction phase, energy and statistical-based features are extracted from the vertical, horizontal, diagonal, and approximate matrices of 2D-DWT. Six features, namely the arithmetic mean, entropy, standard deviation, skewness, kurtosis, and energy, are applied to these four matrices. Moreover, the column vector, which is the maximum value of the columns of any matrices, is expressed with MCV, and the row vector, which is the maximum value of the rows of them, is expressed with MRV. Additionally, six properties, namely arithmetic mean, standard deviation, and entropy of both MCV and MRV, were applied to the four matrices. The same process is also repeated for the second level of decomposition in 2D-DWT. At the end of this phase, a total of 96 features are extracted.
At the feature selection with the FPA-SVM phase, the most suitable ones among the normalized features for each wavelet family are selected with the help of the wrapper approach consisting of the FPA and SVM algorithms. In addition, the fitness function, which takes into account both the number of features used in the model and the model performance, is determined in order to keep the model complexity and computation cost to a minimum level.
In the evaluation of model performance phase, CNN, SVM, and KNN classifier performances are measured with the help of performance metrics and the model with the highest performance is determined.

3. Experimental Results and Discussion

In this section, studies showing the classification efficiency of the proposed hybrid model for real-time classification of plant leaf diseases on a dataset consisting of leaf images of healthy and diseased apple, grape, and tomato plants are presented. The proposed model includes the 2D-DWT signal processing method based on the “sym7” wavelet family, the wrapper approach consisting of FPA and SVM, and the CNN classifier. The efficiency of the proposed optimization algorithm is also compared with the PSO algorithm. The performance of the proposed CNN classifier for the hybrid model is compared with the performances of the SVM and KNN classification algorithms, and its effect on the performance of the plant leaf disease classification model is examined.

3.1. Dataset

In this study, a data set consisting of leaf images of apple, grape, and tomato plants with a size of 256 × 256 pixels was used. The image dataset consisting of apple, grape, and tomato plant diseases used in the study was randomly divided into two independent datasets, 80% and 20%, respectively, for the training and validation phases.
The dataset consisting of apple plant leaf images consists of healthy, black rot, cedar rust, and scab disease classes seen in Figure 3. In the training phase of the models in the study, a total of 1100 apple plant leaf images, 275 images for each class, were used, and in the validation phase, a total of 220 apple plant leaf images, 55 images for each class, were used. The dataset consisting of grape plant leaf images consists of healthy, black rot, black measles, and leaf blight disease classes shown in Figure 4. In the training phase of the models in the study, a total of 1680 grape plant leaf images, 420 images for each class, were used, and in the validation phase, a total of 336 grape plant leaf images, 84 images for each class, were used. The dataset consisting of tomato plant leaf images consists of healthy, bacterial spot, late blight, and yellow leaf curl disease classes presented in Figure 5. In the training phase of the models in the study, a total of 1700 tomato plant leaf images, 425 images for each class, were used, and in the validation phase, a total of 340 tomato plant leaf images, 85 images from each class, were used.
There are 12 classes in total in the study, and the number of plant leaf images belonging to these 12 classes is 5376, of which 4480 are training data and 896 are validation data. Apple, grape, and tomato plant leaf disease images used in the training and validation phases of the model have been taken from the open-source Plant Village dataset [5].
The testing phase of the models used in the study was carried out on the images obtained with the help of the camera on the UAV. In the study, the GoPro Hero4 camera, which can shoot 30 frames per second, was used. The images obtained with the UAV have a resolution of 1920 × 1080 . These images were resized on the embedded system and presented to the model with 256 × 256 resolution. The stabilization of the camera was provided with the help of the gimbal system. In order to cope with environmental conditions, such as blurring, threshold values were used for images taken with the UAV. These threshold values decide whether the image received by the UAV will be processed or not.
The testing phase of the models was done with data as much as with the number of images used in the validation phase of the models. Real-time tests of the aforementioned models were carried out in various provinces of Turkey in August 2022. Apple leaf disease tests were carried out in apple orchards in Eskişehir, grape leaf disease tests were carried out in vineyards in Manisa, and tomato leaf disease tests were carried out in tomato fields in Antalya province. During the time intervals when the tests were carried out, the average temperature in Eskişehir was 27 , the average temperature in Manisa was 30 , and the average temperature in Antalya was 32 . Average humidity values were measured at 62%, 73%, and 78% for Eskişehir, Manisa, and Antalya, respectively. Considering the seasonality of some diseases, the tests in the study were carried out in natural agricultural areas instead of greenhouses in various provinces of Turkey.

3.2. Applying 2D-DWT with Wavelet Families

The features of plant leaf diseases were extracted by 2D-DWT using the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families for filter lengths shown in Table 2. A two-level decomposition was applied to the original image matrix of 256 × 256 pixels. As a result of this decomposition, vertical, horizontal, diagonal, and approximation image matrices were obtained for each wavelet family. In the first level of decomposition, four image matrices with 130 × 130 pixel size and four image matrices with 67 × 67 pixel size in the second level of decomposition, a total of eight image matrices were obtained.

3.3. Extraction of Statistical and Entropy-Based Features

Plant leaf disease is diagnosed by tests such as enzyme-linked immunosorbent assay (ELISA) and polymerase chain reaction (PCR) performed in the laboratory environment, as well as visual inspections made by experienced individuals, be they a botanist or farmer. While this is the right approach, it is a costly and highly labor-intensive process as it requires the installation of laboratory equipment. However, these traditional methods based on experience and laboratory testing are not suitable for real-time detection of plant leaf disease as they are time-consuming and allow expert error under heavy workload. Since microscopic evaluation and diagnostic experiments such as ELISA and PCR do not allow real-time detection of plant leaf disease, the distinguishing features of diseases have been extracted. Extracted features define the characteristic structure of plant leaf diseases, unlike the features obtained by visual image analysis techniques.
Statistical and entropy-based features were used in the study to classify plant leaf diseases with high accuracy. A total of 96 features were extracted, and each feature was labeled as seen in Table 3.
Statistical and entropy-based features given in Table 4 were applied to the image matrices, I = { I V 1 , I H 1 ,   I D 1 ,   I A 1 ,   I V 2 ,   I H 2 ,   I D 2 ,   I A 2 } , obtained from first and second-level decomposition for each wavelet family. Half of the extracted features consist of the mean, standard deviation, and entropy values of MRV, which is the maximum row vector of the I j i ( x , y ) image, and MCV, which is the maximum column vector of the I j i ( x , y ) image. Thus, the sensitivity of the proposed plant leaf disease classification model to the changes in the rows and columns of the I j i ( x , y ) image matrix was increased [41].

3.4. Feature Selection with FPA-SVM Method

In order to keep both the classification performance high and the computational complexity to a minimum level, it is necessary to determine the subset consisting of the most suitable features out of the 96 features extracted from eight image matrices obtained by applying the 2D-DWT technique. FPA, which mimics the reproduction process of flowering plants, selects the most suitable features that keep the SVM classifier performance high among the normalized features for each wavelet family with the help of the fitness function, which is defined as:
f i t n e s s = μ × γ R ( ε C ) + σ × ( c a r d i n a l i t y   o f   t h e   s e l e c t e d   s u b s e t t o t a l   n u m b e r   o f   f e a t u r e s )
where γ R ( ε C ) indicates the classification error rate of classifier. μ and σ represent the significance of classification quality and subset length, respectively. Here, μ [ 0 , 1 ] and σ = ( 1 μ ) . The fitness function considers both the number of features used in the model and model performance to keep model complexity and computational cost to minimum levels.
Due to its high performance for the selected feature subset and fast response time, the wrapper approach in the study is composed of the SVM algorithm based on statistical learning theory. The effectiveness of the proposed optimization algorithm for the wrapper approach to feature selection is emphasized by comparing it with the PSO algorithm. The parameter values of the FPA-SVM and PSO-SVM wrapper approaches used in the study are presented in Table 5. It should be noted that in the models created for each wavelet family, feature selection is performed according to the parameter values in Table 5.
The total selection rates for the best feature groups belonging to the biorthogonal, Coiflet, Daubechies, Fejer–Korovkin, and symlet wavelet families after 50 iterations in FPA were given as percentages in the heat map in Figure 6. It is seen that the highest number of features were extracted from I H 1 , I D 1 , and I A 1 for biorthogonal; I V 1 , I V 2 , and I H 2 for Coiflets; I V 1 , I D 1 , and I A 1 for Daubechies; I V 1 , I D 1 , and I H 2 for Fejer–Korovkin; and I V 1 , I D 1 , and I H 2 for symlets in eight image matrices obtained from 2D-DWT analysis.
The F 4 feature, one of the statistical and entropy-based features, was used at least 10% in all wavelet families. In addition to this feature, F 5 for biorthogonal, F 2 for Coiflets, F 8 for Daubechies, F 5 for Fejer–Korovkin, and F 11 for symlets were selected at least 10%. It can be seen from Figure 6 that these two features specified for each wavelet family constitute at least 25% of the features selected for the plant leaf disease classification model. On the other hand, in the proposed plant leaf disease classification model, it is noted that the features F 6 and F 7 for biorthogonal; F 1 , F 3 , and F 12 for Coiflets; F 6 and F 9 for Daubechies; F 6 and F 12 for Fejer–Korovkin; and F 1 and F 6 for symlets were used the least.

3.5. Evaluation of Plant Leaf Disease Classification Models and Discussion

All models mentioned in the study of plant leaf disease images were performed on a personal computer with an Intel Core i7–10875H processor, an 8 GB NVIDIA RTX 3070 graphics card, and 16 GB of RAM. All codes for the models were compiled by MATLAB 2021b. The models created in the study were tested on 896 pieces of data, including 12 plant leaf disease classes, taken with the camera on the UAV. All models created were run 50 times, and the performance of the models was computed as the mean and standard deviation. In the study, the effects of both optimization and classifier algorithms on model performance were examined. The performances of the CNN, KNN, and SVM classification models created with the features selected by both the FPA-SVM and PSO-SVM wrapper approaches were measured in terms of accuracy metrics, and the results are presented in Table 6 and Table 7, respectively.
Twenty-three feature subsets were selected from the features generated by the 2D-DWT method based on the “sym7” wavelet family by the FPA-SVM wrapper approach. The performances of CNN, SVM, and KNN classifiers created with 23 selected features have been measured at 99.55%, 97.54%, and 93.97%, respectively. The plant leaf disease classification model, which includes the 2D-DWT method based on the “sym7” wavelet family, the wrapper approach consisting of FPA and SVM, and the CNN classifier, has been proposed as the model with the best classification performance in the study. For the wrapper approach consisting of FPA and SVM, the best performance among the models created with the SVM and KNN classifier algorithms has been obtained with the features generated by the “sym7” wavelet family. The performance of the CNN classifier has been measured over 95% for wavelet families excluding “coif2” and “db2”. The performance of the SVM classifier has been measured at over 90% for wavelet families excluding “bior1.1” and “db1”. The performance of the KNN classifier has been measured at over 85% for wavelet families excluding “bior1.1”, “bior1.3”, “bior1.5”, “bior2.8”, “db1”, “db2”, “fk4”, and “sym4”. It is seen that the performance of the CNN classifier is superior when compared to the SVM and KNN classifiers.
Twenty-two feature subsets were selected from the features generated by the 2D-DWT method based on the “db6” wavelet family by the PSO-SVM wrapper approach. The performances of CNN, SVM, and KNN classifiers created with 22 selected features have been measured at 94.87%, 90.15%, and 87.00%, respectively. For the wrapper approach consisting of PSO and SVM, the best performance among the models created with the SVM and KNN classifier algorithms has been obtained with the features generated by the “db6” wavelet family. The performance of the CNN classifier has been measured over 90% for wavelet families excluding “bior2.6”. The performance of the SVM classifier has been measured over 85% for wavelet families, excluding “bior1.1”, “bior1.3”, “bior1.5”, “bior2.6”, “coif5”, “db1”, “db8”, “db9”, “fk4”, “fk6”, “fk14”, “sym2”, and “sym7”. The performance of the KNN classifier has been measured over 80% for wavelet families, excluding “bior1.3”, “bior1.5”, “bior2.6”, “bior3.3”, “coif1”, “coif5”, “db3”, “db8”, “db9”, “db10”, “fk4”, “fk6”, “fk8”, “fk14”, “sym2”, and “sym7”. It is seen that the performance of the CNN classifier is superior when compared to the SVM and KNN classifiers.
Considering the model structure and input size of CNN, it can be said that the model complexity of the proposed classification model is low since it has high performance with fewer features selected by the FPA-SVM wrapper approach. When Table 6 and Table 7 are evaluated together, considering both the classification model performance and the number of selected features, it is seen that the proposed FPA technique is successful compared to the PSO algorithm for CNN, SVM, and KNN classifiers.
The CNN, SVM, and KNN classifier model performances, which were obtained as a result of the 50 times training-validation process for the selected features with the FPA-SVM wrapper approach proposed in the study, are compared with each other for the “bior”, “coif”, “db”, “fk”, and “sym” wavelet families in Figure 7.
It is seen that the performances of the models obtained with the CNN classifier are obviously higher than the performances of the models obtained with both SVM and KNN classifiers. It should also be noted that the SVM classifier performance is higher than the KNN classifier performance for the specified wavelet families. The performance of all classifiers created with the selected features with the FPA-SVM wrapper approach proposed for the “sym7” wavelet family has the best performance. It is noteworthy that the proposed FPA optimization technique improves CNN classifier performance by approximately 5% and SVM and KNN classifier performances by approximately 8% compared to the PSO. It is also seen that the performance of all classifiers created with the features selected with the PSO-SVM wrapper approach for the “db6” wavelet family has the best performance. The performances of various classifiers created with the features selected by the FPA-SVM and PSO-SVM wrapper approaches for wavelet families that provide the best classification performance are measured with the precision, recall, and F1 score performance metrics as well as the accuracy metric and are presented in Table 8.
In the study, the best-performing model, which includes the 2D-DWT signal processing method based on the “sym7” wavelet family and the wrapper approach consisting of FPA-SVM and the CNN classifier, was embedded in the NVIDIA Jetson Nano artificial intelligence (AI) application development kit on the UAV, and real-time classification tests have been performed on apple, grape, and tomato plants. The results obtained as a result of the experiments performed on the test images are presented in Figure 8 in the form of a confusion matrix. It is seen that the model run in real time has a very low-level classification error for only two disease classes out of 12. The model predicts 2 out of 84 grape (black measles) test data as grape (black rot) disease classes. Similarly, the model classifies 2 out of 85 tomato (bacterial spot) test data as apple (healthy). The proposed model can classify both plant species and plant diseases in real time with high accuracy.
Moreover, a total of 227 models created, except for the model proposed in the study, were embedded on the NVIDIA Jetson Nano AI application development kit on the UAV, and the performances of all models were measured. When the performances of the models were evaluated in terms of computational time, it was seen that the computational time of the proposed model was 8.1% better than the second-best model in the study. As a result of the experimental studies, it is concluded that our robust hybrid classification model, which is created by selecting the lowest number of features with an optimization algorithm with low computational complexity, can classify plant leaf diseases with precision in real time.

4. Conclusions and Future Work

Detection of plant diseases at an early phase is of great importance for appropriate treatment. Automatic detection and classification of plant leaf diseases instead of visual inspection contributes to increasing agricultural product productivity. In this study, a new robust hybrid classification model based on swarm optimization-supported feature selection, including machine learning and deep learning algorithms, is proposed to classify diseases in apple, grape, and tomato plants with high accuracy in real time. The proposed model includes the 2D-DWT signal processing method based on the “sym7” wavelet family, the wrapper approach consisting of FPA-SVM, and the CNN classifier. In order for the model to be robust, features extracted from various wavelet families, such as biorthogonal, Coiflets, Daubechies, Fejer–Korovkin, and symlets, are used with the 2D-DWT method. The features that keep the classifier performance high for each wavelet family are selected with a wrapper approach consisting of FPA and SVM, so that the computational complexity of the model has been kept to a minimum level. The CNN classification model is created using the lowest number of features that can keep the classification performance high. The CNN classifier with a single layer of classification without a feature extraction layer is used to minimize the complexity of the model and to overcome the model hyperparameter problem. Our proposed model for plant leaf disease classification is embedded in the NVIDIA Jetson Nano AI development kit on the UAV, and its real-time performance is tested. The experimental results obtained prove that the proposed model can detect and classify the specified plant leaf diseases in real time with high accuracy. The proposed classification model with fast and low computational cost will contribute to an increase in agricultural efficiency by classifying plant leaves with high accuracy at an early phase. In future studies, algorithms that can cope with the chaoticity on 2D images will be tested on classifier model structures in order to detect and classify non-specific diseases in similar leaf images.

Author Contributions

Conceptualization, A.A.; methodology, İ.Y. and A.A.; software, İ.Y. and A.A.; writing—original draft preparation, İ.Y. and A.A.; writing—review and editing, A.A.; supervision, A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Savary, S.; Ficke, A.; Aubertot, J.; Hollier, C. Crop losses due to diseases and their implications for global food production losses and food security. Food Secur. 2012, 4, 519–537. [Google Scholar] [CrossRef]
  2. Ons, L.; Bylemans, D.; Thevissen, K.; Cammue, B.P. Combining biocontrol agents with chemical fungicides for integrated plant fungal disease control. Microorganisms 2020, 8, 1930. [Google Scholar] [CrossRef] [PubMed]
  3. Tewari, V.K.; Pareek, C.M.; Lal, G.; Dhruw, L.K.; Singh, N. Image processing based real-time variable-rate chemical spraying system for disease control in paddy crop. Artif. Intell. Agric. 2020, 4, 21–30. [Google Scholar] [CrossRef]
  4. Chadha, S.; Sharma, M.; Sayyed, A. Advances in sensing plant diseases by imaging and machine learning methods for precision crop protection. In Microbial Management of Plant Stresses; Woodhead Publishing: Cambridge, UK, 2021; pp. 157–183. [Google Scholar]
  5. Mohanty, S.P.; Hughes, D.P.; Salathé, M. Using deep learning for image-based plant disease detection. Front. Plant Sci. 2016, 7, 1419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhang, N.; Yang, G.; Pan, Y.; Yang, X.; Chen, L.; Zhao, C. A review of advanced technologies and development for hyperspectral-based plant disease detection in the past three decades. Remote Sens. 2020, 12, 3188. [Google Scholar] [CrossRef]
  7. Barbedo, J.G.A.; Koenigkan, L.V.; Santos, T.T. Identifying multiple plant diseases using digital image processing. Biosyst. Eng. 2016, 147, 104–116. [Google Scholar] [CrossRef]
  8. Barbedo, J.G.A. Impact of dataset size and variety on the effectiveness of deep learning and transfer learning for plant disease classification. Comput. Electron. Agric. 2018, 153, 46–53. [Google Scholar] [CrossRef]
  9. Almadhor, A.; Rauf, H.T.; Lali, M.I.U.; Damaševičius, R.; Alouffi, B.; Alharbi, A. AI-driven framework for recognition of guava plant diseases through machine learning from DSLR camera sensor based high resolution imagery. Sensors 2021, 21, 3830. [Google Scholar] [CrossRef]
  10. Ali, H.; Lali, M.I.; Nawaz, M.Z.; Sharif, M.; Saleem, B.A. Symptom based automated detection of citrus diseases using color histogram and textural descriptors. Comput. Electron. Agric. 2017, 138, 92–104. [Google Scholar] [CrossRef]
  11. Wen, D.M.; Chen, M.X.; Zhao, L.; Ji, T.; Li, M.; Yang, X.T. Use of thermal imaging and Fourier transform infrared spectroscopy for the pre-symptomatic detection of cucumber downy mildew. Eur. J. Plant Pathol. 2019, 155, 405–416. [Google Scholar] [CrossRef]
  12. Singh, V.; Misra, A.K. Detection of plant leaf diseases using image segmentation and soft computing techniques. Inf. Process. Agric. 2017, 4, 41–49. [Google Scholar] [CrossRef] [Green Version]
  13. Saeed, F.; Khan, M.A.; Sharif, M.; Mittal, M.; Goyal, L.M.; Roy, S. Deep neural network features fusion and selection based on PLS regression with an application for crops diseases classification. Appl. Soft Comput. 2021, 103, 107164. [Google Scholar] [CrossRef]
  14. Ozguven, M.M.; Altas, Z. A new approach to detect mildew disease on cucumber (Pseudoperonospora cubensis) leaves with image processing. J. Plant Pathol. 2022, 104, 1–10. [Google Scholar] [CrossRef]
  15. Tan, L.; Lu, J.; Jiang, H. Tomato leaf diseases classification based on leaf images: A comparison between classical machine learning and deep learning methods. AgriEngineering 2021, 3, 542–558. [Google Scholar] [CrossRef]
  16. Kamal, K.C.; Yin, Z.; Wu, M.; Wu, Z. Depthwise separable convolution architectures for plant disease classification. Comput. Electron. Agric. 2019, 165, 104948. [Google Scholar]
  17. Atila, Ü.; Uçar, M.; Akyol, K.; Uçar, E. Plant leaf disease classification using EfficientNet deep learning model. Ecol. Inform. 2021, 61, 101182. [Google Scholar] [CrossRef]
  18. Li, Y.; Nie, J.; Chao, X. Do we really need deep CNN for plant diseases identification. Comput. Electron. Agric. 2020, 178, 105803. [Google Scholar] [CrossRef]
  19. Thangaraj, R.; Anandamurugan, S.; Kaliappan, V.K. Automated tomato leaf disease classification using transfer learning-based deep convolution neural network. J. Plant Dis. Prot. 2021, 128, 73–86. [Google Scholar] [CrossRef]
  20. Too, E.C.; Yujian, L.; Njuki, S.; Yingchun, L. A comparative study of fine-tuning deep learning models for plant disease identification. Comput. Electron. Agric. 2019, 161, 272–279. [Google Scholar] [CrossRef]
  21. Darwish, A.; Ezzat, D.; Hassanien, A.E. An optimized model based on convolutional neural networks and orthogonal learning particle swarm optimization algorithm for plant diseases diagnosis. Swarm Evol. Comput. 2020, 52, 100616. [Google Scholar] [CrossRef]
  22. Zhong, J.; Huang, Y. Time-frequency representation based on an adaptive short-time Fourier transform. IEEE Trans. Signal Process. 2010, 58, 5118–5128. [Google Scholar] [CrossRef]
  23. Li, L.; Cai, H.; Han, H.; Jiang, Q.; Ji, H. Adaptive short-time Fourier transform and synchrosqueezing transform for non-stationary signal separation. Signal Process. 2020, 166, 107231. [Google Scholar] [CrossRef]
  24. Kim, C.H.; Aggarwal, R. Wavelet transforms in power systems. Part 1: General introduction to the wavelet transforms. Power Eng. J. 2000, 14, 81–87. [Google Scholar]
  25. Grossmann, A.; Morlet, J. Decomposition of Hardy functions into square integrable wavelets of constant shape. SIAM J. Math. Anal. 1984, 15, 723–736. [Google Scholar] [CrossRef] [Green Version]
  26. Mallat, S.G. A theory for multiresolution signal decomposition: The wavelet representation. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 11, 674–693. [Google Scholar] [CrossRef] [Green Version]
  27. Keinert, F. Biorthogonal wavelets for fast matrix computations. Appl. Comput. Harmon. Anal. 1994, 1, 147–156. [Google Scholar] [CrossRef] [Green Version]
  28. Monzón, L.; Beylkin, G.; Hereman, W. Compactly supported wavelets based on almost interpolating and nearly linear phase filters (coiflets). Appl. Comput. Harmon. Anal. 1999, 7, 184–210. [Google Scholar] [CrossRef] [Green Version]
  29. Daubechies, I. Ten Lectures on Wavelets; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1992; ISBN 0-89871-274-2. [Google Scholar]
  30. Cohen, A.; Sun, Q. An arithmetic characterization of the conjugate quadrature filters associated to orthonormal wavelet bases. SIAM J. Math. Anal. 1993, 24, 1355–1360. [Google Scholar] [CrossRef]
  31. Daubechies, I. Orthonormal bases of compactly supported wavelets. Commun. Pure Appl. Math. 1988, 41, 909–996. [Google Scholar] [CrossRef] [Green Version]
  32. Daubechies, I. The wavelet transform, time-frequency localization and signal analysis. IEEE Trans. Inf. Theory 1990, 36, 961–1005. [Google Scholar] [CrossRef] [Green Version]
  33. Nielsen, M. On the construction and frequency localization of finite orthogonal quadrature filters. J. Approx. Theory 2001, 108, 36–52. [Google Scholar] [CrossRef] [Green Version]
  34. Mallat, S. A Wavelet Tour of Signal Processing, 2nd ed.; Elsevier: Berkeley, CA, USA, 1999; ISBN 978-0-12-466606-1. [Google Scholar]
  35. Gonzalez, R.C.; Woods, R.E. Digital Image Processing, 4th ed.; Pearson Education: Rotherham, UK, 2002; ISBN 978-0-13-335672-4. [Google Scholar]
  36. Yang, X.S. Flower Pollination Algorithm for Global Optimization. In International Conference on Unconventional Computing and Natural Computation; Springer: Berlin, Germany, 2012; pp. 240–249. [Google Scholar]
  37. Abdel-Basset, M.; El-Shahat, D.; El-Henawy, I.; Sangaiah, A.K. A modified flower pollination algorithm for the multidimensional knapsack problem: Human-centric decision making. Soft Comput. 2018, 22, 4221–4239. [Google Scholar] [CrossRef]
  38. Abdel-Basset, M.; Shawky, L.A. Flower pollination algorithm: A comprehensive review. Artif. Intell. Rev. 2019, 52, 2533–2557. [Google Scholar] [CrossRef]
  39. Gu, J.; Wang, Z.; Kuen, J.; Ma, L.; Shahroudy, A.; Shuai, B.; Liu, T.; Wang, X.; Wang, G.; Cai, J.; et al. Recent advances in convolutional neural networks. Pattern Recognit. 2018, 77, 354–377. [Google Scholar] [CrossRef] [Green Version]
  40. Ma, C.; Zhang, H.H.; Wang, X. Machine learning for big data analytics in plants. Trends Plant Sci. 2014, 19, 798–808. [Google Scholar] [CrossRef]
  41. Karasu, S.; Saraç, Z. The effects on classifier performance of 2D discrete wavelet transform analysis and whale optimization algorithm for recognition of power quality disturbances. Cogn. Syst. Res. 2022, 75, 1–15. [Google Scholar] [CrossRef]
Figure 1. Pollination types and pollinators in the FPA [38].
Figure 1. Pollination types and pollinators in the FPA [38].
Biology 11 01732 g001
Figure 2. Framework of the proposed robust hybrid model to classify plant leaf diseases in real time.
Figure 2. Framework of the proposed robust hybrid model to classify plant leaf diseases in real time.
Biology 11 01732 g002
Figure 3. Leaf image examples belonging to healthy and diseased classes of apple plants.
Figure 3. Leaf image examples belonging to healthy and diseased classes of apple plants.
Biology 11 01732 g003
Figure 4. Leaf image examples belonging to healthy and diseased classes of grape plants.
Figure 4. Leaf image examples belonging to healthy and diseased classes of grape plants.
Biology 11 01732 g004
Figure 5. Leaf image examples belonging to healthy and diseased classes of tomato plants.
Figure 5. Leaf image examples belonging to healthy and diseased classes of tomato plants.
Biology 11 01732 g005
Figure 6. Heat map of selected features of all iterations of FPA for the (a) biorthogonal, (b) Coiflet, (c) Daubechies, (d) Fejer–Korovkin, and (e) symlet wavelet families.
Figure 6. Heat map of selected features of all iterations of FPA for the (a) biorthogonal, (b) Coiflet, (c) Daubechies, (d) Fejer–Korovkin, and (e) symlet wavelet families.
Biology 11 01732 g006
Figure 7. Performances of classifiers used in the study for the (a) biorthogonal, (b) Coiflet, (c) Daubechies, (d) Fejer–Korovkin, and (e) symlet wavelet families.
Figure 7. Performances of classifiers used in the study for the (a) biorthogonal, (b) Coiflet, (c) Daubechies, (d) Fejer–Korovkin, and (e) symlet wavelet families.
Biology 11 01732 g007
Figure 8. Confusion matrix for the proposed apple, grape, and tomato plant leaf disease classification model.
Figure 8. Confusion matrix for the proposed apple, grape, and tomato plant leaf disease classification model.
Biology 11 01732 g008
Table 1. Scaling and mother wavelet functions of wavelets that provide the best classification performance for each wavelet type used in the study.
Table 1. Scaling and mother wavelet functions of wavelets that provide the best classification performance for each wavelet type used in the study.
Wavelet TypesScaling Function
φ ( x , y ) = φ ( x ) φ ( y )
Horizontal Wavelet
ψ H ( x , y ) = ψ ( x ) φ ( y )
Vertical Wavelet
ψ V ( x , y ) = φ ( x ) ψ ( y )
Diagonal Wavelet
ψ D ( x , y ) = ψ ( x ) ψ ( y )
Biorthogonal spline, bior2.4Biology 11 01732 i001Biology 11 01732 i002Biology 11 01732 i003Biology 11 01732 i004
Coiflets, coif1Biology 11 01732 i005Biology 11 01732 i006Biology 11 01732 i007Biology 11 01732 i008
Daubechies, db5Biology 11 01732 i009Biology 11 01732 i010Biology 11 01732 i011Biology 11 01732 i012
Fejer-Korovkin, fk18Biology 11 01732 i013Biology 11 01732 i014Biology 11 01732 i015Biology 11 01732 i016
Symlets, sym7Biology 11 01732 i017Biology 11 01732 i018Biology 11 01732 i019Biology 11 01732 i020
Table 2. Wavelet families and filter lengths used in the study.
Table 2. Wavelet families and filter lengths used in the study.
Wavelet FamilyFilter Length
Biorthogonal(1.) 1, 3, 5, (2.) 2, 4, 6, 8, (3.) 1, 3, 5
Coiflet1, 2, 3, 4, 5
Daubechies1, 2, 3, 4, 5, 6, 7, 8, 9, 10
Fejer–Korovkin4, 6, 8, 14, 18, 22
Symlet2, 3, 4, 5, 6, 7, 8
Table 3. Index labels of features extracted from image matrices.
Table 3. Index labels of features extracted from image matrices.
I V 1 I H 1 I D 1 I A 1 I V 2 I H 2 I D 2 I A 2
F 1 113253749617385
F 2 214263850627486
F 3 315273951637587
F 4 416284052647688
F 5 517294153657789
F 6 618304254667890
F 7 719314355677991
F 8 820324456688092
F 9 921334557698193
F 10 1022344658708294
F 11 1123354759718395
F 12 1224364860728496
Table 4. Statistical and entropy-based features and their equations.
Table 4. Statistical and entropy-based features and their equations.
LabelFeature NameFeature Expression
F 1 Arithmetic mean m e a n = 1 m × n x m y n | I j i ( x , y ) |
F 2 Entropy e n t r o p y = x m y n I j i ( x , y ) l o g | I j i ( x , y ) |
F 3 Standard deviation s t d = 1 m × n x m y n ( | I j i ( x , y ) | m e a n ) 2
F 4 Skewness s k w = 1 m × n x m y n ( | I j i ( x , y ) | m e a n s t d ) 3
F 5 Kurtosis k r t s = 1 m × n x m y n ( | I j i ( x , y ) | m e a n s t d ) 4
F 6 Energy e n e r g y = x m y n ( I j i ( x , y ) ) 2
F 7 MRV mean M R V m e a n = 1 m x m M R V ( x )
F 8 MCV mean M C V m e a n = 1 n y n M C V ( y )
F 9 Standard deviation of MRV M R V s t d = 1 m x m ( M R V ( x ) M R V m e a n ) 2
F 10 Standard deviation of MCV M C V s t d = 1 n y n ( M C V ( y ) M C V m e a n ) 2
F 11 MRV entropy M R V e n t r o p y = x m | M R V ( x ) | l o g | M R V ( x ) |
F 12 MCV entropy M C V e n t r o p y = y n | M C V ( y ) | l o g | M C V ( y ) |
Table 5. Parameter values of wrapper feature selection approaches used in the study.
Table 5. Parameter values of wrapper feature selection approaches used in the study.
ParametersFPA-SVMPSO-SVM
Number of solutions3030
Maximum number of iterations5050
Number of features9696
Threshold0.70.7
Other parametersswitch probability = 0.4
levy component = 1.5
cognitive factor = 2
social factor = 2
inertia weight = 1
Fitness functionmaximization of classifier performance & minimization of the number of selected features
Table 6. The performance of the classifiers for selected features by the FPA-SVM wrapper approach.
Table 6. The performance of the classifiers for selected features by the FPA-SVM wrapper approach.
WaveletsNum.Selected FeaturesAccuracy (%)
CNNSVMKNN
Biorthogonal (bior)1.1243, 4, 5, 14, 15, 16, 23, 24, 26, 33, 36, 45, 47, 48, 55, 56, 65, 75, 80, 85, 87, 88, 89, 9296.90 ± 0.0989.26 ± 0.3683.66 ± 0.04
1.3238, 17, 20, 22, 23, 25, 26, 28, 32, 37, 41, 44, 45, 52, 58, 60, 62, 63, 65, 74, 77, 84, 8596.76 ± 0.1191.05 ± 0.0283.88 ± 0.07
1.5263, 4, 6, 15, 16, 21, 22, 27, 28, 29, 34, 37, 40, 43, 45, 48, 63, 64, 65, 69, 71, 76, 83, 90, 91, 9497.68 ± 0.6991.56 ± 0.6384.22 ± 0.16
2.2174, 6, 28, 29, 34, 35, 41, 47, 50, 52, 55, 61, 62, 64, 65, 76, 8396.45 ± 0.2092.14 ± 0.2791.00 ± 0.96
2.4211, 4, 5, 8, 23, 27, 28, 30, 34, 36, 41, 42, 47, 56, 57, 58, 65, 66, 87, 92, 9398.08 ± 0.0293.53 ± 0.3388.50 ± 0.78
2.6194, 5, 9, 11, 16, 21, 23, 28, 31, 32, 33, 37, 40, 41, 52, 53, 57, 58, 7597.34 ± 0.4793.53 ± 0.1191.12 ± 0.07
2.82110, 17, 19, 25, 28, 29, 33, 35, 37, 41, 47, 66, 67, 69, 70, 81, 82, 83, 84, 88, 9597.50 ± 0.2992.37 ± 0.0783.21 ± 0.04
3.1201, 2, 3, 4, 5, 9, 15, 16, 19, 24, 26, 28, 40, 52, 61, 67, 70, 76, 82, 8997.34 ± 0.1392.90 ± 0.0789.13 ± 0.04
3.3194, 5, 8, 16, 23, 26, 28, 29, 38, 52, 56, 62, 71, 74, 77, 78, 86, 87, 9397.34 ± 0.1393.95 ± 0.8092.21 ± 0.20
3.5202, 4, 14, 16, 20, 21, 22, 28, 36, 37, 40, 41, 43, 46, 47, 51, 52, 58, 62, 7296.90 ± 0.1393.39 ± 0.1390.00 ± 0.60
Coiflets (coif)1244, 5, 8, 11, 14, 15, 16, 19, 28, 34, 35, 49, 50, 52, 58, 61, 62, 64, 68, 70, 84, 85, 89, 9197.77 ± 0.2295.29 ± 0.3689.69 ± 0.07
2142, 4, 12, 28, 34, 41, 50, 52, 59, 62, 64, 76, 83, 9393.86 ± 0.1191.38 ± 0.4789.87 ± 0.42
3182, 4, 5, 6, 7, 8, 26, 29, 43, 45, 52, 56, 62, 64, 70, 76, 79, 9696.32 ± 0.2292.17 ± 0.5489.89 ± 0.27
4237, 8, 17, 18, 21, 26, 28, 34, 35, 37, 42, 45, 49, 56, 57, 62, 65, 67, 80, 86, 88, 90, 9597.57 ± 0.0992.63 ± 0.1185.87 ± 1.27
5194, 5, 18, 26, 27, 30, 36, 50, 51, 52, 59, 69, 70, 71, 72, 77, 83, 87, 8997.54 ± 0.1192.37 ± 0.2787.30 ± 0.47
Daubechies (db)1215, 8, 17, 18, 19, 25, 26, 28, 35, 36, 37, 44, 45, 53, 56, 65, 72, 77, 82, 83, 9594.60 ± 0.5187.48 ± 0.1380.71 ± 0.31
2222, 4, 11, 12, 16, 17, 20, 24, 27, 32, 37, 45, 52, 58, 59, 60, 64, 67, 78, 80, 88, 9696.70 ± 0.4991.27 ± 0.0284.02 ± 0.65
3233, 4, 7, 15, 19, 21, 25, 28, 35, 39, 42, 50, 51, 52, 53, 67, 68, 75, 78, 80, 81, 93, 9697.43 ± 0.2293.04 ± 0.6085.94 ± 0.11
4194, 12, 23, 26, 28, 29, 33, 34, 36, 37, 50, 53, 58, 60, 65, 67, 70, 87, 9597.34 ± 0.0992.59 ± 0.2788.44 ± 0.27
5213, 4, 8, 14, 23, 26, 28, 32, 40, 41, 46, 48, 50, 58, 65, 70, 72, 75, 82, 83, 8798.19 ± 0.4294.11 ± 0.5890.67 ± 0.04
6192, 4, 10, 21, 26, 28, 32, 38, 39, 41, 46, 47, 64, 72, 76, 82, 83, 91, 9395.92 ± 0.1691.38 ± 0.2589.24 ± 0.60
7184, 6, 7, 8, 17, 19, 21, 23, 28, 32, 33, 36, 41, 58, 61, 64, 72, 8595.63 ± 0.2090.29 ± 0.5686.05 ± 0.22
8271, 2, 4, 8, 12, 16, 22, 24, 25, 27, 28, 29, 32, 42, 50, 54, 56, 67, 73, 76, 82, 87, 88, 92, 94, 95, 9697.88 ± 0.3393.62 ± 0.0285.60 ± 0.67
9243, 4, 7, 8, 10, 22, 28, 29, 32, 40, 41, 43, 45, 48, 50, 54, 56, 59, 74, 79, 87, 90, 91, 9497.88 ± 0.2293.93 ± 0.1686.07 ± 0.20
10221, 4, 8, 17, 18, 19, 24, 26, 28, 36, 38, 42, 44, 46, 50, 68, 77, 82, 86, 87, 88, 9097.88 ± 0.1192.30 ± 0.1186.41 ± 0.13
Fejer-Krovkin (fk)4247, 8, 9, 16, 17, 20, 21, 23, 26, 28, 35, 37, 42, 45, 57, 61, 65, 69, 75, 79, 80, 86, 88, 9696.99 ± 0.2290.31 ± 0.4782.19 ± 0.71
6204, 5, 10, 11, 23, 27, 28, 30, 31, 32, 33, 37, 44, 52, 53, 55, 57, 58, 71, 7597.52 ± 0.0993.62 ± 0.1385.69 ± 0.20
8233, 4, 8, 11, 21, 28, 32, 35, 38, 39, 40, 42, 45, 56, 67, 72, 73, 77, 81, 83, 91, 92, 9497.97 ± 0.2093.21 ± 0.4786.58 ± 0.02
14213, 10, 11, 15, 17, 24, 26, 28, 29, 31, 32, 49, 53, 67, 68, 75, 76, 77, 80, 82, 9197.70 ± 0.0493.15 ± 0.1887.17 ± 0.33
18234, 5, 9, 10, 12, 14, 17, 21, 25, 26, 27, 28, 35, 52, 62, 63, 65, 69, 70, 71, 73, 74, 7997.99 ± 0.1194.75 ± 0.1190.33 ± 0.18
22212, 28, 29, 35, 36, 40, 41, 45, 46, 51, 53, 58, 64, 65, 69, 70, 79, 85, 86, 91, 9496.81 ± 0.1690.96 ± 0.5686.03 ± 0.31
Symlets (sym)2252, 4, 11, 13, 17, 22, 24, 28, 29, 35, 38, 64, 65, 69, 71, 72, 75, 76, 79, 80, 81, 87, 90, 92, 9398.28 ± 0.1894.13 ± 0.4087.14 ± 0.13
3212, 4, 6, 7, 8, 9, 19, 25, 28, 32, 35, 46, 52, 59, 65, 71, 75, 81, 84, 85, 8997.21 ± 0.1192.95 ± 0.4785.83 ± 0.45
4203, 5, 11, 15, 19, 26, 27, 28, 35, 50, 60, 68, 69, 70, 74, 75, 76, 85, 87, 9697.21 ± 0.5692.28 ± 0.3184.22 ± 0.29
5243, 4, 6, 11, 12, 14, 19, 26, 28, 31, 32, 39, 46, 48, 52, 47, 65, 66, 67, 83, 86, 88, 90, 9197.68 ± 0.1393.24 ± 0.6085.65 ± 0.18
6204, 11, 27, 28, 35, 37, 38, 39, 53, 54, 56, 58, 69, 70, 71, 72, 77, 83, 86, 9197.23 ± 0.1393.15 ± 0.2985.71 ± 0.11
7231, 4, 5, 8, 16, 21, 22, 26, 28, 34, 36, 40, 44, 45, 46, 47, 50, 52, 57, 80, 82, 87, 9499.55 ± 0.1397.54 ± 0.2593.97 ± 0.04
8234, 13, 16, 28, 30, 35, 48, 50, 52, 63, 65, 67, 68, 69, 73, 76, 79, 82, 88, 90, 92, 93, 9598.06 ± 0.2995.47 ± 0.1690.42 ± 0.58
Table 7. The performance of the classifiers for selected features by the PSO-SVM wrapper approach.
Table 7. The performance of the classifiers for selected features by the PSO-SVM wrapper approach.
WaveletsNum.Selected FeaturesAccuracy (%)
CNNSVMKNN
Biorthogonal (bior)1.1195, 11, 19, 26, 28, 29, 32, 40, 41, 45, 56, 62, 63, 64, 70, 75, 85, 87, 9691.71 ± 0.0484.93 ± 1.0980.40 ± 0.36
1.3213, 4, 5, 9, 11, 14, 15, 19, 21, 25, 27, 39, 46, 62, 63, 74, 76, 77, 84, 91, 9392.03 ± 0.2981.96 ± 0.0274.28 ± 0.18
1.5233, 5, 8, 9, 10, 26, 28, 29, 31, 34, 37, 40, 41, 49, 56, 59, 62, 66, 79, 81, 87, 91, 9292.90 ± 0.2584.64 ± 0.5478.17 ± 0.02
2.2171, 2, 4, 16, 26, 27, 28, 52, 53, 65, 69, 70, 78, 85, 86, 87, 9691.11 ± 0.0986.85 ± 0.7885.31 ± 0.20
2.4183, 4, 5, 6, 17, 24, 26, 27, 28, 36, 40, 47, 53, 69, 76, 79, 80, 8792.41 ± 0.0286.42 ± 0.1385.55 ± 0.04
2.6154, 16, 44, 47, 50, 52, 57, 60, 69, 70, 73, 76, 80, 81, 8988.14 ± 0.1882.34 ± 0.1679.19 ± 0.63
2.81912, 16, 23, 28, 35, 41, 42, 43, 45, 51, 59, 66, 67, 72, 73, 76, 79, 90, 9693.10 ± 0.4587.32 ± 0.0281.49 ± 0.22
3.1201, 9, 12, 14, 15, 16, 20, 24, 28, 29, 31, 32, 34, 57, 63, 77, 87, 88, 92, 9492.96 ± 0.0986.42 ± 0.1382.47 ± 0.20
3.3211, 2, 6, 12, 13, 29, 35, 37, 40, 44, 52, 53, 57, 62, 64, 67, 70, 77, 83, 86, 8992.79 ± 0.0985.29 ± 0.3378.23 ± 0.47
3.5192, 5, 7, 12, 14, 16, 28, 31, 38, 42, 51, 60, 66, 67, 72, 79, 83, 89, 9292.74 ± 0.2087.16 ± 0.6583.21 ± 0.27
Coiflets (coif)1193, 4, 10, 11, 17, 35, 36, 37, 50, 56, 59, 61, 68, 73, 84, 88, 89, 92, 9592.23 ± 0.0985.29 ± 0.2279.71 ± 0.22
2224, 13, 14, 15, 24, 28, 33, 35, 36, 39, 50, 52, 53, 58, 66, 68, 69, 76, 87, 89, 93, 9692.94 ± 0.2787.92 ± 0.1885.58 ± 0.16
3201, 4, 6, 7, 13, 19, 21, 26, 28, 29, 45, 46, 52, 63, 68, 71, 76, 79, 82, 9093.17 ± 0.0489.44 ± 0.2085.64 ± 0.02
4174, 12, 17, 20, 23, 26, 36, 46, 50, 52, 53, 56, 57, 58, 72, 79, 9593.43 ± 0.0189.33 ± 0.0985.71 ± 0.02
5178, 16, 18, 20, 23, 24, 28, 31, 32, 33, 39, 42, 43, 49, 72, 76, 7991.58 ± 0.2784.44 ± 0.1678.83 ± 0.02
Daubechies (db)1164, 5, 11, 17, 28, 35, 37, 39, 44, 51, 52, 56, 65, 88, 91, 9593.70 ± 0.4284.84 ± 0.3380.82 ± 0.13
2198, 16, 17, 20, 27, 28, 35, 36, 37, 48, 52, 57, 60, 63, 64, 65, 80, 87, 8891.71 ± 1.0785.71 ± 0.2081.83 ± 0.11
3194, 7, 8, 24, 26, 28, 32, 34, 35, 36, 46, 53, 54, 58, 61, 76, 83, 87, 9191.40 ± 0.0285.60 ± 0.4778.81 ± 0.56
4234, 12, 17, 22, 23, 26, 35, 40, 41, 42, 44, 50, 51, 54, 55, 60, 62, 63, 65, 67, 72, 76, 9494.21 ± 0.1189.55 ± 0.3683.54 ± 0.71
5214, 8, 17, 22, 23, 24, 26, 28, 31, 40, 43, 58, 64, 65, 67, 72, 73, 80, 81, 88, 8993.79 ± 0.2089.75 ± 0.5685.69 ± 0.18
6244, 6, 16, 17, 19, 21, 22, 26, 28, 29, 32, 35, 38, 41, 43, 53, 75, 76, 79, 80, 82, 86, 90, 9194.87 ± 0.0490.15 ± 0.2987.00 ± 0.18
7195, 7, 10, 12, 14, 15, 21, 26, 28, 37, 40, 50, 57, 59, 65, 83, 89, 92, 9493.28 ± 0.3886.87 ± 0.0981.77 ± 0.25
8174, 14, 17, 18, 19, 24, 26, 28, 29, 32, 34, 50, 62, 65, 66, 73, 8091.56 ± 0.2584.10 ± 0.2978.21 ± 0.16
9185, 22, 26, 28, 35, 36, 37, 39, 44, 46, 47, 56, 62, 77, 86, 88, 92, 9392.65 ± 0.0184.30 ± 0.0279.04 ± 0.33
10225, 9, 11, 12, 14, 15, 16, 26, 36, 37, 40, 46, 47, 53, 55, 60, 68, 70, 80, 91, 92, 9393.72 ± 0.0785.87 ± 0.1378.17 ± 0.31
Fejer-Krovkin (fk)41910, 12, 17, 28, 29, 32, 41, 43, 58, 72, 76, 78, 83, 85, 87, 88, 93, 94, 9691.89 ± 0.4782.63 ± 0.3679.01 ± 0.13
6204, 7, 17, 22, 25, 32, 34, 35, 37, 39, 44, 47, 48, 53, 54, 61, 64, 80, 86, 9090.60 ± 0.3882.50 ± 0.1177.09 ± 0.27
8212, 6, 8, 16, 28, 36, 39, 40, 42, 48, 51, 59, 60, 64, 68, 69, 70, 71, 72, 91, 9392.23 ± 0.2585.93 ± 0.4279.24 ± 0.02
14163, 15, 19, 26, 28, 31, 32, 48, 52, 55, 57, 58, 62, 72, 77, 8590.93 ± 0.0784.24 ± 0.3879.44 ± 0.16
18205, 8, 15, 16, 24, 26, 31, 34, 37, 49, 52, 61, 62, 76, 79, 80, 85, 88, 92, 9492.88 ± 0.1186.47 ± 0.2980.49 ± 1.34
22214, 7, 8, 9, 20, 21, 22, 25, 26, 34, 45, 51, 52, 62, 64, 69, 75, 76, 81, 84, 9592.74 ± 0.3687.29 ± 0.4581.98 ± 0.27
Symlets (sym)2204, 5, 9, 11, 21, 23, 26, 28, 29, 32, 35, 41, 46, 52, 67, 69, 73, 79, 80, 8891.89 ± 0.0283.86 ± 0.0278.41 ± 0.63
3192, 4, 5, 12, 15, 21, 28, 37, 43, 44, 45, 51, 57, 62, 65, 73, 76, 94, 9592.61 ± 0.3887.18 ± 0.3382.74 ± 0.13
4183, 17, 19, 21, 26, 28, 40, 52, 54, 60, 65, 72, 74, 79, 85, 90, 92, 9591.56 ± 0.1385.89 ± 0.1680.75 ± 0.49
52310, 16, 17, 18, 21, 22, 26, 28, 29, 31, 39, 56, 64, 72, 75, 81, 82, 85, 91, 92, 93, 95, 9693.41 ± 0.0987.96 ± 0.3380.96 ± 0.25
6164, 16, 17, 21, 27, 28, 37, 45, 60, 61, 66, 67, 68, 77, 80, 9291.36 ± 0.1686.96 ± 0.1183.28 ± 0.22
7185, 9, 16, 19, 20, 28, 41, 47, 50, 52, 64, 67, 70, 78, 80, 78, 90, 9590.33 ± 0.0283.59 ± 0.0279.50 ± 0.25
8154, 23, 28, 29, 32, 39, 46, 51, 52, 60, 67, 70, 76, 78, 8890.13 ± 0.1886.65 ± 0.2584.04 ± 1.09
Table 8. Comparison of the performance of various classifiers created with the selected features by the proposed FPA-SVM and PSO-SVM wrapper approaches for wavelet families providing the best classification performance.
Table 8. Comparison of the performance of various classifiers created with the selected features by the proposed FPA-SVM and PSO-SVM wrapper approaches for wavelet families providing the best classification performance.
Feature Selection MethodWaveletClassifierAccuracy (%)Precision (%)Recall (%)F1 score (%)
PSO-SVMdb6CNN94.8795.0694.6194.77
SVM90.1590.2189.6289.73
KNN87.0087.6485.8686.26
FPA-SVM
(proposed)
sym7CNN99.5599.5299.6099.56
SVM97.5497.5896.4496.92
KNN93.9794.2593.8293.84
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yağ, İ.; Altan, A. Artificial Intelligence-Based Robust Hybrid Algorithm Design and Implementation for Real-Time Detection of Plant Diseases in Agricultural Environments. Biology 2022, 11, 1732. https://doi.org/10.3390/biology11121732

AMA Style

Yağ İ, Altan A. Artificial Intelligence-Based Robust Hybrid Algorithm Design and Implementation for Real-Time Detection of Plant Diseases in Agricultural Environments. Biology. 2022; 11(12):1732. https://doi.org/10.3390/biology11121732

Chicago/Turabian Style

Yağ, İlayda, and Aytaç Altan. 2022. "Artificial Intelligence-Based Robust Hybrid Algorithm Design and Implementation for Real-Time Detection of Plant Diseases in Agricultural Environments" Biology 11, no. 12: 1732. https://doi.org/10.3390/biology11121732

APA Style

Yağ, İ., & Altan, A. (2022). Artificial Intelligence-Based Robust Hybrid Algorithm Design and Implementation for Real-Time Detection of Plant Diseases in Agricultural Environments. Biology, 11(12), 1732. https://doi.org/10.3390/biology11121732

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