# A Curvelet-SC Recognition Method for Maize Disease

Authors: Jing Luo; Shuze Geng; Chunbo Xiu; Dan Song; Tingting Dong
Journal: Journal of Electrical and Computer Engineering
 (2015)

---

## Abstract

Because the corn vein and noise influence the contour extraction of the maize leaf disease, we put forward a new recognition algorithm based on Curvelet and Shape Context (SC). This method can improve the speed and accuracy of maize leaf disease recognition. Firstly, we use Seeded Regional Growing (SRG) algorithm to segment the maize leaf disease image. Secondly, Curvelet Modulus Correlation (CMC) method is put forward to extract the effective contour of maize leaf disease. Thirdly, we combine CMC with the SC algorithm to obtain the histogram features and then use these features we obtain to calculate the similarities between the template image and the target image. Finally, we adoptn-fold cross-validation algorithm to recognize diseases on maize leaf disease database. Experimental results show that the proposed algorithm can recognize 6 kinds of maize leaf diseases accurately and achieve the accuracy of 94.446%. Meanwhile this algorithm has guiding significance for other diseases recognition to an extent.


---

## Body

## 1. Introduction

In recent years, maize leaf diseases have a great influence on the yield of maize and quality. Therefore, an important problem caused by them on agriculture shows that it is urgent to protect maize from these diseases and recognize the species of those diseases. There are 6 common maize diseases: Rust spots, Leaf blight, Gray leaf spot, Brown patch,Curvularia leaf spot, and Small spot disease. Their images are shown in Figure 1. At present, the extracted disease features include color features [1], texture features, and shape features [2–5]. The color features include RGB color feature, HIS color feature, RGB color histogram, the standard deviation, the mean, the variance, and the correlation of RGB. The color features are easily influenced by external environmental light, which makes the disease recognition difficult. The color features of diseases are suitable for transient diseases segmentation but not suitable for long-term diseases segmentation. The texture features include the inertia value of gray level cooccurrence matrix, norm of the wavelet transform coefficient, chromaticity moments, and color cooccurrence texture. Compared with the color features of diseases, texture features are not easily affected by external light and are more stable for diseases recognition. The shape features include the area, perimeter, roundness, eccentricity, and rectangle degree. The shape features are suited to some conditions where the diseases have obvious shape features and the edge between disease and healthy tissue is clear. The shape features matching and recognition algorithms based on the transform domain are invariant in photometric transformation, convolution transformation, and rotation transformation, which are not affected by the external environment and perform stably [6]. Correspondingly, the disease recognition methods involve Artificial Neural Network (ANN) method [7, 8], Support Vector Machine (SVM) method [9], and Template Match (TM) method [5]. As for ANN method and SVM method, with the increase of the number of samples, the disease recognition model has become more complex, the network structure and the model parameters become more complicated, and the recognition precision and its speed will be affected.Figure 1
Maize leaf diseases images.
(a)
Leaf blight
(b)
Brown spot
(c)
Gray leaf spot
(d)
Rust spot
(e)
Curvularia leaf spot
(f)
Small spotHuang et al. [5] proposed a Principal Component Analysis-Shape Context (PCA-SC) method that combines SC and Principal Component Analysis (PCA). The algorithm transforms the characteristic matrix extracted by SC algorithm into the corresponding covariance matrix firstly. Secondly, PCA algorithm is used to reduce the dimension of the covariance matrix according to the characteristic value order from big to small so as to form a new feature matrix to implement shape matching and object recognition. The algorithm has strong anti-interference ability, high recognition rate, and high recognition accuracy. But it is hard to identify the similar edge.Target texture edge feature can be described more accurately by Curvelet in multiscale and multidirection [9–11]. Based on discussion above, this paper puts forward a maize leaf diseases recognition mode which is based on Curvelet and SC for these 6 common diseases of maize: Rust spots, Leaf blight, Gray leaf spot, Brown patch,Curvularia leaf spot, and Small spot disease. First of all, this paper puts forward CMC method to realize the extraction of the effective contour of maize leaf disease. Then we combine CMC with the SC algorithm [5] to obtain the histogram features of this disease image and use the matching degree between the target image shape histogram and the template image shape histogram to recognize the disease. Finally, we adopt n-fold cross-validation algorithm to recognize diseases on maize leaf disease database. Experimental results show that the algorithm we put forward can recognize those maize leaf diseases accurately. Therefore, this algorithm has guiding significance for maize leaf disease recognition.


