基于运动序列分割的运动捕获数据关键帧提取

更新时间:2023-05-20 18:56:01 阅读量: 实用文档 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

第20卷第6期

计算机辅助设计与图形学学报

voI.20,No.6

2008年6月

JOURNALOFCOMPUTER—AIDEDDESIGN&COMPUTERGRAPHICS

June,2008

基于运动序列分割的运动捕获数据关键帧提取

朱登明1’2’3’

王兆其”

11(中国科学院计算技术研究所虚拟现实技术实验室北京

100190)

21(中国科学院计算技术研究所智能信息处理重点实验室北京100190)3’(中国科学院研究生院北京100049)

(mdzhu@ict.ae.cn)

摘要采用线性时不变系统把高维运动数据映射到低维状态空问;在低维状态空间中,定义了姿态之间的相似性度量;并采用误差平方和准则对时序的低维数据点集进行运动分割,分割点上的运动姿态被定义为关键帧.实验结果表明:该算法能够较好地提取出运动序列中的关键帧,并且这些关键帧能够很好地概括原始运动序列的内容.

关键词运动捕获数据;关键帧;线性时不变系统;运动分割中图法分类号TP391.41

ExtractionofKeyframefromMotion

CaptureDataBased

on

MotionSequence

Segmentation

ZhuDengmin91 2 3’

Wang

Zhaoqil’

”(VirtualRealityLaboratory,Institute

ofComputingTechnology,Chir"seAcademy

ofSciences,BeOing

100190)

2’(KeyLaboratory

ofIntelligentIn,ormationProcessing,Institute

ofComputingTechnology,ChineseAcademyofSciences,Beijing

100190)

University

ofChineseAcademyofSciences-Beijing

100049)

Abstract

Lineartime-invariantsystemis

usedtoderive

an

explicit

mappingbetweenthe

high—

dimensionalmotioncapturedataandthelOW—dimensional

state

variables.AsimilaritymetriciSdefined

to

measurethe

differencebetweendifferentposesinthelow-dimensionalstatespace.andthen

themethodofmeansquared

error

is

employed

to

divide

the

motion

sequence

into

asequence

of

concatenatedsegments.Theposes

at

these

segmentationpointsare

thendefined

as

keyframes.Experimental

results

show

thattheextractedkeyframes

by

our

method

can

give

goodvisual

summarizationoftheoriginalmotionsequence.

Keywordsmotioncapturedata;keyframe;lineartime-invariantsystem;motionsegmentation

运动捕获技术是一种新的动画生成技术‘卜21,它角色动画研究的热点问题,被广泛地应用于电影广实时地检测、记录表演者的肢体在三维空间中的运告、数字娱乐、体育仿真等领域.在运动捕获技术的动轨迹,捕捉表演者的动作,并将其转化为数字化的广泛使用过程中,越来越多的三维运动数据库可供动画师选择.但在基于高频采集到的运动数据里存为直观,效果更加生动逼真.运动捕获技术是近年来

在大量的冗余信息,为了更加有效地把运动捕获

收稿日期:2007—10—23;修回日期:2008—03—05.基金项目:国家自然科学基金(60573162,60533070,60603082)#国家。八六三”高技术研究向为虚拟现实,多模式人机接口、人工智能等.

