CN103116901B - 基于运动特征的人体运动插值计算方法 - Google Patents

基于运动特征的人体运动插值计算方法 Download PDF

Info

Publication number
CN103116901B
CN103116901B CN201310031870.0A CN201310031870A CN103116901B CN 103116901 B CN103116901 B CN 103116901B CN 201310031870 A CN201310031870 A CN 201310031870A CN 103116901 B CN103116901 B CN 103116901B
Authority
CN
China
Prior art keywords
transition
motion
interpolation
data
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310031870.0A
Other languages
English (en)
Other versions
CN103116901A (zh
Inventor
魏小鹏
张强
姚一
周东生
李晴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University
Original Assignee
Dalian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University filed Critical Dalian University
Priority to CN201310031870.0A priority Critical patent/CN103116901B/zh
Publication of CN103116901A publication Critical patent/CN103116901A/zh
Application granted granted Critical
Publication of CN103116901B publication Critical patent/CN103116901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Processing Or Creating Images (AREA)

Abstract

本发明公开了基于运动特征的人体运动插值计算方法,属于图像处理技术领域。本发明的技术方案具体包括:1.对高维人体运动数据进行降维处理。2.提取运动特征。3.计算过渡窗口长度。4.生成过渡动画。通过基于运动特征的人体运动插值窗口长度计算,着重解决在人体运动插值中,过渡窗口选择复杂,通用性差,生成过渡动画跳帧、漏帧的问题。

Description

基于运动特征的人体运动插值计算方法
技术领域
本发明涉及基于运动特征的人体运动插值计算方法,属于图像处理技术领域。
背景技术
近年来,随着计算机软硬件技术的进步,计算机动画技术得到了飞速的发展,计算机动画是指采用图形与图像的处理技术,以实体造型和真实感显示技术为基础,借助于编程或动画制作软件生成一系列的景物画面。其涉及到图像处理技术、运动控制原理、视频技术和艺术等众多领域,以独特的特点逐步成为一个多种学科和技术综合的领域。其中,运动捕捉技术的不断发展,人们可以利用捕捉设备抓取到的数据生成更多样化且复杂的虚拟人运动。但是,由于运动捕捉设备昂贵且受限于捕捉的外部条件,不可能对所需的人体运动数据每次都进行捕捉,运动数据重用技术和基于运动图的人体动画合成技术由此产生。将捕获的人体运动数据根据运动类型分类保存构建运动图,合成新的人体运动数据时,只需遍历运动图并结合插值技术即可由现有的人体运动数据片段合成新的所需人体运动数据。因此,基于运动捕捉数据的运动合成技术研究是计算机动画技术重要的研究领域之一。
自2002年基于运动图的人体运动合成方法产生后,在人体运动数据合成领域,基于运动图的人体运动合成方法已成为目前的主要方法。如附图2所示,运动图的基本方法是将运动捕捉数据按运动类型分段,每帧为一个节点,然后计算每一个节点间的相似性,满足所设阈值的节点间构建为边,最后形成运动图,进行人体运动合成的过程就是在运动图上搜索所需路径的过程。尽管经过多年的发展,但其仍存在很多问题,例如运动图中点、边的定义,过渡点的选取,运过渡动画插值计算等。
发明内容
为了克服上述的不足,本发明提出了基于运动特征的人体运动插值计算方法。本方法通过使用流形学习的方法对高维运动数据降维,并根据降维后的低维特征曲线完成人体运动特征的提取,最后基于提取的运动特征计算过渡窗口长度,确定过渡区间,使用线性插值算法完成过渡动画的构建。通过基于运动特征的人体运动插值窗口长度计算,着重解决在人体运动插值中,过渡窗口选择复杂,通用性差,生成过渡动画跳帧、漏帧的问题。
本发明的技术方案是:通过PCA降维算法对高维人体运动数据进行降维处理,使用降维后的低维数据,绘制低维特征曲线。根据低维特征曲线提取人体运动特征。在运动特征的基础上选取过渡窗口长度,确定过渡区间,最后使用线性插值算法完成过渡动画的构建。
本发明的技术方案具体包括以下步骤:
1.对高维人体运动数据进行降维处理
使用PCA、流形学习算法,对高维人体运动数据进行降维处理,得到原始运动序列的低维流形结构,根据运动类型的不同,绘制与之匹配的低维特征曲线。
2.提取运动特征
根据生成的人体运动数据低维特征曲线,在低维特征曲线转折处就是高维空间中人体运动姿态的改变处,我们可以将低维特征曲线上的这些点提取出来。
3.计算过渡窗口长度
根据所绘制的低维特征曲线,在需要进行插值构建过渡动画的2段运动的跳转点位置,分别向前或向后取离跳转点最近的一个特征曲线转折点,跳转点和求得的特征曲线转折点之间即为过渡窗口,此窗口内包含的帧数即为过渡区间长度。
4.生成过渡动画
将经过基于运动特征提取的过渡窗口内数据作为插值的原始数据集,假设其分别为数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn],其中,Fi表示数据集中的一帧数据,其可表示为:
F(i)={T,R1,R2,…,Rk}
在上面的公式中,T为当前人体的根关节点,R为其余关节点,k为当前人体骨架结构包含的除根关节点外的关节点总数。对数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn]中每一组相对应的帧Fi,对当前帧Fi上相对应的关节点Rk使用线性插值算法得到插值点,使用这些插值点合成一帧新的数据,以此类推,最终完成过渡动画的构建。
本发明原理:在遍历运动图的过程中,不满足所设阈值的跳转点间,需要使用插值技术进行插补动画的构建,插值结果的好坏直接影响到最终合成动画的质量。进行两段运动间插值计算过渡动画时,过渡长度的选择非常重要,选取合适的过渡窗口长度,是进行插值计算的基础和难点。虽然理论上过渡窗口越长,则生成的过渡动画越平滑,但是在实际应用中,过长的过渡窗口在计算插值动画时,会消耗大量的时间,同时,根据运动特征的不同,过长的过渡窗口会使生成的过渡动画与人们的需求相违。同样,过短的过渡窗口,会使生成的插值动画生硬,无法达到平滑过渡的要求。合适的过渡长度可以使合成的结果平滑、自然。
本发明与现有技术相比具有以下优点:
与基于固定过渡窗口长度的方法相比,在走到跑过渡中,右侧小臂和脚腕等关节处,由于过渡窗口小,无法充分进行插值计算,产生了跳帧、漏帧的现象,以致生成的插值动画不够平滑。同样在由跑到跳的过渡中,人体的上半身下行过快,插值生成的插补动画不足,产生的过渡动画不自然。通过本发明所提的基于运动特征提取的插值窗口计算,过渡窗口大小分别为22帧,经过线性插值计算,平滑自然的由走过渡到跑再由跑过渡到跳。由此对比可见本发明所提算法在规律及不规律两类运动类型上,均优于固定窗口长度的插值计算。
与基于测地距离的过渡窗口计算方法相比,在走到跑过渡中,上半身过渡平滑,但脚步出现了滑步现象,以致生成的插值动画不够平滑,需要再进行去除滑步的处理。而在在由跑到跳的过渡中,与采用的固定窗口长度存在的问题类似,人体的上半身下行过快,插值生成的插补动画不足,产生跳帧现像。可见本发明提出的基于运动特征的插值窗口长度计算算法优于基于测地距离的过渡窗口长度计算算法。
附图说明
图1本发明算法流程图。
图2标准运动图计算原理。
图3人体运动数据低维特征曲线示意图。
图4过渡动画计算方法示意图。
图5本发明生成过渡动画结果示意图。
图6基于固定过渡长度生成的过渡动画示意图1。
图7基于固定过渡长度生成的过渡动画示意图2。
图8基于测地距离过渡长度生成的过渡动画示意图1。
图9基于测地距离过渡长度生成的过渡动画示意图2。
具体实施方式
本发明的实施示例是在以本发明技术方案为前提下进行实施的,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述实施例。实例选取一组走到跑再到跳的运动进行验证,运动数据低维特征曲线如附图3所示:表示的是原始m维空间中的运动数据映射到k维空间中后得到的低维特征曲线;过渡动画计算方法如图4所示:在motionA的第i帧处和motionB的第j帧处为跳转点,分别向前、向后取一定长度的过渡窗口;本实施例生成过渡动画结果如图5所示。
具体实施步骤为:
步骤1:提取运动特征。使用PCA流形学习算法,对高维人体运动数据进行降维处理,并绘制3维运动数据特征曲线,此特征曲线即为低维空间对高维人体运动姿态的映射。
(1)对于一个长度为n帧的运动序列计算平均值Mavg
M a v g = 1 n Σ i = 1 n M i
(2)定义一个差值di=Mi–Mavg,i=1,2,...,n,可得到一个矩阵D:
D=[d1,…,dn]
(3)计算协方差矩阵:
C = 1 n Σ i = 1 n d i d i T = 1 T DD T ,
因为C为半正定矩阵,对其进行对角化,得:
C=VΛVT
其中,Λ=diag(λ12,…,λn),VVT=1,假设C的秩为p,那么C有p个非零特征值,记为
λ1≥λ2≥…≥λp>0
我们按照特征值的大小依次排列,取前k个特征值所对应的特征向量v1,v2,…,vk作为我们的投影方向,那么对于任何一帧运动Mi=x,计算yi=xTvw,i=1,2,…,k,我们就可以用[y1,y2,…,yk]T来表示Mi,这样就可以将数据从n维降到k维。
(4)根据降维后得到的低维数据,绘制低维特征曲线。
步骤2:计算过渡窗口长度。根据所绘制的低维特征曲线,在需要进行插值构建过渡动画的2段运动的跳转点位置,分别向前或向后取离跳转点最近的一个特征曲线转折点,跳转点和求得的特征曲线转折点之间即为过渡窗口,此窗口内包含的帧数即为过渡区间长度。如图4所示。
步骤3:生成过渡动画。将经过基于运动特征提取的过渡窗口内数据作为插值的原始数据集,假设其分别为数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn],其中,Fi表示数据集中的一帧数据,其可表示为:
F(i)={T,R1,R2,…,Rk},
在上面的公式中,T为当前人体的根关节点,R为其余关节点,k为当前人体骨架结构包含的除根关节点外的关节点总数。对数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn]中每一组相对应的帧Fi,对当前帧Fi上相对应的关节点Rk使用线性插值算法得到插值点,使用这些插值点合成一帧新的数据,最终完成过渡动画的构建。
附图5为使用本发明所提算法对一组由走到跑再到跳的人体运动数据合成结果。如图5所示:经过降维、流形学习、运动特征提取、过渡窗口长度确定、插值生成过渡动画这一系列计算,最终生成了一段平滑、自然的运动。
与基于固定过渡窗口长度的方法相比,在图6所示的走到跑过渡中,右侧小臂和脚腕等关节处,由于过渡窗口小,无法充分进行插值计算,产生了跳帧、漏帧的现象,以致生成的插值动画不够平滑。同样的,在图7中,在由跑到跳的过渡中,人体的上半身下行过快,插值生成的插补动画不足,产生的过渡动画不自然。与此相对的图4可见,通过本发明所提的基于运动特征提取的插值窗口计算,过渡窗口大小分别为22帧,经过线性插值计算,平滑自然的由走过渡到跑再由跑过渡到跳。由此对比可见本发明所提算法在规律及不规律两类运动类型上,均优于固定窗口长度的插值计算。
与基于测地距离的过渡窗口计算方法相比,在图8所示的走到跑过渡中,上半身过渡平滑,但脚步出现了滑步现象,以致生成的插值动画不够平滑,需要再进行去除滑步的处理。而在图9中,在由跑到跳的过渡中,与采用的固定窗口长度存在的问题类似,人体的上半身下行过快,插值生成的插补动画不足,产生跳帧现像。可见本发明提出的基于运动特征的插值窗口长度计算算法优于基于测地距离的过渡窗口长度计算算法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种基于运动特征的人体运动插值计算方法,其特征在于:包括如下步骤:
步骤一、对高维人体运动数据进行降维处理
(1)对于一个长度为n帧的运动序列计算平均值Mavg
(2)定义一个差值di=Mi–Mavg,i=1,2,...,n,可得到一个矩阵D:
D=[d1,…,dn]
(3)计算协方差矩阵:
因为C为半正定矩阵,对其进行对角化,得:
C=VΛVT
其中,Λ=diag(λ12,…,λn),VVT=1,假设C的秩为p,那么C有p个非零特征值,记为
λ1≥λ2≥…≥λp>0,
按照特征值的大小依次排列,取前k个特征值所对应的特征向量v1,v2,…,vk作为我们的投影方向,那么对于任何一帧运动Mi=x,计算yi=xTvw,w=1,2,…,k,我们就可以用[y1,y2,…,yk]T来表示Mi,这样就可以将数据从n维降到k维;
步骤二、提取运动特征
使用步骤一中降维得到的原始数据在k维空间的低维映射绘制人体运动数据在k维空间下的特征曲线,根据生成的人体运动数据低维特征曲线与原始运动状态对比可以得到在低维特征曲线转折处就是高维空间中人体运动姿态的改变处,低维特征曲线上这些转折点提取出来;
步骤三、计算过渡窗口长度
根据所绘制的低维特征曲线和提取的转折点,在需要进行插值构建过渡动画的2段运动的Mi和Mi+1的跳转点位置,分别向前或向后取离跳转点最近的一个特征曲线转折点,跳转点和求得的特征曲线转折点之间即为过渡窗口,此窗口内包含的帧数即为过渡区间长度;
步骤四、生成过渡动画
将经过基于运动特征提取的过渡窗口内数据作为插值的原始数据集,假设其分别为数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn],其中,Fi表示数据集中的一帧数据,其可表示为:
F(i)={T,R1,R2,…,Rk},
在上面的公式中,T为当前人体的根关节点,R为其余关节点,k为当前人体骨架结构包含的除根关节点外的关节点总数;对数据集A[F1,F2,…,Fn]和数据集B[F1,F2,…,Fn]中每一组相对应的帧Fi,对当前帧Fi上相对应的关节点Rk使用线性插值算法得到插值点,使用这些插值点合成一帧新的数据,最终完成过渡动画的构建。
CN201310031870.0A 2013-01-28 2013-01-28 基于运动特征的人体运动插值计算方法 Active CN103116901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310031870.0A CN103116901B (zh) 2013-01-28 2013-01-28 基于运动特征的人体运动插值计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310031870.0A CN103116901B (zh) 2013-01-28 2013-01-28 基于运动特征的人体运动插值计算方法

