• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Early Detection of Colletotrichum Kahawae Disease in Coffee Cherry Based on Computer Vision Techniques

    2024-01-20 13:02:08RaveenaSelvanarayananSurendranRajendranandYouseefAlotaibi

    Raveena Selvanarayanan ,Surendran Rajendran,★ and Youseef Alotaibi

    1Department of Computer Science and Engineering,Saveetha School of Engineering,Saveetha Institute of Medical and Technical Science,Chennai,602105,India

    2Department of Computer Science,College of Computer and Information Systems,Umm Al-Qura University,Makkah,21955,Saudi Arabia

    ABSTRACT Colletotrichum kahawae (Coffee Berry Disease) spreads through spores that can be carried by wind,rain,and insects affecting coffee plantations,and causes 80%yield losses and poor-quality coffee beans.The deadly disease is hard to control because wind,rain,and insects carry spores.Colombian researchers utilized a deep learning system to identify CBD in coffee cherries at three growth stages and classify photographs of infected and uninfected cherries with 93%accuracy using a random forest method.If the dataset is too small and noisy,the algorithm may not learn data patterns and generate accurate predictions.To overcome the existing challenge,early detection of Colletotrichum Kahawae disease in coffee cherries requires automated processes,prompt recognition,and accurate classifications.The proposed methodology selects CBD image datasets through four different stages for training and testing.XGBoost to train a model on datasets of coffee berries,with each image labeled as healthy or diseased.Once the model is trained,SHAP algorithm to figure out which features were essential for making predictions with the proposed model.Some of these characteristics were the cherry’s colour,whether it had spots or other damage,and how big the Lesions were.Virtual inception is important for classification to virtualize the relationship between the colour of the berry is correlated with the presence of disease.To evaluate the model’s performance and mitigate excess fitting,a 10-fold cross-validation approach is employed.This involves partitioning the dataset into ten subsets,training the model on each subset,and evaluating its performance.In comparison to other contemporary methodologies,the model put forth achieved an accuracy of 98.56%.

    KEYWORDS Computer vision;coffee berry disease;colletotrichum kahawae;XG boost;shapley additive explanations

    1 Introduction

    Coffee is a very significant and profitable crop.It is farmed on a big scale(42,000 hectares between 2015 and 2022)as well as a small one(128,000 hectares),with a total annual yield of around 70,000 tons.Arabica coffee is the most extensively consumed type worldwide.In early 90’s Kenya,75% of coffee crop loss was due to CBD (1920),in Ethiopia,CBD caused an estimated 30% crop loss in(1990),In Brazil,CBD caused an estimated 10% crop loss in the 2000 s.CBD is a serious threat to the coffee industry because it leads to a decrease in coffee quality and flavour,and it is important to take steps to control it else,coffee production has to be cut in half [1].Coffee berry disease,also known as Colletotrichum kahawae is a parasite that affects all stages of coffee cherry growth,but it is most visible as microscopic,dark,sunken,brown patches on the green berry that start small and later enlarge in area,eating the entire fruit and rendering it black.CBD first emerges on coffee berries as tiny,water-soaked sores.These lesions might grow and become buried with time.The berries may also dry out and shrivel.These spots can spread throughout the fruit,resulting in a flavourless,dry product[2].Coffee berry disease(CBD)is a devastating fungal disease that can cause significant yield losses in coffee crops.The disease is caused by the fungus Colletotrichum kahawae,which can infect coffee cherries at any stage of development.However,the susceptibility of coffee cherries to CBD varies depending on the stage of development as Greene cherry,yellow cherry,chartreuse cherry,yellow-red cherry,red cherry,and overripe cherry are all varieties of cherries[3].In the proposed work to identify the Colletotrichum kahawae accurately stages are classified into three categories namely Early stages of development (pinhead to pea stage) In the early stages of development,coffee cherries are more resistant to CBD infection as shown in Fig.1.However,if the infection does occur,it can lead to the formation of small,dark spots on the surface of the cherry.These spots can eventually enlarge and coalesce,causing the cherry to rot and fall off the tree,in middle stages of development (pea to cherry stage).In the middle stages of development,coffee cherries are more susceptible to CBD infection.Infection at this stage can cause the formation of large,dark spots on the surface of the cherry.These spots can eventually spread to the entire cherry,causing it to rot and fall off the tree,and in the late stages of development(cherry to the ripe stage),coffee cherries are less susceptible to CBD infection[4].However,if infection does occur,it can cause the formation of dark,sunken lesions on the surface of the cherry.These lesions can eventually lead to the mummification of the cherry,but the beans inside the cherry are usually not affected.CBD is a parasite affecting all phases of coffee cherry development,but most obviously shows as microscopic,dark,sunken,brown patches on the green berry that begin small and then enlarge in area,eating the entire berry to turn black.These spots can spread throughout the fruit,producing a tasteless and dry product[5].Even though this disease does not kill trees,crop losses may surpass 80%.Water,insects,or movement through the coffee crop often distribute the spores during harvest.Colletotrichum kahawae is a more dangerous fungal disease that damages the fruit from the fourth to the twelfth week after blossoming.Producing cherries with black spots that may or may not mature properly,drying,and remaining connected to the branch covers the cherry,decreasing growth and quality.

    The popular machine-learning algorithm XG Boost combines the decision trees to identify the features and improve the accuracy of quick predictions.XGBoost is trained for classifying coffee berries infected with CBD [6].Image segmentation is done using the XGBoost model trained using CBD-infected and non-CBD-infected coffee cherry pictures.The model identified CBD-infected cherries 80.18%accurately.SHAP(Shapley Additive explanations)can be used to interpret machine learning model predictions and extract features from the segmented images.SHAP can identify and explain how each prediction attribute affects a model’s forecast.SHAP identifies the most important elements by explaining how each factor influences model prediction.SHAP could identify key CBD features early in development.The data may improve CBD detection systems.SHAP reasons for the model’s CBD-infected cherry prediction.SHAP listed predictions and their effects.Cherry colour may matter more than size early on.Cherry size may matter more than colour late in development [7].XGBoost and SHAP are used together can help to improve the accuracy of CBD classification and to gain insights into the importance of different features for classification.Virtual inspection is the computer vision technique used to inspect the agricultural environment.It is used to detect defects or abnormalities in CBD.Virtual inception detects diseases in coffee cherries via deep learning and virtual inspection.The virtual inspection system measures coffee cherry foliage colour,texture,and Shape.This data trains a deep-learning network to identify disease patterns.After training the deep learning model with coffee cherry illnesses in new and processed images from XGBoost and SHAP analysis.Feed the model additional photos and ask it to predict the disease.At different phases of coffee cherry development,a virtual examination can identify CBD.CBD symptoms are mild early on,making it hard to detect using typical methods.However,a virtual examination can discover CBD early on,preventing disease propagation and yield losses.CBD symptoms are more apparent in late development,making it easier to identify using traditional methods.Virtual inspections can validate the diagnosis and severity of the infection[8].Virtual inspection will compare the real-time images and processed images to identify cherry is infected or not.

    Figure 1:Coffee berry diseases

    Detecting coffee berry disease (CBD) is still limited by a lack of high-quality data for training,the difficulty of accurately identifying CBD in coffee plants,and the presence of other diseases with similar symptoms.Detection is also not very robust to environmental changes or disease appearance.This can make real-world usage of these models problematic [9].Researchers,data scientists,coffee growers,and other stakeholders must collaborate to close these research gaps.Together,can create more accurate,robust,interpretable,and economical machine learning and computer vision models to lessen CBD’s influence on coffee harvests.Section 2 covers the existing works from various crop and leaf classification and segmentation using deep learning models.

    2 Related Work

    Numerous studies have employed coffee cherries as a means of identifying various disorders.Below,we present a comprehensive elucidation of the literature review,as depicted in Table 1.A novel deep learning approach has been developed for the detection of potato blight in Solanum tuberosum L.and tomato plants[10].CNN and ResNet-9 process 3990 photos to identify diseased tomato and potato leaf regions.With excellent accuracy and the ability to detect diverse blight illnesses,it is robust to noise and imagines oscillations.Weakness is the lack of computational resources needed to train and deploy the method.Sanida et al.proposed an efficient hybrid CNN classification model for tomato crop disease.The model is based on a combination of convolutional neural networks (CNNs) and inception modules [11].CNN used in the model is a 13-layer CNN with 3×3 convolution kernels.The inception module is a building block that combines multiple convolutional layers with different filter sizes.Strength can achieve high accuracy in tomato crop disease classification,and weakness is requiring a large dataset of labelled images for training.Raghavendran et al.proposed a bibliometric analysis of the research on plant crop disease identification using machine learning (ML) and deep learning(DL)[12].The bibliometric analysis aims to identify the trends,gaps,and challenges in this research area.Strength is used to identify the research topics that have been investigated,and weakness does not provide insights into the quality of the research.

    Merga et al.proposed to examine the genetic variability of biennial-bearing coffee(Coffea arabica L.) germplasm and its effect on selection efficiency [13].Coffee trees produce high and low yields biennially.This feature can limit coffee output by causing growers to lose income and yields.The coffee germplasm was evaluated for biennial bearing over,for three years and data on biennial bearing was analyzed using variance component analysis.The strength identifies the most important genetic factors that contribute to biennial bearing,and the weakness is results may not be generalizable to all coffee germplasm.Sudha et al.proposed an early cashew leaf anthracnose detection using image processing and machine learning[14].The fungal disease anthracnose can severely harm cashew crops.Preventing crop losses requires early disease diagnosis and control.Convolutional neural network (CNN) to develop the early detection system,and train on a dataset of images of healthy and diseased cashew leaves.Strength is non-destructive,which means that the leaves do not need to be damaged to be tested,and weakness may not be able to detect all cases of anthracnose disease.Sadiq et al.developed deep learning for potato leaf disease identification and classification[15].A convolutional neural network(CNN)is used to extract features from potato leaf photos and classify them by illness.Deep learningbased potato leaf disease diagnosis and categorization.The research suggests using a CNN to extract features from potato leaf photos and categorize them by disease.The preprocessed images fine-tuned the VGG16 model to categorize new potato leaf images.Strength is a non-destructive which means that the leaves do not need to be damaged to be tested and weakness may not be able to detect all cases of potato leaf diseases.

    Melese established a mathematical model for the coffee berry borer (CBB) that accounts for temperature and rainfall fluctuation[16].The model then optimizes CBB regulation to minimize coffee production impact.The optimal control strategy is found by minimizing a cost function that represents the economic loss due to the CBB.Strength can be used to develop optimal control strategies for minimizing the impact of the CBB on coffee production,and weakness is optimal control strategy may not be feasible to implement in practice.Milkisa et al.proposed a method for visualizing and classifying coffee diseases using machine learning and image processing techniques [17].CNN to classify coffee leaf images by disease.The CNN is a good image classification algorithm.Strength can be used to visualize coffee diseases and identify the affected areas,and the weakness is it may not be able to detect all cases of coffee diseases.Miranda proposed a random forest technique to categorize Landsat 8 OLI images as CBN-infected or healthy [18].Random forest is a tree-based ensemble technique with high accuracy and robustness.Preprocessed Landsat 8 OLI pictures improved contrast and reduced noise.Random forest algorithm was trained on Landsat 8 OLI pictures tagged CBN-infected or healthy.The trained random forest algorithm classified new Landsat 8 OLI pictures.Strength is a remote sensing method,which means that it can be used to monitor large areas of coffee crops,and weakness is may not be able to detect CBN in the early stages of infection.

    Abuhayi et al.developed a CNN-based coffee disease classification method that concatenates features from different CNN layers [19].VGG16,a CNN model pre-trained on a huge dataset of natural images.Coffee leaf photos tagged with disease categories were used to fine-tune the VGG16 model.The paper concatenates CNN layer characteristics and combines CNN layers to improve classification accuracy.Strength is a relatively inexpensive method to implement,and weakness may not be able to detect all cases of coffee diseases.Zenebe et al.characterized and determine the virulence of a Colletotrichum kahawae isolate from Gidami,Western Ethiopia.The isolate was characterized using morphological,cultural,and molecular approaches.Leaf disc inoculation assays were utilized to assess and isolate pathogenicity.The virulence of the isolate was determined by the percentage of leaves that developed symptoms.Strength provides information about the characteristics and virulence of a Colletotrichum kahawae isolate from Gidami,Western Ethiopia,and weakness did not investigate the factors that influence the virulence of the isolate [20].Section 3 discusses the methodology and algorithm.

    3 Materials and Methods

    The proposed model developed a method for analyzing CBD and outlined how the phases of coffee cherry are categorized into two classes(CBD and non-CBD).The novelty of the research work feature selection and feature extraction.XGBoost is a tree-based and extraction-based classification technique used for ranking the feature.It selects the important features and remove less important ones.SHAP analysis for feature selection to calculate the SHAP values for each feature.Select the features with the highest SHAP values.Evaluate the performance of the model using the selected features.XGBoost and SHAP analysis for coffee berry disease early detection will be a set of features that are most relevant to the prediction of the disease and help to prevent significant damage to coffee crops.

    Research work was implemented based on the untrained SHAP model to identify and categories coffee berry diseases from the image’s dataset.The primary characteristics of the dataset are visible as microscopic,dark,sunken,brown blotches on the green berry that begin small and expand in size,devouring the fruit and turning it black.These features are extracted and targeted to detect the proposed SHAP model for classification and detection.If the trained model does not help,then the untrained model is ideal because it can take advantage of all the unique characteristics of the dataset.Image is captured,preprocessed and trained using the SHAP model to analyze whether the coffee cherry is infected or healthy.This initiative highlights and acknowledges the need for creating a quick,cost-effective,and trustworthy system that enables advancement in agriculture.Using machine learning and deep neural networks,this research will create a SHAP-based XGBoost system for disease detection in coffee berries.This approach was created to help farmers who do not wish to journey far in search of plant diseases and their treatments.Therefore,this model may be used by farmers,scientists,and botanists to quickly and accurately diagnose coffee berry illnesses.As shown in Fig.2,compared the importance of the SHAP feature to the XGBoost feature in the“Discussions”section to understand how each of them analyses each segment.

    Figure 2:Working model for coffee berry disease prediction

    3.1 Dataset and Preprocessing

    The research conducted by gathering CBD data from the Zenodo community’s public cherry coffee fruits at various stages.The datasets comprise 640 photos of cherry coffee fruits in various stages of ripening with a 256×256 resolution.Sample 1 is 130,Sample 2 is 160,Sample 3 is 160,Sample 4 is 78,and Healthy cherry is 112,In the proposed work,the coffee berry disease identification stage is split into four groups based on attributes of size,colour,taste,variety,nutritional value,the coffee beans’position inside the pulp,skin of the coffee cherry,and uses.Growth and Variety stages are Greene cherry and yellow cherry(Early Development Stage),chartreuse cherry,and yellow-red cherry(Middle Development Stage),red cherry,and overripe cherry Immature(Early Development Stage)as shown in Table 2.The specimen portraits are from coffee farms in Karnataka-coorg’s Nilgiris District,Yercaud,and Kodaikanal.It is estimated that there are approximately 250 coffee plants on the farms,and the data acquired shows 475 coffee cherries,270 of which are healthy and 205 of which are unhealthy.These cherry samples were collected at random from 15 Coffee plants.On the same day,coffee cherry samples are gathered and photographs are shot with a camera of 2688×1520,or 4 Megapixels.640 coffee cherry photos have been collected from Zenodo communities and nature.All experiments were conducted in a Jupyter Notebook and Python Libraries utilizing the Windows 11 operating system.Feature selection is achieved using the Recursive Feature Elimination technique which can be used to select features by recursively eliminating features that are not important.Hardware components for GPUs (graphics processing units) and TPUs (tensor processing units) are used in cloud computing platform AWS[30].

    Table 2: Attributes of coffee berry disease

    Preprocessing is implemented in Python Libraries for image preprocessing based computer vision tool,sample coffee berry disease and healthy images are loaded.Images are resized using the default size of 256×256 pixels.Noise is removed by applying a Gaussian filter.They are effective in reducing noise in an image by preserving the edge smoothness,adjusting the scale space,and bell-shaped smoothening without overshooting.Contrast enhancement which converts to grayscale and applies histogram equalization.Linear normalization is a very important stage it changes the input image into a range of pixel values represented from 0 to 255,images are too dark or shaded in the normalization process it becomes clearer.

    where the input image is the value of the input channel,min is the minimum value of the input channel in the dataset,max is the minimum value of the input channel in the dataset,and the output channel is the normalized value of the input channel.if the Input image is 100,the min is 0 and the max is 255,then the result will be 127.5.This is because(100-0)/(255-0)=127.5.Data argumentation is achieved using Random cropping: which involves randomly cropping the images of diseased and healthy coffee cherries,image flipping: which involves flipping the images of diseased and healthy coffee cherries horizontally or vertically,and image rotation: which involves rotating the images of diseased and healthy coffee cherries by a random angle all the following techniques helps to improve the performance by reducing overfitting,make more robust by changing the orientation of the images.The recursive Feature Elimination Algorithm is a feature selection method that works by removing the least important features from the dataset.The procedure begins by determining the relevance of each feature by fitting a model to the entire dataset.The model is re-fit once the least essential feature has been eliminated from the data set.This is done until an adequate number of features has been selected are stored in the selected_features list and returned as the output of the algorithm.

    3.2 Spilt Data

    In research,datasets are collected from Immature Cherry 130,Half-ripe Cherry 320,Overripe Cherry 112,and Dry Cherry 78.Training,Testing,Validating,and Holding out progress are involved in all four groups.The training data set is used to train the model and change the weights,while the validation set is used to fine-tune the model’s hyperparameters and prevent over fitting.A testing set is used to evaluate the trained model’s final performance and can assist in determining whether the model is ready for deployment or requires any tweaks.Holdout is often used to evaluate the model’s performance after all necessary adjustments based on the training,validation,and testing sets.Following that,processing and model development will handle missing values and normalization.The proportion is based on dataset availability and the nature of the coffee berry disease problem[31].Group 1 has a dataset for immature cherries,Group 2 has a dataset for Half-rip cherries consisting of Half-ripe and mature cherries,Group 3 has a dataset for ripe cherries,and Group 4 has a dataset for dry cherries as shown in Table 2.In all the four groups are split into training 60%,Testing 20%,Validation 10%,and Holdout 10%.The purpose of splitting the dataset into four groups is to ensure that the model is not over fitting when training the data,set to select the features,and perform 10 cross folds validation.

    3.3 XGBoost Model Analysis

    An XG Boost is a tree-based model that is used to train and classify coffee berries as either infected or healthy.The algorithm will start by building a simple decision tree.This tree will split the data into four groups,based on the value of one of the features.For illustration,the tree might split the data based on the size of the berry.Most visible as microscopic and dark,sunken with brown patches on the green berry that start small and later enlarge in the area,eating the entire fruit and rendering it black.The algorithm will subsequently evaluate tree performance.If the decision of the tree is inaccurate,the algorithm will build a new tree that fixes the deficiencies.This will be performed until the algorithm converges or a maximum number of trees is achieved.“Sample 1”consists of the immature Cherry 160 dataset split into training 60%,Testing 20%,Validation 10%,and Holdout 10%.Trained using XGBoost parameters,such as maximum depth,learning rate,and the number of estimators as shown in Fig.3.The maximum depth of each tree is a hyperparameter that controls the complexity of the model.A higher depth limit allows the model to learn more complicated patterns in the data but can also result in overfitting.A reduced maximum depth will keep the algorithm from overfitting.

    where max_depth is the maximum depth of each tree in the XGBoost model,n_features are the number of features in the dataset,and d_max is the maximum depth specified by the user as shown in Eq.(2).The maximum depth of a tree is its total amount of levels.A tree with a maximum depth of 1 has one level,and 2 has two levels,etc.A deeper tree can learn more intricate feature-target variable relationships[32,33].

    where learning_rate is the learning rate of the XGBoost model and t is the number of trees in the XGBoost model as shown in Eq.(3).The learning rate regulates model updates per iteration.A faster learning rate will speed up model learning.A number of estimators is a number of trees built in the XGBoost.A higher level of estimators results in high accuracy.The XGBoost technique can train a model to categorize new coffee berries as infected or healthy.Start by setting the maximum depth to 3,learning rate to 0.1,and estimators to 100.This will yield a complicated,accurate model [34].Additionally,evaluated model performance.

    Figure 3:Training and iterations for XGBoost model

    XGBoost model takes as input two arguments (data,target) imported into XGBoost Library then,creates an XGBoost classifier model.It defines the XGBoost model parameters by training the model on the data and target arguments to make predictions on the data argument and returns the predictions.“Sample 2”consists of Half-Ripe Cherry 320 dataset split into training 60%,Testing 20%,Validation 10%,and Holdout 10%.Trained using XGBoost parameters,such as maximum depth,learning rate,and the number of estimators.“Sample 3”consists of Ripe Cherry 112 dataset split into training 60%,Testing 20%,Validation 10%,and Holdout 10%.Trained using XGBoost parameters,such as maximum depth,learning rate,and the number of estimators.“Sample 4”consists of Overripe Cherry 78 dataset split into training 60%,Testing 20%,Validation 10%,and Holdout 10%.Trained using XGBoost parameters,such as maximum depth,learning rate,and the number of estimators[35,36].If the model is inaccurate,increase the maximum depth or estimators.Decrease the learning rate slightly.The algorithm will then evaluate the tree’s performance.If the tree is not accurate enough,the algorithm will construct a new tree that attempts to repair the prior tree’s flaws.This method will be repeated until the algorithm converges or a certain number of trees are attained.The final model will be a composite of all of the trees that were constructed.The model will decide which tree to employ based on the characteristics of the new coffee berry.The tree that is employed will then predict if the berry is diseased or healthy.

    3.4 SHAP Analysis

    SHAP (Shapley Additive explanations) analysis is used for coffee berry disease analysis to gain interpretability and insights into the predictions made by machine learning models providing a quantitative measure of feature importance[37]is maintained by calculating each feature’s contribution to the prediction of a specific instance.if the SHAP analysis shows that the size of the coffee berry is the most important feature,this means that the model is more likely to classify a berry as infected if it is larger than average.SHAP analysis can also be used to identify outliers.

    An outlier is a case in which the model predicts wrong.SHAP analysis can assist in identifying the features that are most different for the outlier,which can aid in understanding why the model made the incorrect prediction and improve the model by removing outliers by focusing on important features as shown in Eq.(4).The SHAP value for a CBD-specific feature and instance is calculated as the difference between the average prediction when the feature is included and the average prediction when the feature is excluded,considering all possible combinations of features.

    whereΦ(p,x)represents the SHAP value of feature j for instance x,m(S)is the average prediction for a subset of features S,m(S ∪{p})is the average prediction for the subset S with the feature j included and v is the total number of features.The formula calculates the difference in predictions between subsets that include the feature p and subsets that do not include the feature p.It sums up the differences over all possible combinations of features,considering all possible subsets S that exclude the feature p.The first line of code includes the SHAP library.The Explainer object is created using the second line.The SHAP values will be computed using this object.The SHAP values are computed in the third line.The SHAP values are returned by the fourth line.The fifth line is responsible for data loading.The target variable is loaded on the sixth line.The Explainer object is created using the model,data,and target in the seventh line.The eighth line employs the Shap.plots.bar()function to plot the SHAP values[38].

    3.5 Virtual Inspection and Diagnosis

    The trained XGBoost model is deployed in a virtual inspection system.Processed SHAP analysis images of coffee cherries from four different samples are captured using Colour based segmentation.Colour-based segmentation is a procedure for dividing the image into different colour subsections.Region-Growing is a technique for restoring together pixels of similar colour.This can be done by thresholding the image to identify pixels that are a certain colour,such as brown or black.Threshold value can be adjusted depending of the dataset.However,starting point is to use a threshold value is greater than or equal to 0.5 The region growing technique can then be used to group these pixels to form regions that represent the infected cherries [39,40].The result can be accomplished by beginning with a cherry pixel and then repeatedly expanding the zone to include all nearby pixels with comparable colours which identifies and predicts whether the cherries are healthy or affected by the disease as illustrated in Fig.4.This technique can be used to detect coffee berry disease in coffee cherries by identifying the characteristic colours of infected cherries.The initial stage in colour-based segmentation is to extract the image’s colour information using a variety of approaches,including RGB,HSV,and YCbCr colour schemes.Region-Growing can be effective in detecting coffee berry disease if the disease symptoms exhibit a cohesive region within the coffee berry.It can handle variations in lighting conditions and image quality and is capable of capturing local dependencies within the region[41].The contribution of these techniques to virtual inspection is that they can be used to develop more accurate and efficient methods for identifying and classifying objects in images.This can be used to improve the quality of inspection and to reduce the cost of inspection.

    3.6 Result Analysis

    The coffee berry disease prediction results using XGBoost-based SHAP and virtual inspection were promising.The model was able to achieve an accuracy of 98.56%,as shown in Table 3,which is a significant improvement over previous methods.The model also identified the most important features for predicting coffee berry disease,which can be used to develop targeted interventions.The SHAP values,which are a measure of the importance of each feature in the model,showed that the most important features for predicting coffee berry disease were the cherry colour,the cherry Shape,and the presence of lesions.These features are all visible through virtual inspection,which means that the model can be used to predict coffee berry disease without the need for physical sampling and identify coffee plants at risk of infection,which can help reduce the spread of the disease and protect coffee production.Section 4 deals with the tenfold cross-validation with the existing algorithm.Section 4 shows the evaluation setup using XGBoost and SHAP by comparing the performance evaluation.

    Figure 4:Colour-based segmentation with region growing technique

    Table 3: Performance evaluation for the XGBoost

    4 Result and Discussion

    Implementing XGBoost and SHAP analysis to compare the proposed model’s performance about,concerning its accuracy and losses,respectively.

    4.1 Evaluation Setup

    Python 3.6.5 simulates the suggested model on PC i5-8600k,GeForce 1050Ti 4 GB,16 GB RAM,250 GB SSD,and 1 TB HDD.The parameter settings are learning rate:0.01,dropout:0.5,batch size:5,epochs: 45.The suggested model was created using Keres Python framework and library version 2.7.The tests used input image sizes from 32×32 by 3 to 256×256 by 3.The proposed strategy appeared to outperform others for a 224×224×3 image size.The suggested model is evaluated using True Positive,True Negative,False Negative,and False Positive metrics using the XGBoost model and SHAP analysis.

    4.2 Performance Evaluation of the XGBoost

    True Positive(TP):the percentage of cases that were both diseased and accurately expected to be diseased as indicated in Eqs.(5)-(9).False Positive(FP):the percentage of cases that were truly healthy but were misdiagnosed as ill.False Negative(FN):the percentage of cases that were genuinely ill but were misdiagnosed as healthy.True Negative(TN):the percentage of cases that were both healthy and accurately anticipated to be healthy.

    As a consequence,Table 3 demonstrates that the proposed method classified the coffee cherry photos into five categories—Healthy,Sample 1,Sample 2,Sample 3,and Sample 4 with an overall accuracy of 80.18%.The suggested method has a precision rate of 83.03%,a recall rate of 87.29%,a specificity rate of 88.62%,and an F1-Score of 89.42%.

    4.3 Performance Evaluation of the SHAP Analysis over XGBoost

    Samples generated from performance evaluation for the XGBoost are processed to achieve accuracy.Additionally,utilizing the XAI tool,the developers employed the SHAP method to thoroughly describe the algorithm.The SHAP calculates the value of characteristics based on their residual effect 60 using Shapley values [42].The SHAP assessment can be performed on all specimens to discover which qualities provide the most value.Furthermore,the SHAP can be utilized as just one component to explain what factors contributed to the model’s ability to predict a particular category and how much they impacted the model’s accuracy as illustrated in Table 4.

    Table 4: Performance evaluation for the SHAP analysis

    The proposed approach properly identified 489 images out of 640,i.e.,85,170,62,44,and 112 for Sample 1,Sample 2,Sample 3,Sample 4,and Healthy samples,respectively.In Sample 1,0 photos were classed as FN(Healthy),whereas only 4 photos in Sample 2 were classified as FP.This is because,since that Sample 1 coffee cherry has features with Healthy images(Sample 3).Furthermore,Sample 3 and Sample 4 had the lowest gaps between coffee cherry photos,with Sample 3 having three FN(Healthy) coffee cherry images and Sample 4 having only one FP (Sample 3) image.5 photographs of healthy coffee cherries round up the Sample 1 classifications.Suggested method exhibits accuracy rates of 85.412%,90.564%,92.646%,92.426%,and 90.862%for Sample 1,Sample 2,Sample 3,Sample 4,and Healthy samples,respectively.

    4.4 Performance Evaluation for 10-Fold Cross-Validation

    Performance evaluation involves partitioning the dataset into ten distinct folds.The model is subsequently trained on nine folds and evaluated on the remaining fold.The aforementioned procedure is iterated in a total of ten cross folds,and subsequently,the outcomes are computed as an average[43].The model would be trained on 90%of the dataset and tested on the remaining 10%.This process would be repeated 10 times,and the results would be averaged from the Fig.5.

    Figure 5:Coffee berry disease prediction

    50 healthy coffee cherry images out of 448,48 of which were Sample 1,were misclassified.10 photographs were mislabeled as dry cherry due to minute variations between Sample 1 and Healthy shots.Similar to this,445,420,and 430 of 500 coffee cherry Samples 2,3,and 4 were correctly identified.Thus,the suggested method separates cherry images into five categories.The grade-wise accuracy ratings are 91.89%,92.23%,94.57%,95.97%,and 97.23%.The suggested model outperforms than existing model with an accuracy of 98.56%in Table 5.

    Table 5: Performance evaluation for the 10-fold cross-validation

    4.5 Comparison with the Existing Methodologies

    Early on the Keras library is employed to conduct a comparative analysis of the performance of various Machine Learning algorithms,namely SVM Linear Kernel,Random Forest,Naive Bayes,DenseNet 121,and ResNet-50.These algorithms are executed,guided,and evaluated using the gathered dataset.The evaluation of the accuracy of several procedures was conducted.The graph presented for SVM Linear Kernel in Fig.6,Random Forest in Fig.7,Naive Bayes in Fig.8,DenseNet 121 in Fig.9 and ResNet-50 in Fig.10.Figures illustrates the relationship between accuracy and the number of epochs for the individual model.The performance of various machine learning models in classifying coffee berries as infected or healthy was evaluated.Specifically,a Support Vector Machine (SVM) with a linear kernel achieved an accuracy of 85.12%.Similarly,a Random Forest model achieved an accuracy of 84.02%,while a Naive Bayes model achieved an accuracy of 81.21%.Additionally,the DenseNet 121 model achieved an accuracy of 85.12%,and the ResNet-50 model achieved an accuracy of 85.12%.These results are depicted in Fig.11.

    Figure 6:SVM linear kernel

    Figure 7:Random forest

    Figure 8:Na?ve bayes

    Figure 9:DenseNet 121

    Figure 10:ResNet-50

    Figure 11:Comparison between other models

    Fig.12 specifies the loss throughout training and testing;as the number of iterations grows,the losses decrease.In the process of model training,it is expected that the training loss will exhibit a decreasing trend.Additionally,in the absence of overfitting,the testing loss should demonstrate a level of similarity to the training loss.In epoch 1 training loss is 1.8 and the testing loss is 1.9,in epoch 2 training loss is 1.7 and testing loss is 1.8,in epoch 3 training loss is 1.6 and testing loss is 1.7,in epoch 4 training loss is 1.5 and testing loss is 1.6,The loss is greater at first and decreases as the number of repeats grows.When 18 iterations were completed,the loss was 0.9%,and when 26 iterations were completed,the loss was less than 0.3%.

    Figure 12:Training and testing loss

    Fig.13 also depicts a full confusion matrix for the classification task performed by the proposed model,the Boosted SHAP Analysis.

    Figure 13:Sample(A),Sample(B),Sample(C),Sample,(D)Confusion matrix for testing set

    Fig.14 shows that the proposed method exceeds its competitors in many parameters.In comparison to other contemporary methodologies,the model put forth achieved an accuracy of 98.56%.

    Figure 14:Performance evaluation

    5 Conclusion

    This research paper presents a unique and quality machine learning-based model for classifying and identifying coffee cherry disease based on the stages of cherry development.Coffee berry disease(CBD) is a devastating disease that can destroy coffee crops and reduce the quality and quantity of coffee production.In the proposed work,cherry datasets are divided into three stages and 6 varieties(healthy,Sample 1,Sample 2,Sample 3,and Sample 4) coffee cherry dataset.XGBoost Model is utilized for image segmentation and classification based on the colour of the cherry,presence of spots on the cherry,size of the spots,shape of the spots,and location of the spots on the cherry.SHAP model is utilized to highlight the important features and extract features from the segmented images.These images are processed using the trained model,which predicts whether the cherries are healthy or affected by the disease.The model is trained by combining XGBoost and SHAP analysis to improve the quality and accuracy of the performance.Virtual inspection will compare the real-time images and processed images to identify cherry is infected or not.The performance of the current approaches and the suggested strategy is evaluated using the three Stages of Coffee Cherry data set,with an accuracy of 98.56%.The proposed method exceeds its competitors.Future research will concentrate on prescribing fertilizers for plant disease control in early diagnosis and prevention of spread.Excessive fertilizer use will damage plants and destroy soil quality.

    Acknowledgement:The authors would like to express their gratitude to the Deanship for Research&Innovation,Ministry of Education in Saudi Arabia for providing financial support for this research study.

    Funding Statement:The current research study has received financial support from the Deanship for Research &Innovation,Ministry of Education in Saudi Arabia,under the Auspices of Project Number:IFP22UQU4281768DSR122.

    Author Contributions:Conceptualization,Y.A.,and R.S.;methodology,S.R.;software,R.S.;validation,Y.A.,S.R.,and R.S.;formal analysis,Y.A.;investigation,S.R.;resources,Y.A.;data curation,R.S.;writing—original draft preparation,S.R.;writing—review and editing,Y.A.;visualization,Y.A.;supervision,R.S.;project administration,S.R.;funding acquisition,Y.A.All authors have read and agreed to the published version of the manuscript.

    Availability of Data and Materials:The dataset utilized and analyzed in our research is publicly accessible to the public cherry coffee beans from different stages in the Zenodo communities(https://zenodo.org/record/4914786),(https://www.nature.com/articles/s41598-019-38775-5/figures/1).The coding system along with additional data are accessible upon adequate request from the initial and coauthor authors.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    日本黄色视频三级网站网址| 国内少妇人妻偷人精品xxx网站 | 成人午夜高清在线视频| av片东京热男人的天堂| 亚洲中文日韩欧美视频| 波多野结衣高清无吗| 一进一出好大好爽视频| 99riav亚洲国产免费| 精品久久蜜臀av无| 色老头精品视频在线观看| 欧美日本亚洲视频在线播放| 亚洲九九香蕉| 精品久久蜜臀av无| 欧美乱码精品一区二区三区| 亚洲乱码一区二区免费版| АⅤ资源中文在线天堂| 最近在线观看免费完整版| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看 | 精品久久久久久久人妻蜜臀av| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区视频了| 婷婷六月久久综合丁香| 最近最新免费中文字幕在线| 日本三级黄在线观看| 欧美日本视频| 最好的美女福利视频网| 中文字幕熟女人妻在线| 久久精品国产99精品国产亚洲性色| 国产激情欧美一区二区| 久久欧美精品欧美久久欧美| 叶爱在线成人免费视频播放| 亚洲欧美精品综合久久99| 无限看片的www在线观看| 99国产极品粉嫩在线观看| 十八禁网站免费在线| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 成人av一区二区三区在线看| 九九热线精品视视频播放| 午夜精品久久久久久毛片777| 久久久久久久久免费视频了| 久久久久久人人人人人| 亚洲性夜色夜夜综合| 久久久精品国产亚洲av高清涩受| 波多野结衣巨乳人妻| 日本免费一区二区三区高清不卡| 国产精品精品国产色婷婷| 精品无人区乱码1区二区| 一区二区三区国产精品乱码| 999久久久国产精品视频| 国产一区二区在线av高清观看| 好男人在线观看高清免费视频| 日本一本二区三区精品| 我的老师免费观看完整版| 国产欧美日韩精品亚洲av| 亚洲色图 男人天堂 中文字幕| 99热这里只有是精品50| 少妇粗大呻吟视频| 国产精品自产拍在线观看55亚洲| 婷婷丁香在线五月| 午夜两性在线视频| www.999成人在线观看| 午夜福利视频1000在线观看| 黄色a级毛片大全视频| 亚洲一区二区三区色噜噜| 欧美黑人欧美精品刺激| 日本一区二区免费在线视频| 麻豆一二三区av精品| 又粗又爽又猛毛片免费看| 久久香蕉国产精品| 又爽又黄无遮挡网站| 亚洲自偷自拍图片 自拍| 少妇人妻一区二区三区视频| 欧美成人一区二区免费高清观看 | 亚洲中文字幕日韩| 亚洲国产欧美人成| 99热6这里只有精品| 亚洲中文日韩欧美视频| 日韩 欧美 亚洲 中文字幕| 两个人视频免费观看高清| 特大巨黑吊av在线直播| 午夜两性在线视频| 午夜福利免费观看在线| 在线看三级毛片| 欧美日韩国产亚洲二区| 黄色女人牲交| 在线观看午夜福利视频| 啦啦啦观看免费观看视频高清| 麻豆久久精品国产亚洲av| 长腿黑丝高跟| 女人被狂操c到高潮| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 欧美乱妇无乱码| 啦啦啦韩国在线观看视频| 国产爱豆传媒在线观看 | 亚洲av五月六月丁香网| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 中文资源天堂在线| 国产真人三级小视频在线观看| 免费人成视频x8x8入口观看| 怎么达到女性高潮| 久久久国产欧美日韩av| 国产成人影院久久av| 久久久久精品国产欧美久久久| 亚洲成人免费电影在线观看| 精品久久久久久久毛片微露脸| 亚洲免费av在线视频| 99国产精品一区二区三区| 色噜噜av男人的天堂激情| 欧美成人一区二区免费高清观看 | 日韩欧美 国产精品| 国产三级黄色录像| 亚洲av电影不卡..在线观看| 国产一级毛片七仙女欲春2| 国产午夜精品久久久久久| 99热6这里只有精品| 欧美乱色亚洲激情| 宅男免费午夜| 中文字幕精品亚洲无线码一区| 欧美日本视频| 少妇被粗大的猛进出69影院| 狂野欧美白嫩少妇大欣赏| 国产午夜精品论理片| 国产一区二区在线观看日韩 | 国产91精品成人一区二区三区| 精品欧美国产一区二区三| 天堂av国产一区二区熟女人妻 | 18禁黄网站禁片午夜丰满| 韩国av一区二区三区四区| 国产在线观看jvid| 毛片女人毛片| 99国产精品一区二区三区| 久久婷婷人人爽人人干人人爱| 给我免费播放毛片高清在线观看| 中文资源天堂在线| 国产av一区二区精品久久| 男人的好看免费观看在线视频 | 黄频高清免费视频| 久久久久精品国产欧美久久久| 999精品在线视频| 久久久国产精品麻豆| 可以免费在线观看a视频的电影网站| 亚洲午夜精品一区,二区,三区| 美女黄网站色视频| 国产视频内射| 777久久人妻少妇嫩草av网站| 亚洲18禁久久av| 日韩有码中文字幕| 在线观看66精品国产| 男女床上黄色一级片免费看| 首页视频小说图片口味搜索| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 婷婷丁香在线五月| 精品国产超薄肉色丝袜足j| 波多野结衣高清作品| 禁无遮挡网站| 精品高清国产在线一区| 99久久精品热视频| 国产精品美女特级片免费视频播放器 | 男女那种视频在线观看| 亚洲自偷自拍图片 自拍| 999久久久精品免费观看国产| www.自偷自拍.com| 五月玫瑰六月丁香| 香蕉国产在线看| 国产人伦9x9x在线观看| 免费无遮挡裸体视频| 伊人久久大香线蕉亚洲五| 亚洲国产欧美人成| 在线观看日韩欧美| 成人18禁高潮啪啪吃奶动态图| 日日夜夜操网爽| 亚洲一区中文字幕在线| 国产99白浆流出| 精品久久久久久,| 国产蜜桃级精品一区二区三区| videosex国产| 麻豆成人午夜福利视频| 欧美一级毛片孕妇| 脱女人内裤的视频| 亚洲一区二区三区色噜噜| 香蕉av资源在线| 一个人观看的视频www高清免费观看 | 久久久久久人人人人人| 巨乳人妻的诱惑在线观看| 极品教师在线免费播放| 操出白浆在线播放| e午夜精品久久久久久久| 免费无遮挡裸体视频| 免费电影在线观看免费观看| 在线看三级毛片| 亚洲专区国产一区二区| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久久久免费视频| 午夜福利18| 男插女下体视频免费在线播放| 91九色精品人成在线观看| 丰满的人妻完整版| 国产v大片淫在线免费观看| 久久精品91蜜桃| 我要搜黄色片| 一级毛片精品| 丰满人妻一区二区三区视频av | 久久人妻福利社区极品人妻图片| 一级毛片女人18水好多| 欧美在线黄色| 亚洲人成伊人成综合网2020| 我的老师免费观看完整版| а√天堂www在线а√下载| 啦啦啦观看免费观看视频高清| 18禁黄网站禁片免费观看直播| 午夜精品在线福利| 变态另类成人亚洲欧美熟女| 18禁美女被吸乳视频| 怎么达到女性高潮| 久久精品aⅴ一区二区三区四区| 国产av又大| 欧美性长视频在线观看| 国产黄a三级三级三级人| 成人亚洲精品av一区二区| 国产高清视频在线播放一区| 午夜免费激情av| 国产伦在线观看视频一区| 99热只有精品国产| 狂野欧美激情性xxxx| 悠悠久久av| 久久久久久久精品吃奶| 国产成+人综合+亚洲专区| 88av欧美| 级片在线观看| 男女下面进入的视频免费午夜| 日本一区二区免费在线视频| 妹子高潮喷水视频| 99久久无色码亚洲精品果冻| 国产成+人综合+亚洲专区| 亚洲国产欧洲综合997久久,| 欧美中文综合在线视频| 搡老岳熟女国产| 在线观看日韩欧美| 999久久久精品免费观看国产| av视频在线观看入口| 久久 成人 亚洲| 99久久精品热视频| 欧美日韩亚洲国产一区二区在线观看| 女人高潮潮喷娇喘18禁视频| 日韩有码中文字幕| 亚洲成人免费电影在线观看| a级毛片在线看网站| 亚洲国产日韩欧美精品在线观看 | 久久这里只有精品中国| 狂野欧美激情性xxxx| 国产精品综合久久久久久久免费| 禁无遮挡网站| 小说图片视频综合网站| 一个人免费在线观看电影 | 亚洲专区中文字幕在线| 欧美乱妇无乱码| 亚洲av五月六月丁香网| 最近视频中文字幕2019在线8| 久久热在线av| 91麻豆精品激情在线观看国产| 又爽又黄无遮挡网站| 久久久久久人人人人人| 真人一进一出gif抽搐免费| 白带黄色成豆腐渣| av视频在线观看入口| 在线国产一区二区在线| 久久久久久久午夜电影| 国产激情欧美一区二区| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 亚洲精品在线观看二区| 国产麻豆成人av免费视频| 一区二区三区激情视频| 在线观看午夜福利视频| 国产精品98久久久久久宅男小说| cao死你这个sao货| 99国产精品一区二区三区| 亚洲av日韩精品久久久久久密| 两个人看的免费小视频| 亚洲片人在线观看| 午夜福利成人在线免费观看| 久久久国产欧美日韩av| 欧美一区二区国产精品久久精品 | 啦啦啦免费观看视频1| 久久这里只有精品19| 亚洲精品中文字幕在线视频| 成人手机av| 国产乱人伦免费视频| 欧美日韩乱码在线| 国产伦一二天堂av在线观看| 精品久久久久久久毛片微露脸| 久久久精品大字幕| 免费电影在线观看免费观看| 午夜久久久久精精品| 一a级毛片在线观看| 欧美黄色淫秽网站| 亚洲真实伦在线观看| 久99久视频精品免费| 精品国产超薄肉色丝袜足j| a级毛片在线看网站| 欧美精品啪啪一区二区三区| 国产精品自产拍在线观看55亚洲| 久久久久久久久中文| 一二三四社区在线视频社区8| www.熟女人妻精品国产| 久久亚洲精品不卡| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 久久精品国产99精品国产亚洲性色| 波多野结衣高清无吗| 国产真实乱freesex| 亚洲中文字幕日韩| 午夜福利视频1000在线观看| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看 | 麻豆一二三区av精品| 日本精品一区二区三区蜜桃| 精品久久久久久久久久免费视频| 国产成年人精品一区二区| 国产高清videossex| av在线播放免费不卡| 亚洲一区高清亚洲精品| 91字幕亚洲| 午夜精品久久久久久毛片777| 欧美三级亚洲精品| 一本一本综合久久| 亚洲精品中文字幕在线视频| 国产精华一区二区三区| 亚洲精品在线美女| 久久九九热精品免费| 精品不卡国产一区二区三区| 日韩大码丰满熟妇| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 黄片大片在线免费观看| 精品国产超薄肉色丝袜足j| 日韩高清综合在线| 亚洲熟女毛片儿| 一边摸一边抽搐一进一小说| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 成年版毛片免费区| 1024视频免费在线观看| 午夜两性在线视频| 欧美成狂野欧美在线观看| 日本黄色视频三级网站网址| 真人一进一出gif抽搐免费| 亚洲精品在线美女| 亚洲一卡2卡3卡4卡5卡精品中文| 老熟妇仑乱视频hdxx| 国产真实乱freesex| 亚洲成人久久性| 麻豆国产av国片精品| 精品国产超薄肉色丝袜足j| 男人的好看免费观看在线视频 | 在线看三级毛片| 日韩欧美国产一区二区入口| 午夜激情福利司机影院| 免费看a级黄色片| av片东京热男人的天堂| 97人妻精品一区二区三区麻豆| www国产在线视频色| 美女大奶头视频| 久久草成人影院| 国产在线观看jvid| 国产成人啪精品午夜网站| 亚洲片人在线观看| 午夜久久久久精精品| 久久天堂一区二区三区四区| 精品一区二区三区av网在线观看| 亚洲精品粉嫩美女一区| 99久久99久久久精品蜜桃| 亚洲性夜色夜夜综合| 国产熟女xx| 在线观看美女被高潮喷水网站 | 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 我要搜黄色片| 老熟妇乱子伦视频在线观看| 亚洲精华国产精华精| 黄色成人免费大全| 青草久久国产| 欧美一区二区精品小视频在线| 国产成人aa在线观看| 97碰自拍视频| 又爽又黄无遮挡网站| 精华霜和精华液先用哪个| 精品欧美一区二区三区在线| 久久婷婷人人爽人人干人人爱| 久久久久精品国产欧美久久久| 久久天堂一区二区三区四区| 精品国产乱子伦一区二区三区| 亚洲熟妇熟女久久| 亚洲av成人一区二区三| 国产黄色小视频在线观看| 久久久久国产精品人妻aⅴ院| 啪啪无遮挡十八禁网站| 一区二区三区激情视频| 操出白浆在线播放| 午夜影院日韩av| 国产v大片淫在线免费观看| 男人舔女人下体高潮全视频| 亚洲欧洲精品一区二区精品久久久| 亚洲,欧美精品.| 精品国产乱子伦一区二区三区| 国产午夜精品久久久久久| 看黄色毛片网站| 精品国产乱码久久久久久男人| 日韩国内少妇激情av| 免费搜索国产男女视频| 精品国产亚洲在线| 亚洲人与动物交配视频| 美女大奶头视频| 日韩欧美在线二视频| 亚洲,欧美精品.| 欧美最黄视频在线播放免费| a级毛片a级免费在线| 国产成+人综合+亚洲专区| 欧美成人性av电影在线观看| 国产视频一区二区在线看| 香蕉国产在线看| 久久精品人妻少妇| 18禁裸乳无遮挡免费网站照片| 波多野结衣巨乳人妻| 叶爱在线成人免费视频播放| 白带黄色成豆腐渣| 两性夫妻黄色片| 天天躁夜夜躁狠狠躁躁| 欧美黄色淫秽网站| 久久精品国产亚洲av香蕉五月| xxx96com| videosex国产| 国内精品一区二区在线观看| 国产精品久久视频播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲成a人片在线一区二区| 国产真人三级小视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 深夜精品福利| 午夜日韩欧美国产| 日日干狠狠操夜夜爽| 麻豆成人av在线观看| 真人一进一出gif抽搐免费| www.自偷自拍.com| 午夜福利成人在线免费观看| 亚洲男人天堂网一区| 看片在线看免费视频| 我要搜黄色片| 国产精品亚洲一级av第二区| 国产精品美女特级片免费视频播放器 | 国产精品香港三级国产av潘金莲| 黑人欧美特级aaaaaa片| 精品久久久久久久毛片微露脸| 国产精品乱码一区二三区的特点| 亚洲成人精品中文字幕电影| 欧美国产日韩亚洲一区| 变态另类丝袜制服| 一本大道久久a久久精品| 欧美黑人巨大hd| 国产99白浆流出| 精品无人区乱码1区二区| 欧美在线黄色| 久久久久久久久免费视频了| 成在线人永久免费视频| 黄色成人免费大全| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播| 成人亚洲精品av一区二区| 日韩欧美三级三区| 欧美高清成人免费视频www| 国产成人精品久久二区二区91| 亚洲人成网站在线播放欧美日韩| 18禁黄网站禁片午夜丰满| 久久国产乱子伦精品免费另类| 久久伊人香网站| 天堂动漫精品| 国产午夜精品论理片| 中亚洲国语对白在线视频| 悠悠久久av| 又爽又黄无遮挡网站| 国产欧美日韩一区二区三| 精品国产乱子伦一区二区三区| or卡值多少钱| av免费在线观看网站| 成人av在线播放网站| 亚洲美女黄片视频| 亚洲天堂国产精品一区在线| 久久久国产精品麻豆| 日韩精品中文字幕看吧| 国产亚洲av嫩草精品影院| 99国产极品粉嫩在线观看| 国产成人啪精品午夜网站| 长腿黑丝高跟| 天天一区二区日本电影三级| 午夜精品在线福利| 亚洲国产精品999在线| 中文在线观看免费www的网站 | 成年版毛片免费区| 成人av在线播放网站| av国产免费在线观看| 国产在线精品亚洲第一网站| 亚洲狠狠婷婷综合久久图片| 色综合婷婷激情| av在线天堂中文字幕| 又黄又粗又硬又大视频| 18禁美女被吸乳视频| 一二三四在线观看免费中文在| 欧美乱妇无乱码| 久久久久久免费高清国产稀缺| 久9热在线精品视频| 久久久久久免费高清国产稀缺| 好男人在线观看高清免费视频| 国产一区二区三区视频了| 精品熟女少妇八av免费久了| 日韩 欧美 亚洲 中文字幕| 啦啦啦观看免费观看视频高清| 日本免费a在线| 午夜两性在线视频| 国产69精品久久久久777片 | 欧美色视频一区免费| 精品久久久久久,| 国产精品一区二区免费欧美| 免费观看人在逋| 亚洲精品在线美女| 成年女人毛片免费观看观看9| 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久| 毛片女人毛片| 97人妻精品一区二区三区麻豆| 淫妇啪啪啪对白视频| 性色av乱码一区二区三区2| 国产一级毛片七仙女欲春2| 国产精品亚洲av一区麻豆| 免费在线观看完整版高清| 国产伦在线观看视频一区| 日本一二三区视频观看| 国产成人系列免费观看| 看免费av毛片| 亚洲黑人精品在线| 嫩草影视91久久| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 全区人妻精品视频| 亚洲一区高清亚洲精品| 亚洲国产欧洲综合997久久,| 国内精品一区二区在线观看| 亚洲成人久久爱视频| 国产主播在线观看一区二区| 中文字幕av在线有码专区| av超薄肉色丝袜交足视频| 国产三级中文精品| 亚洲欧洲精品一区二区精品久久久| 亚洲性夜色夜夜综合| 久久精品91蜜桃| 中出人妻视频一区二区| 变态另类丝袜制服| 777久久人妻少妇嫩草av网站| 午夜激情福利司机影院| 免费电影在线观看免费观看| 亚洲自偷自拍图片 自拍| 99在线视频只有这里精品首页| 男女下面进入的视频免费午夜| 99国产精品99久久久久| 日韩大码丰满熟妇| 欧美av亚洲av综合av国产av| 99精品在免费线老司机午夜| a级毛片在线看网站| 精品人妻1区二区| 九色国产91popny在线| 亚洲欧美一区二区三区黑人| 国产精品综合久久久久久久免费| 日本三级黄在线观看| 欧美性猛交╳xxx乱大交人| 国产99白浆流出| 欧美极品一区二区三区四区| 亚洲av日韩精品久久久久久密| 国产视频一区二区在线看| 久久精品国产99精品国产亚洲性色| АⅤ资源中文在线天堂| 国产精品1区2区在线观看.| 日本撒尿小便嘘嘘汇集6| 久久久久久久久久黄片| 亚洲一区中文字幕在线| 非洲黑人性xxxx精品又粗又长| 丰满人妻一区二区三区视频av | 美女扒开内裤让男人捅视频| 九九热线精品视视频播放| 熟女电影av网| 免费看日本二区| 国产精品九九99| 精品少妇一区二区三区视频日本电影| 首页视频小说图片口味搜索| 国产午夜精品久久久久久| 女同久久另类99精品国产91| 老司机午夜福利在线观看视频| 久久精品综合一区二区三区| 国产精品久久久久久久电影 | 夜夜夜夜夜久久久久| 免费高清视频大片| 一二三四在线观看免费中文在| 女人被狂操c到高潮| 亚洲精品国产精品久久久不卡| 国产亚洲av高清不卡| 91麻豆av在线| a级毛片a级免费在线| 日本黄大片高清| 在线播放国产精品三级| 黄频高清免费视频| 亚洲全国av大片|