钉(Graduate“抽象运动”,以驱动虚拟人模型,使动画制作过程更发展计划(2006AA012336,2007AA012320);中国科学院“科技助残行动计划”项目(KGCX2一Yw一610).朱登明,男,1973年生。博士研究生,副研究员,CCF会员,主要研究方向为虚拟现实、人体动画合成等.王兆其,男,1966年生,博士,研究员,博士生导师,CCF高级会员,主要研究方

788计算机辅助设计与图形学学报2008焦

数据用于人体运动合成,以及更加方便地存储、传输或浏览运动数据库的内容,需要提取出表示运动数据内容的关键帧.运动数据有其内在的特征,首先,运动数据是高维向量,和视频、图像文件一样,它们对存储空间要求很大.利用提取的关键帧,可以有效地节省存储空间.其次,采用人工的方法对运动数据库进行标注,容易造成主观理解不统一的问题.而较好的关键帧提取方法可以有效地提取出表示运动序列内容的关键帧,通过浏览运动序列的关键帧,用户可以快速地了解运动数据库的内容.

由于人体运动数据是一种特殊的时序运动信号,并且表示人体运动的数据都是高维向量,不同关节之间运动的耦合性很强,如人在正常行走时,人的手臂摆动和腿部运动就存在关联.因此怎样分析蕴含在原始运动数据内的运动特征,并自动提取出表示运动序列内容的关键帧,已成为一个研究的难点.已有的运动数据的关键帧提取方法主要分为两大类:一类是等间隔采样的关键帧提取方法.对原始运动序列进行等间隔重采样,但在原始运动过程中,由于动作节奏或频率并不一样,均匀采样的方法会导致在动作缓慢部分过采样、在动作变化剧烈的地方欠采样,从而丢失了原始运动的细节信息;另一类方法就是自适应关键帧提取方法.根据运动数据的内容特征,自动提取能表示运动序列的关键帧.自适应关键帧提取方法主要分为基于聚类和曲线简化的关键帧提取方法.

Zhuang等[31采用无监督聚类的方法提取视频关键帧,用于视频内容分析.Liu等[41采用聚类方法提取运动数据的关键帧,将每类的首帧作为关键帧,并在此基础之上构建运动索引树.Loy等[5]同样采用聚类的方法提取图像序列的关键帧,但其是以每类的中心作为关键帧.由于聚类方法能够很好地把相似的样本归为一类,因此利用这类方法提取的关键帧可以较好地表示原始样本的内容;但聚类方法很少考虑样本之间的时序相关性,容易导致对运动序列的分析失真.

Lim等[6]首先提出了一种基于曲线简化‘71的运动数据的关键帧提取方法,将每帧运动数据作为高维空间中的一个数据点,这条曲线的凹、凸点就相当于运动序列中的关键帧;该方法的缺点是忽略了运动序列的时序相关性和不同关节运动之间的耦合性,采用简单的欧氏距离度量不同数据帧之间的相

似性,结果不能反映真实的差异.沈军行等嘲也采用类似的方法提取运动数据的关键帧,但他们采用基于四元数距离的方法度量数据帧之间的相关性.杨涛等[9]提出了一种基于分层曲线简化的方法提取运动数据的关键帧,引入骨骼夹角作为运动特征,并以此确定候选关键帧,最后采用分层曲线简化算法精选候选关键帧,以获得最终关键帧集合.Assa等[10]采用多维尺度分析的方法,把高维运动数据映射到低维欧氏空间,再采用曲线简化技术提取运动数据的关键帧.同样的曲线简化技术也被用于视频序列的关键帧提取[-11-12].

本文采用的关键帧提取方法类似于聚类方法,但与之又有不同:首先,本文采用线性时不变系统,把原始高维运动数据映射到低维状态空间,该方法的优点是在降维时充分地考虑了运动序列之间的时序性,同时抑制了原始高维数据中的噪声;其次,在低维状态空间中,定义姿态之间的相似性度量;最后,采用误差平方和准则对时序的低维数据点集进行运动分割,段与段之间的分割点反映了原始运动序列中姿态变化最为剧烈的关键点,也就是能表示运动序列内容的关键帧.实验结果表明,采用本文方法提取的关键帧,可以很好地概括运动序列的内容.

1运动数据分析

1.1运动数据表示

在表示复杂的人体运动时,通常都把人体模型简化为具有16个关节、51个自由度的骨架模型.简化的骨架模型可以看成是由关节连接的骨骼段的层次结构组成,每一个骨骼都有一个局部坐标系与之相连.其中,根关节有6个自由度,3个表示根关节在世界坐标系中的位置,3个表示根关节在世界坐的3个自由度.在该简化的骨架模型中,子关节的运动受父关节的影响进行相应的运动,而子关节又可以根据自已的特征进行运动,其父关节的运动不受帧运动数据m(£)可以表示为m(£)=(p。(£),q,(£),标系中的朝向,而其他关节只有相对于父关节旋转影响.本文采用单位四元数来表示关节的旋转,其优点是可以避免欧拉角表示带来的奇异点.因此,每一q2(£),…,qM(£));其中,M表示关节数目,p1(£)表示根关节在世界坐标系中的位置,q。(f)∈S3表示根关节在世界坐标中的朝向,q。(£)表示第7/个关节相

6期朱登明等:基于运动序列分割的运动捕获数据关键帧提取

789

对于父关节的旋转.在实际应用中,通常采用对数映射把四元数从S3空间映射到R3空间[13。4。,因此,定

义向量

只一(pl(£),log(ql(f)),log(qz(£)),… log(qM(f)))

(1)

1.2运动数据降维

线性时不变系统[15]是实际系统的一种理想化模型,是对实际系统经过近似化和工程化处理后所导出的一类理想化系统.很多研究者采用线性时不变系统对采集到的人体运动数据进行建模分析,用于动画的生成与合成[16-t7].由于线性时不变系统能够有效地提取出高维运动数据里面蕴含的低维结构,因此本文采用该系统对运动数据进行降维分析.线性系统模型能被状态空间方程描述为[1刀

j“2触t+帆

(2)

【Y。一Cr;+H

其中,毗和',。是高斯白噪声,服从高斯分布:职~N(0,Q),n~N(O,R),Q和R是协方差矩阵;Y。是捕获得到的运动数据集,维度为P,其表示见式(1),为已知参数;毛是低维的隐含状态变量,维度为k,且忌《户;A为状态转移矩阵,dim[A]=k×k;c是输出矩阵,dimEC]----p×k.式(2)中的未知参数{A,C,R,Q,而}可以通过EM算法求得[1引.通过式(2),可以把高维的运动数据集{y,,y:,…,弘}映射到低维

状态空间中的数据点集{工,,勋,…,款).这些低维的

数据点集保持了原始运动数据之间的时序相关性,反映了蕴涵在原始高维数据集里面的低维结构,揭示了原始运动数据的内在结构属性.

2关键帧提取

2.1运动姿态的相似性度■

本文将在低维隐含状态空间中度量人体姿态的相似性.考虑到运动数据是时序的运动序列,在度量姿态的相似性时,还需考虑它的运动学特征.定义姿态Y;与y』之间的距离为

d(yf,Yj)一d(xi,而)一d,(置,毛)+弘 d,(xi,Xj)

(3)

其中,肛是权重参数;

dp(毛,而)=

(工g一却)2,

fllt

d。(置,而)=

(h一%)2,

f11

峨=(工H1一置)/h,h是采样时间间隔,它反映了姿态的运动学特征.

2.2运动序列分割

对时序的数据序列进行分段分析是一种被广泛

应用于模式识别、语音分析、信息检索中的非常有效的技术,它将相似的样本尽可能地归为一段,不相似的样本尽量划分到不同的段中.由于低维状态变量{工l,x:,…,墨}反映了原始高维运动数据{Y,,y。,…,Y。}的内在结构属性,为了能更好地提取出运动序列中的关键帧,本文采用误差平方和准则,在低维状态空间中对状态变量进行分段分析,选取段与段之间的分割点作为原始运动序列中的关键点.由于这些关键点反映了运动序列中姿态变化最为剧烈的点,本文中这些关键点上的运动姿态就被定义为关

键帧.

设原始运动序列的样本数目为T帧,则低维状态空间中隐含状态变量的表示为{工。,工z,…,工T}.本文采用类似于文献[6,9—10]的方法,原始运动序列的首帧和尾帧被默认为运动序列的关键帧.设要提取的关键帧的数目为N,则原始样本序列被分割为N一1段,采用二层动态规划算法可以求出样本序列的最优分割点.由于捕获到的运动数据都有其内在的动力学属性,为了刻画运动的动力学特征,需要用户预先给定运动片断的最小帧长,设每段的最小

帧长为z曲,则N一1≤l产I.设第s段的起始点为

L‘rainJ

第b帧,结束点为第e帧,则第s段的样本均值为

只一F晶善xt

“)

结合式(3)和式(4),定义

D,(6,P)=芝:d(工f,t)

(5)

则可以递归定义第咒段结束于第P帧的最优路径得分

D。(已)=

min

[D。(6,P)+D,l(6—1)]

(H一”‘lmin<6<o

(6)

其中e一6≥z商。.基于式(6),可以采用动态规划算法求出运动序列的全局最优分割点.算法步骤如下:

Stepl.初始化.Do(O)=O,西。(O)=oo,l≤n≤N--1.Step2.从e=l曲到T,百1(P>=移1(1,#),Zmjn≤£≤t

Step3.递推从咒一2到N一1,使

D2(P)一min'rain(K‘I-D2(6,£)+石l(6—1)],

21.i.≤e≤丁;

百3(e)=min[西3(6,e)+D2(6—1)],

790

计算机辅助设计与图形学学报

3l。i。≤P≤丁f;

DN-l(P)一

min[声N-1(6,P)+DⅣ-2(6—1)],‘Ⅳ一Z)‘lmm<Kt

的运动类型来测试本文算法的有效性,并和基于均匀采样、基于瞎线简化‘争71以及将分段算法直接用于原始运动捕获数据的关键帧提取方法进行了实验

比较.

~一l

(N一2)lml。≤e≤丁.

Step4.最优解.{e1,P2,…,eN—l}一argminDⅣ一1(丁);

实验1.翻筋斗.选择运动序列为280帧的翻筋斗运动作为实验样本,采用线性时不变系统,把运动数据从51维的高维运动空间映射到10维的低维状态空间.设要提取的关键帧数目为6帧,l耐。=15,则运动序列分割的计算时间为t=8.2min,实验结果

如图1所示.

其中{e,,e:,…,e。一。}表示该运动序列的全局最优分割点,且ew一。=T.这些最优分割点加上运动序列的首帧就构成了该运动序列的关键帧.

3实验结果与讨论

本文采用Vicon运动跟踪设备捕获真实的人体运动数据,运动数据的采样频率为60Hz,每帧运动数据的自由度为51.所有的实验都是在IntelPIV

3.4GHzCPU,1

实验2.跳跃运动.跳跃运动的序列长度为313帧,低维状态空间的维度为9维,要提取的关键帧数目为7帧,l。。一15,运动序列分割的计算时间为£=

10.1

GB内存的PC机上进行的.

min,实验结果如图2所示.

实验中,选择了翻筋斗、跳跃和跑步这3种不同

a基于均匀采样算法b基于曲线简化算法

c将分段算法直接用于原始捕获数据

d本文算法

图1从翻筋斗运动序列中提取6帧关键帧

a基于均匀采样算法b基于曲线简化算法

c将分段算法直接用于原始捕获数据d本文算法

图2从跳跃运动序列中提取7帧关键帧

实验3.跑步运动.跑步运动的序列长度为250帧,低维状态空间的维度为8维,要提取的关键帧数

目为11帧,z。i。=10,运动序列分割的计算时间为

£=6.3

min,实验结果如图3所示.

6期朱登明等:基于运动序列分割的运动捕获数据关键帧提取

791

a基于均匀采样算法b基于曲线简化算法

c将分段算法直接用于原始捕获数据d本文算法

图3从跑步运动序列中提取11帧关键帧

从图1~3中可以看出,基于均匀采样算法提取运动序列的关键帧容易在运动缓慢的部分过采样,在运动剧烈的部分欠采样;而基于曲线简化和将分段算法直接用于原始运动捕获数据的关键帧提取方法,是直接在原始运动捕获数据的空间中度量运动的相似性,由于没有考虑人体不同关节之间的运动存在耦合性,导致在度量姿态之间的相似性时,其结果不能反映真实姿态之间的差异.基于曲线简化算法只考虑了静态帧与帧之间的距离,没有考虑运动序列在时序上的相关性.本文算法在低维特征空间中度量姿态之间的相似性,采用线性系统模型降维,可以在低维空间中保持原始捕获数据之间的时序相关性,同时可以有效地抑制原始运动捕获数据中的噪声.本文算法能够更好地提取出运动序列的关键帧,并且这些关键帧能够更好地表示原始运动序列

的内容.

归进行的,需要搜索的空间范围较大、比较耗时.因此,本文算法对运动数据的规模有一定的限制,运动序列的长度在1500帧以下比较合适.

到目前为止,判断关键帧提取结果的好坏并没有一个统一的客观度量标准,主要是从视觉效果来判断,带有一定的主观性.从实验结果看,与基于均匀采样、曲线简化、将分段算法直接用于原始捕获数据的关键帧提取方法相比,本文方法能更好地提取出运动序列的关键帧,更好地表示原始运动序列的内容.基于曲线简化算法很难控制要提取的关键帧数目,需要反复地尝试不同的误差阈值以满足需求,这个过程相当耗时;而本文算法可以很容易地控制关键帧的提取数目.

[12

考文献

Dasgupta

A,NakamuraY.Makingfeasiblewalkingmotion

capture

结论

ofhumanoidrobotsfromhumanmotionProceedingsof

IEEE

data[c]//

on

InternationalConferenceRobotics

andAutomation,Detroit,MI。1999。2:1044-1049

本文提出了一种新的关键帧提取方法,首先采用线性时不变系统把高维的运动数据映射到低维的隐含状态空间,这些低维的状态变量保持了原始高维运动数据之间的时序相关性,它反映了原始高维运动数据内在的结构属性,同时线性系统模型能够有效地抑制原始捕获数据中的噪声;其次在低维状态空间中,采用误差平方和准则对时序的低维状态变量进行运动分割,段与段之间的分割点反映了运动序列中姿态变化最为剧烈的关键点,本文中这些关键点被定义为运动序列的关键帧,它能很好地表示原始运动序列的内容.由于本文采用动态规划算法对运动序列进行自动分割,而动态规划算法是递

E2]Delaney

of

B.Onthetrailoftheshadowwomanlthemystery

motion

capture口].IEEE

14—19

ComputerGraphics

and

Applications,1998,18(5)l

[3]Zhuang

T,RuiY,HuangTS,eta1.Adaptivekeyframe

extractionusingunsupervised

on

clustering[C]//Proceedingof

ImageProcessing(ICIP98)。

InternationalConference

Chicago,IL,1998,1:866-870

[4]LiuF,Zhuang

motion

index

T,Wu

F,eta1.3Dmotionretrievalwith

Vision

tree[J].

Computerand.Image

Understanding,2003。92(2/3):265—284

[5]DoyG。SullivanJ,CarlssonS.Pose—basedclusteringin

action

sequences[c]//Proceedings

Workshop

on

ofthe1stIEEE

in

3D

International

Higher—LevelKnowledge

ModelingandMotion

Analysis。Nice,2003:66-72

792

计算机辅助设计与图形学学报

2008年

[6J

Lira

IS,Thalmann

D.Key—postureextractionoutofhuman

motiondatabycurve

simplification[C]//Proceedings

ofthe23rd

Annual

International

Conferenceof

the

IEEE

Engineering

in

Medicineand

Biology

Society,Istanbul,

2001,2:1167-1169

[7]Ramer

U.An

iterative

procedure

for

the

polygonal

approximationof

plane

curves[J].Computer

Graphics

and

ImageProcessing.1972,1(3)I

244-256

[8]Shenjunxing,Sun

Shouqian,Pan

Yunhe.Key—frame

extractionfrom

motion

capture

data[J].Journal

of

Computer—AidedDesign&Computer

Graphics。2004.16

(5):719-723(inChinese)

(沈军行,孙守迁,潘云鹤.从运动捕获数据中提取关键帧[J].计算机辅助设计与图形学学报,2004,16(5):719—

723)

[9]YangTao,Xiao

Jun,Wu

Fei,eta1.Extractionofkeyframe

ofmotion

capture

databased

on

layered

curve

simplification

EJ].Journal

of

Computer—AidedDesign&Computer

Graphics,2006,18(11):1691—1697(inChinese)(杨涛,肖俊,吴飞,等.基于分层曲线简化的运动捕获数据关键帧提取[J].计算机辅助设计与图形学学报,

2006,18(11):1691-1697)

[10]AssaJ,CaspiY,Cohen—Or

D.Action

synopsis:poseselection

and

illustration[J].

ACM

Transactions

on

Graphics,2005,24(3):667—676

[11]DeMenthon

D,Kobla

V,DoermannD.Videosummarization

by

curve

simplification

It]//Proceedings

of

the

6th

ACM

InternationalConfereneeon

Multimedia,Bristol,1998:2il一

218

[12]ZhaoL,QiW,Li

Sz,etal.Key—frameextractionand

shot

retrievalusingnearest

feature

line(NFL)[c]//Proceedings

of

ACM

Workshopson

Multimedia,Los

Angeles,CA,

2000:217-220

[13]Grassia

S.Practicalparameterization

ofrotationsusing

the

exponential

map[J].Journal

of

Graphics

Tools,1998,3

(3):29-48

[14]LeeJ,ChaiJ,Reitsma

PSA,eta1.Interaetivecontrolof

avatars

animated

with

human

motion

data[J].ACM

Transactions

on

Graphics,2002,21(3):4911-500

[15]Ljung

L.Systems

identification--theory

forthe

user[M].

2nd

ed.Beijing:Tsinghua

UniversityPress,2002:97-104

[16]LiY,Wang

S,ShumH

Y.Motion

texture:a

two—level

statistical

model

for

character

motion

synthesis[J].ACM

Transactions

on

Graphics,2002。21(3):465-472

[17]Hsu

E,Pulli

K,Popovi6J.Style

translation

forhuman

motion[J].ACM

Transactions

on

Graphics,2005,24(3)2

1082-1089

[18]Ghahramani

z,HintonGE.Parameterestimation

forlineardynamical

systems[R].Toronto#University

of

Toronto。

1996

基于运动序列分割的运动捕获数据关键帧提取

作者:作者单位:

朱登明, 王兆其, Zhu Dengming, Wang Zhaoqi

朱登明,Zhu Dengming(中国科学院计算技术研究所虚拟现实技术实验室,北京100190;中国科学院计算技术研究所智能信息处理重点实验室,北京100190;中国科学院研究生院,北京100049), 王兆其,Wang Zhaoqi(中国科学院计算技术研究所虚拟现实技术实验室,北京,100190)

计算机辅助设计与图形学学报

JOURNAL OF COMPUTER-AIDED DESIGN & COMPUTER GRAPHICS2008,20(6)2次

刊名:英文刊名:年,卷(期):被引用次数:

参考文献(18条)

1.Lee J;Chai J;Reitsma P S A Interactive control of avatars animated with human motion data[外文期刊] 2002(03)

2.Grassia F S Practical parameterization of rotations using the exponential map[外文期刊] 1998(03)3.Zhao L;Qi W;Li S Z Key-frame extraction and shot retrieval using nearest feature line (NFL) 20004.杨涛;肖俊;吴飞 基于分层曲线简化的运动捕获数据关键帧提取[期刊论文]-计算机辅助设计与图形学学报2006(11)

5.沈军行;孙守迁;潘云鹤 从运动捕获数据中提取关键帧[期刊论文]-计算机辅助设计与图形学学报 2004(05)6.Ramer U An iterative procedure for the polygonal approximation of plane curves[外文期刊] 1972(03)7.Lim I S;Thalmann D Key-posture extraction out of human motion data by curve simplification 20018.Loy G;Sullivan J;Carlsson S Pose-based clustering in action sequences 2003

9.Liu F;Zhuang Y T;Wu F 3D motion retrieval with motion index tree[外文期刊] 2003(2-3)10.Zhuang Y T;Rui Y;Huang T S Adaptive key frame extraction using unsupervised clustering 199811.Delaney B On the trail of the shadow woman:the mystery of motion capture[外文期刊] 1998(05)12.Ghahramani Z;Hinton G E Parameter estimation for linear dynamical systems 199613.Hsu E;Pulli K;Popovic J Style translation for human motion[外文期刊] 2005(03)

14.Li Y;Wang T S;Shum H Y Motion texture:a two-level statistical model for character motionsynthesis[外文期刊] 2002(03)

15.Ljung L Systems identification-theory for the user 2002

16.DeMenthon D;Kobla V;Doermann D Video summarization by curve simplification 1998

17.Assa J;Caspi Y;Cohen-Or D Action synopsis:pose selection and illustration[外文期刊] 2005(03)18.Dasgupta A;Nakamura Y Making feasible walking motion of humanoid robots from human motion capturedata 1999

引证文献(3条)

1.刘贤梅.郝爱民.赵丹 基于混合遗传算法的人体运动捕获数据关键帧提取[期刊论文]-模式识别与人工智能2011(5)

2.夏利民.邓克捷 基于自分裂竞争学习算法的关键帧提取[期刊论文]-计算机工程与应用 2011(2)

3.刘更代.潘志庚.程熙.李灵.张明敏 人体运动合成中的机器学习技术综述[期刊论文]-计算机辅助设计与图形学学报 2010(9)

本文链接:/Periodical_jsjfzsjytxxxb200806018.aspx

本文来源:https://www.bwwdw.com/article/kh44.html

Top