Publications (2)

Publication Number Publication Date
CN103116901A CN103116901A (zh) 2013-05-22
CN103116901B true CN103116901B (zh) 2016-03-30

Family

ID=48415265

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310031870.0A Active CN103116901B (zh) 2013-01-28 2013-01-28 基于运动特征的人体运动插值计算方法

Country Status (1)

Country Link
CN (1) CN103116901B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109068053B (zh) * 2018-07-27 2020-12-04 香港乐蜜有限公司 图像特效展示方法、装置和电子设备
CN109166409B (zh) * 2018-10-10 2021-02-12 长沙千博信息技术有限公司 一种手语转换方法及装置
CN113658300B (zh) * 2021-08-18 2023-05-30 北京百度网讯科技有限公司 动画播放方法、装置、电子设备以及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122391A (zh) * 2010-12-13 2011-07-13 中国人民解放军国防科学技术大学 一种运动捕获数据自动分割方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122391A (zh) * 2010-12-13 2011-07-13 中国人民解放军国防科学技术大学 一种运动捕获数据自动分割方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于PCA和多约简SVM的多级说话人辨识;刘雪燕;《计算机应用》;20080131;全文 *
基于运动捕捉数据的人体运动编辑与合成技术研究;杜会军;《中国优秀硕士学位论文全文数据库.信息科技辑》;20100228;论文第26、30、35、51页 *
基于运动捕获数据的人体运动编辑技术研究;李婷;《中国优秀硕士学位论文全文数据库.信息科技辑》;20100531;论文第18-24页 *