## 2. Maize Leaf Disease Recognition Method

### 2.1. Disease Image Graying

In order to improve the speed of maize leaf disease recognition, we need to transform the color disease image to gray disease image. According to the sensitivity of human eyes to the image, we adopt formula (1) for image transformation. Consider(1)fi,j=0.3Ri,j+0.59Gi,j+0.11Bi,j.


### 2.2. Disease Image Segmentation

This paper adopts the SRG [12–14] to segment the maize leaf disease image. The basic principles of SRG algorithm are as follows.Step 1.
According to the relationship between the image pixels, the setS of seeds is selected.Step 2.
The setS is marked as M (M1,M2,M3,…,Mi,…,Mn). Seeds of the group are placed on the stack.Step 3.
The thresholdTi of growth center is calculated with Otsu, and R is used for the termination condition of region growing.Step 4.
Bottom pointer and top pointer of stack are compared, and the seeds at the bottom of stack are selected as the seeds of region growing when bottom pointer is less than top pointer.Step 5.
The points are added into the seedsM when threshold is less than termination condition (Ti<R) in the four neighborhoods of the seed points.Step 6.
Steps 4 and 5 are repeated untilTi>R.Maize leaf disease image is shown in Figure2, and the gray image of maize leaf disease is shown in Figure 3. The segmentation result of Figure 3 is shown in Figure 4.Figure 2
Maize leaf disease original image.Figure 3
Gray maize leaf disease image.Figure 4
Segmentation result of maize leaf disease.


### 2.3. Contour Extraction of Maize Leaf Disease Using Curvelet Modulus Correlation

Compared with wavelet, target texture feature can be described more accurately by Curvelet in multiscale and multidirection. In order to reduce the interference of the contour extraction of disease image, the CMC method which is used for obtaining the disease contour is presented. The method is based on the amplitude of the target edge that is large at adjoining scales, while the amplitude of the noise edge can attenuate rapidly at adjoining scales. According to the size of the maize leaf disease, we divide the image into 5 scales. The correlated coefficient between the third scale and the second scale is described as(2)corS,Ni,j=BS,Ni,j×BS-1,N′i′,j′,where(3)S=3,N=1,K,P3,N′=N2,N′=1,…,P2,i′=i2,i′=1,…,16or18,i=1,…,32or34;j′=j=22.B(N,S) is the Curvelet coefficient at direction N and scale S. P3 is the direction number of the third scale. P2 is the direction number of the second scale.The concrete implementation steps are as follows.Step 1.
The correlated modulus between the third scale and the second scale is calculated.Step 2.
Formula (4) is adopted to detect the edge points of the third scale:(4)corS,2N-1i,j-corS,2Ni,j>Ts,s=3,where Ts is the threshold.Step 3.
The following rules are used to detect the edge points at the second scale.
If there are 3 or more than 3 correlation coefficient matrices belonging to the edge, the corresponding matrix is the edge. Otherwise, it is noise.Step 4.
The following rules are used to detect the edge points of the first scale, the fourth scale, and the fifth scale:(5)BS,N>Ps,s=1,4,5,where Ps is the threshold value.Step 5.
The edge points of each scale are reconstituted to accomplish the edge detection.Figure6 is the detection result of Figure 5. This result shows that CMC method can extract the disease contour clearly. It is good for subsequent disease feature extraction.Figure 5
Segmentation result of Gray leaf spot.Figure 6
Detection result of Figure5 by CMC.


### 2.4. Feature Extraction of Maize Leaf Disease

