融合深度图和三维模型的人体运动捕捉方法

1.融合深度图和三维模型的人体运动捕捉方法，其特征在于：
由以下步骤实现：
步骤一：采集人体动作的深度信息，去除运动目标背景，获取完整的人体动作深度信息；
步骤二：将提取的人体动作深度信息信息转化成人体的三维点云信息，对人体进行三维重建，获得人体动作的三维模型；
步骤三：重复步骤一和步骤二，基于大量采集的人体动作深度信息，建立人体动作的三维模型数据库M＝{Y1,Y2,…，Yn},M表示三维模型数据库，Yn表示第n个三维模型；
步骤四：根据人体结构及人体动作的三维模型构建出人体动作骨架，建立人体动作的骨架数据库G＝{S1,S2,…Sn},G表示人体骨架数据库，其中Sn表示第n个人体骨架数据，其中人体动作的骨架数据与三维模型数据一一对应；
步骤五：提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配，按相似性排序输出人体动作骨架，相似性距离值最小的骨架作为最优骨架，以最优骨架作为运动捕获结果；
步骤一中，基于采集的人体动作深度信息，去除运动目标背景，获取完整的人体动作深度信息的具体步骤为：
(1)采集的一个人体动作深度图用F(x,y,dp)表示，其中，x、y分别为像素坐标系下的横坐标和纵坐标，dp为深度信息；假设基于深度信息分割背景区域和目标区域的阈值为：
T＝(maxDepthVulue+minDepthVulue)/2
其中,maxDepthVulue和minDepthVulue分别为图像深度值的最大和最小值，将T记录在T0中，根据阈值T，将F(x,y,dp)分为背景区域和目标区域，分别求出背景区域和目标区域的平均深度值du1和du2；
(2)重新计算T＝(du1+du2)/2，判断T与T0是否相等，如果不相等，那么将T记录在T0里，重复上述步骤，直到T＝T0成立，终止算法；把最后得到的T作为最佳阈值对F(x,y,dp)进行分割，去除背景，得到完整的人体动作深度信息d0；
步骤二中，将提取的人体动作深度信息转化成人体的三维点云信息，对人体进行三维重建，获得人体动作的三维模型的具体步骤为：
(1)将提取出的深度信息d0归一化，假设深度值的个数为N，计算出深度值中的最大值max d0(k)和最小值min d0(k)，归一化后的深度值为：
z(n)＝(d0(n)-min d0(k))/(max d0(k)-min d0(k))k，n＝1，2，...，N(2)令Z＝z(n),世界坐标系是以摄像机为原点的，经过标定后摄像机可以看做是理想的成像模型，根据简单的相似三角形变换原理：
从中可以计算得到世界坐标系的X、Y值，其中，x、y分别为像素坐标系的横坐标和纵坐标，f为摄像机的焦距；最终可以得到人体动作的三维点云信息(X,Y,Z)；
步骤四中，根据人体结构及人体动作的三维模型构建出人体动作骨架的具体步骤为：
(1)在深度图中将身体躯干近似视为四边形，记为Q；
Q上方的两个顶点就是肩关节点a1和a2的位置；
颈关节点b为a1、a2连线的中点；
沿b点向上，最上端就是头部节点c的位置；
Q下方的两个顶点为髋关节d1和d2的位置；
胯关节点e在d1、d2的中点位置；
(2)确定手关节点和肘关节位置从a1、a2处开始进行搜索，若手臂伸直，根据上下手臂长度比例确定肘关节点f1、f2和手关节点g1、g2的位置；若手臂弯曲，在肘关节点的位置有拐点，拐点的位置为f1、f2，沿f1、f2继续搜索至终点为手关节点g1、g2的位置；
同理，根据手关节点和肘关节位置的确定方法可以确定出膝关节点h1、h2和脚关节点i1、i2的位置；
(3)将各个关节点坐标按人体结构顺序用直线相连接，就得到了人体动作骨架；
步骤五中，提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配的具体步骤为：
(1)将三维模型表示的人体动作用层次聚类算法进行聚类：
①假设一个参数m表示最终聚类簇的数目，将每个输入的数据作为单独的一个数据簇D0，求出与每个D相邻的最近的簇Dx，数据簇间的距离可以看做是数据簇中心间的距离；
②将距离最近的两个数据簇Dp和Dq进行合并，产生的新的数据簇Dn，再计算Dn与其它簇的距离，如果当前数据簇的数目大于m，那么就转到②继续进行数据簇的合并，否则算法结束，得到聚类结果D＝{u1,u2,…,um}，其中u表示每个类的值；
(2)计算待识别人体动作序列各个动作帧与三维模型数据库中的人体动作序列各个动作帧之间的距离：
①假设待识别的人体动作和数据库中任意一个三维模型聚类结果分别表示为：
Di＝{u′1，u′2，…，u′m}
Dj＝{v1，v2，…，vm}
从聚类树的根节点开始，进行深度优先搜索直至叶节点；
②在遍历搜索的过程中对Di和Dj中各个类的对象计算距离，将各个类的距离求和，得到两颗聚类树的欧式距离：
(3)用正则时间规整的方法计算出最优匹配路径，得到最优匹配序列：
①假设待识别的三维模型动作序列为 三维模型数据库中的
一组动作序列为 m和n分别表示两个动作序列的长度，两个序列
中各个动作帧之间的距离用上述方法可以计算得到；
②因为m和n可能不相等，所以其对应的匹配关系可能有很多种，利用正则时间规整的方法，计算出序列间的最小距离J，同时获得最佳匹配路径Px和Py：
这里，Wx＝W(px)∈{0,1}m×l和Wy＝W(py)∈{0,1}n×l是两个二进制的规整矩阵，对每一步的t∈{1:l}有wpt,t＝1，其他情况时wpt,t＝0，其中l是正则时间规整算法自动选择的需要匹配的动作帧的数目，l≥max(m,n)；φ(·)是正则项：
Px和Py表示所有可能的匹配路径，需要满足的约束条件为：
Ψ＝{{Px,Py}|Px∈(1:m)l,Py∈(1:n)l
边界：
单调性：
连续性：
和 是线性变化矩阵，d≤min(dx,dy)，需满足的约束条件为：
其中，λ∈[0,1]；
③重复以上步骤，利用正则时间规整方法，计算三维模型数据库中每一组动作序列Y分别与X的最小距离J及最佳匹配路径Px和Py；然后求出各个J的最小值和其对应的Px和Py，从而获得对应的最佳匹配序列。

融合深度图和三维模型的人体运动捕捉方法\n技术领域\n[0001] 本发明属于多媒体信息检索技术领域，具体涉及一种融合深度图和三维模型的人体运动捕捉方法。\n背景技术\n[0002] 人体运动捕获技术是多媒体信息检索领域的热点问题，特别是在影视动画、游戏等的发展中，具有广泛应用前景，国内外许多研究机构正致力于这一方向的研究。近年来，随着运动捕获技术的快速发展，三维影视动画、游戏、新一代人交互等领域的兴起，非常多复杂和逼真的人体动作需要快速的被捕获应用，需要一种快速有效的方法进行人体运动的捕获。目前已提出的光学式的运动捕获方法，主要基于计算机视觉原理，通过对目标上特定光点的监视和跟踪来完成运动捕捉任务。但上述的运动捕获方法有几点不足：\n[0003] (1)光学式的运动捕获在使用时需要在表演者身上粘贴标记点，并需要表演者穿上特制的表演用服装，交互起来并不方便。当运动复杂时，不同部位的标记点有可能发生混淆、遮挡，产生错误结果，这时需要人工干预后处理过程。\n[0004] (2)虽然它可以捕捉实时运动，但后处理(包括标记点的识别、跟踪、空间坐标的计算)的工作量较大，对于表演场地的光照、反射情况有一定的要求，装置定标也较为烦琐。\n发明内容\n[0005] 本发明的目的是提供一种融合深度图和三维模型的人体运动捕捉方法，能有效克服现有运动捕获方法中动作幅度受限制、捕捉效果失真度高和误差较大的技术缺陷。\n[0006] 本发明所采用的技术方案是：\n[0007] 融合深度图和三维模型的人体运动捕捉方法，其特征在于：\n[0008] 由以下步骤实现：\n[0009] 步骤一：采集人体动作的深度信息，去除运动目标背景，获取完整的人体动作深度信息；\n[0010] 步骤二：将提取的人体动作深度信息信息转化成人体的三维点云信息，对人体进行三维重建，获得人体动作的三维模型；\n[0011] 步骤三：重复步骤一和步骤二，基于大量采集的人体动作深度信息，建立人体动作的三维模型数据库M＝{Y1,Y2,…，Yn},M表示三维模型数据库，Yn表示第n个三维模型；\n[0012] 步骤四：根据人体结构及人体动作的三维模型构建出人体动作骨架，建立人体动作的骨架数据库G＝{S1,S2,…Sn},G表示人体骨架数据库，其中Sn表示第n个人体骨架数据，其中人体动作的骨架数据与三维模型数据一一对应；\n[0013] 步骤五：提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配，按相似性排序输出人体动作骨架，相似性距离值最小的骨架作为最优骨架，以最优骨架作为运动捕获结果；\n[0014] 步骤一中，基于采集的人体动作深度信息，去除运动目标背景，获取完整的人体动作深度信息的具体步骤为：\n[0015] (1)采集的一个人体动作深度图用F(x,y,dp)表示，其中，x、y分别为像素坐标系下的横坐标和纵坐标，dp为深度信息；假设基于深度信息分割背景区域和目标区域的阈值为：\n[0016] T＝(maxDepthVulue+minDepthVulue)/2\n[0017] 其中,maxDepthVulue和minDepthVulue分别为图像深度值的最大和最小值，将T记录在T0中，根据阈值T，将F(x,y,dp)分为背景区域和目标区域，分别求出背景区域和目标区域的平均深度值du1和du2；\n[0018] (2)重新计算T＝(du1+du2)/2，判断T与T0是否相等，如果不相等，那么将T记录在T0里，重复上述步骤，直到T＝T0成立，终止算法；把最后得到的T作为最佳阈值对F(x,y,dp)进行分割，去除背景，得到完整的人体动作深度信息d0；\n[0019] 步骤二中，将提取的人体动作深度信息转化成人体的三维点云信息，对人体进行三维重建，获得人体动作的三维模型的具体步骤为：\n[0020] (1)将提取出的深度信息d0归一化，假设深度值的个数为N，计算出深度值中的最大值max d0(k)和最小值min d0(k)，归一化后的深度值为：\n[0021] z(n)＝(d0(n)-min d0(k))/(max d0(k)-min d0(k))k，n＝1，2，...，N[0022] (2)令Z＝z(n),世界坐标系是以摄像机为原点的，经过标定后摄像机可以看做是理想的成像模型，根据简单的相似三角形变换原理：\n[0023]\n[0024]\n[0025] 从中可以计算得到世界坐标系的X、Y值，其中，x、y分别为像素坐标系的横坐标和纵坐标，f为摄像机的焦距；最终可以得到人体动作的三维点云信息(X,Y,Z)；\n[0026] 步骤四中，根据人体结构及人体动作的三维模型构建出人体动作骨架的具体步骤为：\n[0027] (1)在深度图中将身体躯干近似视为四边形，记为Q；\n[0028] Q上方的两个顶点就是肩关节点a1和a2的位置；\n[0029] 颈关节点b为a1、a2连线的中点；\n[0030] 沿b点向上，最上端就是头部节点c的位置；\n[0031] Q下方的两个顶点为髋关节d1和d2的位置；\n[0032] 胯关节点e在d1、d2的中点位置；\n[0033] (2)确定手关节点和肘关节位置从a1、a2处开始进行搜索，若手臂伸直，根据上下手臂长度比例确定肘关节点f1、f2和手关节点g1、g2的位置；若手臂弯曲，在肘关节点的位置有拐点，拐点的位置为f1、f2，沿f1、f2继续搜索至终点为手关节点g1、g2的位置；\n[0034] 同理，根据手关节点和肘关节位置的确定方法可以确定出膝关节点h1、h2和脚关节点i1、i2的位置；\n[0035] (3)将各个关节点坐标按人体结构顺序用直线相连接，就得到了人体动作骨架；\n[0036] 步骤五中，提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配的具体步骤为：\n[0037] (1)将三维模型表示的人体动作用层次聚类算法进行聚类：\n[0038] ①假设一个参数m表示最终聚类簇的数目，将每个输入的数据作为单独的一个数据簇D0，求出与每个D相邻的最近的簇Dx，数据簇间的距离可以看做是数据簇中心间的距离；\n[0039] ②将距离最近的两个数据簇Dp和Dq进行合并，产生的新的数据簇Dn，再计算Dn与其它簇的距离，如果当前数据簇的数目大于m，那么就转到②继续进行数据簇的合并，否则算法结束，得到聚类结果D＝{u1,u2,…,um}，其中u表示每个类的值；\n[0040] (2)计算待识别人体动作序列各个动作帧与三维模型数据库中的人体动作序列各个动作帧之间的距离：\n[0041] ①假设待识别的人体动作和数据库中任意一个三维模型聚类结果分别表示为：\n[0042] Di＝{u′1，u′2，…，u′m}\n[0043] Dj＝{v1，v2，…，vm}\n[0044] 从聚类树的根节点开始，进行深度优先搜索直至叶节点；\n[0045] ②在遍历搜索的过程中对Di和Dj中各个类的对象计算距离，将各个类的距离求和，得到两颗聚类树的欧式距离：\n[0046]\n[0047] (3)用正则时间规整的方法计算出最优匹配路径，得到最优匹配序列：\n[0048] ①假设待识别的三维模型动作序列为 三维模型数据库\n中的一组动作序列为 m和n分别表示两个动作序列的长度，两个\n序列中各个动作帧之间的距离用上述方法可以计算得到；\n[0049] ②因为m和n可能不相等，所以其对应的匹配关系可能有很多种，利用正则时间规整的方法，计算出序列间的最小距离J，同时获得最佳匹配路径Px和Py：\n[0050]\n[0051] 这里，Wx＝W(px)∈{0,1}m×l和Wy＝W(py)∈{0,1}n×l是两个二进制的规整矩阵，对每一步的t∈{1:l}有 其他情况时 其中l是正则时间规整算法自动选择的需要匹配的动作帧的数目，l≥max(m,n)；φ(·)是正则项：\n[0052]\n[0053] Px和Py表示所有可能的匹配路径，需要满足的约束条件为：\n[0054] Ψ＝{{Px,Py}|Px∈(1:m)l,Py∈(1:n)l\n[0055] 边界：\n[0056] 单调性：\n[0057] 连续性：\n[0058] 和 是线性变化矩阵，d≤min(dx,dy)，需满足的约束条件\n为：\n[0059]\n[0060]\n[0061] 其中，λ∈[0,1]；\n[0062] ③重复以上步骤，利用正则时间规整方法，计算三维模型数据库中每一组动作序列Y分别与X的最小距离J及最佳匹配路径Px和Py；然后求出各个J的最小值和其对应的Px和Py，从而获得对应的最佳匹配序列。\n[0063] 本发明具有以下优点：\n[0064] (1)本发明所涉及的方法使用相关设备采集深度信息，表演者无需在人体上安装传感器或添加标记点，活动的幅度大，对于表演场地要求不多，方便且容易实现。\n[0065] (2)本发明针对以往无标记点的运动捕获技术主要是基于视频方法，拍摄出二维图像序列达到运动捕获的目的，本发明融合了深度信息，解决了在没有深度信息的二维图像中身体的自遮挡导致身体部位信息缺失的问题，保证了信息获取的完整性和可靠性。\n[0066] (3)本发明使用层次聚类法，数据簇间距离和规则的相似度容易定义，将数据量较大的三维点云数据缩小为几个类，以便于后续利用深度优先搜索方法计算两帧三维模型的欧式距离，减少了计算量。\n[0067] (4)本发明在运动序列匹配上使用正则时间规整(Canonical Time Warping,CTW)的方法，该方法具有算法鲁棒的优点，并且即使测试序列模式与数据库中参考序列模式的时间尺度不能完全一致，它仍能较好地完成测试序列与参考序列之间的模式匹配。CTW是基于动态规划的方法，步骤五第(3)步的②中公式是CTW方法的目标函数，目标函数中加有正则项和两个二进制规整矩阵，为目标函数提供了唯一解的可能性，提高了两个序列匹配的精度，使匹配的时间大幅降，保证了运动捕获的速度和精度。\n附图说明\n[0068] 图1为本发明方法的总体流程图。\n[0069] 图2为本发明步骤一的具体流程框图。\n[0070] 图3为本发明步骤五的具体流程框图。\n具体实施方式\n[0071] 下面结合具体实施方式对本发明进行详细的说明。\n[0072] 本发明所涉及的融合深度图和三维模型的人体运动捕捉方法，由以下步骤实现：\n[0073] 步骤一：应用相关设备采集人体动作的深度信息，去除运动目标背景，获取完整的人体动作深度信息。\n[0074] 基于采集的人体动作深度信息，去除运动目标背景，获取完整的人体动作深度信息的具体步骤为：\n[0075] (1)采集的一个人体动作深度图用F(x,y,dp)表示，其中，x、y分别为像素坐标系下的横坐标和纵坐标，dp为深度信息；假设基于深度信息分割背景区域和目标区域的阈值为：\n[0076] T＝(maxDepthVulue+minDepthVulue)/2\n[0077] 其中,maxDepthVulue和minDepthVulue分别为图像深度值的最大和最小值，将T记录在T0中，根据阈值T，将F(x,y,dp)分为背景区域和目标区域，分别求出背景区域和目标区域的平均深度值du1和du2；\n[0078] (2)重新计算T＝(du1+du2)/2，判断T与T0是否相等，如果不相等，那么将T记录在T0里，重复上述步骤，直到T＝T0成立，终止算法；把最后得到的T作为最佳阈值对F(x,y,dp)进行分割，去除背景，得到完整的人体动作深度信息d0。\n[0079] 步骤二：将提取的人体动作深度信息转化成人体的三维点云信息，对人体进行三维重建，获得人体动作的三维模型，具体步骤为：\n[0080] (1)将提取出的深度信息d0归一化，假设深度值的个数为N，计算出深度值中的最大值最大值max d0(k)和最小值min d0(k)，归一化后的深度值为：\n[0081] z(n)＝(d0(n)-mind0(k))/(maxd0(k)-mind0(k))k，n＝1，2，...，N[0082] (2)令Z＝z(n)，世界坐标系是以摄像机为原点的，经过标定后摄像机可以看做是理想的成像模型，根据简单的相似三角形变换原理：\n[0083]\n[0084]\n[0085] 从中可以计算得到世界坐标系的X、Y值，其中，x、y分别为像素坐标系的横坐标和纵坐标，f为摄像机的焦距。最终可以得到人体动作的三维点云信息(X,Y,Z)。\n[0086] 步骤三：重复步骤一和步骤二，基于大量采集的尽可能多的人体动作深度信息，建立人体动作的三维模型数据库M＝{Y1,Y2,…，Yn},M表示三维模型数据库，Yn表示第n个三维模型。\n[0087] 步骤四：根据人体结构及人体动作的三维模型构建出人体动作骨架，建立人体动作的骨架数据库G＝{S1,S2,…Sn},G表示人体骨架数据库，其中Sn表示第n个人体骨架数据，其中人体动作的骨架数据与三维模型数据一一对应。\n[0088] 根据人体结构及人体动作的三维模型构建出人体动作骨架的具体步骤为：\n[0089] (1)在深度图中将身体躯干近似视为四边形，记为Q；\n[0090] Q上方的两个顶点就是肩关节点a1和a2的位置；\n[0091] 颈关节点b为a1、a2连线的中点；\n[0092] 沿b点向上，最上端就是头部节点c的位置；\n[0093] Q下方的两个顶点为髋关节d1和d2的位置；\n[0094] 胯关节点e在d1、d2的中点位置；\n[0095] (2)确定手关节点和肘关节位置从a1、a2处开始进行搜索，若手臂伸直，根据上下手臂长度比例确定肘关节点f1、f2和手关节点g1、g2的位置；若手臂弯曲，在肘关节点的位置有拐点，拐点的位置为f1、f2，沿f1、f2继续搜索至终点为手关节点g1、g2的位置；\n[0096] 同理，根据手关节点和肘关节位置的确定方法可以确定出膝关节点h1、h2和脚关节点i1、i2的位置；\n[0097] (3)将各个关节点坐标按人体结构顺序用直线相连接，就得到了人体动作骨架。\n[0098] 步骤五：提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配，按相似性排序输出人体动作骨架，相似性距离值最小的骨架作为最优骨架，以最优骨架作为运动捕获结果。\n[0099] 提取待识别人体动作的深度信息，基于深度信息构建待识别人体动作的三维模型，然后与三维模型数据库中的人体动作进行相似性匹配的具体步骤为：\n[0100] (1)将三维模型表示的人体动作用层次聚类算法(THA)进行聚类：\n[0101] ①假设一个参数m表示最终聚类簇的数目，将每个输入的数据作为单独的一个数据簇D0，求出与每个D相邻的最近的簇Dx，数据簇间的距离可以看做是数据簇中心间的距离；\n[0102] ②将距离最近的两个数据簇Dp和Dq进行合并，产生的新的数据簇Dn，再计算Dn与其它簇的距离，如果当前数据簇的数目大于m，那么就转到②继续进行数据簇的合并，否则算法结束，得到聚类结果D＝{u1,u2,…,um}，其中u表示每个类的值；\n[0103] (2)计算待识别人体动作序列各个动作帧与三维模型数据库中的人体动作序列各个动作帧之间的距离：\n[0104] ①假设待识别的人体动作和数据库中任意一个三维模型聚类结果分别表示为：\n[0105] Di＝{u′1，u′2，…，u′m}\n[0106] Dj＝{v1，v2，…，vm}\n[0107] 从聚类树的根节点开始，进行深度优先搜索直至叶节点；\n[0108] ②在遍历搜索的过程中对Di和Dj中各个类的对象计算距离，将各个类的距离求和，得到两颗聚类树的欧式距离：\n[0109]\n[0110] (3)用正则时间规整(Canonical Time Warping,CTW)的方法计算出最优匹配路径，得到最优匹配序列：\n[0111] ①假设待识别的三维模型动作序列为 三维模型数据库\n中的一组动作序列为 m和n分别表示两个动作序列的长度，两个\n序列中各个动作帧之间的距离用上述方法可以计算得到；\n[0112] ②因为m和n可能不相等，所以其对应的匹配关系可能有很多种，利用正则时间规整的方法，计算出序列间的最小距离J，同时获得最佳匹配路径Px和Py：\n[0113]\n[0114] 这里，Wx＝W(px)∈{0,1}m×l和Wy＝W(py)∈{0,1}n×l是两个二进制的规整矩阵，对每一步的t∈{1:l}有 其他情况时 其中l是正则时间规整算法自动选择的需要匹配的动作帧的数目，l≥max(m,n)。φ(·)是正则项：\n[0115]\n[0116] Px和Py表示所有可能的匹配路径，需要满足的约束条件为：\n[0117] Ψ＝{{Px,Py}|Px∈(1:m)l,Py∈(1:n)l\n[0118] 边界：\n[0119] 单调性：\n[0120] 连续性：\n[0121] 和 是线性变化矩阵，需满足的约\n束条件为：\n[0122]\n[0123]\n[0124] 其中，λ∈[0,1]；\n[0125] ③重复以上步骤，利用正则时间规整方法，计算三维模型数据库中每一组动作序列Y分别与X的最小距离J及最佳匹配路径Px和Py。然后求出各个J的最小值和其对应的Px和Py，从而获得对应的最佳匹配序列。\n[0126] 本发明的内容不限于实施例所列举，本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换，均为本发明的权利要求所涵盖。