一种缺失交通数据的填充设备、装置及方法与流程

文档序号:18755346发布日期:2019-09-24 22:12阅读:241来源:国知局
一种缺失交通数据的填充设备、装置及方法与流程

本申请涉及数据填充领域,特别涉及一种缺失交通数据的填充设备、装置、方法及可读存储介质。



背景技术:

由于高速路上检测仪器硬件的损坏或软件运作失误等原因,在观察高速公路交通数据的过程中发现,缺失数据是非常常见的。数据缺失会影响到后面的数据分析步骤,导致无法对高速公路上实际的交通状态进行合理分析处理,因而难以对该高速路的交通情况作出准确预测以及报警。

可见,对高速公路的缺失交通数据进行处理是高速公路智能交通数据质量优化控制系统中至关重要的一步,然而现有的缺失数据填充方案的准确性较低,导致即使在填充缺失数据后也无法提升数据分析结果的可靠性。



技术实现要素:

本申请的目的是提供一种缺失交通数据的填充设备、装置、方法及可读存储介质,用以解决现有的缺失数据填充方案的准确性较低,导致即使在填充缺失数据后也无法提升数据分析结果的可靠性的问题。具体方案如下:

第一方面,本申请提供了一种缺失交通数据的填充设备,包括:

存储器:用于存储计算机程序;

处理器:用于执行所述计算机程序,以实现以下步骤:

获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

优选的,所述处理器还用于:

获取预先设置的最大偏差相似性准则算法的算法参数的初始取值范围;将所述算法参数作为骨干粒子群算法中骨干粒子的位置向量,利用骨干粒子群算法对所述骨干粒子的位置向量进行优化,得到所述算法参数在所述初始取值范围上的最优取值。

优选的,所述处理器具体用于:

初始化骨干粒子群算法中骨干粒子的位置向量;

对所述骨干粒子的位置向量进行更新;

判断更新后的位置向量是否发生越界行为;

若没有发生,则确定所述更新后的位置向量的适应度数值;

根据所述适应度值,分别对所述骨干粒子的最优位置向量和整个骨干粒子群的最优位置向量进行更新;

在达到预设终止条件时,确定最优位置向量。

优选的,所述处理器具体用于:

根据knn填充算法,确定所述聚类簇中缺失交通数据与所述聚类簇中各个完整交通数据之间的欧氏距离;

确定最小的预设数量的欧式距离所对应的目标完整交通数据;

根据所述目标完整交通数据的加权平均值,对所述聚类簇中的缺失交通数据进行填充。

第二方面,本申请还提供了一种缺失交通数据的填充装置,包括:

数据获取模块:用于获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;

聚类模块:用于根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;

数据填充模块:用于利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

优选的,缺失交通数据的填充装置还包括:

初始范围设置模块:用于设置最大偏差相似性准则算法的算法参数的初始取值范围;

最优取值确定模块:用于将所述算法参数作为骨干粒子群算法中骨干粒子的位置向量,利用骨干粒子群算法对所述骨干粒子的位置向量进行优化,得到所述算法参数在所述初始取值范围上的最优取值。

优选的,所述数据填充模块包括:

欧式距离确定单元:用于根据knn填充算法,确定所述聚类簇中缺失交通数据与所述聚类簇中各个完整交通数据之间的欧氏距离;

目标完整交通数据确定单元:用于确定最小的预设数量的欧式距离所对应的目标完整交通数据;

数据填充单元:用于根据所述目标完整交通数据的加权平均值,对所述聚类簇中的缺失交通数据进行填充。

第三方面,本申请还提供了一种缺失交通数据的填充方法,包括:

获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;

根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;

利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

第四方面,本申请还提供了一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时用于实现以下步骤:

获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

本申请所提供的一种缺失交通数据的填充设备、装置、方法及可读存储介质,能够获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;根据算法参数,利用最大偏差相似性准则算法对缺失交通数据集进行聚类,得到聚类簇集合;最终利用knn填充算法分别对聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

可见,该方案考虑到交通数据特点,基于最大偏差相似性准则算法与knn填充算法对缺失交通数据进行填充,一方面,针对knn填充算法易产生最邻近噪点的问题,使用最大偏差相似性准则来针对缺失交通数据进行聚类,另一方面,针对最大偏差相似性准则算法的参数选取困难的问题,利用骨干粒子群算法对其参数优化。最终实现了准确填充缺失交通数据的目的,保证了填充后的交通数据的分析结果的可靠性。

