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

    Effective Frameworks Based on Infinite Mixture Model for Real-World Applications

    2022-08-24 12:56:18NorahSalehAlghamdiSamiBourouisandNizarBouguila
    Computers Materials&Continua 2022年7期

    Norah Saleh Alghamdi, Sami Bourouisand Nizar Bouguila

    1Department of Computer Sciences, College of Computer and Information Sciences,Princess Nourah bint Abdulrahman University, P.O.Box 84428, Riyadh 11671, Saudi Arabia

    2College of Computers and Information Technology, Taif University, Taif, 21944, Saudi Arabia

    3The Concordia Institute for Information Systems Engineering (CIISE), Concordia University,Montreal, QC H3G 1T7, Canada

    Abstract: Interest in automated data classification and identification systems has increased over the past years in conjunction with the high demand for artificial intelligence and security applications.In particular, recognizing human activities with accurate results have become a topic of high interest.Although the current tools have reached remarkable successes, it is still a challenging problem due to various uncontrolled environments and conditions.In this paper two statistical frameworks based on nonparametric hierarchical Bayesian models and Gamma distribution are proposed to solve some realworld applications.In particular, two nonparametric hierarchical Bayesian models based on Dirichlet process and Pitman-Yor process are developed.These models are then applied to address the problem of modelling grouped data where observations are organized into groups and these groups are statistically linked by sharing mixture components.The choice of the Gamma mixtures is motivated by its flexibility for modelling heavy-tailed distributions.In addition, deploying the Dirichlet process prior is justified by its advantage of automatically finding the right number of components and providing nice properties.Moreover, a learning step via variational Bayesian setting is presented in a flexible way.The priors over the parameters are selected appropriately and the posteriors are approximated effectively in a closed form.Experimental results based on a real-life applications that concerns texture classification and human actions recognition show the capabilities and effectiveness of the proposed framework.

    Keywords: Infinite Gamma mixture model; variational Bayes; hierarchical Dirichlet process; Pitman-Yor process; texture classification; human action recognition

    1 Introduction and Literature Review

    Data clustering has been the subject of wide research to the present days [1-4].The goal of clustering is to group the observed data into separate subgroups, such as the data in each subgroup shows some similarities to each other.Various clustering approaches have been developed in the past and some of them are based on distance metrics (such as the K-means and SOM algorithms).It is noteworthy that these algorithms are sensitive to noise, are not flexible when dealing with incomplete data, and may not definitely capture heterogeneity inherent in complex datasets.Another important alternative is the model-based clustering approach [5].It is effective for modelling the structure of clusters since data are supposed to be produced by some distributions.In particular, finite mixture models (such as Gaussian mixture model) have shown to be effective (in terms of discovering complex patterns and grouping them into similar clusters) for many computer vision and pattern recognition applications [6].Nevertheless, Gaussian model cannot fit well complex non-Gaussian shapes.To cope with the disadvantages related to conventional Gaussian assumption, many contributions have occurred to develop other flexible mixture models and were interested to non-Gaussian behavior of real datasets [7,8].For instance,Gamma (GaMM) mixtures have demonstrated to offer high flexibility and ease of use than Gaussian [9] for many image processing and pattern recognition problems.This is due their compact analytical form which is able to cover long-tailed distributions and to approximate data with outliers.This mixture has been used with success for a range of interesting problems [7,10-12]especially when dealing with proportional data.An illustration of how mixture models can be applied for data modelling, classification and recognition is given in Fig.1.

    Figure 1: Block diagram of application of mixture model for data classification and recognition

    Unfortunately, the inadequacy of finite mixture model has been apparent when selecting the appropriate number of mixture components.In other word, model selection (i.e., model’s complexity)is one of the difficult problems within finite mixture models.The crucial problem of“how many groups in the dataset?”still remains of great interest for various data mining fields since determining inappropriate number of clusters may conduct to poor generalization capability.This problem can be solved by considering an infinite number of components via nonparametric Bayesian methods[13,14], principally within“Dirichlet process (DP)”[15].Indeed, a Dirichlet Process (DP) is a parameterized stochastic process characterized by a base distribution and it can be defined as a probability distribution over discrete distributions.Numerous studies have been devoted to infinite mixture models which are emerged to cope with the challenging problem of model selection.

    Two sound hierarchical Bayesian alternatives to the conventional DP named as hierarchical Dirichlet process (HDP) and hierarchical Pitman-Yor process (HPYP) have exposed encouraging results especially when dealing with modeling grouped data [16,17].Indeed,within hierarchical models,mixture components can be shared through data groups (i.e., parameters are shared among groups).In such cases, it is possible to make a Bayesian hierarchy on the DP and the base distribution of the DP is distributed according to another DP.

    Thus, the main contributions of this manuscript are first to extend our previous works about the Gamma mixture by investigating two efficient nonparametric hierarchical Bayesian models based on both Dirichlet and Pitman-Yor processes mixtures of Gamma distributions.Indeed, in order to reach enhanced modeling performance, we consider Gamma distribution which is able to cover longtailed distributions and to approximate accurately visual vectors.Another critical issue when dealing with mixture models is model parameters estimation.Accordingly, we propose to develop an effective variational Bayes learning algorithm to estimate the parameters of the implemented models.It is noteworthy to indicate that the complexity of variational inference-based algorithm still remains less than Markov Chain Monte Carlo-Based Bayesian inference and leads to faster convergence.Finally,the implemented hierarchical Bayesian models and the variational inference approaches are validated via challenging real-life problems namely human activity recognition and texture categorization.

    This manuscript is organized as follows.In Section 2 we introduce the hierarchical DP and PYP mixtures of Gamma distributions which are based on stick-breaking construction.In Section 3, we describe the details of our variational Bayes learning framework.Section 4 reports the obtained results, which are based on two challenging applications, to verify the merits and effectiveness of our framework, and Section 5 is devoted to conclude the manuscript.

    2 Model Specification

    In this section, we start by briefly presenting finite Gamma mixture model and then we present our nonparametric frameworks based on hierarchical Dirichlet and Pitman-Yor processes mixtures.

    2.1 Finite Gamma Mixture Model

    If aD-dimensional random vector= (Y1,...,YD) is distributed according to a multidimensional Gamma distribution, then its probability density function (pdf) is defined as:

    2.2 Hierarchical Dirichlet Process Mixture Model

    The hierarchical Dirichlet process (HDP) is an an effective nonparametric Bayesian method to modelling grouped data, which allows the mixture models to share components.Here, observed data are arranged into groups (i.e., mixture model) that we want to make them statistically linked.HDP It is built on the Dirichlet process (DP) as well described in [17] for each group data.It is noteworthy that the DP has acquired popularity in machine learning to handle nonparametric problems [18].The DP was presented as a prior on probability distributions and this makes it extremely appropriate for specifying infinite mixture models thanks to the use of the stick-breaking process [19].In the case of hierarchical Dirichlet process (HDP), the DPs for all groups share a base distribution which is itself distributed according to a Dirichlet process.Let’s assume that we have a grouped data set Y separated intoMgroups, such that each group is associated with a DPGj, thus the HDP takes part an indexed set of DPsGjthat share a global (or base) distributionG0which is itself distributed as a DP with base distributionHand concentration parameter γ:

    Here, the hierarchical Dirichlet process is represented using the stick-breaking construction[19,20].The global measureG0is distributed according toDP(γ,H) and it can be expressed as

    where δΩkrepresents an atom concentrated at Ωk, and {Ωk} is a set of independent random variables drawn fromH.The variable {ψk} are denoted as stick-breaking that verify= 1.AsG0is defined as the base distribution of the DPGjand has the stick-breaking representation as shown in Eq.(3),thenGjincludes all the atoms Ωkwith distinct weights (by following the property of Dirichlet process[18]).On the other side, we carry out another stick-breaking process to construct each group-level DPGjaccording to [21] as

    where {πjt} is a set of stick-breaking weights which shall be positive and sum to one and δωjtare grouplevel Dirac delta atoms at ωjt.As ωjt(group-level atom) is distributed according toG0, thus it will take onto Ωk(base-level atoms) with probability ψk.

    Next, we introduce a latent indicatorWjtkas an indicator variable, such thatWjtk∈{0,1} (in order to indicate which group-level atom maps to).Wjtk= 1 if ωjt(group-level atom) maps to the Ωk(globallevel atom) that is indexed byk; otherwise,Wjtk= 0.Accordingly, we can have ωjt=.By this way, there is no need to keep a representation for ωjt.The indicator variable= (Wjt1,Wjt2,...) is distributed as:

    Given that ψ is a function of ψ′according to Eq.(4), it is possible to rewrite the indicator variablep() as:

    According to Eq.(4) the stick lengths ψ′are drawn from a specific Beta distribution and their realization is determined as

    To complete the description of the HDP mixture model, given a grouped observation (data) Y,we associate each data pointwith a factor(hereiindexes the data in each groupj) such thatand= (θj1,θj2,...) are distributed according toF(θji) andGj, respectively.In this case the likelihood function can be written as:

    whereF(θji) is the probability distribution ofYjigiven θji.The base distributionHprovides the prior distribution for θji.This setting (i.e., hierarchical Dirichlet process (HDP) mixture model) plays an important role and ensures that each group is associated with a mixture model, and the components of the mixture are shared across different groups.

    As each θjiis distributed according toGj(see Eq.(9)), it takes the value ωjtwith probability πjt.We then introduce another indicator variableZjit∈{0,1} for θjias

    That is, the indicatorZjitis used to indicate which component θjibelongs to.In particular,Zjitis equal to 1 if θjiis associated with componentt(also maps to the group-level atom ωjt); else,Zjit= 0.Thus, we can write.As ωjtmaps to Ωk, we then can write.

    Based on the stick-breaking construction of the DP (see Eq.(5)), we notice thatis a function of, therefore the previous equation becomes

    Finally, according to Eq.(5), the prior distribution of π′is a Beta and it is given as follows:

    2.3 Hierarchical Pitman-Yor Process Mixture Model

    The Pitman-Yor process (PYP) [22] is a two-parameter extension to the DP (i.e., is a generalization of DP) that permits modelling heavier-tailed distributions.It can be applied to build hierarchical models.It offers a sophisticated way to cluster data such that the number of clusters is unknown.It is characterized by an additional discount parameter γain addition to the concentration parameter γb, that satisfying the 0<γa<1,γb>-γa.Similar to DP, the sample drawn from PYP also associated a probability measureH[23].Here, a hierarchical Pitman-Yor process (HPYP) is introduced where the base measure for a PYP is itself a draw from a PYP.Specifically, HPYP defines the global-level measureG0and group-level distributionGj(that is the indexedGjshares a same baseG0which itself follows a PYP).This behavior makes the HPYP especially suitable for complex visual data modeling and classification.We can use the HPYP to cluster data by applying the stick-breaking construction that defines the base measure as follows:

    where {Λk} is a set of independent random variables drawn fromHand δΛkis an atom (probability mass) at Λk.The random variables ηkrepresent the stick-breaking weights that satisfying= 1.The stick-breaking representation for the group-level PYPGjis expressed as:

    where {pjt} are the stick-breaking weights.ψjtis the atom of second-level PYP that is distributed according toG0.Then, a global-level indicator variablesIand a group-level indicator variablesCare introduced.Here,Cis used to map θjito group-level atom ψjtand the indicatorIis used to map the atom θjito base-level atom Λk.

    2.4 Hierarchical Infinite Gamma Mixture Model

    In this subsection, we introduce two hierarchical infinite mixture models with Gamma distributions.In this case, each vector= (Yji1,...,YjiD),fromthe grouped data, is drawn froma hierarchical infinite Gamma mixture model.Then, the two likelihood functions of these hierarchical models given the unknown parameters of Gamma and latent variables can be expressed as follows:

    Next, we have to place conjugate distributions over the unknown parameters α and β.As α and β are positive, then it is convenient that they follow Gamma distributions G(.).Thus, we have

    3 Model Learning via Variational Bayes

    Variational inference [3,24] is a well-defined method deterministic approximation method that is used in order to approximate posterior probability via an optimization process.In this section,we propose to develop a variational learning framework of our hierarchical infinite Gamma mixture models.Here,Θ=represents both unknown and the latent variables.Our objective is to estimate a suitable approximationq(Θ) for the true posterior distributionp(Θ|Y) via a process of maximizing the lower bound of lnp(Y) given as

    In particular, we adopt one of the most successfully variational inference techniques namely the factorial approximation (or mean fields approximation) [3], which is able to offer effective updates.Thus, we apply this method to fully factorizeq(Θ) of on HDPGaM and HPYPGaM mixtures into disjoint factors.Then, we apply a truncation method as previously applied in [20] to truncate the variational approximations into global truncation levelKand group truncation levelTas follow:

    where the truncation levelsKandTwill be optimized over the learning procedure.The approximated posterior distribution is then factorized as

    For a specific variational factorqs(Θs), the general equation of the optimal solution is expressed as [3]:

    where〈.〉i/=sdenotes an expectation with respect to all the distributions ofqi(Θi) except fori=s.The parametric forms for the variational posteriors (for each factor) are determined on the basis of Eq.(26) as

    where the corresponding hyperparameters in the above equations can be calculated as a similar way in[25,26].The complete variational Bayes inference algorithms of both HDPGaM and HPYPGaM are summarized in Algorithm 1 and Algorithm 2, respectively.

    Algorithm 1: HDPGaM: Proposed Hierarchical Dirichlet process gamma mixture algorithmimages/BZ_1121_265_950_1515_1565.png

    Algorithm 2: HPYPGaM: Proposed Hierarchical Pitman-Yor process gamma mixture algorithmimages/BZ_1121_265_1730_1515_2347.png

    4 Experimental Results

    The principal purpose of the experiment section is to investigate the performance of the developed two frameworks based on HDP mixture and HPYP mixture model with Gamma distributions.Hence,we propose to compare them with other statistical models using two challenging applications: Texture categorization and human action categorization.In all these experiments, the global truncation levelKand the group level truncation levelTare both initialized to 120 and 60, respectively.For HDP mixture, we set the hyperparameters of the stick lengths γ and λ as (0.25, 0.25).The parameters of HPYP mixture γa,γb,βaand βbare initialized to (0.5, 0.25, 0.5, 0.25).The hyperparameters of Gamma base distribution are initialized by sampling from priors.

    4.1 Texture Classification

    In thisworkwe are primarily motivated by the problem of modeling and classifying texture images.Contrary to natural images which include certain objects and structures, texture images are very special case of images that do not include a well-defined shape.Texture pattern is one of the most important elements in visual multimedia content analysis.It forms the basis for solving complex machine learning and computer vision tasks.In particular, texture classification supports a wide range of applications, including information retrieval, image categorization [27-30], image segmentation [27,31,32], material classification [33], facial expression recognition [29], and object detection [28,34].The goal here to classify texture images using the two hierarchical infinite mixtures and also by incorporating three different representations (to extract relevant features from images) from the literature, namely Local Binary Pattern (LBP) [35], Local Binary Pattern (LBP) [35], scale-invariant feature transform (SIFT) [36], and dense micro-block difference (DMD) [37].A deep review for these methods is outside the scope of current work.Instead, we focus on some powerful feature extraction methods that have shown interesting state-of-the-art results.

    4.1.1 Methodology

    For this application, we start by extracting features from input images and then we model them using the proposed HDPGaM and HPYPGaM.Each imageIjis considered as a group and is related to a infinite mixture modelGj.Next, every vectorYjiofIjis assumed generated fromGj, whereGjrepresents visual words.The next step is to generate a global vocabulary to share it among all groups viaG0(global infinite model).It is noteworthy that the the building of the visual vocabulary, here, is part of our hierarchical models and therefore, the size of the vocabulary (i.e., number of components)is is inferred automatically from the data thanks to the characteristic of nonparametric Bayesian models.Regarding SIFT features, the bag-of-visual-words model is adopted here to calculate the histogram of visual words from each input image.Regarding the set of DMD descriptors, these are obtained after extracting DMD features and then encoding them though th Fisher vector method as proposed in[37].The resulting descriptors are able to attain good discrimination thanks to their invariance with respect to scale, resolution, and orientation.Finally, each image is represented by a multidimensional vector of high-order statistics1The Matlab code for the features is available at http://www.cs.tut.fi/ mehta/texturedmd..

    4.1.2 Dataset and Results

    We conducted our experiments of texture classification using the proposed hierarchical HDP Gamma mixture (referred to as HDPGaM) and HPYP Gamma mixture (referred to as HPYPGaM)on three publicly available databases.The first one namely UIUCTex [38] contains 25 texture classes and each class has 40 images.The second dataset namely UMD [39] contains 25 textures classes containing each one 40 images.The third dataset namely KTH-TIPS [33] includes 10 classes and each one contains 81 images.Some sample texture images from each class and each dataset are shown in Fig.2.We use 10-fold cross-validation technique to partition these databases and to study the performance.In addition, the evaluation process and the obtained results are based on 30 runs.

    In order to quantify the performance of the proposed frameworks (HDPGaM and HPYPGaM),we proceed by evaluating and comparing the obtained results with seven other methods namely infinite mixture of Gaussian distribution (inGM), infinite mixture of generalized Gaussian distribution (inGGM), infinite mixture of multivariate generalized Gaussian distribution (inMGGM),Hierarchical Dirichlet Process mixture of Gaussian distribution (HDPGM), hierarchical Pitman-Yor process mixture of Gaussian distribution (HPYPGM), Hierarchical Dirichlet Process mixture of generalized Gaussian distribution (HDPGGM), and hierarchical Pitman-Yor process mixture of generalized Gaussian distribution (HPYPGGM).

    Figure 2: Texture samples in different categories for different datasets (a) KTH-TIPS, (b) UIUCTex,and (c) UMD dataset

    We run all methods 30 times and calculate the average classification accuracy which are depicted in Tabs.1-3 respectively.According to these results, we can notice that HDPGaM and HPYPGaM have the highest achieved accuracy for the three databases in terms of the texture classification accuracy rate.It is noted that when comparing these results by considering the Student’s t-test, the differences in performance are statistically significant between our frameworks and the rest of methods.In particular, results indicate the benefits of our proposed models in terms of texture modeling and classification capabilities which surpass those obtained by HDPGM, HDPGGM,HPYPGM, andHPYPGGM.By contrast, the worst performance is obtained within the infinite Gaussian mixture models.It should be noted that the proposed frameworks outperform the other methods and that the three adopted feature extraction methods (SIFT, LBP and DMD).Thus, these results confirm the merit of the proposed methods.Due to the effectiveness of DMD descriptor for describing and modelling observed texture images, we also find that DMDachieves better accuracy compared with both SIFT and LBP.It shows the merits of DMD which is able to consider all possible fine details images at different resolutions.We also note that with HPYP mixture we can reach better results compared to HDP mixture and this is for all tested distributions.This can be explained by the fact that HPYP mixture model has better generalization capability and better capacity to model heavier-tailed distribution (PYP prior can lead to better modeling ability).

    Table 3: The average accuracy results of texture classification using different algorithms for the 3 texture-datasets using DMD features

    Table 1: The average accuracy results of texture classification using different algorithms for the 3 texture-datasets using SIFT features

    Table 2: The average accuracy results of texture classification using different algorithms for the 3 texture-datasets using LBP features

    Table 3: Continued

    4.2 Human Actions Categorization

    Visual multimedia recognition has been a challenging research topic which could attract many applications such as actions recognition [40,41], image categorization [42,43], biomedical image recognition [44], and facial expressions [29,30].In this work, we are focusing on a particular problem that has received a lot of attention namely Human actions recognition (HAR) through sequence of videos.Indeed, the intention of recognizing activities is to identify and analyze various human actions.At present, HAR is one the hot computer vision topics not only in research but also in industries where automatic identification of any activity can be useful, for instance, for monitoring, healthcare,robotics, and security-based applications [45].Recognizing manually activities is very challenging and time consuming.This issue has been addressed and so various tools have been implemented such as in [40,45-47].However, precise recognition of actions is still required using advanced and efficient algorithms in order to deal with complex situations such as noise, occlusions, and lighting.

    We perform here the recognition of Human activities using the proposed frameworks HDPGaM and HPYPGaM.Our methodology is outlined as following: First, we extract and normalize SIFT3D descriptor [40] from observed images.These features are then quantized as visual words via bag-ofwords (BOW) model and K-means algorithm [48].Then, these features are quantized as visual words via K-means algorithm [48].Then, a probabilistic Latent Semantic Analysis (pLSA) [49] is adopted to construct ad-dimensional vector.In particular, each imageIjis considered as a“group”and is associated with an infinite mixture modelGj.Thus, we suppose that Each SIFT3D feature vector is drawn from the infinite mixture modelGjand“visual words”denote mixture components ofGj.

    On the other hand, a global vocabulary is generated and shared between all groups via the global-modelG0of the proposed hierarchical models.This setting agrees with the purpose behind the hierarchical process mixture model.It is also noted that the building of the visual vocabulary is part of the hierarchical process mixture models and this step is not carried out separately via the k-means algorithm as many other approaches do.It is also noted that the generation of the visual vocabulary is part of our hierarchical process mixture models and this process is not performed separately via kmeans as many other approaches did.In fact, it is due to the property of the nonparametric Bayesian model that the number of components in the global-level mixture model can be deduced from the data.We conducted our experiments of actions recognition using a publicly dataset known as KTH human action dataset [46].This database contains 2391 sequences of different actions grouped into 6 classes.It also represents four scenarios (outdoors (s1), outdoors with scale variation (s2), outdoors with different clothes (s3) and indoors (s4)).Some scenarios from this dataset are given in Fig.3.We randomly divided this dataset into 2 subsets to train the developed frameworks and to evaluate its robustness.

    Our purpose through this application is to show the advantages of investigating our proposed hierarchical models HDPGaM and HPYPGaM over other conventional hierarchical mixtures and other methods from the state of the art.Therefore, we focused first on evaluating the performance of HDPGaM and HPYPGaM over Hierarchical Dirichlet Process mixture of Gaussian distribution(HDPGM), hierarchical Pitman-Yor process mixture of Gaussian distribution (HPYPGM), Hierarchical Dirichlet Process mixture of generalized Gaussian distribution (HDPGGM), and hierarchical Pitman-Yor process mixture of generalized Gaussian distribution (HPYPGGM).It is noted that we learned all the implemented models using variational Bayes.The average recognition performances of our frameworks and models based on HDP mixture and HPYP mixture are depicted in Tab.4.

    Figure 3: Sample frames of the KTH dataset actions with different scenarios

    Table 4: Average recognition performance (%) obtained using our frameworks and other models based on HDP mixture and HPYP mixture for KTH database

    As we can see in this table, the proposed frameworks were able to offer the highest recognition rates (82.27% for HPYPGaM and 82.13% for HDPGaM) among all tested models.For different runs, we have p-values<0.05 and therefore, the differences in accuracy between our frameworks and other models are statistically significant according to Student’s t-test.Next, we compared our models against other mixture models (here finite Gaussian mixture (GMM) and finite generalized Gaussian mixture (GGMM) and methods from the literature.The obtained results are given in Tab.5.

    Table 5: Average recognition performance (%) obtained using our frameworks and other methods from the literature for KTH database

    Accordingly, we can observe that models again are able to provide higher discrimination rate than the other methods.Clearly, these results confirm the effectiveness of our frameworks for activities modeling and recognition compared to other conventional Dirichlet and Pitman-Yor processes based on Gaussian distribution.Another remark is that our model HPYPGaM outperforms our second model HDPGaM for this specific application and this demonstrates the advantages of using hierarchical Pitman-Yor process over Dirichlet process which is flexible enough to be used for such recognition problem.

    5 Conclusions

    In this paper two non-parametric Bayesian frameworks based on both hierarchical Dirichlet and Pitman-Yor processes and Gamma distribution are proposed.The Gamma distribution is considered because of its flexibility for semi-bounded data modelling.Both frameworks are learned using variational inference which has certain advantages such as easy assessment of convergence and easy optimization by offering a trade-off between frequentist techniques and MCMC-based ones.An important property of our approach is that it does not need the specification of the number of mixture components in advance.We carried out experiments on texture categorization and human action recognition to demonstrate the performance of our models which can be used further for a variety of other computer vision and pattern recognition applications.

    Acknowledgement:The authors would like to thank Taif University Researchers Supporting Project number(TURSP-2020/26),Taif University,Taif,Saudi Arabia.They would like also to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R40), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

    Funding Statement:The authors would like to thank Taif University Researchers Supporting Project number(TURSP-2020/26),Taif University,Taif,Saudi Arabia.They would like also to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R40),Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

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

    纯流量卡能插随身wifi吗| 午夜福利影视在线免费观看| 日本撒尿小便嘘嘘汇集6| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩综合在线一区二区| 午夜福利一区二区在线看| 50天的宝宝边吃奶边哭怎么回事| 国产精品 国内视频| 考比视频在线观看| 成人av一区二区三区在线看| 午夜日韩欧美国产| 免费久久久久久久精品成人欧美视频| 桃红色精品国产亚洲av| 女人久久www免费人成看片| 97在线人人人人妻| 又紧又爽又黄一区二区| 久久久久久久久免费视频了| 十分钟在线观看高清视频www| 欧美日韩亚洲综合一区二区三区_| 欧美黑人精品巨大| 三上悠亚av全集在线观看| 我要看黄色一级片免费的| 国产高清国产精品国产三级| 亚洲熟女精品中文字幕| 亚洲av电影在线进入| av天堂在线播放| 国产精品99久久99久久久不卡| 久久人妻av系列| 日韩一卡2卡3卡4卡2021年| 久久久久久久久久久久大奶| 岛国毛片在线播放| 女人高潮潮喷娇喘18禁视频| 老熟妇乱子伦视频在线观看| 午夜免费鲁丝| 久久久久国产一级毛片高清牌| 丝袜人妻中文字幕| 天天添夜夜摸| 一进一出好大好爽视频| 99国产极品粉嫩在线观看| 又紧又爽又黄一区二区| 一本久久精品| 91老司机精品| 精品午夜福利视频在线观看一区 | 啦啦啦 在线观看视频| av欧美777| 大型黄色视频在线免费观看| 少妇 在线观看| 午夜福利视频精品| 正在播放国产对白刺激| 国产成人av教育| 久久免费观看电影| 最近最新免费中文字幕在线| 亚洲九九香蕉| 免费少妇av软件| 国产精品熟女久久久久浪| 国产精品一区二区免费欧美| 一本色道久久久久久精品综合| 91国产中文字幕| 男女下面插进去视频免费观看| 99国产精品免费福利视频| 天天影视国产精品| 中文字幕制服av| 欧美日韩国产mv在线观看视频| 久久ye,这里只有精品| av片东京热男人的天堂| 男女下面插进去视频免费观看| 亚洲精品美女久久久久99蜜臀| 999久久久国产精品视频| 丝袜美足系列| 在线播放国产精品三级| 成人三级做爰电影| 人妻久久中文字幕网| 无人区码免费观看不卡 | 欧美精品av麻豆av| 大码成人一级视频| 国产欧美日韩综合在线一区二区| 桃红色精品国产亚洲av| 电影成人av| 亚洲七黄色美女视频| 亚洲性夜色夜夜综合| 日本黄色日本黄色录像| 久久天堂一区二区三区四区| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠躁躁| 超碰97精品在线观看| 水蜜桃什么品种好| 我要看黄色一级片免费的| 美国免费a级毛片| 欧美日韩亚洲国产一区二区在线观看 | 久久亚洲真实| av有码第一页| 黄色视频不卡| 捣出白浆h1v1| 亚洲av美国av| 亚洲视频免费观看视频| 搡老岳熟女国产| 亚洲国产成人一精品久久久| 国产精品亚洲一级av第二区| 曰老女人黄片| 在线观看舔阴道视频| 久久99热这里只频精品6学生| www.精华液| 免费av中文字幕在线| 亚洲综合色网址| 国产成人欧美在线观看 | 国产熟女午夜一区二区三区| 成人手机av| 国产av又大| 91国产中文字幕| tocl精华| 99热国产这里只有精品6| 中文字幕人妻熟女乱码| 香蕉国产在线看| 免费看十八禁软件| 久久天躁狠狠躁夜夜2o2o| 亚洲精品中文字幕一二三四区 | 亚洲精品国产一区二区精华液| 亚洲avbb在线观看| 高清av免费在线| 国产精品 欧美亚洲| 欧美精品av麻豆av| 精品第一国产精品| 中文字幕av电影在线播放| 国产成人免费无遮挡视频| 脱女人内裤的视频| 亚洲精品av麻豆狂野| 久久国产亚洲av麻豆专区| 国产一区二区 视频在线| 午夜免费鲁丝| 亚洲自偷自拍图片 自拍| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自偷自拍图片 自拍| 老司机亚洲免费影院| 免费在线观看视频国产中文字幕亚洲| 精品国产国语对白av| 一区二区三区乱码不卡18| 不卡av一区二区三区| 国产亚洲精品一区二区www | 搡老乐熟女国产| 搡老乐熟女国产| 国产精品国产高清国产av | 男男h啪啪无遮挡| 中文字幕av电影在线播放| 久久国产亚洲av麻豆专区| 色播在线永久视频| 精品国产国语对白av| 亚洲欧美一区二区三区久久| 咕卡用的链子| 国产精品久久久久成人av| 夜夜骑夜夜射夜夜干| 无限看片的www在线观看| 精品久久久精品久久久| 女人爽到高潮嗷嗷叫在线视频| 黑人操中国人逼视频| 搡老熟女国产l中国老女人| 久久精品亚洲精品国产色婷小说| 夜夜骑夜夜射夜夜干| 欧美 亚洲 国产 日韩一| 亚洲色图av天堂| av免费在线观看网站| 亚洲成人免费电影在线观看| 蜜桃国产av成人99| 天天躁夜夜躁狠狠躁躁| 人人妻,人人澡人人爽秒播| 午夜免费成人在线视频| 午夜福利视频在线观看免费| 国产xxxxx性猛交| 精品一品国产午夜福利视频| 考比视频在线观看| 亚洲av美国av| 19禁男女啪啪无遮挡网站| 亚洲精品粉嫩美女一区| 美女高潮到喷水免费观看| 久久精品人人爽人人爽视色| 亚洲男人天堂网一区| 欧美午夜高清在线| 人妻一区二区av| 精品乱码久久久久久99久播| a级片在线免费高清观看视频| 国产淫语在线视频| 三上悠亚av全集在线观看| 亚洲全国av大片| 欧美亚洲日本最大视频资源| 19禁男女啪啪无遮挡网站| 国产在线视频一区二区| 精品久久久久久久毛片微露脸| 亚洲性夜色夜夜综合| 精品国产一区二区三区久久久樱花| 久久久久久久精品吃奶| 免费看十八禁软件| 日本av免费视频播放| 国产精品美女特级片免费视频播放器 | 9色porny在线观看| 看免费av毛片| 午夜福利视频精品| 亚洲av片天天在线观看| 亚洲午夜理论影院| 欧美一级毛片孕妇| 欧美黄色淫秽网站| av超薄肉色丝袜交足视频| 18在线观看网站| 香蕉国产在线看| 别揉我奶头~嗯~啊~动态视频| 欧美乱妇无乱码| 麻豆国产av国片精品| 国产区一区二久久| 日本wwww免费看| 久久精品国产99精品国产亚洲性色 | 免费女性裸体啪啪无遮挡网站| 久久久国产精品麻豆| 国产欧美亚洲国产| 国产精品.久久久| 俄罗斯特黄特色一大片| 日韩欧美免费精品| 国产成人精品无人区| 亚洲av国产av综合av卡| 制服诱惑二区| 无人区码免费观看不卡 | 91九色精品人成在线观看| 色婷婷久久久亚洲欧美| 一二三四在线观看免费中文在| 侵犯人妻中文字幕一二三四区| 天堂俺去俺来也www色官网| 久久久久久久大尺度免费视频| 久久久久久久国产电影| 欧美av亚洲av综合av国产av| 80岁老熟妇乱子伦牲交| 19禁男女啪啪无遮挡网站| 18禁美女被吸乳视频| 亚洲人成电影观看| 国产97色在线日韩免费| 制服人妻中文乱码| 国产精品免费视频内射| 男女午夜视频在线观看| 波多野结衣av一区二区av| 1024香蕉在线观看| 国产av精品麻豆| 欧美亚洲 丝袜 人妻 在线| 久久久精品国产亚洲av高清涩受| 熟女少妇亚洲综合色aaa.| av国产精品久久久久影院| 午夜福利影视在线免费观看| 亚洲欧美日韩另类电影网站| 精品亚洲乱码少妇综合久久| 久久精品国产亚洲av香蕉五月 | 女人高潮潮喷娇喘18禁视频| av在线播放免费不卡| 亚洲av日韩精品久久久久久密| cao死你这个sao货| 欧美日韩一级在线毛片| 精品久久蜜臀av无| 亚洲精品久久午夜乱码| 亚洲专区字幕在线| 超碰97精品在线观看| 亚洲少妇的诱惑av| 国产成人精品久久二区二区免费| 国产老妇伦熟女老妇高清| 久久久久久人人人人人| 国产高清激情床上av| 亚洲精品国产精品久久久不卡| 精品国产乱码久久久久久男人| 超色免费av| 首页视频小说图片口味搜索| 国产精品一区二区免费欧美| 超碰成人久久| 在线亚洲精品国产二区图片欧美| 一个人免费看片子| 日本黄色日本黄色录像| 亚洲黑人精品在线| 首页视频小说图片口味搜索| 50天的宝宝边吃奶边哭怎么回事| www日本在线高清视频| 12—13女人毛片做爰片一| 精品一品国产午夜福利视频| 成年人免费黄色播放视频| 成人免费观看视频高清| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区 | 欧美激情久久久久久爽电影 | 女同久久另类99精品国产91| 新久久久久国产一级毛片| 国产一卡二卡三卡精品| 极品教师在线免费播放| 高清毛片免费观看视频网站 | 99精品在免费线老司机午夜| 国产精品一区二区精品视频观看| 乱人伦中国视频| 国产不卡av网站在线观看| 久久久久精品人妻al黑| 欧美在线一区亚洲| 午夜久久久在线观看| 在线 av 中文字幕| 丁香欧美五月| 免费观看人在逋| 国产成人av教育| 日日摸夜夜添夜夜添小说| 美女视频免费永久观看网站| 黄片小视频在线播放| 超色免费av| 我的亚洲天堂| 国产一区二区三区视频了| 中文字幕人妻丝袜一区二区| 青草久久国产| 黄片播放在线免费| 老汉色∧v一级毛片| 性高湖久久久久久久久免费观看| 久久久国产成人免费| 成人影院久久| 一级a爱视频在线免费观看| 人妻 亚洲 视频| 成人手机av| 国产日韩一区二区三区精品不卡| 黄网站色视频无遮挡免费观看| 久久av网站| 亚洲色图 男人天堂 中文字幕| 搡老乐熟女国产| 女同久久另类99精品国产91| 国产av一区二区精品久久| 精品福利永久在线观看| 亚洲国产av影院在线观看| 脱女人内裤的视频| 欧美日韩亚洲高清精品| 天堂动漫精品| 大香蕉久久成人网| 久久人人爽av亚洲精品天堂| 一边摸一边做爽爽视频免费| 亚洲少妇的诱惑av| 免费观看人在逋| 日韩大片免费观看网站| 咕卡用的链子| 乱人伦中国视频| 肉色欧美久久久久久久蜜桃| 国产成人一区二区三区免费视频网站| 一级毛片精品| 国产欧美日韩一区二区三区在线| 高清视频免费观看一区二区| 欧美亚洲 丝袜 人妻 在线| 国产精品 欧美亚洲| 国产精品一区二区精品视频观看| 久热这里只有精品99| 久久免费观看电影| 久久青草综合色| 日韩制服丝袜自拍偷拍| 亚洲精品粉嫩美女一区| 欧美日韩福利视频一区二区| 脱女人内裤的视频| 欧美 日韩 精品 国产| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 一本久久精品| 国产成人系列免费观看| 性色av乱码一区二区三区2| 老司机福利观看| 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 99国产精品一区二区蜜桃av | 超色免费av| 国产日韩一区二区三区精品不卡| 女同久久另类99精品国产91| av视频免费观看在线观看| 母亲3免费完整高清在线观看| 欧美黄色淫秽网站| 欧美日韩国产mv在线观看视频| 国产深夜福利视频在线观看| 久久精品人人爽人人爽视色| 国产国语露脸激情在线看| 精品少妇久久久久久888优播| 亚洲精品国产区一区二| 99在线人妻在线中文字幕 | 亚洲熟女精品中文字幕| 亚洲精品av麻豆狂野| 好男人电影高清在线观看| 亚洲性夜色夜夜综合| 国产精品亚洲av一区麻豆| 女性被躁到高潮视频| 亚洲精品久久成人aⅴ小说| 色精品久久人妻99蜜桃| 国产在线精品亚洲第一网站| 飞空精品影院首页| 亚洲性夜色夜夜综合| 国产成人欧美在线观看 | 精品少妇久久久久久888优播| 在线亚洲精品国产二区图片欧美| 国产成人精品在线电影| 久久久久视频综合| 国产精品98久久久久久宅男小说| 久久av网站| 三级毛片av免费| 在线十欧美十亚洲十日本专区| 午夜福利乱码中文字幕| 国产色视频综合| 丰满少妇做爰视频| 天天躁狠狠躁夜夜躁狠狠躁| 大型av网站在线播放| 在线播放国产精品三级| 在线观看www视频免费| av有码第一页| bbb黄色大片| 久久精品亚洲熟妇少妇任你| 一进一出好大好爽视频| 成年人黄色毛片网站| 一区二区三区国产精品乱码| 97在线人人人人妻| 在线看a的网站| 真人做人爱边吃奶动态| 免费少妇av软件| 久9热在线精品视频| 国产精品国产高清国产av | 丝袜美腿诱惑在线| 黄色a级毛片大全视频| 日本欧美视频一区| 日本a在线网址| av网站在线播放免费| 国产精品免费视频内射| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看 | 久久中文字幕一级| 亚洲国产中文字幕在线视频| 丝瓜视频免费看黄片| 亚洲七黄色美女视频| 正在播放国产对白刺激| 狠狠精品人妻久久久久久综合| 老司机影院毛片| 欧美成人免费av一区二区三区 | 国产精品亚洲av一区麻豆| 国产精品一区二区精品视频观看| 日韩欧美三级三区| 黄片大片在线免费观看| 中文字幕色久视频| 一区二区三区乱码不卡18| 悠悠久久av| 久久精品人人爽人人爽视色| 老司机福利观看| 建设人人有责人人尽责人人享有的| 成人免费观看视频高清| 成年女人毛片免费观看观看9 | 精品午夜福利视频在线观看一区 | 国产xxxxx性猛交| 国产色视频综合| 9191精品国产免费久久| 国产精品麻豆人妻色哟哟久久| 久久精品91无色码中文字幕| 亚洲色图综合在线观看| 日本欧美视频一区| 桃红色精品国产亚洲av| 色综合婷婷激情| avwww免费| 麻豆国产av国片精品| 国产亚洲精品一区二区www | 丝袜喷水一区| 丝袜人妻中文字幕| 热re99久久国产66热| 国产精品香港三级国产av潘金莲| 亚洲色图 男人天堂 中文字幕| 欧美乱码精品一区二区三区| 考比视频在线观看| 日本一区二区免费在线视频| 在线 av 中文字幕| 三上悠亚av全集在线观看| 精品人妻在线不人妻| 久久免费观看电影| 国产成人免费无遮挡视频| 精品一区二区三卡| 亚洲熟女精品中文字幕| 无限看片的www在线观看| 亚洲第一青青草原| e午夜精品久久久久久久| 午夜福利影视在线免费观看| 蜜桃在线观看..| 老司机午夜十八禁免费视频| 老司机在亚洲福利影院| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 成人av一区二区三区在线看| 高清黄色对白视频在线免费看| 看免费av毛片| 日韩视频在线欧美| 后天国语完整版免费观看| 超碰成人久久| 精品人妻熟女毛片av久久网站| bbb黄色大片| 一级黄色大片毛片| 国产av精品麻豆| 最近最新中文字幕大全电影3 | 亚洲久久久国产精品| 成年女人毛片免费观看观看9 | 久久久精品94久久精品| 韩国精品一区二区三区| 久久精品91无色码中文字幕| 9191精品国产免费久久| kizo精华| 狂野欧美激情性xxxx| 黄色 视频免费看| 69av精品久久久久久 | 午夜福利一区二区在线看| 菩萨蛮人人尽说江南好唐韦庄| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 免费在线观看日本一区| h视频一区二区三区| 极品人妻少妇av视频| 亚洲中文字幕日韩| 亚洲国产精品一区二区三区在线| 久久国产亚洲av麻豆专区| 国产一区二区三区综合在线观看| 18禁裸乳无遮挡动漫免费视频| 欧美日韩黄片免| 欧美在线一区亚洲| 欧美日韩国产mv在线观看视频| 免费在线观看完整版高清| 老司机午夜福利在线观看视频 | 国产精品久久久人人做人人爽| e午夜精品久久久久久久| 国产一区二区三区视频了| 人成视频在线观看免费观看| 亚洲av片天天在线观看| 高清视频免费观看一区二区| 午夜激情久久久久久久| 国产精品影院久久| 男人操女人黄网站| 午夜激情av网站| 精品国产乱码久久久久久男人| av国产精品久久久久影院| 欧美午夜高清在线| 欧美成人免费av一区二区三区 | 亚洲 欧美一区二区三区| 他把我摸到了高潮在线观看 | 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人不卡在线观看播放网| 一区在线观看完整版| 超碰97精品在线观看| 汤姆久久久久久久影院中文字幕| 午夜激情av网站| 大型av网站在线播放| 夜夜骑夜夜射夜夜干| 极品教师在线免费播放| 精品少妇内射三级| 十八禁网站网址无遮挡| 国产一区二区 视频在线| 国产精品国产av在线观看| 亚洲avbb在线观看| 日韩精品免费视频一区二区三区| 亚洲免费av在线视频| av线在线观看网站| 男人操女人黄网站| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看| 欧美激情久久久久久爽电影 | 精品久久久久久久毛片微露脸| 免费在线观看日本一区| 国产精品免费一区二区三区在线 | 日本a在线网址| 国产成人精品久久二区二区免费| 久久av网站| 国产片内射在线| 国产主播在线观看一区二区| 人妻久久中文字幕网| 亚洲第一av免费看| 久久久久久亚洲精品国产蜜桃av| 伊人久久大香线蕉亚洲五| 日本一区二区免费在线视频| 久久婷婷成人综合色麻豆| 91国产中文字幕| 黑人巨大精品欧美一区二区mp4| 大型av网站在线播放| 国产亚洲精品一区二区www | 精品少妇黑人巨大在线播放| 制服人妻中文乱码| 国产一卡二卡三卡精品| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三| 搡老乐熟女国产| 日韩欧美一区二区三区在线观看 | 少妇裸体淫交视频免费看高清 | 老司机在亚洲福利影院| 日韩大码丰满熟妇| 欧美亚洲日本最大视频资源| 99热国产这里只有精品6| 亚洲国产欧美日韩在线播放| 午夜精品久久久久久毛片777| 成人av一区二区三区在线看| 最黄视频免费看| 亚洲国产成人一精品久久久| tocl精华| 久久人人97超碰香蕉20202| 18禁黄网站禁片午夜丰满| 午夜精品国产一区二区电影| 亚洲免费av在线视频| 国产xxxxx性猛交| 搡老乐熟女国产| 国产午夜精品久久久久久| 一级毛片女人18水好多| 国产精品一区二区在线观看99| 黄色a级毛片大全视频| 欧美黑人欧美精品刺激| 亚洲欧美日韩高清在线视频 | 午夜成年电影在线免费观看| av视频免费观看在线观看| av天堂久久9| 久久久水蜜桃国产精品网| 亚洲欧美一区二区三区久久| 中文字幕人妻丝袜制服| 女性被躁到高潮视频| 亚洲国产欧美一区二区综合| 亚洲欧美激情在线| 欧美人与性动交α欧美精品济南到| 自线自在国产av| 成人国语在线视频| 国产av一区二区精品久久| 久久影院123| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 美女高潮到喷水免费观看| 欧美乱码精品一区二区三区|