The SC method presented by Belongie [15] is an operator which is used for describing the shape. Its core idea is point set matching.The flow of maize leaf disease recognition based on SC is as follows.Step 1.
The obtained edge detecting point setx,y is transformed into log-polar (log⁡r,θ).
The transformation relation between rectangular coordinates and log-polar is as follows:(6)r=x2+y2;ρ=log⁡x2+y2;θ=arctanyx.Step 2.
The points are calculated in log-polar.
The obtained result in the polar coordinate is transformed into 12 directions, and each direction is 30 degrees.log2r is the step, and each direction is divided into 5 portions. So the whole space is divided into 60 bins. For a fixed point of the contour called pi, its property is described by the number of pixels in each bin. Thus, every feature point can be described by a row vector of 60 dimensions (characteristic attributes). If the number of the edge points is n, the disease feature can be described by a matrix (n×60). In this paper n equals 100.
Figure7 shows the log-polar coordinate.Figure 7
Schematic diagram of log-polar coordinate.Step 3.
The similarity between the template image and the target image is calculated.The similarity between the pointPi of the model image and the point Qj of the target image is calculated by(7)Cij=CPi,Qj=12∑k=1Khik-hjk2hik+hjk,where hik is the corresponding histogram of point Pi, hjk is the corresponding histogram of point Qj, and K=60.IfCij is smaller, the point Pi and the point Qj are more matching.Thus, formula (8) is used to calculate the similarity between the model image and the target image. Consider(8)C=∑i=1n∑j=1mCij/mn.If the number of target disease image contours isn and the number of model disease image contours is m, we can get n×m matching. In this paper, the threshold T1 is set to 50. If the matched-degree C is less than T1, the model image and the target image are matched; otherwise, the model image and the target image are not matched.Figure8 shows a feature point of 60 characteristic attributes in the log-polar coordinates.Figure 8
60 characteristic attributes of a feature point.Figure9 shows 100-feature extraction of the Gray leaf spot disease (Figure 6). Each feature point has 60 attributes, where x represents the number of feature points and y represents 60 bins. z represents the number of pixels of each bin.Figure 9
Feature points and characteristic attributes.The Leaf blight and Gray leaf spot segmentation and feature extraction are shown in Figures10 and 11.Figure 10
Leaf blight and Gray leaf spot segmentation results.
(a)
(b)Figure 11
Feature extraction results of Gray leaf spot and Leaf blight disease.Figure12 shows the matching results of Figures 10(a) and 10(b).Figure 12
The final matching results.From Figure12, the matching results are 31% (31%<50%); the algorithm can distinguish between the Leaf blight and the Gray leaf spot disease.


## 2.1. Disease Image Graying

In order to improve the speed of maize leaf disease recognition, we need to transform the color disease image to gray disease image. According to the sensitivity of human eyes to the image, we adopt formula (1) for image transformation. Consider(1)fi,j=0.3Ri,j+0.59Gi,j+0.11Bi,j.


## 2.2. Disease Image Segmentation

This paper adopts the SRG [12–14] to segment the maize leaf disease image. The basic principles of SRG algorithm are as follows.Step 1.
According to the relationship between the image pixels, the setS of seeds is selected.Step 2.
The setS is marked as M (M1,M2,M3,…,Mi,…,Mn). Seeds of the group are placed on the stack.Step 3.
The thresholdTi of growth center is calculated with Otsu, and R is used for the termination condition of region growing.Step 4.
Bottom pointer and top pointer of stack are compared, and the seeds at the bottom of stack are selected as the seeds of region growing when bottom pointer is less than top pointer.Step 5.
The points are added into the seedsM when threshold is less than termination condition (Ti<R) in the four neighborhoods of the seed points.Step 6.
Steps 4 and 5 are repeated untilTi>R.Maize leaf disease image is shown in Figure2, and the gray image of maize leaf disease is shown in Figure 3. The segmentation result of Figure 3 is shown in Figure 4.Figure 2
Maize leaf disease original image.Figure 3
Gray maize leaf disease image.Figure 4
Segmentation result of maize leaf disease.


## 2.3. Contour Extraction of Maize Leaf Disease Using Curvelet Modulus Correlation

Compared with wavelet, target texture feature can be described more accurately by Curvelet in multiscale and multidirection. In order to reduce the interference of the contour extraction of disease image, the CMC method which is used for obtaining the disease contour is presented. The method is based on the amplitude of the target edge that is large at adjoining scales, while the amplitude of the noise edge can attenuate rapidly at adjoining scales. According to the size of the maize leaf disease, we divide the image into 5 scales. The correlated coefficient between the third scale and the second scale is described as(2)corS,Ni,j=BS,Ni,j×BS-1,N′i′,j′,where(3)S=3,N=1,K,P3,N′=N2,N′=1,…,P2,i′=i2,i′=1,…,16or18,i=1,…,32or34;j′=j=22.B(N,S) is the Curvelet coefficient at direction N and scale S. P3 is the direction number of the third scale. P2 is the direction number of the second scale.The concrete implementation steps are as follows.Step 1.
The correlated modulus between the third scale and the second scale is calculated.Step 2.
Formula (4) is adopted to detect the edge points of the third scale:(4)corS,2N-1i,j-corS,2Ni,j>Ts,s=3,where Ts is the threshold.Step 3.
The following rules are used to detect the edge points at the second scale.
If there are 3 or more than 3 correlation coefficient matrices belonging to the edge, the corresponding matrix is the edge. Otherwise, it is noise.Step 4.
The following rules are used to detect the edge points of the first scale, the fourth scale, and the fifth scale:(5)BS,N>Ps,s=1,4,5,where Ps is the threshold value.Step 5.
The edge points of each scale are reconstituted to accomplish the edge detection.Figure6 is the detection result of Figure 5. This result shows that CMC method can extract the disease contour clearly. It is good for subsequent disease feature extraction.Figure 5
Segmentation result of Gray leaf spot.Figure 6
Detection result of Figure5 by CMC.