附图说明

为了更清楚的说明本申请实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。

图1为一种缺失交通数据的填充设备实施例一的结构示意图;

图2为一种缺失交通数据的填充设备实施例二的填充过程示意图;

图3为本申请提供的对比实验的实验结果示意图;

图4为一种缺失交通数据的填充设备装置实施例的功能框图;

图5为一种缺失交通数据的填充方法的实现流程图。

具体实施方式

为了使本技术领域的人员更好地理解本申请方案,下面结合附图和具体实施方式对本申请作进一步的详细说明。显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。

目前,采集到的交通数据往往存在缺失数据,缺失数据为后续数据分析过程带来诸多干扰,导致数据分析结果失真。然而,现有的数据填充方法无法准确填充缺失数据,无法提升缺失交通数据分析结果的可靠性。针对该问题,本申请提供一种缺失交通数据的填充设备、装置、方法及可读存储介质,利用基于骨干粒子群算法优化的最大偏差相似性准则算法对缺失交通数据进行聚类,并利用knn填充算法分别对各个聚类簇中的缺失交通数据进行填充,实现了准确填充缺失交通数据的目的,提升了缺失交通数据的分析结果的可靠性。

下面对本申请提供的一种缺失交通数据的填充设备实施例一进行介绍,参见图1,实施例一包括:

存储器100:用于存储计算机程序;

处理器200:用于执行所述计算机程序,以实现以下步骤:

获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

上述缺失交通数据集包括缺失交通数据,此外还包括完整交通数据,具体形式为时间序列数据,其中缺失交通数据是指在特定时间未按计划从特定主体获得特定变量的交通数据。

在本实施例中,缺失交通数据的填充设备可以是pc(personalcomputer,个人电脑),也可以是服务器、平板电脑、网络节点等。

存储器100至少包括一种类型的可读存储介质,具体可以包括以下任意一项或多项可读存储介质:闪存、硬盘、多媒体卡、卡型存储器、磁性存储器、磁盘、光盘。在一些实施例中,存储器100可以是缺失交通数据的填充设备的内部存储单元,例如缺失交通数据的填充设备的硬盘;也可以是缺失交通数据的填充设备的外部存储设备,例如缺失交通数据的填充设备上的插接式硬盘,智能存储卡(smartmediacard,smc),安全数字(securedigital,sd)卡,闪存卡(flashcard)等。进一步地,存储器100还可以既包括缺失交通数据的填充设备的内部存储单元,也包括外部存储设备。存储器100不仅可以用于存储安装于缺失交通数据的填充设备的应用软件及各类数据,例如计算机程序的代码等,还可以用于暂时地存储已经输出或者将要输出的数据。

处理器200在一些实施例中可以是一中央处理器(centralprocessingunit,cpu)、控制器、微控制器、微处理器或其他数据处理芯片,用于运行存储器100中存储的程序代码或处理数据,例如执行计算机程序等。

存储器100与处理器200之间的总线可以是外设部件互连标准总线或扩展工业标准结构总线等,该总线具体可以分为地址总线、数据总线、控制总线。为便于表示,图1中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。

具体的,在利用骨干粒子群算法优化最大偏差相似性准则算法的算法参数的过程中,处理器200用于:

获取预先设置的最大偏差相似性准则算法的算法参数的初始取值范围;将所述算法参数作为骨干粒子群算法中骨干粒子的位置向量,利用骨干粒子群算法对所述骨干粒子的位置向量进行优化,得到所述算法参数在所述初始取值范围上的最优取值。

作为一种具体的实施方式,在利用骨干粒子群算法对骨干粒子的位置向量进行优化的过程中,处理器200用于:

初始化骨干粒子群算法中骨干粒子的位置向量;

对所述骨干粒子的位置向量进行更新;

判断更新后的位置向量是否发生越界行为;

若没有发生,则确定所述更新后的位置向量的适应度数值;

根据所述适应度值,分别对所述骨干粒子的最优位置向量和整个骨干粒子群的最优位置向量进行更新;

在达到预设终止条件时,确定最优位置向量。

具体的,在利用knn填充算法分别对聚类簇集合中聚类簇的缺失交通数据进行填充的过程中,处理器200用于:

根据knn填充算法,确定所述聚类簇中缺失交通数据与所述聚类簇中各个完整交通数据之间的欧氏距离;

确定最小的预设数量的欧式距离所对应的目标完整交通数据;

根据所述目标完整交通数据的加权平均值,对所述聚类簇中的缺失交通数据进行填充。

本实施例所提供一种缺失交通数据的填充设备,考虑到交通数据特点,基于最大偏差相似性准则算法与knn填充算法对缺失交通数据进行填充,一方面,针对knn填充算法易产生最邻近噪点的问题,使用最大偏差相似性准则来针对缺失交通数据进行聚类,另一方面,针对最大偏差相似性准则算法的参数选取困难的问题,利用骨干粒子群算法对其参数优化。最终实现了准确填充缺失交通数据的目的,保证了填充后的交通数据的分析结果的可靠性。

下面开始详细介绍本申请提供的一种缺失交通数据的填充设备实施例二,实施例二基于上述实施例一实现,并在实施例一的基础上进行了一定程度上的拓展。

实施例二中缺失交通数据的填充设备主要包括存储器和处理器,本实施例主要对缺失交通数据的填充设备的工作流程进行介绍,参见图2,工作流程具体包括:

步骤s21、获取缺失交通数据集;

具体的,上述缺失交通数据集可以为高速公路的缺失交通数据集,根据缺失交通数据集构建n×m的时间序列数据ai=(ai1,ai2,...,aim),其中i=1,2,...,n,aij(j=1,2,...,m)为交通数据记录。

步骤s22、执行骨干粒子群算法的初始化操作;

初始化操作主要包括参数设置和种群初始化两个过程,下面分别就这两个过程进行介绍:

参数设置:设置粒子数目为n;粒子群算法迭代最大次数tmax;最大偏差相似性准则算法的参数α,β,γ各自的取值范围分别为[radown,raup],[rβdown,rβup]与[rγdown,rγup]。其中,radown为α的下界,取值范围为[0.3,0.8];raup为α的上界,取值范围为[0.9,1.6];rβdown为β的下界,取值范围为(0,1);rβup为β的上界,取值范围为[0.3,1);rγdown为γ的下界,取值范围为(0,0.1];rγup为γ的上界,取值范围为[0.1,0.5]。作为一种具体的实施方式,本实施例设置骨干粒子群种群数目为n=30,最大迭代次数tmax=1000,radown=0.7,raup=0.9,rβdown=m-1,rβup=1-α,rγdown=0.05,rγup=0.25,m=96。

初始化种群:对每个粒子定义位置为x,将最大偏差相似性准则算法的参数(α,β,γ)作为骨干粒子的位置向量,对每个i,第i只磷虾的位置为xi=(xi1,xi2,xi3),其中xi1∈[radown,raup],xi2∈[rβdown,rβup],xi3∈[rγdown,rγup]。

步骤s23、执行骨干粒子群算法,得到最优粒子位置,以作为最大偏差相似性准则算法的最优参数;

步骤s24、根据上述最优参数,利用最大偏差相似性准则算法对缺失交通数据集进行聚类,得到聚类簇集合;

如图2所示,聚类簇集合包括n个聚类簇。

步骤s25、对于聚类簇集合中的各个聚类簇,根据knn填充算法对所述聚类簇中的缺失交通数据进行填充;

步骤s26、得到完整交通数据集。

作为一种具体的实施方式,本实施例的步骤s23具体包括以下步骤:

步骤s231:根据骨干粒子群每个个体的位置,计算适应度函数的函数值;

步骤s232:更新骨干粒子群中个体位置;

具体的,更新个体粒子位置公式如式(1)所示:

其中,pbij(t)为i粒子在当前迭代的最优位置的j维度的量,gbj(t)是当前迭代最优解的位置向量的j维度的量。

步骤s233、判断粒子个体是否有越界行为,若有则重新初始化,否则进入s234;

步骤s234、重新计算全部粒子个体的适应度值;

步骤s235、更新所有粒子个体的pb,并将粒子群全局历史最优位置赋给gb;

更新公式如式(2)所示:

步骤s236、判断是否达到终止条件,如果达到进入步骤s237,否则跳转到s232;

步骤s237、输出最优粒子位置。

