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

    Automated Grading of Breast Cancer Histopathology Images Using Multilayered Autoencoder

    2022-08-24 03:29:48ShakraMehakUsmanAshrafRabiaZafarAhmedAlghamdiAhmedAlfakeehFawazAlasseryHabibHamamandMuhammadShafiq
    Computers Materials&Continua 2022年5期

    Shakra Mehak,M.Usman Ashraf,Rabia Zafar,Ahmed M.Alghamdi,Ahmed S.Alfakeeh,Fawaz Alassery,Habib Hamam and Muhammad Shafiq

    1Knowledge Unit of System and Technology,University of Management&Technology,Sialkot Campus,51310,Pakistan

    2Department of Computer Science,University of Management and Technology,Sialkot Campus,51310,Pakistan

    3Faculty of Computing and Technology,University of Engineering&Technology,Narowal,51610,Pakistan

    4College of Computer Science and Engineering,University of Jeddah,21493,Saudi Arabia

    5Faculty of Computing and Information Technology,King Abdulaziz University,Jeddah,21589,Saudi Arabia

    6Department of Computer Engineering,College of Computers and Information Technology,Taif University,Taif,Saudi Arabia

    7Faculty of Engineering,Moncton University,NB,E1A3E9,Canada

    8Department of Information and Communication Engineering,Yeungnam University,Gyeongsan,38541,Korea

    Abstract: Breast cancer(BC)is the most widely recognized cancer in women worldwide.By 2018,627,000 women had died of breast cancer(World Health Organization Report 2018).To diagnose BC, the evaluation of tumours is achieved by analysis of histological specimens.At present, the Nottingham Bloom Richardson framework is the least expensive approach used to grade BC aggressiveness.Pathologists contemplate three elements,1.mitotic count,2.gland formation, and 3.nuclear atypia, which is a laborious process that witness’s variations in expert’s opinions.Recently,some algorithms have been proposed for the detection of mitotic cells,but nuclear atypia in breast cancer histopathology has not received much consideration.Nuclear atypia analysis is performed not only to grade BC but also to provide critical information in the discrimination of normal breast,non-invasive breast(usual ductal hyperplasia,atypical ductal hyperplasia)and pre-invasive breast(ductal carcinoma in situ) and invasive breast lesions.We proposed a deep-stacked multi-layer autoencoder ensemble with a softmax layer for the feature extraction and classification process.The classification results show the value of the multilayer autoencoder model in the evaluation of nuclear polymorphisms.The proposed method has indicated promising results, making them more fit in breast cancer grading.

    Keywords: Breast cancer;nuclear atypia;autoencoder

    1 Introduction

    Breast cancer(BC)alludes to a hostile tumour that has been created from nuclei in the breast.BC can originate in three areas:(1)lobules,which are the milk drain organs;(2)ducts,which are paths that empty the drain out of the lobules to the nipples; or (3) stromal tissues, which incorporate the fatty stringy connective tissues of the breast.Among them,the first two are the most common classes of BC[1].Diagnosis with higher precision is desirable to provide treatment efficiently.Regardless of extensive development in the diagnosis and management of breast cancer,it is still the second leading cause of high mortality[2].However,it is necessary to diagnose cancer in the early phases and exact regions to avoid severity.Pathologists use different medical imaging modalities,such as mammography studies,ultrasound scanning, computerized scanning and magnetic resonance imaging, to obtain images of the breast for later analysis[3].The analysis was performed to obtain the aggressiveness of cancerous cells.The aggressiveness of cells is categorically explained by using a grading system that has impacts on treatment.Cancer grading determines how the tumour nuclei look under microscopic tests[4,5].It is not the same as measuring the tumour size and its aggression and metastasis.

    The Nottingham grading system(NGS)proposed by Elston and Ellis is a universal grading system prescribed by the WHO to evaluate BC aggressiveness by pathologists worldwide [6].It compares the appearance of BC tissues with the appearance of normal breast tissues, as shown in Tab.1.It consists of three grades(I),(II),and(III)achieved from the expansion of these three measures:mitotic count,gland formation and nuclear atypia.These three criteria are valued as 1,2,or 3(Tab.1)[6].The evaluation of the BC grade is semi-quantitative.The computation of the three parameters provides tumour aggressiveness.The minimum conceivable score is 3(1+1+1),and the maximum conceivable score is 9(3+3+3).Patients with a score of 8 or 9 are assigned Grade-III.Grade-II alludes to scores of 6 or 7,while Grade-I alludes to scores of 3,4 or 5.High-grade nuclear atypia shows a high deviation in cell shape.Low-grade nuclear atypia,for the most part,has round nuclei[5].

    Table 1:Nottingham grading system by Elston and Ellis

    The rest of this paper is structured as follows.Section 2 outline the latest research on automated image analysis in histopathology and contribution of the paper.Section 3 describes the proposed methodology.Section 4 presents the results and discussion.In the last section, we have drawn our conclusion.

    2 Automated Image Analysis in Histopathology

    The microscopic study of cell morphology plays a vital role for medical experts in decision making[7].Automated nuclei segmentation and classification is a repetitive activity and challenging for pathological images.Due to rapid development in digital pathology,various cancer diagnosis and grading systems have been proposed,including brain image analysis[8–10],cervix[11],lungs[12],liver[13]and breast[5,14–18].In this regard,we can find more systems in[19–22].

    BC image analysis can be applied to handle numerous pathology jobs,such as mitosis detection[18,23],tissue segmentation[24],histological classification or cancer grading[15].Automatic analysis is commonly performed by using hematoxylin and eosin (H&E)-stained slides and deep learning techniques to enhance model performance[25].Breast cancer diagnosis from histopathological images always remains the benchmark in clinical pathology [20–24].Histopathology is the investigation of microscopic variations or abnormalities in tissues.Distinctive types of nuclei in BC histological images call attention to the difficulties in their recognition [26].The images that are obtained through traditional radiology,such as mammograms,ultrasound and MRI,are greyscale and do not cover interesting regions, whereas histopathology images have better quality in diagnosing various diseases, including an extensive variety of cancers, due to distinctive features achieved through the staining process in biopsy tests[27].Cancer detection and grading are performed on histopathological images that are obtained by different high-resolution scanners,such as Aperio XT and Hamamatsu NanoZoomer.Cancer cell identification and grading are achieved at different magnifications, such as 10×, 20× and 40×, because of the tissue structure morphology [28].The obtained whole slide images are exposed to pathologists,who examined the images manually to segment regions of interest for assessing critical information for classification[29].The proper classification of breast lesions can provide critically important diagnostic information to prevent under-and over-treatment and properly guide patient treatment.

    The major problem with the current manual diagnosis system is inter-and intra-observer variability and reproducibility problems.Computer-aided diagnosis(CAD)is a developing interdisciplinary area that consolidatescomponents of machine learning and digital image processing with medicinal information[21].Novel CAD techniques for analysing breast tissue will enable pathologists to have more accurate and reproducible diagnoses, leading to improvements in the ability of pathologists to diagnose and prognosis breast cancer.Several studies exist about breast cancer detection and classification using deep learning.Sadad et al.[4]proposed the methodology to classify pre-segmented breast lumps as cancerous and non-cancerous in X-ray images.They used a convolutional neural network and pre-trained models, combined data augmentation and transfer learning to address overfitting problems.The authors only used limited pre-trained networks on a digital database for screening mammography (DDSM), which might be inadequate to generalize the results of their study.Huynh et al., [19] presented a computer-aided system that extracted 219 breast lesions from breast mammographic images.They used CNN as a feature extractor and support vector machine for classification purposes.Although classification accuracy improved, due to the small number of training data,it might experience overfitting[16].Khan et al.,[21]proposed a framework intended to score cancerous tumour tissues based on ER and PR hormone receptors by using image processing techniques.This system scores ER/PR pervasiveness more accurately than conventional methods that usually assist pathologists in exploring histopathological images.The results are recorded in a database to facilitate further analysis.Doyle et al., [30] segregated poor quality from high-grade BC tumour slides by applying a mix of surface-based components for recognizing tumours from non-cancer,and afterward, extended design elements to recognize low from high evaluations.However, this method depends on the precision of the position of nuclei.In their dataset, the nuclei were commented on physically.However, Abirami et al., [31] proposed an artificial neural network-based classifier to detect micro-classifications in every area of breast cancer mammograms.They used discrete wavelet transform(DWT)for feature extraction and neural networks such as multi-layer perceptron and radial basis function for the classification of BC mammographic images.

    Evaluation of the histopathological images of nuclear polymorphism BC is a qualitative method of treatment for assessment and prognosis.In 2020, a study is conducted for mitotic cell detection by using deep learning techniques instead of image processing techniques, as those techniques are costly and not up to the mark.For detection of mitotic cell Faster-RCNN has been used, but this generates results that has unacceptable rate of false positives.These false positives are generated because of small size of mitotic cell or it is similar looking to normal cell.To reduce false positives post processing is applied.In post processing two deep learning method RESNET 50 and DENSENET-201 has been applied that results in refined set of mitotic cell candidates, which are further refined by using statistical parameters such as shape,texture etc.This technique is tested on two open source databases ICPR 2012,ICPR 2014 and for generalization testing on TUPAC16.On all these datasets the proposed methodology performs well[32].In another study,authors proposed two deep learning models for semantic segmentation and classification to increase robustness.The findings indicate that by combining traditional image processing methodologies and deep learning approaches with hand-crafted features,the efficiency of traditional image processing methodologies and deep learning approaches can be significantly improved[33].

    Different from the existing techniques,the proposed method in our study selects only critical cell nuclei features and utilizes them to grade NA into I,II and III accordingly.The results comprise the prediction of nuclear atypia,which will not merely help pathologists in diagnosis but also bring robust and productive changes in the field of medical imaging.

    To identify histological grades in breast cancer images, we developed an automated grading of nuclei using unsupervised feature extraction.Our major contributions in this paper are as follows:

    ·The proposed system combines a patch-based automated feature extraction method with controlled decision fusion for feature extraction and classification.

    ·An autoencoder-based model is proposed that automatically isolates discriminative patches from high-resolution images for data training.

    ·Our model attains several state-of-the-art outcomes classifying whole slide images(WSIs)into cancer grades and is closely related to inter-observer agreement among pathologists.

    ·Experimental evidence shows that training deep learning models on high-resolution WSI blocks can produce promising results,making them more suitable for breast cancer grading.

    3 Materials and Methods

    3.1 Dataset Description

    In this research, a publicly available nuclear atypia dataset containing Grade-I, II and III BC anomalies by the ICPR (International Conference on Pattern Recognition), is used to assess the effectiveness of the proposed model.The dataset was generated by using 300 manually tagged whole 50 high power field(HPF)slide images of size 1539×1376.These images were scanned by an Aperio scanner at×20 magnification.A solitary WSI of 20×magnification tissue sample comprises almost 1010 multi-channel, representing a highly complex multi-scale structure [34].An HPF is a territory of 0.262 mm2,which is a surface comparable to that of a magnifying lens field breadth of 0.58 mm.The Aperio scanner has a resolution of 0.2273 μm horizontal and 0.222753 μm vertical per pixel.For each frame,two files are maintained.One contains a list of marked scores by each pathologist,and the other is the final score of nuclear atypia decided on majority opinions.The score for nuclear atypia requires a wide zone to have the ability to assess the shape and size of an expansive populace of nuclei.For this purpose,the pathologists worked at×20 magnification,and the final score was given by two or three pathologists.The dataset has been split into training and validation sets.The training and validation sets contain 70 percent and 30 percent of total dataset respectively.Each set contains WSIs for three grades,19%of training dataset is grade-1%,40%is grade-II and remaining is grade-III.In validation set,grade-I,II and III are 16%,42%and 42%respectively.The data used in our experiments also accessible through[35].

    3.2 Data Augmentation

    For deep learning models,the fundamental object is the dataset.Deep learning algorithms truly perform well when the dataset is large[36].This is mainly accurate in the medical imaging field,where access to data is highly protected due to patient secrecy concerns [37].To overcome this matter, we used data augmentation to increase the dataset.The main purpose of data augmentation is generating genuine-looking training data by applying a transformation to an example without altering its label[38–40].Some non-categorical image transformation techniques are used to produce new examples from the training data for image classification,such as cropping,flipping,shifting,and colour jittering,adding Gaussian noise and rotating images at different angles.To balance our dataset,we applied data augmentation of the WSI at different angles(45°,90°,120°and 270°).Algorithms 1 was used for data augmentation.

    3.3 Patch Extraction for Whole Slide Nuclei Image Classification

    Patch extraction is normally used to represent larger scan.Deep neural networks,though,operate on smaller image dimensions as stated in[41]that input size of image may not be larger than 350 by 350 pixels.To reduce the significant image size,input space,computational time and system specifications,WSIs of size 1539×1376 were divided into 128 by 128 size patches,20 manually annotated patches by expert,were used in our study as patch with higher grade will represent WSI[42].As far as number of patches concerns,it is observed that deep neural networks with larger input size may need large number of training parameters i.e.,neurons.This may cause difficulty for neural network during training[41].Every patch consisted of associated features that were representative of the WSI.All patches were assigned a unique identifier and grade according to the Nottingham grade system by a pathologist.Our patch extraction approach is elaborated in Algorithm-2.The overall architecture of the proposed methodology is illustrated in Fig.1.

    Algorithm 1:Data Augmentation Input–Whole slide image(WSI)Output–Augmented whole slide images Step 1:Load WSI I.Step 2:Create image data augmentation generator A (f (i, j), degrees) by a given number of degrees from 0 to 360.Step 3:Generate augmented images.Step 4:Stop.

    Algorithm 2:Generation of patches Input–Whole slide image(WSI)Output–Patches of size n×n of the input image Step 1:Let I be the epidemic patch count and N be the total number of patches.Step 2:Generate patches Step 2.1:Convolve window on the image.Step 2.2:If patch size decreases from n×n,go to step 4.Step 2.3:Extract patch.Step 3:if I=N go to step 4,otherwise go to step 2.Step 4:Stop.

    Figure 1:Framework for nuclei grading based on the Nottingham grading system

    3.4 Feature Learning with Two-Layered Autoencoders

    The architecture of the autoencoder mainly consists of two blocks of layers[34].The first layer is the input layer,which resides as an encoder that is intended to perform the encoding of input data into a lower-dimensional latent compression;a hidden layer can be seen as a new feature representation of input data.The other block contains an output layer that is trained to reconstruct the input from the hidden representation[43,44].The autoencoder simultaneously reduces the divergence between feed data and its reconstruction.The output feature set,which is obtained after applying the transformation function,of every connected layer is the input of successive layers[45].Conventionally,the autoencoder takes input and finds weights that are evaluated on the basis of a threshold.If that threshold matches,data training is stopped; otherwise, weights are again calculated to meet that required criterion.The working architecture of the autoencoder is shown in Fig.2.The autoencoder finds the optimal parameters P= (W,B,B′)by using a hidden layer that executes a transformation function,Eq.(1).

    X= {x1,x2,x3...xn} patches of dimension 128×128 are mapped to dx vectors.The encoder maps the dx vector to an n-dimensional feature vector, i.e., h = ft(x) ∈Vdh.Each neuron in the encoder has a vector of weights that will react to a particular nuclear structure.Later, the features of image patches are passed to the decoder.The decoder tries to change the map to recreate the original input for validation.The high-level feature vector of each patch x, i.e., h ∈V represent cell nuclei atypia.The n feature vector of n patches (X) along with n target values signified asy (m) ∈ {I, II, III},i.e.,{h (m),are then given as input to the softmax layer,which reduces the cost function and acts as a classifier.The list of notations used in modelling is given in Tab.2.We used an autoencoder combined with a softmax layer to shape a deep network.The softmax function can be traditionally used to characterize a categorical distribution that works on the basis of probability function Eq.(2)over different likely results.

    Figure 2:Feature learning process from BC histopathology

    Table 2:Summary of the notations

    In proposed framework,the autoencoder has three hidden layers,the output of the hidden layer is given as input to the next layer.The first hidden layer generates 3000 features which were given as input to the second layer that filter out 300 features from given feature set.Later on,these 300 features maps into 50 features by third layer.The final 50 features are used in classification as shown in Fig.2.

    3.5 Classification

    3.5.1 Support Vector Machine

    The support vector machine(SVM)is a supervised machine learning classifier that combines the influence of conventional theoretical statistical approaches[46].SVMs are commonly used classifiers in various machine learning-based healthcare areas,such as medical imaging[47]and bioinformatics[48].SVM incorporates the capability to overcome the issue of overfitting by using the idea of error minimization.The SVM classifier looks for the ideal hyperplane,which is in the middle of the classes[49].This ideal unscrambling hyperactive plane has numerous refined measurable qualities.SVM additionally learns a non-linear decision function that captures extreme edge hyperplanes by a kernel in high-dimensional vector space that is mapped from input space [50].Typically, SVM addresses noisy data problems and is usually used for binary classification, but it can be specifically stretched out for multi-class classification.In this study,SVM was used for paired characterization of aggressive tumours of Grades I,II and III,as shown in Fig.3.For classification,we designed three parameters for each grade as:xj=(0 0 1),i.e.,Grade-I,b)yj=(0 1 0),i.e.,Grade-II,c)zj=(1 0 0),i.e.,Grade-III.

    Figure 3:Multi-class support vector machine using three hyperplanes,classifying Grades-I,II and III

    3.5.2 Decision Tree

    Numerous machine learning(ML)algorithms are available,among which j48 and random forest are the most frequently used decision trees for classification [51].A decision tree is a method that classifies data as an n-array hierarchy.The initial step,known as tree building,generates a decision tree model through training data and a function, i.e., entropy for obtaining the latest knowledge trends,given in Eq.(3).The other step is tree pruning,in which the insignificant branches do not contribute.To classify,from the root hub,the property estimation of the entered records is tested until a leaf hub is achieved for its respective class.The purpose of tree pruning is to decrease the background distortion or noise in the training data [10,52,53].In this study, random forest (Fig.4) and J48 decision tree(Fig.5)classifiers are used.The random forest(RF)algorithm classifies a given problem by executing the combination of a few decision trees where each tree depends on the value of a distinctly sampled random vector.The distinction of such ensembles of classifiers is that their tree-based components are developed from a specific measure of irregularity,and random forest is characterized as a non-specific standard of the randomized gathering of decision trees.In addition to randomization, achieved by developing trees using bootstrapping, a second layer of randomization is presented at the hub level when growing the tree.RF predicts the output by mapping test features to rules that are generated by using information gain,i.e.,Eq.(4).Finally,it considers the highly voted predicted target as the final prediction from random forest.

    wherepnis the probability that i occurs in the system.A is the target class,and X is the attribute.

    Figure 4:RF architecture that generates N trees to classify nuclei classification

    The J48 decision tree for classification purposes was also used.J48 uses entropy to calculate the homogeneity of the sample data.If the entropy is zero, then the sample data are completely homogeneous, and if the sample data are divided unequally, then it has entropy one.The supplementary features of J48 are detecting missing values,finding persistent characteristic value ranges,determining principles and decision tree pruning.This algorithm delivers the principles for the prediction of the objective variable.In different algorithms, the classification is performed recursively until each sheet is pure; that is, the classification of the data must be as precise as might be expected in the circumstances[54].

    Figure 5:The J48 decision tree shows that the features with the highest information gain classify them into their perspective classes

    3.5.3 Multi-Layer Perceptron Network

    Multi-layer perceptron (MLP) is a supervised neural network consisting of many processing parameters,i.e.,weights that are linked with each other.MLPs are neural systems with an input layer comprised of entities,at least one computation hidden layer that uses the default activation function and an output layer comprising computation nodes.An optimization class is used to minimize the given loss function [55].All attributes of the system are standardized before being given as input.In the input layer,the feature vector is given to the source node that generates the input signal.The input signal moves onward on a path via layers and produces an output signal.As a final point,the computational process terminates on the output layer,which produces results as shown in Fig.6.

    Figure 6:Multi-layer perceptron with N hidden Neurons

    3.5.4 Performance Metrics for Classification

    The quantitative performance of our framework with different models is shown in Tabs.3–7, where we computed our results using metrics in Eqs.(5)–(8).The performance of automatic grading on the basis of nuclei is quantified in terms of the F-measure, precision, and recall(sensitivity)[56,57].

    X = correctly identified instances, Y = incorrectly identified instances, N = number of identified instances that do not align with manual labels.

    Table 3:Results for SVM

    Table 4:Results for RF

    Table 5:Results for J48

    Table 6:Results for MLP

    Table 7:Results for softmax

    4 Results and Discussion

    The proposed framework is evaluated on a publicly available nuclei dataset.This dataset consists of 30 HPF whole slide images at×20 magnification.The study is implemented on a desktop computer(16 GB RAM with Intel Core(TM)3.2 GHz processor)and Intel graphic processing unit.For software implementation,we performed experiments using MATLAB 2016a.

    Deep learning models require a large dataset for effective learning, we augmented the data and extracted them into 128×128 equal size patches to reduce the computational complexity.Our input patch contains significant A101_01= 128×128 =16,384 pixels characterized as learning features.Each patch contains malignant and benign areas of different grades.For feature learning,we used an unsupervised deep learning-based autoencoder model that contains hidden layers along with encoding and decoding layers.The features extracted from a single patch are pixel vectors.This vector is s a portion vector of pixel forces with degree(128×128×3)×1.Thus,there are Sx=S0=128×128×3 input units in the input layer.The first,second and third hidden layers have h(i)=3000,j(i)=300 and K(i)=50,where i=1,2,3,...,n hidden units,respectively.Finally,the softmax layer classifies these 50 neurons into 3 grades accordingly.

    Along with automatic classification with autoencoders,we used various other classifiers because softmax layer classification does not show promising results for this problem.The overall accuracy that was achieved using the softmax function was 63%.The extracted features are given as input to four different machine learning classifiers, i.e.; SVM, RF, J48 and MLP.For SVM, a normalized feature vector is given as input to the RBF kernel with cost 0.1 and precision threshold 0.0005.Second,we use RF for classification with attribute selection and noise removal performed by information gain and pruning,respectively.Along with RF,we also classify grades by considering another decision tree algorithm known as J48.Finally,we use a multi-layer neural network that has 50 hidden neurons and sigmoids as the activation function.The results are compiled at patch-level,however the final label of WSIs can be predicated by considering grade that is assigned to maximum patches of respective WSI.All four classifiers show better results than softmax,as shown in Tabs.3–7.

    Figure 7:Performance of proposed classifies:(a)recall;(b)F-measure,(c)precision;(d)accuracy

    In the classification problem, the sizes of the training and testing data are constrained.We used the N-fold cross-validation technique, where N=10.Various N-fold cross-validations attempt distinctive things with a comparable learning design,and the data index afterwards conveys differing results because of the effect of sporadic assortment in selecting the folds.Stratification reduces the assortment, yet it most likely does not abstain from it totally.We used the N-fold cross-validation method for estimating a tuning parameter t(subset size).This strategy reserves a certain quantity of data for training and testing.N-fold cross-validation is used to predict the error rate of feature learning techniques.We divided our data set into N equal parts,for each N=1,2,3,...,10 fit the model with t parameter to the other N-1 parts,givingα-N(t)and computed its error in predicting the Nth part.Since it was a hectic and tedious job to have a specialist or pathologist that physically recognized every single nucleus in each of these patches(to give a ground truth to quantitative assessment),we compared the results of classifiers with ground truth details provided along with the dataset.The metrics used to evaluate the nuclei scoring of each method include sensitivity,precision,F-measure and accuracy.The performance measures of all classifiers are shown in Fig.7 including(a)recall;(b)F-measure,(c)precision;(d)accuracy.

    5 Conclusions

    The core aim of this study was to explore challenges in the automated nuclei grading framework for H&E images.The qualitative analysis of BC histopathology images was achieved by performing patch extraction,feature learning,and classification of nuclei.The outcome of the proposed methodology was comparable with pathologists’observations given along the dataset.This technique can conceivably help build a computerized grading system for BC.A significant number of studies have been conducted in the field of computer-aided diagnostic tools,emphasizing nuclei segmentation and classification in various image modalities.Nevertheless,there are some open areas to delve into,which may have unique challenges that should be covered in future research.This study was designed to focus specifically on feature engineering using multi-layered autoencoder.Later on in part b of this study,we will test our extracted features using more advance classifiers.Furthermore,we intended to compute the BC grade by combining mitosis and nuclei regions.We also plan to conduct an investigation of the relationship between different grading parameters to evaluate the morphology of the breast in BC.

    Acknowledgement:We deeply acknowledge Taif University for supporting this study through Taif University Researchers Supporting Project Number(TURSP-2020/150),Taif University,Taif,Saudi Arabia.

    Funding Statement:This work was supported by Taif University(in Taif,Saudi Arabia)through the Researchers Supporting Project Number(TURSP-2020/150).

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

    精品高清国产在线一区| 丰满的人妻完整版| 亚洲av熟女| 很黄的视频免费| 淫妇啪啪啪对白视频| 亚洲色图av天堂| 免费看a级黄色片| 美女高潮到喷水免费观看| www.www免费av| 91麻豆av在线| 一本综合久久免费| 午夜影院日韩av| 91成年电影在线观看| 一级作爱视频免费观看| 国产精品电影一区二区三区| 国产99白浆流出| 午夜免费观看网址| 日韩大码丰满熟妇| 国产在线精品亚洲第一网站| 亚洲三区欧美一区| 在线观看一区二区三区| 日本 欧美在线| 色综合欧美亚洲国产小说| 午夜日韩欧美国产| 无遮挡黄片免费观看| 午夜福利免费观看在线| 极品人妻少妇av视频| 国产在线精品亚洲第一网站| 女人爽到高潮嗷嗷叫在线视频| 成人欧美大片| 日本一区二区免费在线视频| 国产人伦9x9x在线观看| 久久久久国产精品人妻aⅴ院| 午夜精品国产一区二区电影| 日韩 欧美 亚洲 中文字幕| 老汉色av国产亚洲站长工具| 啦啦啦韩国在线观看视频| 亚洲最大成人中文| ponron亚洲| 国产高清激情床上av| 色综合站精品国产| 亚洲aⅴ乱码一区二区在线播放 | 国语自产精品视频在线第100页| 最好的美女福利视频网| 最近最新中文字幕大全免费视频| 高清在线国产一区| 亚洲精品在线美女| 欧美日本中文国产一区发布| bbb黄色大片| 国产成人免费无遮挡视频| 国产成人欧美| 久久精品人人爽人人爽视色| 日本vs欧美在线观看视频| 亚洲国产日韩欧美精品在线观看 | 久久中文字幕一级| 美女扒开内裤让男人捅视频| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 波多野结衣一区麻豆| 国内精品久久久久精免费| 999久久久精品免费观看国产| 亚洲av美国av| 久99久视频精品免费| 久久国产亚洲av麻豆专区| 一区二区三区高清视频在线| 成人亚洲精品一区在线观看| 成人三级做爰电影| 一级a爱视频在线免费观看| 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看 | 琪琪午夜伦伦电影理论片6080| 日韩有码中文字幕| 欧美乱码精品一区二区三区| 99re在线观看精品视频| www.www免费av| 久久人人爽av亚洲精品天堂| 黄色 视频免费看| 19禁男女啪啪无遮挡网站| 日韩有码中文字幕| 又黄又粗又硬又大视频| 男女做爰动态图高潮gif福利片 | 一区二区三区国产精品乱码| 中文字幕av电影在线播放| 狠狠狠狠99中文字幕| 久久亚洲真实| 高清毛片免费观看视频网站| 欧美中文综合在线视频| 国产成人精品久久二区二区免费| 免费无遮挡裸体视频| 久久国产精品男人的天堂亚洲| 久久久水蜜桃国产精品网| 18美女黄网站色大片免费观看| 性欧美人与动物交配| 久久久久久久久免费视频了| 在线十欧美十亚洲十日本专区| 日韩免费av在线播放| 国产主播在线观看一区二区| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 色综合站精品国产| 十八禁人妻一区二区| 99久久精品国产亚洲精品| aaaaa片日本免费| av天堂久久9| 欧美乱码精品一区二区三区| 一边摸一边做爽爽视频免费| 午夜两性在线视频| 女人被躁到高潮嗷嗷叫费观| 亚洲九九香蕉| √禁漫天堂资源中文www| 久久国产精品影院| 成人免费观看视频高清| 精品日产1卡2卡| 欧美日韩精品网址| 国产精品二区激情视频| av片东京热男人的天堂| 亚洲色图av天堂| 亚洲,欧美精品.| 在线观看日韩欧美| 欧美性长视频在线观看| 精品国产美女av久久久久小说| 桃红色精品国产亚洲av| 精品高清国产在线一区| 国产又爽黄色视频| 欧美成人一区二区免费高清观看 | 欧美日韩亚洲综合一区二区三区_| e午夜精品久久久久久久| 99国产精品一区二区三区| avwww免费| 男人的好看免费观看在线视频 | 亚洲天堂国产精品一区在线| 黑人巨大精品欧美一区二区蜜桃| 91av网站免费观看| 中文字幕最新亚洲高清| 此物有八面人人有两片| av网站免费在线观看视频| 级片在线观看| 自线自在国产av| 日韩欧美免费精品| 日韩欧美国产在线观看| 神马国产精品三级电影在线观看 | 黑丝袜美女国产一区| 少妇的丰满在线观看| 一二三四在线观看免费中文在| www.999成人在线观看| 99在线视频只有这里精品首页| 中文字幕高清在线视频| 亚洲片人在线观看| av福利片在线| 老司机在亚洲福利影院| 亚洲av日韩精品久久久久久密| 老司机午夜福利在线观看视频| 麻豆成人av在线观看| 十八禁网站免费在线| 午夜老司机福利片| 麻豆成人av在线观看| 天堂影院成人在线观看| 国产欧美日韩一区二区三区在线| 大型av网站在线播放| 最近最新中文字幕大全电影3 | 啦啦啦免费观看视频1| 麻豆av在线久日| 黄色a级毛片大全视频| 伊人久久大香线蕉亚洲五| 制服丝袜大香蕉在线| 成人18禁在线播放| 涩涩av久久男人的天堂| 欧美在线一区亚洲| 高潮久久久久久久久久久不卡| 此物有八面人人有两片| 1024视频免费在线观看| 国产亚洲精品综合一区在线观看 | 村上凉子中文字幕在线| 国产成人影院久久av| 久久 成人 亚洲| 欧美在线黄色| 后天国语完整版免费观看| 国产熟女xx| 婷婷六月久久综合丁香| 可以在线观看的亚洲视频| 99国产精品一区二区三区| 亚洲第一青青草原| 欧美黄色淫秽网站| 悠悠久久av| 精品免费久久久久久久清纯| 一级毛片女人18水好多| 久久精品国产亚洲av香蕉五月| 精品欧美国产一区二区三| 无限看片的www在线观看| 亚洲成人国产一区在线观看| 国产单亲对白刺激| 精品国产亚洲在线| 一夜夜www| 亚洲情色 制服丝袜| 午夜影院日韩av| 午夜免费激情av| 国产亚洲精品av在线| 欧美黑人精品巨大| 性少妇av在线| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 亚洲情色 制服丝袜| 国产亚洲精品第一综合不卡| 国产av又大| 老司机深夜福利视频在线观看| 午夜福利一区二区在线看| 在线观看免费午夜福利视频| 亚洲专区国产一区二区| 精品国产乱码久久久久久男人| 精品第一国产精品| 亚洲自拍偷在线| 精品免费久久久久久久清纯| 精品午夜福利视频在线观看一区| 天堂影院成人在线观看| 在线观看免费视频日本深夜| 黑人巨大精品欧美一区二区蜜桃| 国产精品,欧美在线| 欧美在线一区亚洲| 大码成人一级视频| 亚洲国产欧美网| 久久久久久人人人人人| www.999成人在线观看| 日韩高清综合在线| 午夜成年电影在线免费观看| 久久精品亚洲精品国产色婷小说| 国产男靠女视频免费网站| av有码第一页| 真人一进一出gif抽搐免费| 国产在线精品亚洲第一网站| 久久国产精品影院| 日本精品一区二区三区蜜桃| 丝袜人妻中文字幕| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 在线国产一区二区在线| 国产男靠女视频免费网站| 久久久久久久精品吃奶| 国产片内射在线| 少妇 在线观看| 日本五十路高清| 亚洲成人精品中文字幕电影| 日本 欧美在线| 久久人妻av系列| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 国产高清激情床上av| 一级,二级,三级黄色视频| 亚洲欧美日韩高清在线视频| 婷婷精品国产亚洲av在线| 亚洲无线在线观看| 日韩一卡2卡3卡4卡2021年| 久久 成人 亚洲| 国产99白浆流出| 99热只有精品国产| 99riav亚洲国产免费| 久久久久精品国产欧美久久久| 欧美不卡视频在线免费观看 | 无人区码免费观看不卡| 在线观看舔阴道视频| 色在线成人网| 亚洲欧美一区二区三区黑人| 欧美日韩亚洲国产一区二区在线观看| 欧美激情久久久久久爽电影 | 欧美日本视频| 久久久久久大精品| 欧美最黄视频在线播放免费| 中文字幕人妻丝袜一区二区| 成人18禁高潮啪啪吃奶动态图| 欧美在线黄色| 欧美一级毛片孕妇| 人人妻人人澡欧美一区二区 | 久久国产精品男人的天堂亚洲| 女警被强在线播放| 国产成人欧美| 日韩欧美免费精品| 亚洲天堂国产精品一区在线| 99re在线观看精品视频| 老司机午夜福利在线观看视频| 嫩草影视91久久| 欧美乱妇无乱码| 国产精品99久久99久久久不卡| 亚洲九九香蕉| 好看av亚洲va欧美ⅴa在| 亚洲精品国产精品久久久不卡| 国产片内射在线| 中文字幕人成人乱码亚洲影| 91麻豆av在线| 欧美av亚洲av综合av国产av| 久久久久久亚洲精品国产蜜桃av| 亚洲全国av大片| 国产亚洲av嫩草精品影院| 国产乱人伦免费视频| 少妇裸体淫交视频免费看高清 | 黄色片一级片一级黄色片| 动漫黄色视频在线观看| 亚洲五月婷婷丁香| 一级a爱视频在线免费观看| 日本a在线网址| 9191精品国产免费久久| 久久久精品欧美日韩精品| 国内精品久久久久精免费| 亚洲视频免费观看视频| 中文字幕最新亚洲高清| 18禁国产床啪视频网站| 久久影院123| 婷婷六月久久综合丁香| svipshipincom国产片| 啦啦啦免费观看视频1| 国产亚洲欧美98| 久久人人97超碰香蕉20202| 欧美色视频一区免费| 亚洲欧美激情在线| 大香蕉久久成人网| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av片天天在线观看| 一个人免费在线观看的高清视频| 色播在线永久视频| 露出奶头的视频| 久久青草综合色| 黄片小视频在线播放| 成人18禁在线播放| 午夜免费观看网址| 国产日韩一区二区三区精品不卡| 99热只有精品国产| 欧美一级a爱片免费观看看 | 99re在线观看精品视频| 校园春色视频在线观看| 亚洲美女黄片视频| 女人精品久久久久毛片| 免费看a级黄色片| 午夜老司机福利片| 久久久久亚洲av毛片大全| 一进一出抽搐动态| 韩国精品一区二区三区| 欧美色欧美亚洲另类二区 | 国产精品久久电影中文字幕| 在线观看一区二区三区| 亚洲国产精品久久男人天堂| av欧美777| 老司机福利观看| 村上凉子中文字幕在线| 精品欧美一区二区三区在线| 国产精品香港三级国产av潘金莲| 亚洲国产精品久久男人天堂| 午夜免费观看网址| 国产成人av激情在线播放| 美女高潮喷水抽搐中文字幕| 看免费av毛片| 久久久国产精品麻豆| 男男h啪啪无遮挡| av有码第一页| 91在线观看av| 欧美日韩乱码在线| 99久久99久久久精品蜜桃| 精品一品国产午夜福利视频| 又大又爽又粗| 中国美女看黄片| 亚洲在线自拍视频| 久久国产精品男人的天堂亚洲| 久久久久久久精品吃奶| 亚洲在线自拍视频| 村上凉子中文字幕在线| 1024视频免费在线观看| 亚洲国产日韩欧美精品在线观看 | 亚洲最大成人中文| 久久精品影院6| 老司机福利观看| 日本免费一区二区三区高清不卡 | 亚洲色图av天堂| 久久精品人人爽人人爽视色| 日韩欧美国产在线观看| 日本五十路高清| 制服诱惑二区| 亚洲午夜精品一区,二区,三区| 久久久久久大精品| 久久国产精品人妻蜜桃| 日韩成人在线观看一区二区三区| 香蕉国产在线看| 亚洲色图av天堂| 国产av精品麻豆| 精品第一国产精品| 亚洲国产精品999在线| 国产极品粉嫩免费观看在线| 国产一区二区三区视频了| 1024视频免费在线观看| 色老头精品视频在线观看| 俄罗斯特黄特色一大片| 欧美最黄视频在线播放免费| 日韩免费av在线播放| 精品国产一区二区三区四区第35| x7x7x7水蜜桃| 国产精品久久久久久亚洲av鲁大| 成年人黄色毛片网站| 成年版毛片免费区| 日韩视频一区二区在线观看| 久久午夜综合久久蜜桃| 国产在线观看jvid| 精品熟女少妇八av免费久了| 国产乱人伦免费视频| 亚洲狠狠婷婷综合久久图片| 久久人妻av系列| 一a级毛片在线观看| 国产成人av教育| 欧美成人午夜精品| 男女之事视频高清在线观看| 国产激情欧美一区二区| 少妇被粗大的猛进出69影院| 亚洲成a人片在线一区二区| 色哟哟哟哟哟哟| 最近最新免费中文字幕在线| 91九色精品人成在线观看| 久久精品影院6| av天堂久久9| 国产又色又爽无遮挡免费看| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜| 俄罗斯特黄特色一大片| 久久人妻av系列| 无遮挡黄片免费观看| 一边摸一边抽搐一进一小说| 精品久久久精品久久久| 黄频高清免费视频| 少妇 在线观看| 亚洲成av人片免费观看| 亚洲av成人不卡在线观看播放网| 亚洲成a人片在线一区二区| 久久国产精品男人的天堂亚洲| 法律面前人人平等表现在哪些方面| 两个人免费观看高清视频| 亚洲欧洲精品一区二区精品久久久| 一区福利在线观看| 最近最新中文字幕大全电影3 | 免费av毛片视频| 欧美一级a爱片免费观看看 | 国产午夜精品久久久久久| 亚洲av片天天在线观看| 99国产精品一区二区蜜桃av| 亚洲国产高清在线一区二区三 | 亚洲欧美精品综合久久99| 日本 av在线| www.熟女人妻精品国产| 宅男免费午夜| 国产精品野战在线观看| 免费高清在线观看日韩| 两性午夜刺激爽爽歪歪视频在线观看 | 精品人妻1区二区| 女人精品久久久久毛片| 男女做爰动态图高潮gif福利片 | 国产免费av片在线观看野外av| 最近最新中文字幕大全电影3 | 一级毛片高清免费大全| 日日干狠狠操夜夜爽| 在线观看免费视频日本深夜| 午夜精品在线福利| 久久香蕉国产精品| 色播在线永久视频| 亚洲国产精品成人综合色| 黄网站色视频无遮挡免费观看| 一区二区三区精品91| 亚洲 国产 在线| 老司机午夜福利在线观看视频| 午夜福利一区二区在线看| 午夜免费鲁丝| 91av网站免费观看| 一边摸一边做爽爽视频免费| 无人区码免费观看不卡| 在线观看一区二区三区| 怎么达到女性高潮| 欧美日韩乱码在线| 亚洲三区欧美一区| 性少妇av在线| 日韩欧美一区二区三区在线观看| 黄色毛片三级朝国网站| 大型av网站在线播放| 精品国产亚洲在线| 两人在一起打扑克的视频| 中文字幕人成人乱码亚洲影| 国产av一区在线观看免费| 婷婷六月久久综合丁香| 欧美中文综合在线视频| 亚洲一区二区三区色噜噜| 亚洲成人久久性| 在线观看日韩欧美| 成人国产一区最新在线观看| 久久久国产成人免费| 亚洲色图av天堂| 中国美女看黄片| 亚洲熟女毛片儿| 91国产中文字幕| 黄色 视频免费看| 国产国语露脸激情在线看| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 美女午夜性视频免费| 午夜成年电影在线免费观看| 国产亚洲精品第一综合不卡| 女人被狂操c到高潮| 涩涩av久久男人的天堂| 国内精品久久久久久久电影| 日本五十路高清| 亚洲成a人片在线一区二区| 国产av在哪里看| 99久久精品国产亚洲精品| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 成人国产综合亚洲| 亚洲性夜色夜夜综合| 日日爽夜夜爽网站| 亚洲av美国av| 久久热在线av| 波多野结衣av一区二区av| 大码成人一级视频| 亚洲第一av免费看| 亚洲人成伊人成综合网2020| 悠悠久久av| 午夜福利视频1000在线观看 | 亚洲熟妇中文字幕五十中出| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人精品巨大| 在线av久久热| 亚洲精品国产区一区二| 免费女性裸体啪啪无遮挡网站| 美国免费a级毛片| 亚洲精品国产精品久久久不卡| 午夜免费鲁丝| 亚洲成人国产一区在线观看| 91麻豆av在线| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 纯流量卡能插随身wifi吗| 国产色视频综合| 在线观看www视频免费| 制服诱惑二区| 国产一区二区三区视频了| 国产精品 国内视频| 亚洲专区字幕在线| 天堂动漫精品| 久久久久久久久免费视频了| 国产成人系列免费观看| 91麻豆av在线| 午夜福利欧美成人| 欧美日韩亚洲国产一区二区在线观看| 最好的美女福利视频网| 国产成人一区二区三区免费视频网站| 在线观看66精品国产| 免费看十八禁软件| 日韩精品中文字幕看吧| 久久久国产精品麻豆| 国产精品野战在线观看| 久久久久久国产a免费观看| 亚洲情色 制服丝袜| 国产高清激情床上av| 久久性视频一级片| 一二三四在线观看免费中文在| 1024视频免费在线观看| 又黄又粗又硬又大视频| 男女之事视频高清在线观看| 一本综合久久免费| 日韩欧美免费精品| 亚洲一区二区三区不卡视频| 成人国产一区最新在线观看| 国产欧美日韩一区二区精品| 亚洲av成人一区二区三| 亚洲专区国产一区二区| 国产成年人精品一区二区| 国内精品久久久久久久电影| 女人高潮潮喷娇喘18禁视频| 亚洲欧美激情综合另类| 乱人伦中国视频| 真人一进一出gif抽搐免费| 亚洲av片天天在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 嫩草影视91久久| 国产单亲对白刺激| 亚洲久久久国产精品| 不卡av一区二区三区| 日韩有码中文字幕| 亚洲视频免费观看视频| 动漫黄色视频在线观看| 夜夜夜夜夜久久久久| 18禁国产床啪视频网站| 天天一区二区日本电影三级 | 亚洲美女黄片视频| 亚洲欧美日韩高清在线视频| 黄频高清免费视频| 18禁观看日本| 成人18禁在线播放| 午夜亚洲福利在线播放| 日韩免费av在线播放| 亚洲天堂国产精品一区在线| 精品免费久久久久久久清纯| 国产欧美日韩一区二区三区在线| 伦理电影免费视频| 亚洲av第一区精品v没综合| 男人舔女人的私密视频| 国产aⅴ精品一区二区三区波| 日本a在线网址| 成人国产一区最新在线观看| 久久热在线av| 久久人妻福利社区极品人妻图片| 少妇粗大呻吟视频| 久久久久精品国产欧美久久久| 久久人妻熟女aⅴ| 好男人电影高清在线观看| 中文字幕最新亚洲高清| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| 亚洲天堂国产精品一区在线| 两个人视频免费观看高清| 韩国精品一区二区三区| 我的亚洲天堂| 色精品久久人妻99蜜桃|