## 2.4. Feature Extraction of Maize Leaf Disease

The SC method presented by Belongie [15] is an operator which is used for describing the shape. Its core idea is point set matching.The flow of maize leaf disease recognition based on SC is as follows.Step 1.
The obtained edge detecting point setx,y is transformed into log-polar (log⁡r,θ).
The transformation relation between rectangular coordinates and log-polar is as follows:(6)r=x2+y2;ρ=log⁡x2+y2;θ=arctanyx.Step 2.
The points are calculated in log-polar.
The obtained result in the polar coordinate is transformed into 12 directions, and each direction is 30 degrees.log2r is the step, and each direction is divided into 5 portions. So the whole space is divided into 60 bins. For a fixed point of the contour called pi, its property is described by the number of pixels in each bin. Thus, every feature point can be described by a row vector of 60 dimensions (characteristic attributes). If the number of the edge points is n, the disease feature can be described by a matrix (n×60). In this paper n equals 100.
Figure7 shows the log-polar coordinate.Figure 7
Schematic diagram of log-polar coordinate.Step 3.
The similarity between the template image and the target image is calculated.The similarity between the pointPi of the model image and the point Qj of the target image is calculated by(7)Cij=CPi,Qj=12∑k=1Khik-hjk2hik+hjk,where hik is the corresponding histogram of point Pi, hjk is the corresponding histogram of point Qj, and K=60.IfCij is smaller, the point Pi and the point Qj are more matching.Thus, formula (8) is used to calculate the similarity between the model image and the target image. Consider(8)C=∑i=1n∑j=1mCij/mn.If the number of target disease image contours isn and the number of model disease image contours is m, we can get n×m matching. In this paper, the threshold T1 is set to 50. If the matched-degree C is less than T1, the model image and the target image are matched; otherwise, the model image and the target image are not matched.Figure8 shows a feature point of 60 characteristic attributes in the log-polar coordinates.Figure 8
60 characteristic attributes of a feature point.Figure9 shows 100-feature extraction of the Gray leaf spot disease (Figure 6). Each feature point has 60 attributes, where x represents the number of feature points and y represents 60 bins. z represents the number of pixels of each bin.Figure 9
Feature points and characteristic attributes.The Leaf blight and Gray leaf spot segmentation and feature extraction are shown in Figures10 and 11.Figure 10
Leaf blight and Gray leaf spot segmentation results.
(a)
(b)Figure 11
Feature extraction results of Gray leaf spot and Leaf blight disease.Figure12 shows the matching results of Figures 10(a) and 10(b).Figure 12
The final matching results.From Figure12, the matching results are 31% (31%<50%); the algorithm can distinguish between the Leaf blight and the Gray leaf spot disease.


## 3. Experimental Results and Analysis

Experimental environment is as follows: the operating system is windows 8, Pentium (R), quad core, 8 G internal storage, and the software is matlabR2011b. The experimental images are from the maize leaf disease database established by Institute of Crop Science, Chinese Academy of Agricultural Sciences.


### 3.1. Experimental Subjects

744 images are selected to verify the accuracy of the proposed algorithm. Every disease responds to 124 images, and there are 6 kinds of diseases. There are four templates for each type of disease and 120 test samples. Figure13 shows four templates of Brown spot.Figure 13
Brown spot samples.


### 3.2. Experimental Results