根据前文描述,可以理解的是,本实施例对于骨干粒子每一维度的界限如下:

xminj≤xij(t+1)≤xmaxj(3)

其中xminj及xmaxj分别为全局j维度上的下界和上界。

作为一种具体的实施方式,本实施例的步骤s24具体包括以下步骤:

步骤s241、根据最优粒子位置,确定最大偏差相似性准则算法精度的最优参数;

步骤s242、将交通数据曲线集合赋给f集合,得f={x1,x2,...,xn},将交通数据曲线数量赋给nf,得nf=n,并根据最大偏差相似性准则算法算出n0以及m0;

本实施例中,设最大偏差相似性准则算法中聚类结果数据集为d,还未被分类的数据曲线集合为f,f集合里面的数据曲线数量是nf,以及xi作为核心的相似数据集合为s(xi),其中数据曲线数量n(xi)=|s(xi)|,s(xi)集合里其余的曲线跟核心xi总的度量距离d(xi),初始化:nf=0,d(xi)=0。

最大偏差相似性准则内容如下:

(1)设xi和xj所相应时间点之间的绝对差值为sijk,其中i,j=1,2,...,n,k=1,2,...,m,绝对差值公式如式(7)所示:

sijk=|xik-xjk|(7)

(2)设所有小于等于γ的sijk个数为nij,即nij={sijk≤γ},同时nij被称为xi跟xj之间的相似时点数;相反设mij是能够连续满足sijk>γ这一条件的sijk最多的个数,同时称之为xi跟xj的最大连续偏离时点数,其中i,j的取值范围为[1,n];k∈{1,2,…,m};γ是提前设置的一个常数值,并满足0≤γ≤1,是度量xi与xj对应的时间点之间的数据相似性的阈值。当sijk大于γ时,可确认xik跟xjk是不相似的,相反则表示它们是相似的;关于mij的计算公式如式(8)所示:

以高速公路交通数据xi作为参考对象,根据公式(7)来算出xj和参考对象之间的mij以及nij,最大偏差相似性准则如下两条:

(a)设相似度为n0,并让n0=[α×m],其中α是预先设置的常数,取值范围:0≤α≤1-m-1,使得nij≥n0;

(b)设偏离度为m0,并让m0=[β×m],其中β是预先设置的常数,取值范围:0≤β≤1-α,使得mij≤m0;

若满足最大偏差相似性准则,可以确认xj跟xi是相似的。

步骤s243、对全部i,j=1,2,...,n且i≤j,根据最大偏差相似性准则算法算出nij与mij,并根据式(9)得到xi和xj之间的度量距离:

具体的距离公式为:

步骤s244、对全部i,j=1,2,...,n且i≤j,将xi当做中心点,同时把算得的nij和mij分别跟n0、m0对比,对于能够达到最大偏差相似性准则的数据曲线放进s(xi)集合里面,并让s(xi)=s(xi)∪{xj},同样的使s(xj)=s(xj)∪{xi},再让n(xi)=|s(xi)|,通过式(10)得到d(xi):

步骤s245、判断nf是否等于0,若是则跳转到s249,反之则执行s246;

步骤s246、在f集合中根据n(xi)值进行降序排序,找到最大n(xi)值对应的对象,同时设y1,y2,...,yt是符合要求的所有对象,yi∈f,i=1,2,...,t,并求数据曲线z0,要求d(z0)=min{d(yi)},且z0∈{y1,y2,...,yt},如果存在不止一个这样的数据曲线z0,要从中选一个。

步骤s247、更新最大偏差相似性准则算法聚类结果数据集d,即d=d∪{s(z0)},并更新未被分类的数据曲线集合为f以及其包含的数据曲线数量nf,即f=f-s(z0),nf=nf-n(xi)。

步骤s248、计算未聚类数据曲线的n值以及距离d,即选取所有z∈f,让s(z)=s(z)-s(z0),得到n(z)=|s(z)|,同时通过式(10)算出d(z),然后跳转到s234;

步骤s249、输出最后的最大偏差相似性准则算法聚类结果数据集d。

作为一种具体的实施方式,本实施例的步骤s25具体包括以下步骤:

s251、对于各个缺失数据,分别判断其属于s24确定的聚类簇集合中哪个聚类簇,根据缺失数据所在簇内的数据构建完整记录的数据矩阵af=(af1,af2,...,afm),设缺失数据的数量为t,缺失目标为ai(i=1,2,...,t);

s252、计算出同一簇内的每个正常数据记录与该目标ai的欧式距离,放到集合di;

s253、对欧式距离集合di里的值进行排序筛选,选取前k个最小的欧式距离值所对应的高速路交通数据作为缺失目标ai的k邻近值;

s254、得到k个邻近值后,计算这k个值的加权平均值fi,即该缺失数据ai的填补值,令ai=fi,并令i=i+1;

s255、确定i和t的大小关系,若i﹥t执行s256,否则返回s252;

s256、在所有聚类簇的缺失数据被全部填充完毕,结束算法。

本实施例中,加权平均值f的计算公式如下:

其中e为常数,d为欧式距离,xj为k个邻近值的其中一个完整数据。

可见,本实施例提供了一种缺失交通数据的填充设备,上文主要对该设备执行的方案进行了描述,根据上文可知,该方案提供了基于最大偏差相似性准则与knn的缺失数据填充算法来对缺失交通数据进行填充的方案。具体的,针对最大偏差相似性准则算法的参数选取困难问题,该方案基于骨干粒子群算法对其参数优化;针对knn填充算法会产生最邻近噪点的问题,该方案使用参数优化后的最大偏差相似性准则算法来针对缺失交通数据中缺失属性对完整交通数据进行聚类,并利用knn填充算法对聚类得到各个聚类簇进行填充。显著提升了缺失交通数据的填充方案的准确性和实时性。

为证明本实施例的实施效果,本申请还进行了仿真实验,对本申请中基于参数优化后的最大偏差相似性准则算法的knn填充算法、基于传统的最大偏差相似性准则算法的knn填充算法、knn填充算法三种方案进行了对比。如图3所示,在缺失比例由0到50%的区间中华,本申请中基于参数优化后的最大偏差相似性准则算法的knn填充算法的实施效果均优于其他两种算法,缺失交通数据填补的rmse值更小,效果更优。

下面对本申请实施例提供的一种缺失交通数据的填充装置进行介绍,如图4所示,该装置包括:

数据获取模块401:用于获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;

聚类模块402:用于根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;

数据填充模块403:用于利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

作为一种具体的实施方式,该装置还包括:

初始范围设置模块:用于设置最大偏差相似性准则算法的算法参数的初始取值范围;

最优取值确定模块:用于将所述算法参数作为骨干粒子群算法中骨干粒子的位置向量,利用骨干粒子群算法对所述骨干粒子的位置向量进行优化,得到所述算法参数在所述初始取值范围上的最优取值。

作为一种具体的实施方式,所述数据填充模块403包括:

欧式距离确定单元:用于根据knn填充算法,确定所述聚类簇中缺失交通数据与所述聚类簇中各个完整交通数据之间的欧氏距离;

目标完整交通数据确定单元:用于确定最小的预设数量的欧式距离所对应的目标完整交通数据;

数据填充单元:用于根据所述目标完整交通数据的加权平均值,对所述聚类簇中的缺失交通数据进行填充。

本实施例的缺失交通数据的填充装置与前述缺失交通数据的填充设备实施例的实施方案相对应,具体的,本实施例的缺失交通数据的填充装置可以为一组功能模块架构,例如app程序,用于实现前述缺失交通数据的填充设备的工作流程。

此外,本申请还提供了一种缺失交通数据的填充方法,如图5所示,包括:

步骤s501、获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;

步骤s502、根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;

步骤s503、利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

最后,本申请还提供了一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时用于实现以下步骤:

获取待填充的缺失交通数据集,并获取利用骨干粒子群算法优化得到的最大偏差相似性准则算法的算法参数;根据所述算法参数,利用最大偏差相似性准则算法对所述缺失交通数据集进行聚类,得到聚类簇集合;利用knn填充算法分别对所述聚类簇集合中聚类簇的缺失交通数据进行填充,得到完整交通数据集。

本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。

结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(ram)、内存、只读存储器(rom)、电可编程rom、电可擦除可编程rom、寄存器、硬盘、可移动磁盘、cd-rom、或技术领域内所公知的任意其它形式的存储介质中。

以上对本申请所提供的方案进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方法及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本申请的限制。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1