Also Published As

Publication number Publication date
CN103116901A (zh) 2013-05-22

Similar Documents

Publication Publication Date Title
CN103003846B (zh) 关节区域显示装置、关节区域检测装置、关节区域归属度计算装置、关节状区域归属度计算装置以及关节区域显示方法
Bera et al. Efficient trajectory extraction and parameter learning for data-driven crowd simulation.
CN111553968B (zh) 一种三维人体重构动画的方法
Saini et al. A review on particle swarm optimization algorithm and its variants to human motion tracking
CN100440257C (zh) 一种虚拟人群运动的三维可视化方法
CN105512621A (zh) 一种基于Kinect的羽毛球动作指导系统
CN110427799A (zh) 基于生成对抗网络的人手深度图像数据增强方法
CN110473266A (zh) 一种基于姿态指导的保留源场景人物动作视频生成方法
CN114821640A (zh) 基于多流多尺度膨胀时空图卷积网络的骨架动作识别方法
CN106709933B (zh) 一种基于非监督学习的运动估计方法
CN110473284A (zh) 一种基于深度学习的运动物体三维模型重建方法
CN110458046A (zh) 一种基于关节点提取的人体运动轨迹分析方法
CN103116901B (zh) 基于运动特征的人体运动插值计算方法
CN102999923A (zh) 基于自适应阈值的运动捕捉数据的关键帧提取方法
CN111462274A (zh) 一种基于smpl模型的人体图像合成方法及系统
CN114170353A (zh) 一种基于神经网络的多条件控制的舞蹈生成方法及系统
CN105740815B (zh) 一种基于深度递归分层条件随机场的人体行为识别方法
CN103021009B (zh) 基于非线性流形学习的运动图过渡点选取方法
CN116740240A (zh) 款式多样的实时服装动画生成方法
Wang et al. Robust face tracking via collaboration of generic and specific models
Zhou et al. Hierarchical learning recurrent neural networks for 3D motion synthesis
CN103699214A (zh) 一种基于三维自然手势的三维跟踪和交互方法
Ammar et al. Comparative Study of latest CNN based Optical Flow Estimation
CN112949419A (zh) 一种基于肢体层次结构的动作识别方法
Guo et al. 3D object detection and tracking based on streaming data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant