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

    Learning Dual-Domain Calibration and Distance-Driven Correlation Filter:A Probabilistic Perspective for UAV Tracking

    2024-01-12 03:47:50TaiyuYanYuxinCaoGuoxiaXuXiaoranZhaoHuZhuandLizhenDeng
    Computers Materials&Continua 2023年12期

    Taiyu Yan ,Yuxin Cao ,Guoxia Xu ,Xiaoran Zhao ,Hu Zhu and Lizhen Deng

    1Jiangsu Province Key Lab on Image Processing and Image Communication,Nanjing University of Posts and Telecommunications,Nanjing,China

    2School of Computer Science,Qufu Normal University,Qufu,China

    3National Engineering Research Center of Communication and Network Technology,Nanjing University of Posts and Telecommunications,Nanjing,China

    ABSTRACT Unmanned Aerial Vehicle (UAV) tracking has been possible because of the growth of intelligent information technology in smart cities,making it simple to gather data at any time by dynamically monitoring events,people,the environment,and other aspects in the city.The traditional filter creates a model to address the boundary effect and time filter degradation issues in UAV tracking operations.But these methods ignore the loss of data integrity terms since they are overly dependent on numerous explicit previous regularization terms.In light of the aforementioned issues,this work suggests a dual-domain Jensen-Shannon divergence correlation filter(DJSCF)model address the probability-based distance measuring issue in the event of filter degradation.The two-domain weighting matrix and JS divergence constraint are combined to lessen the impact of sample imbalance and distortion.Two new tracking models that are based on the perspectives of the actual probability filter distribution and observation probability filter distribution are proposed to translate the statistical distance in the online tracking model into response fitting.The model is roughly transformed into a linear equality constraint issue in the iterative solution,which is then solved by the alternate direction multiplier method(ADMM).The usefulness and superiority of the suggested strategy have been shown by a vast number of experimental findings.

    KEYWORDS Dual-domain weighting;distance measure;correlation filter;ADMM

    1 Introduction

    Smart cities[1]are increasingly emerging as the primary outcome of urbanization and information technology integration as social intelligence advances.Six-generation (6G) technology [2] supports intelligent interconnection of human-machine objects and can generate potential application scenarios[3]such as full regional coverage by making full use of full spectrum resources.However,smart cities lack the support of spatial geographic information,and the application scenarios usually require terminal devices to provide more processing power and are more sensitive to transmission delays.

    However,current devices typically have poor processing power,and the centralized architecture of traditional cloud computing networks may result in poor transmission latency [4].With the development of the Internet of Things(IoT)[5–10],information data is extracted in a timely manner using big data technologies,and edge computing servers are deployed at the edge of the network through mobile edge computing (MEC) [11–13].This can effectively solve the problems of sensitive transmission delay and insufficient computing power of mobile terminals [14],and facilitate the provision of sufficient information for smart city construction[15–18].In MEC networks,UAVs are highly flexible and operable and can dock space networks with ground networks to collaboratively build three-dimensional networks[19].UAVs are able to dynamically adjust their deployment position in the network by tracking and adaptively responding to changes in the network environment [20].UAVs’visual tracking has been used for aerial refueling[21],aerial aircraft tracking[22],and target tracking[23,24].Although there have been many successful research results in visual tracking,it still remains a difficult point to achieve high robustness and accuracy visual tracking with the influence of environmental factors such as occlusion,target in-plane/out-plane rotation,fast motion,blurring,lighting,and deformation,etc.

    Recently,discriminative correlation filters (DCFs) [25] based tracking methods have received widespread attention [26–28].However,these methods have caused the boundary effect problem with adverse effects on tracking performance while using the circular matrix which is calculated in the frequency domain to improve the model efficiency.Therefore,if it is affected by aberrance in challenging scenes.Compared with ordinary target tracking,the targets in the UAV domain are often very small and the existing DCF methods are definitely sensitive to handle the small targets in this case [29,30].The existing DCF methods estimate the center position of the target,and the tracking model will continue to degenerate in the sequences[28].Thus,if the target is small and similar to the background,UAV tracking is more prone to template drift and model degradation during the learning process.

    A variety of visual tracking methods based on DCF have been proposed,and most of the successful strategies are applied to solve the above model degradation problem.For example,spatially regularized discriminative correlation filters (SRDCF) [31] and spatial-temporal regularized correlation filters (STRCF) [32] used spatial regularization and space-time regularization to solve boundary effects.In addition,STRCF mitigates the degradation of time filtering by introducing penalty terms in the context of SRDCF,which makes the tracker subject to more noise interference and reduces tracking reliability.And adaptive spatially-regularized correlation filters (ASRCF) [33]used the adaptive spatial regularization method to further improve the estimation ability of object perception and obtain more reliable filter coefficients.Many DCFs based trackers have suppressed the boundary effects by expanding the search area to obtain better tracking effects.The background-aware correlation filter(BACF)[34]used a cropping matrix to extract patches densely from the background and expand the search area at a lower computational cost.Group feature selection method for discriminative correlation filters(GFS-DCF)[35]advocates channel selection in DCF-based tracking and used an effective low-rank approximation [36].The filter constraint here can smoothly realize the temporal information and adaptively integrate historical information to achieve group feature selection across channels and spatial dimensions.In order to get accurate tracking,a large margin object tracking method with circulant feature maps(LMCF)[37]and an attentional correlation filter network(ACFN)[38],tried to suppress the possible abnormalities through some measures.To solve the aberrance problem of response maps,the aberrance repressed correlation filter(ARCF)[39]added a regular term on the basis of BACF which can suppress the mutation of the response maps and used Euclidean distance to define the direct response approximation between the two response maps.

    In addition,in order to deal with the problem of boundary discontinuity and sample contamination caused by the cosine window,Li et al.[40] proposed a spatial regularization method to remove the cosine window and further introduced binary and Gaussian shape mask functions to eliminate the boundary problem.

    For UAV tracking,there are also some special strategies in the DCF framework.Semantic-aware real-time correlation tracking framework (SARCT) [41] introduced a semantic-aware DCF-based framework for UAV tracking by combing the semantic segmentation module and sharing features for the object detection module and semantic segmentation.Also,the part-based mechanism was also combined with BACF in [42] and the Bayesian inference approach was proposed to achieve a coarse-to-fine strategy.AutoTrack[43]proposed a new method for online automatic adaptive learning of spatiotemporal regularization terms,which uses spatial local response map variants as spatial regularization,so that DCF focuses on the learning of the trusted part of the object,and determines the update processing of the filter according to the global response map variants.Overall,these existing methods share some commonalities.In particular,they often heavily rely on a presumed fixed spatial or temporal filter learning and fixed response map approximation.

    In this paper,UAV tracking is sufficiently far away from the object in traditional tracking scenarios that the small target is fundamentally limited by less knowledge used in online learning and a more unbalanced sample problem.While many existing works focus on the sparsity[44]in spatial smooth or temporal smooth,the direct coupling of the several regularization terms used in several works like[34,41–43]is not a good choice,and the online filter learning will be obfuscated by redundancy.Therefore,it is interesting to ask the problem concerned in this work:what can we measure in online learning:probabilistic response fitting or direct response approximation?

    From the perspective of actual and observed probability filter distribution in small targets for UAV tracking,the probabilistic response fitting in temporal smooth is more robust than the direct response approximation.Although Euclidean distance is useful,it is limited by the largest-scaled feature that would dominate the others.It is hard to judge the “best” similarity measure or a best-performing measure here.In general,we are the first work to transfer the statistical distance to correlation response fitting in the online tracking model.It also could shed light on the performance and behavior of measures in deep learning works like [45].Although the insights of this work are the most similar to our work,the main truth is that the distribution assumption is different.

    Furthermore,a novel unsupervised method called maximally divergent intervals (MDI) [46]searched for continuous intervals of time and region in a space containing anomalous events instead of analyzing the data point by point and attempted to use KL divergence and unbiased KL divergence to achieve better detection results.Not only that,but MDI also introduced JS divergence to demonstrate the superiority of the model and improve its scalability of the model.Similarly,Jiang et al.in[47]also used KL divergence to search for abnormal blocks in higher-order tensors.At the beginning of signal processing,the KL and JS divergence was used to fit the matching degree of the two distribution signals.Inspired by this,we consider the similarity search between the response maps of the filters as the curve fitting problem between the adjacent frame filters.It should be noted that the original Euclidean distance is not applicable because it is prone to filter under-fitting learning.Actually,the under-fitting problem had recently been mentioned in feature embedding problems like [48].While it seems to be rather complicated to estimate a probabilistic response fitting by the above KL or JS procedure,the main benefit of estimating the mean and the covariance lies in the inherent regulation properties in our correlation filter response term.

    While the KL divergence and JS divergence can achieve a good fitting effect.In addition,since the fitting object in the KL divergence is the response map of the adjacent frame filter to the image block[39],the KL divergence term can ensure that the filter has the highest similarity in the response map in the time domain.Therefore,our model uses a temporal-response calibrated filter based on the KL divergence.Since the KL divergence does not measure the distance between the two distributions in space,the KL divergence can be considered as a measure of the information loss of the response maps in the filter in our algorithm.Then we can find the target with the highest response value in the time domain according to the principle of function loss minimization,which is considered our tracking target.Experiments have shown that trackers based on temporal-response calibrated filters have greatly improved accuracy.

    Based on the above-mentioned analysis,we propose a new model named AKLCF.Specifically,we use KL divergence fitting to mitigate the under-fitting of the Euclidean distance.In addition,to ensure the validity of the model,we add model confidence terms to avoid the variations of the model that occurred during the optimization process,which greatly improves the stability of the model.The algorithm flow chart of AKLCF is shown in Fig.1a.Specifically,the multi-channel features are extracted from the previous frame of an image.We use KL divergence to calculate the similarity between response maps,and then build the learning model of AKLCF.At the same time,we further use the JS divergence to investigate the impact of the similarity measurement method on the tracking effect of related filters and propose the AJSCF model.The algorithm flow chart of AJSCF is shown in Fig.1b.Notice that AJSCF is different from AKLCF in divergence.Compared with KL divergence,JS divergence solves the problem of asymmetry of KL divergence and avoids the“gradient disappearance”problem in the process of KL divergence fitting,thereby improving the tracking effect of the model.

    Figure 1:The algorithm flowcharts of aberrance Kullback-Leibler(KL)divergence correlation filter(AKLCF) and aberrance Jensen-Shannon (JS) divergence correlation filter (AJSCF).The structure of the two flowcharts is the same and the main difference is that the former uses KL divergence to achieve similarity fitting,while the latter uses JS divergence

    The UAV tracking will be further influenced by the imbalance of positive and negative samples during training DCF.Although kernelized correlation filter(KCF)[49]gave continuous weighting to samples with different offsets by continuous labeling,and could effectively detect target samples of interest to a certain extent and improve the accuracy of visual tracking,it does not solve the adverse effects of imbalance of positive and negative samples on visual tracking.Therefore,inspired by the collaborative distribution alignment (CDA) model [50],we propose a dual-domain Jensen-Shannon divergence correlation filter(DJSCF)model,which introduces dual-domain weighting into the data items of the loss function,uses a weighting method to suppress the influence of the imbalance of positive and negative samples and combines JS divergence to improve the accuracy of target tracking.The algorithm flow chart of DJSCF is shown in Fig.2.DJSCF adds the sorting and selection of features on the basis of AJSCF,which is realized by dual-domain weighting,and effectively alleviates the problem of imbalance between samples.

    Figure 2:The flow chart of dual-domain Jensen-Shannon divergence correlation filter(DJSCF)

    The main contributions of this work are shown as follows:

    ? To address the under-fitting of the response approximation in the Euclidean distance domain,we suggest the AKLCF model,which makes use of KL divergence fitting.The distortion brought on by noise can be efficiently reduced by the filter based on temporal response calibration.To further assure the model’s validity,we also include model confidence terms.

    ? We introduce JS divergence to propose the AJSCF model,which validates the impact of measurement methods on tracking accuracy and enhances the richness and verifiability of article content,in order to address the asymmetry of KL divergence and the disappearance of gradient.

    ? In order to address the issue of sample imbalance,we present a DJSCF model with JS divergence fitting that incorporates two-domain weighting due to two diagonally weighted matrices in the sparse data item.Several tests were carried out to demonstrate the excellence of our DJSCF.

    2 Preliminaries

    2.1 Distance Measurement

    Target tracking mainly uses the tracking algorithm to select the one with the highest similarity to the target object as the target of the next frame,and then obtains the motion trajectory of the target object in the entire video sequence.Therefore,how accurately represent the similarity between the candidate block and the target block is an important factor that determines the accuracy of target tracking.As the research of target tracking and the difficulty of tracking tasks continue to increase,the distance similarity measurement methods used by tracking algorithms are also continuously improved.In image processing and analysis theory,common similarity measurement methods include Euclidean distance,block distance,checkerboard distance,weighted distance,Bart Charlie coefficient,Harsdorf distance,etc.Among them,Euclidean distance is the most applied and simplest.

    The most commonly used mathematical expression of the Euclidean distance is the norm,which is mainly divided into four norms.Among them,the target tracking algorithm often uses theL1norm and theL2norm to extract sample sparsity characteristics and measure similarity between samples.With the complex changes in the target tracking scene,the Euclidean distance cannot meet the performance requirements of target tracking.Researchers try to apply similarity measurement algorithms in the signal domain to target trackings,such as KL divergence and JS divergence.

    When the probability of an event is p(x),its information quantity is -log(p(x)),and the expectation of information quantity is entropy.Relative entropy is also known as KL divergence.If we have two separate probability distributions P(x) and Q(x) for the same random variable x,we can use KL divergence to measure the difference between the two distributions.In machine learning,P is often used to represent the true distribution of the sample,and Q is used to represent the distribution predicted by the model.Then the KL divergence can calculate the difference between the two distributions,that is the loss value.

    From the KL divergence formula,we can see that the closer the distribution of Q is to P (the more fitting the Q distribution is to P),the smaller the divergence value is,that is,the smaller the loss value is.Because the logarithmic function is convex,the value of KL divergence is non-negative.KL divergence is sometimes called KL distance,but it does not satisfy the property of distance.Because KL divergence is not symmetric and KL divergence does not satisfy triangle inequality.In machine learning,KL divergence can be used to evaluate the gap between label and predictions.In the process of optimization,we only need to pay attention to the cross entropy.

    In probability statistics,JS divergence has the same ability to measure the similarity of two probability distributions as KL divergence.Its calculation method is based on KL divergence and inherits the non-negative nature of KL divergence.However,there is one important difference: JS divergence has symmetry.The formula of JS divergence is as follows:set two probability distributions as P and Q,M=0.5×(P+Q).

    If the two distributions P and Q are far apart and completely do not overlap,then the KL divergence value is meaningless,while the JS divergence value is a constant log 2.This has a great impact on the deep learning algorithm based on gradient descent,which means that the gradient is 0,that is,the gradient disappears.

    2.2 Dual-Domain Weighting

    Due to the differences in background,light,sensor resolution,and other aspects,the data collected by UAVs in different scenes may be in different data domains.There are different numbers of positive samples and negative samples in each domain,and the number of negative samples(background)is far greater than the number of positive samples(target)in single-target tracking.The problem of sample imbalance will lead to the over-fitting of a large sample proportion,and the prediction will tend to classify with more samples.Therefore,there is a strong need to align the number of samples between different domains and train the data to find samples suitable for the target domain.Inspired by the trilateral weighted sparse coding scheme(TWSC)in[51]and the dual-domain collaboration in[50],we use the method of two-field weighting to increase the weight of samples,that is,to add a larger weight to the samples with a smaller number of samples,and to add a smaller weight to the samples with a larger number of samples.

    In the model,we introduce two diagonal matrices,waandwb,as weighted matrices.wais a block diagonal matrix,with each block having the same diagonal element to describe sample attributes in different channels.wbis used to describe the sample variance in the corresponding patch.The formula can be expressed as:

    wherexis the positive sample,yis the negative sample.Our method also introduces the idea of dualdomain weighting,and the supporting materials show the specific processing operations ofwaandwb.In other words,two weight matrices are introduced into the data fidelity term of the sparse coding frame to adaptively represent the statistics of positive and negative samples in each patch of the multichannel,which can make better use of the sparse prior of the natural image.

    3 The Dual-Domain Collaboration and Distances Correlation Filters

    In this section,we first introduce the AKLCF and explain the role of KL divergence in the model tracking process.We further make the derivation of the optimization process in detail.Then,we introduce and derive DJSCF in detail,and introduce how dual-domain weighting and JS divergence can jointly improve the accuracy and robustness of target tracking.Considering that our model is based on the similarity of response maps,Fig.2 shows the response effect graphs of the three methods.The trackers learn both positive samples (red sample) of objects and negative samples (green sample) in the background.

    3.1 Aberrance Kullback-Leibler Divergence Correlation Filter(AKLCF)Model

    Then,we replace the third term of Eq.(4)with the functionR(Mk|Mk-1;xk-1)to obtain the loss function of AKLCF:

    As discussed above,in the formula with Eq.(5) add of is replaced by the Kullback-Leibler divergence term.In addition,we add a model confidence term to ensure the reliability of the model.

    Then we improve the loss functions by minimizing the following unconstraint objective:

    To better reflect the effect of the KL divergence term,we first analyze it in detail to obtain the specific expression of the KL divergence term.Since the KL divergence is a measure of the asymmetry of the difference between the two probability distributions.In our model,the KL divergence is used to measure the difference between the convolution matrix of the current frame and the previous frame,that is:

    For convenience,we assume that the value vectors of filter and image blocks have multivariate Gaussian distributions,then the convolution matrix also has a Gaussian distribution,i.e.,M~N(M0,φ),whereM0is the mean value vector of the convolution matrix distribution andφrepresents the covariance matrix of distribution.In the previous section,the functionMi=Bxd?f dindicates the convolution prediction of the class labelf iand the samplexiini-thdimension with the given multivariate Gaussian distribution.For better understanding,we regard each valuefiin thei-thdimension as the models knowledge about the corresponding featurexiand regard the diagonal entry of covariance matrixφi,ias the confidence of featureMi.The smaller the value ofφi,iis,the more confident the learner is in the mean weight valueMi.Expect diagonal values,other covariance termsφi,jcan be took as the correlations between two weight valueMiandMjfor image featureandIn this way,we model parameter confidence for each element of along all channels with a diagonal Gaussian distribution.Thus,we construct the kullback-leibler divergence term based on the divergence between empirical distribution and probability distribution.As a result,we get theRklas following:

    We substitute the specific expression of the model confidence term and Eq.(8)into Eq.(6)to get:

    In general,the first term of the formula is data fidelity,the second is the space penalty,and the third and fourth terms are our improved contents,namely the model confidence term and the KL divergence term.

    3.2 The Dual-Domain Jensen-Shannon Divergence Correlation Filter(DJSCF)Model

    DJSCF adds dual-domain weighted matrices and JS divergence in the response map which need to be considered in the optimization process,therefore the construction of the loss function in DJSCF tracker is described in detail below:

    Similar to AKLCF,we use JS divergence and model confidence to improve the third term of the objective function,and get the objective function as follows:

    It equals to:

    The AJSCF model is a version that only replaces the KL divergence of AKLCF with JS divergence.The comparison between AJSCF and AKLCF shows the superiority of JS divergence in similarity fitting.At the same time,the latter DJSCF is obtained by adding dual-domain weighting based on AJSCF.And comparing AJSCF with DJSCF can verify the effect of dual-domain weighting on tracking performance.

    3.3 ADMM Decomposition of AKLCF and DJSCF

    AKLCF:The model in Eq.(9) is a convex problem and can be minimized with the ADMM algorithm to obtain the global optimal solution.Since Eq.(9) has a convolution calculation in the time domain,Pasival’s theorem is used to facilitate the calculation by converting the problem into the frequency domain.Then we get the objective function in the frequency domain as follows:

    According to Eq.(8)and Pasival’s theorem,the KL divergence term in the frequency domain is expressed as follows:

    Substitute Eq.(17)into Eq.(16),and using Pasival’s theorem,then we get the specific optimization function.

    whereμis a penalty factor.In the current frame,since the response map in the previous frame is already generated,can be treated as a constant signal,which can simplify the further calculation.The ADMM algorithm is used to divide the optimization function into several sub-problems as follows:

    DJSCF:Similar to AKLCF,according to Pasival’s theorem,we get:

    We get the optimization formulation:

    According to the principle of JS divergence and Eq.(8),we get:

    Then the optimization solution of DJSCF is implemented by ADMM and is divided into the following subproblems:

    3.4 Solution of Sub-Problems of Two Models

    From Eqs.(19) and (23),it is easy to find that except for the subproblemthe solutions of other subproblems in the two models are the same.Thus,we talk about the solution integration of the two models and only discuss them separately on subproblem

    In AKLCF model,the augmented Lagrangian form is expressed as follows:

    Using the ADMM algorithm,two sub-problems can be obtained,namely:

    Each smaller problem can be efficiently calculated and the solution is presented below:

    The solutions of the weighted matricesWaorWbare as follows:

    where the parameterτis calculated as 0.6.To make better use of the constraints and selection of the weighting matrix,we use the sorting method to sort and assign the matrix elements calculated by Eq.(38),thereby getting new matrices as the weighted matrices for the next iteration.The specific processing steps are as follows.

    The core idea of matrix update is to set non-diagonal elements to 0 and reorder the diagonal elements to assign values [52,53].The specific method is to arrange the elements on the diagonal from small to large and replace the original value with the sort position number of each element.For example,if thej-thelement of the diagonal is ranked at theN-thposition,the value of the element becomesW(j)=N,and then use the following formula to obtain thej-thnew element:

    whereαis the weight parameter and the value ofiis a or b.Finally,we convert all the elements on the diagonal to get a new matrixWaorWb.

    Update of Lagrangian Parameter:

    The Lagrangian parameter is updated according to the following equation:

    where the subscriptiandi+1 denote thei-thand thei+1-thiteration,respectively.

    Updating of appearance model:

    The appearance model updating takes an important influence on visual tracking.The appearance of tracked objects often changes in the tracking process under the effect of posture,scale,occlusion,and rotation.Therefore,we apply an online update strategy to update the appearance model and improve the robustness of our method.The appearance model at frame t can be updated as follows:

    whereηis the learning rate of the appearance model.

    3.5 Complexity Analysis

    In each iterative calculation of subproblem f,the transformation of Fast Fourier Transform(FFT)and Inverse Fast Fourier Transform(IFFT)needs to be performed,then the computational complexity is O(DMNlog(MN),and if the number of iterations isK,the total computational complexity of the model isO(KDMN(log(MN)+3)).Considering the DFT and inverse DFT into account,the whole complexity of solvingfis O(DMNlog(MN) to take.As forg,the computational complexity of “KL” isO(DMN),and then the computational complexity ofgisO(DMN),which is same asφandξ.Hence,if the maximum number of iterations isNI,the whole cost of the AKLCF model isO(DMNlog(MN)+3(DMN)Ni).In addition,the DJSCF is improved based on the AKLCF by adding the dual-domain constraint and the JS divergence,thus its robustness and performance of tracks are better than the AKLCF with higher computational complexity.As for the iterative calculation,DJSCF mainly becomes more complicated to solveg.The JS divergence is equivalent to calculating the average of the two KL divergences,thus the computational complexity ofgis alsoO(3DMN),and the calculation amount ofWisO(DMN).The calculation offandξis consistent with that of AKLCF,thus the final calculation complexity of DJSCF reachesO(DMNlog(MN)+6(DMN)Ni).The fps of DJSCF,AJSCF,and AKLCF are 2.6149,4.9097,and 5.9699,respectively in UAV123.

    4 Experiments and Results

    To demonstrate the superiority and effectiveness of our proposed models,we compare them with several state-of-the-art trackers with the HOG and deep features.We chose the 19-layer deep convolutional neural network VGG19[54](Visual Geometry Group)to extract the deep features in our experiments.In order to better explore the robust performance of the models,we conduct comparative experiments on different data sets.The value ofαis set to 0.71,the ADMM iteration is set to 5 and the learning rateηis 0.0192.All experiments of all 21 trackers are carried out by MATLAB R2017a on a computer with an i7-8700 K processor(3.7 GHz),48 GB RAM,and NVIDIA Quadro P2000 GPU.

    4.1 Experiment Datasets and Baseline Methods

    We evaluated the performance of AKLCF,AJSCF,DJSCF,and other trackers on four benchmark datasets,including UAV123@10fps [55],DTB70 [29],UAVDT-S [30],UAVDT-M [30] and VisDrone2018-test-dev[56]datasets.UAV123 contains 123 video sequences,which is the most commonly used and most comprehensive data set for UAV tracking.DTB70 contains 70 video sequences.UAVDT-S and UAVDT-M are made up of 50 challenging videos respectively.It is worth noting that the UAVDT-M data set is used to detect the performance of multi-target tracking,and we modified it to make it also able to judge the performance of our single-target tracker.The VisDrone2018-test-dev dataset has 35 challenging videos.All evaluation criteria are according to the original protocol defined in three benchmarks,respectively[29,30,55].We measured the performance of the tracking methods with one-pass evaluation(OPE)[57].The evaluation method of OPE makes use of marking the real target frame in the first frame of the video sequence,no longer marks the target in the next frame number of images,and carries out a complete tracking of the entire video sequence.The common indicators of OPE are success plots and precision plots.The average pixel error averages the detection error of all frames.If the tracking fails due to accidental factors,it will have a certain impact on the result of this indicator,and the size of the impact is related to the number of frames of accidental factors.When the overlap rate of a frame is greater than the set threshold,take the horizontal axis as the threshold of the overlap threshold and the vertical axis as the success rate to obtain the success plots.The precesion plots are similar to the success plots.The horizontal axis is the threshold,and the vertical axis is the ratio of the number of frames below a certain threshold to the total number of frames.The threshold setting condition is the distance between the detection center and the real target center.When the distance is less than a certain threshold,it is considered to meet the precision.The results are compared with 15 state-of-the-art trackers with both HOG feature-based trackers(i.e.,ARCF-HC[39],LADCF[52],STRCF[32],SRDCF[31],DSST[58],SAMF[59],KCF[49],GFSDCF[35],AutoTrack[43],ECO-HC(with gray-scale)[60],Histograms of Oriented Gradients(HOG-LR)[61]) and deep-based trackers (i.e.,ECO [60],CCOT [62],ASRCF [33],MDNet [63],ADNet [64],CFNet[65],CREST[66],MCPF[67],SiamFC[68],and HDT[69]).

    4.2 Quantitative Analysis of Tracking Dataset

    First of all,we compared the overall performance of AKLCF,AJSCF,and DJSCF with other state-of-the-art hand-crafted trackers on UAV123@10fps[55],DTB70[29],UAVDT-S[30],UAVDTM[30]and VisDrone2018-test-dev[56]datasets.The comparison results are shown in Fig.3.From the resulting graph,it can be seen that each comparison method performs best on the dataset UAVDTM.Due to the limitations of feature extraction ability and computational process,traditional HOGbased methods perform slightly worse than correlation filters and deep learning-based methods.A class of methods based on discriminant correlation filters,combined with spatial adaptation or structured support vector machines,have good computational efficiency and positioning performance and perform well on various datasets.However,the method based on deep learning pays more attention to accuracy and therefore has greater computational complexity.On five different datasets,our method has an accuracy rate that exceeds the lower-performing method by nearly 25%.Even in this individual case on the UAVDT-M dataset,our method is only 0.1%less accurate than the firstplace AutoTrack.Therefore,experimental results comparing multiple datasets show that our method has strong competitiveness in tracking target detection.

    4.3 Quantitative Analysis of Variants and Deep Feature Configuration

    We carried out comparative studies in several representative scenarios,such as background clutter,illumination changes,partial occlusion,size changes,and viewpoint shifts,to reflect the tracking performance of DJSCF in various contexts.Fig.4 also displays the attribute-based performance of each tracker,with each group representing a distinct attribute.In general,a tracker that uses deep learning outperforms other trackers by a large margin.We discovered that all trackers had the lowest scores when it came to background clutter problems,additionally,some trackers struggle with scale change characteristics.Since particle filter frameworks can more effectively explore state space,we can find that particle filter-based approaches(i.e.,HOG-LR)are superior to local sliding window-based methods(i.e.,DSST and KCF)for qualities relevant to motion models.

    Figure 3:Precision and success plots of proposed methods as well as other HOG feature-based trackers on(a)UAV123@10fps(first column),(b)DTB70(second column),(c)UAVDT-S(third column),(d)UAVDT-M (fourth column) and (e) VisDrone2018-test-dev (fifth column).Precision and AUC are marked in the precision plots and success plots,respectively

    Figure 4:Attributes of proposed methods as well as other HOG feature-based trackers including(a)background clutter (first column),(b) illumination variation (second column),(c) partial occlusion(third column),(d)scale variation(fourth column),and(e)viewpoint change(fifth column).Precision and area under the curve(AUC)are marked in the precision plots and success plots,respectively

    To further accurately assess the performance of our AKLCF,AJSCF,and DJSCF,we compared the trackers on the UAVDT-S dataset with deep features.DJSCF,AJSCF,and AKLCF took the top three positions on both the precision and success plots,as shown in Fig.5.Due to the introduction of a weighting matrix to balance the number of positive and negative samples in DJSCF,the effect is better than AJSCF,with values greater than AJSCF by about 1.5%.In terms of precision plots,DJSCF,AJSCF,and AKLCF outperform ECO(0.702),which has good performance on numerous datasets,by 8.83%,6.27%,and 4.84%,and MDNet(0.725),which was rated fourth.The DJSCF score is significantly greater than that of other trackers when it comes to success plots.In comparison to MDNet,which is in fourth position,DJSCF is 3.45% higher,while ECO,which is in fifth place,is 6.43%higher.

    Figure 5:Comparison between DJSCF,AJSCF,and AKLCF trackers and different state-of-the-art deep-based trackers.The value of average precision and average success rate is calculated by OPE results from dataset UAVDT-S

    To better demonstrate the superiority of our model in tracking performance,we selected 4 sets of video sequences in the DTB70 dataset to show the performance of the proposed method.As shown in Fig.6,the different methods in each frame of the image The target tracking results are represented by different colored frames,which are ARCF-HC (red),LADCF (green),AKLCF (blue),DJSCF(black),and STRCF(pink).Because the target of UAV tracking is small,it is difficult to maintain a good tracking effect.For example,in Figs.6a,6d,the LADCF has a serious tracking offset,resulting in tracking failure,and STRCF also has the problem of expanding the tracking area in the Bike video.In short,our methods have always maintained a good tracking effect,and the tracking frame always surrounds the tracking target.

    4.4 Ablation Analysis

    To study the influence of parameter setting on the tracking effect,we carried out a corrosion analysis on the main parameters of the weighting matrix and obtained the parameter setting with the best tracking effect.Because there are two weighting matrices,namely two variable parametersai,i=1when calculating the weighting matrixWa,andi=2when the weighting matrix isWb.We use the control variable method to discuss the tracking performance of our method under different parameters.As shown in Fig.7,we conducted two experiments on the UAV123 dataset.In Fig.7a,we first definea2=0.018and then take differenta1values.In Fig.7b,we definea1=0.02,and then take different values ofa2.As can be seen from Fig.7.(a)whena2=0.018,the tracking effect ofa1=0.02is the best.It can be seen from Fig.7.(b)that whena1=0.02,a2=0.018has the best tracking effect.Therefore,these two values are taken as the final tracking results of the DJSCF algorithm.

    4.5 Robustness Analysis

    To reflect the superior robustness of DJSCF,we conducted comparative experiments between DJSCF,ARCF,and AKLCF based on the motorbike video sequence in the DTB70 dataset.As shown in Fig.8,the robustness of the tracker is mainly reflected by the overlap rate of the target frame estimated by the model and the initial target frame in different frames of images and the overlap value shown on the ordinate.Among them,the blue line represents the ARCF tracker,and the red and green line represents the AKLCF and DJSCF,respectively.It can be seen from Fig.8 that the overlap rate of DJSCF in different frames is higher than the other two,especially between 31 and 51 frames.This shows that although ARCF or AKLCF are affected by a complex environment and reduce the tracking effect,DJSCF can still maintain a high tracking accuracy and demonstrate good robustness.

    Figure 6:Comparison of tracking performance on the DTB70 dataset

    Figure 7:Comparison of tracking results based on different parameters

    Figure 8:Comparison of robustness on the motorbike video sequence in the DTB70 dataset

    5 Conclusion

    To solve the distortion problem of correlation filter caused by boundary effect and background noise,this paper proposes a distorted KL dispersion correlation filter model by introducing KL dispersion,which improves the Euclidean distance of the original model and reduces the tracking offset caused by distortion caused by noise in the filter.In addition,a filter model of JS divergence correlation is proposed based on the two-domain weighted constraint.A large number of experimental results show that the proposed method can improve the accuracy of UAV video target tracking with less speed loss.However,due to the complexity of the DJSCF algorithm,the real-time requirements will be limited.In future work,adding convolution or depth network to the model may further improve the precision and success rate.

    Acknowledgement:Resources and computing environment provided by Nanjing University of Posts and Telecommunications,Nanjing,China,thanks to the support of the school.

    Funding Statement:This work is supported by the National Natural Science Foundation of China under Grant 62072256 and Natural Science Foundation of Nanjing University of Posts and Telecommunications(Grant Nos.NY221057 and NY220003).

    Author Contributions:The authors confirm contribution to the paper as follows: study conception and design: Taiyu Yan,Guoxia Xu and Xiaoran Zhao;data collection: Yuxin Cao and Xiaoran Zhao;analysis and interpretation of results: Taiyu Yan,Guoxia Xu and Hu Zhu;draft manuscript preparation:Taiyu Yan,Yuxin Cao and Lizhen Deng.All authors reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:All code,dataset images,and other accessible materials in this article can be obtained through the email address of the corresponding author.

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

    国产乱人视频| 午夜精品久久久久久毛片777| 久99久视频精品免费| 日本三级黄在线观看| 国产午夜精品论理片| 久久精品91蜜桃| 少妇丰满av| 午夜福利在线观看免费完整高清在 | 国产不卡一卡二| 亚洲av日韩精品久久久久久密| 久久九九热精品免费| www.精华液| 黑人操中国人逼视频| 久久久久久九九精品二区国产| 91麻豆精品激情在线观看国产| 国内毛片毛片毛片毛片毛片| 99精品久久久久人妻精品| 一本精品99久久精品77| 国语自产精品视频在线第100页| 精品久久久久久,| 最近视频中文字幕2019在线8| 午夜精品在线福利| 一级毛片精品| h日本视频在线播放| 成人三级黄色视频| 国产精品女同一区二区软件 | 999久久久国产精品视频| 国产精品一及| 美女cb高潮喷水在线观看 | 99精品欧美一区二区三区四区| а√天堂www在线а√下载| 一级作爱视频免费观看| 精品国产乱子伦一区二区三区| 中文亚洲av片在线观看爽| 久久精品aⅴ一区二区三区四区| av视频在线观看入口| 成年女人毛片免费观看观看9| 韩国av一区二区三区四区| 男人舔奶头视频| 国产三级黄色录像| 一级作爱视频免费观看| 国产激情偷乱视频一区二区| 欧美日韩黄片免| 男女午夜视频在线观看| 亚洲中文字幕日韩| 日韩欧美国产在线观看| 99热这里只有是精品50| 亚洲精品一卡2卡三卡4卡5卡| 日本成人三级电影网站| 免费在线观看影片大全网站| 后天国语完整版免费观看| 免费在线观看视频国产中文字幕亚洲| 91在线观看av| av福利片在线观看| 久久久久久久久中文| 韩国av一区二区三区四区| 岛国在线观看网站| 欧美日韩瑟瑟在线播放| 9191精品国产免费久久| 亚洲av中文字字幕乱码综合| 特级一级黄色大片| 999精品在线视频| 国产精品一及| 可以在线观看的亚洲视频| 日日夜夜操网爽| 日韩欧美 国产精品| 国产成人福利小说| 国产真实乱freesex| 一级毛片女人18水好多| 男人舔女人的私密视频| 51午夜福利影视在线观看| 99久久综合精品五月天人人| 又黄又爽又免费观看的视频| 成人午夜高清在线视频| 最新在线观看一区二区三区| 亚洲最大成人中文| 丰满人妻一区二区三区视频av | 国产精品自产拍在线观看55亚洲| 又粗又爽又猛毛片免费看| 免费高清视频大片| 免费在线观看亚洲国产| 欧美日本视频| 最新中文字幕久久久久 | 好看av亚洲va欧美ⅴa在| 国产男靠女视频免费网站| 欧美性猛交╳xxx乱大交人| 日韩有码中文字幕| 无限看片的www在线观看| 亚洲美女黄片视频| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式 | 久久久水蜜桃国产精品网| 在线观看免费午夜福利视频| 国产精品99久久久久久久久| 香蕉av资源在线| 嫁个100分男人电影在线观看| 国产在线精品亚洲第一网站| 午夜福利在线在线| 欧美激情在线99| 高清在线国产一区| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 亚洲人成伊人成综合网2020| 亚洲自偷自拍图片 自拍| 国产一区在线观看成人免费| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 18禁观看日本| av视频在线观看入口| 日韩欧美在线二视频| 日本成人三级电影网站| 亚洲av电影在线进入| 成人鲁丝片一二三区免费| 欧美乱色亚洲激情| 国产精品av视频在线免费观看| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| av中文乱码字幕在线| 观看美女的网站| 99热这里只有精品一区 | 51午夜福利影视在线观看| 法律面前人人平等表现在哪些方面| 在线观看66精品国产| 国产私拍福利视频在线观看| 在线免费观看不下载黄p国产 | 一边摸一边抽搐一进一小说| 午夜精品一区二区三区免费看| 12—13女人毛片做爰片一| 久久久久久国产a免费观看| 国产97色在线日韩免费| 小说图片视频综合网站| 噜噜噜噜噜久久久久久91| 99久国产av精品| 国产精品1区2区在线观看.| 1024香蕉在线观看| 亚洲人成伊人成综合网2020| 久久久久国产精品人妻aⅴ院| 久久婷婷人人爽人人干人人爱| 男插女下体视频免费在线播放| 国产伦精品一区二区三区四那| 人人妻人人澡欧美一区二区| 免费人成视频x8x8入口观看| 天堂av国产一区二区熟女人妻| 色在线成人网| 真人一进一出gif抽搐免费| 99国产精品99久久久久| 久久这里只有精品中国| 国产麻豆成人av免费视频| 国产免费男女视频| 国产高清视频在线观看网站| 成年女人毛片免费观看观看9| av天堂中文字幕网| 亚洲精品美女久久av网站| 久久久国产成人免费| 久久久精品大字幕| 亚洲在线观看片| 国产精品99久久久久久久久| 亚洲色图av天堂| 久久精品91蜜桃| 久久亚洲精品不卡| 亚洲无线在线观看| 欧美成人一区二区免费高清观看 | 熟妇人妻久久中文字幕3abv| 日本成人三级电影网站| 悠悠久久av| 巨乳人妻的诱惑在线观看| 欧美黄色片欧美黄色片| 国产一区二区在线av高清观看| 成在线人永久免费视频| 欧美成人一区二区免费高清观看 | 97碰自拍视频| 国产精品久久久人人做人人爽| 欧美一区二区国产精品久久精品| 亚洲精品中文字幕一二三四区| 中文在线观看免费www的网站| 国产乱人视频| 国产精品 国内视频| 国产激情欧美一区二区| 一区福利在线观看| 99在线视频只有这里精品首页| 91九色精品人成在线观看| 热99在线观看视频| 欧美日韩一级在线毛片| 手机成人av网站| 亚洲国产欧洲综合997久久,| 亚洲av电影不卡..在线观看| 97碰自拍视频| 给我免费播放毛片高清在线观看| 岛国在线观看网站| 一个人免费在线观看的高清视频| svipshipincom国产片| 12—13女人毛片做爰片一| 亚洲av第一区精品v没综合| 中文在线观看免费www的网站| 日韩欧美国产一区二区入口| 国产亚洲av嫩草精品影院| 韩国av一区二区三区四区| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 手机成人av网站| 国产一区二区激情短视频| 免费看日本二区| 日本黄色视频三级网站网址| 亚洲av熟女| 国产在线精品亚洲第一网站| 亚洲国产看品久久| 国产麻豆成人av免费视频| 网址你懂的国产日韩在线| 九色成人免费人妻av| 免费看日本二区| 国产亚洲精品久久久久久毛片| 亚洲av日韩精品久久久久久密| 欧美日韩亚洲国产一区二区在线观看| 999精品在线视频| 最近最新中文字幕大全免费视频| 精品一区二区三区av网在线观看| 丁香欧美五月| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 噜噜噜噜噜久久久久久91| 久久伊人香网站| 1000部很黄的大片| 国产精品 国内视频| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| 99热6这里只有精品| 国产精品国产高清国产av| 哪里可以看免费的av片| 老司机福利观看| 国内精品美女久久久久久| 叶爱在线成人免费视频播放| 中文字幕久久专区| 欧美午夜高清在线| 老司机午夜福利在线观看视频| 1000部很黄的大片| 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 精品国产三级普通话版| 中出人妻视频一区二区| 好男人电影高清在线观看| 老鸭窝网址在线观看| 在线国产一区二区在线| 老汉色av国产亚洲站长工具| 亚洲天堂国产精品一区在线| 一区福利在线观看| 精品久久蜜臀av无| 国产男靠女视频免费网站| 香蕉国产在线看| 激情在线观看视频在线高清| 免费大片18禁| 人妻丰满熟妇av一区二区三区| www日本黄色视频网| 成人av在线播放网站| 白带黄色成豆腐渣| 久久久久久大精品| 久久午夜综合久久蜜桃| 久久精品国产亚洲av香蕉五月| 亚洲人成伊人成综合网2020| 日韩免费av在线播放| 免费看十八禁软件| 亚洲va日本ⅴa欧美va伊人久久| 哪里可以看免费的av片| 国产精品一区二区免费欧美| 午夜福利欧美成人| 美女高潮的动态| 欧美激情久久久久久爽电影| 国产私拍福利视频在线观看| 法律面前人人平等表现在哪些方面| 成年女人看的毛片在线观看| 亚洲av日韩精品久久久久久密| 丰满人妻熟妇乱又伦精品不卡| 亚洲自拍偷在线| 亚洲国产色片| 999精品在线视频| 中亚洲国语对白在线视频| 在线观看舔阴道视频| 中文字幕人妻丝袜一区二区| 成人18禁在线播放| avwww免费| 在线看三级毛片| 色哟哟哟哟哟哟| 夜夜爽天天搞| 国产精品国产高清国产av| 久久久精品欧美日韩精品| 亚洲在线自拍视频| 91麻豆精品激情在线观看国产| 亚洲欧美精品综合久久99| 日韩欧美国产一区二区入口| 亚洲人与动物交配视频| 精品电影一区二区在线| av女优亚洲男人天堂 | 成人av在线播放网站| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 波多野结衣巨乳人妻| 久久久久久久久久黄片| 午夜久久久久精精品| 啦啦啦免费观看视频1| 国产日本99.免费观看| 一区二区三区高清视频在线| tocl精华| 在线观看免费视频日本深夜| 免费搜索国产男女视频| 日韩精品中文字幕看吧| 69av精品久久久久久| 久久精品国产综合久久久| 国产成人福利小说| 成人欧美大片| 九色成人免费人妻av| 国产高清三级在线| 精品国产美女av久久久久小说| www国产在线视频色| 亚洲精品在线观看二区| 无人区码免费观看不卡| 日本黄大片高清| 国内毛片毛片毛片毛片毛片| x7x7x7水蜜桃| 搡老熟女国产l中国老女人| 男女下面进入的视频免费午夜| 免费观看的影片在线观看| 一区二区三区国产精品乱码| 2021天堂中文幕一二区在线观| 久久久久久大精品| 精品无人区乱码1区二区| 亚洲自偷自拍图片 自拍| 丁香欧美五月| 国产成人影院久久av| 国产探花在线观看一区二区| 男人舔奶头视频| 日本a在线网址| 成人亚洲精品av一区二区| 欧美激情久久久久久爽电影| 99久久精品热视频| 亚洲五月婷婷丁香| 国产伦精品一区二区三区四那| 亚洲精品在线观看二区| 国产伦人伦偷精品视频| 亚洲国产精品久久男人天堂| 国产一区二区三区视频了| 长腿黑丝高跟| aaaaa片日本免费| 日本熟妇午夜| 久久久成人免费电影| 国产美女午夜福利| 亚洲国产精品合色在线| 中文字幕熟女人妻在线| 91九色精品人成在线观看| 噜噜噜噜噜久久久久久91| 看黄色毛片网站| 亚洲中文av在线| 国产av麻豆久久久久久久| 99在线人妻在线中文字幕| 久久久久久久久免费视频了| 五月玫瑰六月丁香| 老司机午夜十八禁免费视频| 在线观看免费视频日本深夜| 欧美一区二区国产精品久久精品| 中出人妻视频一区二区| 我的老师免费观看完整版| www.www免费av| 亚洲自偷自拍图片 自拍| 九色成人免费人妻av| 九色国产91popny在线| 日韩 欧美 亚洲 中文字幕| 在线免费观看不下载黄p国产 | 久久久成人免费电影| 亚洲成a人片在线一区二区| 成年女人看的毛片在线观看| 最近最新免费中文字幕在线| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 国产视频内射| 一区二区三区激情视频| 我要搜黄色片| 啦啦啦韩国在线观看视频| 国内少妇人妻偷人精品xxx网站 | 一级毛片精品| 色综合亚洲欧美另类图片| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 久久久国产精品麻豆| www.熟女人妻精品国产| 日韩人妻高清精品专区| 成人午夜高清在线视频| 精品久久久久久久久久免费视频| 悠悠久久av| 亚洲av第一区精品v没综合| 天堂网av新在线| 亚洲精品一区av在线观看| 18美女黄网站色大片免费观看| 国产激情偷乱视频一区二区| 亚洲国产精品999在线| 一本一本综合久久| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 午夜精品久久久久久毛片777| 欧美激情久久久久久爽电影| 国产成人福利小说| 美女扒开内裤让男人捅视频| 两个人的视频大全免费| 最近最新免费中文字幕在线| 丁香六月欧美| 在线观看一区二区三区| aaaaa片日本免费| 12—13女人毛片做爰片一| 国产高清有码在线观看视频| 国产人伦9x9x在线观看| 午夜视频精品福利| 亚洲av第一区精品v没综合| 免费在线观看日本一区| 亚洲真实伦在线观看| 成人无遮挡网站| 波多野结衣巨乳人妻| 亚洲片人在线观看| 黄色片一级片一级黄色片| 亚洲av五月六月丁香网| 亚洲中文av在线| 中文在线观看免费www的网站| 色在线成人网| 国产麻豆成人av免费视频| 在线a可以看的网站| 综合色av麻豆| 99热这里只有精品一区 | 国产成年人精品一区二区| 久久草成人影院| 12—13女人毛片做爰片一| 岛国在线免费视频观看| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线| 亚洲aⅴ乱码一区二区在线播放| 真人做人爱边吃奶动态| 国产伦精品一区二区三区视频9 | 欧美激情在线99| 亚洲在线自拍视频| 女人高潮潮喷娇喘18禁视频| a在线观看视频网站| 99久久99久久久精品蜜桃| 午夜a级毛片| 国产麻豆成人av免费视频| 成人永久免费在线观看视频| 国产精品免费一区二区三区在线| 18禁黄网站禁片午夜丰满| 成人国产综合亚洲| 国产欧美日韩精品亚洲av| 桃色一区二区三区在线观看| 国产精品1区2区在线观看.| 国产伦人伦偷精品视频| 国产高清videossex| 成人av一区二区三区在线看| 听说在线观看完整版免费高清| 精品久久久久久久久久久久久| 久久精品亚洲精品国产色婷小说| 久久久成人免费电影| 国产亚洲精品久久久com| 国产综合懂色| 最好的美女福利视频网| 99精品久久久久人妻精品| 婷婷精品国产亚洲av| 亚洲九九香蕉| 国产成人av激情在线播放| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 日本成人三级电影网站| 国产高清有码在线观看视频| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 99久国产av精品| 国产爱豆传媒在线观看| 亚洲精品久久国产高清桃花| 欧美绝顶高潮抽搐喷水| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 丰满的人妻完整版| 99久久综合精品五月天人人| www.自偷自拍.com| 十八禁人妻一区二区| 熟妇人妻久久中文字幕3abv| 欧美性猛交黑人性爽| 日韩欧美国产在线观看| 搡老岳熟女国产| 国内精品久久久久久久电影| 一区二区三区激情视频| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久人妻蜜臀av| 国产蜜桃级精品一区二区三区| 首页视频小说图片口味搜索| 日韩人妻高清精品专区| 久99久视频精品免费| 国产毛片a区久久久久| 免费看美女性在线毛片视频| 看片在线看免费视频| 美女高潮喷水抽搐中文字幕| 99久久综合精品五月天人人| 成人av一区二区三区在线看| 天天躁日日操中文字幕| 少妇熟女aⅴ在线视频| 嫩草影院精品99| 99热这里只有是精品50| 国产精品一及| 观看免费一级毛片| 99热只有精品国产| 真实男女啪啪啪动态图| 亚洲男人的天堂狠狠| 国产三级在线视频| 色视频www国产| 午夜福利18| 欧美黄色淫秽网站| 2021天堂中文幕一二区在线观| 99久久精品国产亚洲精品| 日韩精品青青久久久久久| 热99re8久久精品国产| 久久久久亚洲av毛片大全| 三级毛片av免费| 亚洲精品在线观看二区| e午夜精品久久久久久久| 午夜福利18| 久久草成人影院| 午夜福利欧美成人| 国产成人啪精品午夜网站| 精品国产三级普通话版| 日本在线视频免费播放| 1024香蕉在线观看| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 成人一区二区视频在线观看| 19禁男女啪啪无遮挡网站| 国产欧美日韩精品一区二区| 国产精品日韩av在线免费观看| 小蜜桃在线观看免费完整版高清| 噜噜噜噜噜久久久久久91| 欧美日韩国产亚洲二区| 性色av乱码一区二区三区2| 免费在线观看视频国产中文字幕亚洲| 亚洲av第一区精品v没综合| 日韩精品中文字幕看吧| 国产精品久久视频播放| 一区二区三区激情视频| 在线观看日韩欧美| 久久精品人妻少妇| 国产美女午夜福利| 九九久久精品国产亚洲av麻豆 | 一级a爱片免费观看的视频| 国产成人av激情在线播放| 性色av乱码一区二区三区2| 国内精品美女久久久久久| 老司机午夜福利在线观看视频| 天堂动漫精品| 欧美日韩瑟瑟在线播放| 国产欧美日韩一区二区精品| www日本黄色视频网| 国产精品电影一区二区三区| 啪啪无遮挡十八禁网站| 亚洲专区中文字幕在线| 99久久成人亚洲精品观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产精品合色在线| 国产美女午夜福利| 国产精品久久久久久精品电影| 亚洲国产日韩欧美精品在线观看 | 两个人看的免费小视频| 一个人免费在线观看电影 | 国产亚洲欧美98| 亚洲av免费在线观看| 欧美乱色亚洲激情| 欧美色视频一区免费| 男人舔女人下体高潮全视频| 国产欧美日韩精品一区二区| 国产精品久久视频播放| 午夜影院日韩av| 国产成人av激情在线播放| 岛国在线免费视频观看| 国产乱人视频| 国产69精品久久久久777片 | 日本在线视频免费播放| 国产精品一区二区三区四区久久| 日日摸夜夜添夜夜添小说| 久久中文字幕一级| 国产精品乱码一区二三区的特点| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩综合久久久久久 | 国产真实乱freesex| 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产91精品成人一区二区三区| 久久热在线av| 琪琪午夜伦伦电影理论片6080| 成人永久免费在线观看视频| 国产精品 国内视频| 动漫黄色视频在线观看| 国产麻豆成人av免费视频| 亚洲国产精品999在线| 黄色丝袜av网址大全| 九色成人免费人妻av| 成人永久免费在线观看视频| 最近最新免费中文字幕在线| 免费在线观看影片大全网站| 桃色一区二区三区在线观看| 国产精品av视频在线免费观看| 亚洲欧美日韩高清在线视频| 亚洲性夜色夜夜综合| 黄色日韩在线| 精品欧美国产一区二区三| 国产精品久久久久久亚洲av鲁大| 色尼玛亚洲综合影院| 免费观看人在逋| 12—13女人毛片做爰片一| 99久久成人亚洲精品观看| 精品免费久久久久久久清纯| 又粗又爽又猛毛片免费看| 日本一本二区三区精品| 国产精品久久久久久人妻精品电影| 人人妻,人人澡人人爽秒播| 日本一本二区三区精品| 99久久无色码亚洲精品果冻| 人人妻人人澡欧美一区二区|