In the experiment, 6-fold cross validation is used to train and test the proposed algorithm. The 744 plant disease images are divided into 6 independent subsets, and each subset contains 124 plant disease images. Each subset has 4 template images and 120 training samples. The training and testing process is repeated 6 times, and finally the 6-time test results are averaged. The recognition rate between the same kinds of diseases and different kinds of diseases is calculated, respectively.Figure14 shows the results between different kinds of diseases. The Brown spot is selected as template. The remaining 5 kinds of diseases are used as test samples.Figure 14
Experiment results between different kinds of diseases.Figure15 shows results between the same diseases. The Leaf blight, Gray leaf spot, Rust spot and Curvularia leaf spot are selected as templates. The corresponding test samples are used as test samples. For the same diseases between experiments.Figure 15
Experiment results between the same diseases.Figure16 shows recognition rate between the same diseases corresponding to 4 templates.Figure 16
Recognition rates between the same diseases.Figure17 shows recognition rate among the different classes diseases corresponding to 4 templates.Figure 17
Recognition rate among the different classes diseases.From the above, the results are shown in Table1.Table 1
Maize leaf disease recognition accuracy (%).
Recognition rates between the same diseases
Recognition rate among the different classes diseases
The total recognition rate
92.225%
96.667%
94.446%From the experiment results, recognition rates between the same diseases and recognition rate among the different classes diseases reach above 90%. In the same diseases recognition, the recognition rate of Brown patch is the highest (97.385%), and the Small spot disease has the lowest recognition rate (88.673%). In the different classes diseases recognition, Brown patch has the highest recognition rate (98.658%), and Gray leaf spot has the lowest recognition rate (93.434%). Some edges extracted from Gray leaf spot andCurvularia leaf spot are similar, some edges extracted from Small spot and Gray leaf spot are also similar, and some edges extracted from Brown patch,Curvularia leaf spot, and Gray leaf spot are also similar. These similarities lead to the error recognition.


## 3.1. Experimental Subjects

744 images are selected to verify the accuracy of the proposed algorithm. Every disease responds to 124 images, and there are 6 kinds of diseases. There are four templates for each type of disease and 120 test samples. Figure13 shows four templates of Brown spot.Figure 13
Brown spot samples.


## 3.2. Experimental Results

In the experiment, 6-fold cross validation is used to train and test the proposed algorithm. The 744 plant disease images are divided into 6 independent subsets, and each subset contains 124 plant disease images. Each subset has 4 template images and 120 training samples. The training and testing process is repeated 6 times, and finally the 6-time test results are averaged. The recognition rate between the same kinds of diseases and different kinds of diseases is calculated, respectively.Figure14 shows the results between different kinds of diseases. The Brown spot is selected as template. The remaining 5 kinds of diseases are used as test samples.Figure 14
Experiment results between different kinds of diseases.Figure15 shows results between the same diseases. The Leaf blight, Gray leaf spot, Rust spot and Curvularia leaf spot are selected as templates. The corresponding test samples are used as test samples. For the same diseases between experiments.Figure 15
Experiment results between the same diseases.Figure16 shows recognition rate between the same diseases corresponding to 4 templates.Figure 16
Recognition rates between the same diseases.Figure17 shows recognition rate among the different classes diseases corresponding to 4 templates.Figure 17
Recognition rate among the different classes diseases.From the above, the results are shown in Table1.Table 1
Maize leaf disease recognition accuracy (%).
Recognition rates between the same diseases
Recognition rate among the different classes diseases
The total recognition rate
92.225%
96.667%
94.446%From the experiment results, recognition rates between the same diseases and recognition rate among the different classes diseases reach above 90%. In the same diseases recognition, the recognition rate of Brown patch is the highest (97.385%), and the Small spot disease has the lowest recognition rate (88.673%). In the different classes diseases recognition, Brown patch has the highest recognition rate (98.658%), and Gray leaf spot has the lowest recognition rate (93.434%). Some edges extracted from Gray leaf spot andCurvularia leaf spot are similar, some edges extracted from Small spot and Gray leaf spot are also similar, and some edges extracted from Brown patch,Curvularia leaf spot, and Gray leaf spot are also similar. These similarities lead to the error recognition.


## 4. Conclusion

Because the corn vein and noise influence the contour extraction of the maize leaf disease, we put forward a new maize leaf disease recognition algorithm using the Curvelet technique and SC. This method can improve the speed of maize leaf disease recognition. The Curvelet-SC descriptor was put forward based on SC algorithm and the principle of this algorithm was elaborated in the paper. It combined Curvelet and SC for shape matching and texture recognition.With the experimental analysis, the average accuracy of Curvelet-SC descriptor for maize leaf disease recognition is up to 94.446%. Thus, Curvelet-SC shape descriptor is an efficient method for disease recognition. The proposed algorithm can recognize the maize leaf disease, and it also has guiding significance for other disease recognition to an extent.



---

*Source: 164547-2015-04-05.xml*