CN108629831B - 基于参数化人体模板和惯性测量的三维人体重建方法及系统 - Google Patents

基于参数化人体模板和惯性测量的三维人体重建方法及系统 Download PDF

Info

Publication number
CN108629831B
CN108629831B CN201810317074.6A CN201810317074A CN108629831B CN 108629831 B CN108629831 B CN 108629831B CN 201810317074 A CN201810317074 A CN 201810317074A CN 108629831 B CN108629831 B CN 108629831B
Authority
CN
China
Prior art keywords
human body
depth
model
vertex
dimensional
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.)
Expired - Fee Related
Application number
CN201810317074.6A
Other languages
English (en)
Other versions
CN108629831A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201810317074.6A priority Critical patent/CN108629831B/zh
Publication of CN108629831A publication Critical patent/CN108629831A/zh
Application granted granted Critical
Publication of CN108629831B publication Critical patent/CN108629831B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30196Human being; Person

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Processing Or Creating Images (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于参数化人体模板和惯性测量的三维人体重建方法及系统,其中,方法包括:对人体进行深度图拍摄以得到单张深度图像;读取人体四肢上绑定的惯性测量单元的转向测量读数;获取三维点云和重建模型顶点及参数话人体模型顶点之间的匹配点对;根据匹配点对和转向测量读数得到能量函数和人体姿态参数;对能量函数进行求解,重建模型和参数化人体模板与三维点云进行对齐,以得到对其后的模型;优化转换关系;通过深度图更新和补全对其后的模型,并进一步优化参数化人体模板的人体形态参数。该方法求解准确鲁棒,简单易行,运行速度为实时,拥有广阔的应用前景,可以在PC机或工作站等硬件系统上快速实现。

Description

基于参数化人体模板和惯性测量的三维人体重建方法及系统
技术领域
本发明涉及计算机视觉和计算机图形学技术领域,特别涉及一种基于参数化人体模板和惯性测量的三维人体重建方法及系统。
背景技术
人体三维重建是计算机图形学和计算机视觉领域的重点问题。高质量的人体三维模型在影视娱乐、人口数据统计分析等领域有着广泛的应用前景和重要的应用价值。
然而,相关技术中,高质量人体三维模型的获取通常依靠价格昂贵的激光扫描仪或者多相机阵列系统来实现,虽然精度较高,但是也显著存在着一些缺点:第一,扫描过程中要求被人保持绝对静止,微小的移动就会导致扫描结果存在明显的误差;第二,造假昂贵,很难普及到普通民众日常生活中。第三,速度慢,重建一个三维人体模型需要较长的时间;第四,无法处理快速的运动,当人体运动速度较快时,往往由于运动跟踪的不准确导致模型重建出现错误。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种基于参数化人体模板和惯性测量的三维人体重建方法,该方法求解准确鲁棒,简单易行,运行速度快。
本发明的另一个目的在于提出一种基于参数化人体模板和惯性测量的三维人体重建系统。
为达到上述目的,本发明一方面实施例提出了一种基于参数化人体模板和惯性测量的三维人体重建方法,包括以下步骤:对人体进行深度图拍摄以得到单张深度图像;读取人体四肢上绑定的惯性测量单元的转向测量读数;将所述单张深度图像变换为三维点云,以获取所述三维点云和重建模型顶点及参数话人体模型顶点之间的匹配点对;根据所述匹配点对和所述转向测量读数得到能量函数,以共同求解所述重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数;对所述能量函数进行求解,并根据求解结果将所述重建模型和所述参数化人体模板与所述三维点云进行对齐,以得到对其后的模型;根据所述求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系;以及通过深度图更新和补全所述对其后的模型,并进一步优化参数化人体模板的人体形态参数。
本发明实施例的基于参数化人体模板和惯性测量的三维人体重建方法,通过在人体四肢上绑定多个惯性测量单元,同时利用深度相机对人体进行拍摄,从而获得深度图像和转向测量值作为系统输入信息,并基于这些信息完成对动态人体进行实时三维重建的功能。该方法求解准确鲁棒,简单易行,运行速度为实时,拥有广阔的应用前景,可以在PC机或工作站等硬件系统上快速实现。
另外,根据本发明上述实施例的基于参数化人体模板和惯性测量的三维人体重建方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数。
进一步地,在本发明的一个实施例中,所述将所述单张深度图像变换为三维点云,进一步包括:获取深度相机的内参矩阵;根据所述内参矩阵将所述单张深度图像投影到三维空间中,以生成所述三维点云。
进一步地,在本发明的一个实施例中,所述单张深度图像的投影公式为:
Figure BDA0001624135050000021
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000022
为深度相机内参矩阵。
进一步地,在本发明的一个实施例中,求解骨骼的变形矩阵的公式为:
Figure BDA0001624135050000023
其中,K是骨骼
Figure BDA0001624135050000024
的父骨骼(关节)集合,θk表示骨骼k求解出来的转角,
Figure BDA0001624135050000025
则为已知的转轴,
Figure BDA0001624135050000026
为将求解出来的骨骼转动转换为变形矩阵。
进一步地,在本发明的一个实施例中,所述每一个顶点的计算公式为:
Figure BDA0001624135050000027
Figure BDA0001624135050000028
其中,
Figure BDA0001624135050000029
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure BDA00016241350500000210
为所述变形矩阵的旋转部分,
Figure BDA0001624135050000031
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分。
进一步地,在本发明的一个实施例中,所述转换关系的计算公式为:
Figure BDA0001624135050000032
Figure BDA0001624135050000033
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure BDA0001624135050000034
Figure BDA0001624135050000035
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure BDA0001624135050000036
Figure BDA0001624135050000037
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure BDA0001624135050000038
的插值权重。
为达到上述目的,本发明另一方面实施例提出了一种基于参数化人体模板和惯性测量的三维人体重建系统,包括:深度相机,用于对人体进行深度图拍摄以得到单张深度图像;多个惯性测量单元,所述多个惯性测量单元绑定于人体四肢上,以读取转向测量读数;匹配模块,用于将所述单张深度图像变换为三维点云,以获取所述三维点云和重建模型顶点及参数话人体模型顶点之间的匹配点对;运动结算模块,用于根据所述匹配点对和所述转向测量读数得到能量函数,以共同求解所述重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数;坐标转换关系优化模块,用于对所述能量函数进行求解,并根据求解结果将所述重建模型和所述参数化人体模板与所述三维点云进行对齐,以得到对其后的模型,并根据所述求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系;以及模型更新模块,用于通过深度图更新和补全所述对其后的模型,并进一步优化参数化人体模板的人体形态参数。
本发明实施例的基于参数化人体模板和惯性测量的三维人体重建系统,通过在人体四肢上绑定多个惯性测量单元,同时利用深度相机对人体进行拍摄,从而获得深度图像和转向测量值作为系统输入信息,并基于这些信息完成对动态人体进行实时三维重建的功能。该系统具有设备简单,方便部署和可扩展强,所需的输入信息非常容易采集,并且可以实时的获得人体的动态三维模型,此外还能跟踪人体的快速大幅度肢体运动,且求解准确鲁棒,简单易行,运行速度为实时,拥有广阔的应用前景。
另外,根据本发明上述实施例的基于参数化人体模板和惯性测量的三维人体重建系统还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数;
所述单张深度图像的投影公式为:
Figure BDA0001624135050000041
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000042
为深度相机内参矩阵;
求解骨骼的变形矩阵的公式为:
Figure BDA0001624135050000043
其中,K是骨骼
Figure BDA0001624135050000044
的父骨骼(关节)集合,θk表示骨骼k求解出来的转角,
Figure BDA0001624135050000045
则为已知的转轴,
Figure BDA0001624135050000046
为将求解出来的骨骼转动转换为变形矩阵;
所述每一个顶点的计算公式为:
Figure BDA0001624135050000047
Figure BDA0001624135050000048
其中,
Figure BDA0001624135050000049
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure BDA00016241350500000410
为所述变形矩阵的旋转部分,
Figure BDA00016241350500000411
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分;
所述转换关系的计算公式为:
Figure BDA00016241350500000412
Figure BDA00016241350500000413
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure BDA00016241350500000414
Figure BDA00016241350500000415
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure BDA00016241350500000416
Figure BDA00016241350500000417
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure BDA00016241350500000418
的插值权重。
进一步地,在本发明的一个实施例中,所述匹配模块还用于获取深度相机的内参矩阵,并根据所述内参矩阵将所述单张深度图像投影到三维空间中,以生成所述三维点云。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的基于参数化人体模板和惯性测量的三维人体重建方法流程图;
图2为根据本发明一个实施例的基于参数化人体模板和惯性测量的三维人体重建系统的结构示意图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的基于参数化人体模板和惯性测量的三维人体重建方法及系统,首先将参照附图描述根据本发明实施例提出的基于参数化人体模板和惯性测量的三维人体重建方法。
图1为根据本发明一个实施例的基于参数化人体模板和惯性测量的三维人体重建方法流程图。
如图1所示,该基于参数化人体模板和惯性测量的三维人体重建方法包括以下步骤:
在步骤S101中,对人体进行深度图拍摄以得到单张深度图像。
具体而言,使用深度相机对动态人体进行拍摄,获得连续的单张深度图像序列。
在步骤S102中,读取人体四肢上绑定的惯性测量单元的转向测量读数。
具体而言,通过无线通信获取人体四肢上绑定的惯性测量单元在各个时刻的转向测量读数。
在步骤S103中,将单张深度图像变换为三维点云,以获取三维点云和重建模型顶点及参数话人体模型顶点之间的匹配点对。
在本发明的一个实施例中,将深度图像投影到三维空间中变换为一组三维点云包括:
获取深度相机的内参矩阵;
根据内参矩阵将深度图投影到三维空间中变换为一组三维点云。其中,变换的公式为:
Figure BDA0001624135050000051
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000052
为深度相机内参矩阵;
在获取匹配点对方面,使用相机投影公式将三维模型的顶点投影到深度图像上以获得匹配点对。
单张深度图像的投影公式为:
Figure BDA0001624135050000061
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000062
为深度相机内参矩阵。
在步骤S104中,根据匹配点对和转向测量读数得到能量函数,以共同求解重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数。
具体而言,能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数。
进一步地,在本发明的一个实施例中:
Figure BDA0001624135050000063
Figure BDA0001624135050000064
Figure BDA0001624135050000065
Figure BDA0001624135050000066
其中,惯性单元数据项EIMU要求求解得到的骨骼旋转量与测量数值尽可能相近,深度图数据项Edata保证经过非刚性运动驱动后的重建模型和经过人体姿态驱动的参数化人体模型与从深度图获得的三维点云尽可能的对齐;局部刚性运动约束项Ereg可以在使模型整体受局部刚性约束运动的同时保证较大幅度的合理的非刚性运动也能被很好的解算出来,从而使模型更精确的与三维点云对齐;人体参数化人体模型姿态和非刚性运动一致性约束项Ebind用于保证解算出来的人体参数化人体模型姿态和非刚性运动尽可能的一致,从而可以保证最终解算出来的非刚性运动即符合人体骨架动力学模型,又充分的与从深度图中获得的三维点云对齐;人体姿态运动正则项Eprior使用高斯混合模型来约束人体姿态的正确性,不正常的解算姿态会导致该项能量较大,直至姿态解算正确。
在步骤S105中,对能量函数进行求解,并根据求解结果将重建模型和参数化人体模板与三维点云进行对齐,以得到对其后的模型。
具体而言,单张深度图像的投影公式为:
Figure BDA0001624135050000067
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000068
为深度相机内参矩阵。
进一步地,在本发明的一个实施例中,共同求解重建模型上每一个顶点的非刚性运动位置变换参数和人体参数化人体模型姿态参数。最终求解获得的信息为每一个三维模型顶点的变换矩阵和人体参数化人体模型姿态参数,即每个骨骼的单独的变换矩阵。为了实现快速线性求解的要求,本发明实施例的方法对利用指数映射方法对变形方程做如下近似:
Figure BDA0001624135050000071
其中,
Figure BDA0001624135050000072
为截至上一帧的模型顶点vi的累积变换矩阵,为已知量;I为四维单位阵;
其中,
Figure BDA0001624135050000073
Figure BDA0001624135050000074
即上一帧变换后的模型顶点,则经过变换有:
Figure BDA0001624135050000075
对于每个顶点,要求解的未知参数即为六维变换参数x=(v1,v2,v3,wx,wy,wz)T。骨骼运动的线性化方式与非刚性运动相同。
在步骤S106中,根据求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系。
具体而言,求解骨骼的变形矩阵的公式为:
Figure BDA0001624135050000076
其中,K是骨骼
Figure BDA0001624135050000077
的父骨骼(关节)集合,θk表示骨骼k求解出来的转角,
Figure BDA0001624135050000078
则为已知的转轴,
Figure BDA0001624135050000079
为将求解出来的骨骼转动转换为变形矩阵。
进一步地,每一个顶点的计算公式为:
Figure BDA00016241350500000710
Figure BDA00016241350500000711
其中,
Figure BDA00016241350500000712
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure BDA00016241350500000713
为该变形矩阵的旋转部分,
Figure BDA00016241350500000714
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分。
进一步地,转换关系的计算公式为:
Figure BDA00016241350500000715
Figure BDA0001624135050000081
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure BDA0001624135050000082
Figure BDA0001624135050000083
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure BDA0001624135050000084
Figure BDA0001624135050000085
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure BDA0001624135050000086
的插值权重。
在本发明的一个实施例中,利用累积的旋转来优化当前的转换关系,通过要求从第一帧到当前帧的骨骼转动求解结果与惯性测量单元测量值尽可能保持一致,来求解更加精确的转换关系。计算公式为:
Figure BDA0001624135050000087
Figure BDA0001624135050000088
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure BDA0001624135050000089
Figure BDA00016241350500000810
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure BDA00016241350500000811
Figure BDA00016241350500000812
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,用于融合两个旋转矩阵,w和ω分别为RI2C
Figure BDA00016241350500000813
的插值权重,通过将原有的转换关系RI2C
Figure BDA00016241350500000814
加权融合。
在步骤S107中,通过深度图更新和补全对其后的模型,并进一步优化参数化人体模板的人体形态参数。
进一步地,在本发明的一个实施例中,使用深度图像对对齐后的三维模型进行更新和补全,将新获得的深度信息融合到三维模型中,更新三维模型表面顶点位置或为三维模型增加新的顶点,使其更符合当前深度图像的表达。
本发明实施例的基于参数化人体模板和惯性测量的三维人体重建方法,通过在人体四肢上绑定多个惯性测量单元,同时利用深度相机对人体进行拍摄,从而获得深度图像和转向测量值作为系统输入信息,并基于这些信息完成对动态人体进行实时三维重建的功能。该方法求解准确鲁棒,简单易行,运行速度为实时,拥有广阔的应用前景,可以在PC机或工作站等硬件系统上快速实现。
其次参照附图描述根据本发明实施例提出的基于参数化人体模板和惯性测量的三维人体重建系统。
图2是本发明一个实施例的基于参数化人体模板和惯性测量的三维人体重建系统的结构示意图。
如图2所示,该基于参数化人体模板和惯性测量的三维人体重建系统10包括:深度相机100、多个惯性测量单元200、匹配模块300、运动结算模块400、坐标转换关系优化模块500和模型更新模块。
其中,深度相机100,用于对人体进行深度图拍摄以得到单张深度图像。多个惯性测量单元200,多个惯性测量单元绑定于人体四肢上,以读取转向测量读数。匹配模块300,用于将单张深度图像变换为三维点云,以获取三维点云和重建模型顶点及参数话人体模型顶点之间的匹配点对。运动结算模块400,用于根据匹配点对和转向测量读数得到能量函数,以共同求解重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数。坐标转换关系优化模块500,用于对能量函数进行求解,并根据求解结果将重建模型和参数化人体模板与三维点云进行对齐,以得到对其后的模型,并根据求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系。模型更新模块600,用于通过深度图更新和补全对其后的模型,并进一步优化参数化人体模板的人体形态参数。该基于参数化人体模板和惯性测量的三维人体重建系统10
进一步地,在本发明的一个实施例中,能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数;
单张深度图像的投影公式为:
Figure BDA0001624135050000091
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure BDA0001624135050000092
为深度相机内参矩阵;
求解骨骼的变形矩阵的公式为:
Figure BDA0001624135050000093
其中,K是骨骼
Figure BDA0001624135050000094
的父骨骼(关节)集合,θk表示骨骼k求解出来的转角,
Figure BDA0001624135050000095
则为已知的转轴,
Figure BDA0001624135050000096
为将求解出来的骨骼转动转换为变形矩阵;
每一个顶点的计算公式为:
Figure BDA0001624135050000097
Figure BDA0001624135050000098
其中
Figure BDA0001624135050000099
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure BDA00016241350500000910
为该变形矩阵的旋转部分;
Figure BDA00016241350500000911
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分;
转换关系的计算公式为:
Figure BDA0001624135050000101
Figure BDA0001624135050000102
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure BDA0001624135050000103
Figure BDA0001624135050000104
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure BDA0001624135050000105
Figure BDA0001624135050000106
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure BDA0001624135050000107
的插值权重。
进一步地,在本发明的一个实施例中,匹配模块300还用于获取深度相机的内参矩阵,并根据内参矩阵将单张深度图像投影到三维空间中,以生成三维点云。
需要说明的是,前述对基于参数化人体模板和惯性测量的三维人体重建方法实施例的解释说明也适用于该实施例的系统,此处不再赘述。
本发明实施例的基于参数化人体模板和惯性测量的三维人体重建系统,通过在人体四肢上绑定多个惯性测量单元,同时利用深度相机对人体进行拍摄,从而获得深度图像和转向测量值作为系统输入信息,并基于这些信息完成对动态人体进行实时三维重建的功能。该系统具有设备简单,方便部署和可扩展强,所需的输入信息非常容易采集,并且可以实时的获得人体的动态三维模型,此外还能跟踪人体的快速大幅度肢体运动,且求解准确鲁棒,简单易行,运行速度为实时,拥有广阔的应用前景。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (8)

1.一种基于参数化人体模板和惯性测量的三维人体重建方法,其特征在于,包括以下步骤:
对人体进行深度图拍摄以得到单张深度图像;
读取人体四肢上绑定的惯性测量单元的转向测量读数;
将所述单张深度图像变换为三维点云,以获取所述三维点云和重建模型顶点及参数化人体模型顶点之间的匹配点对;
根据所述匹配点对和所述转向测量读数得到能量函数,以共同求解所述重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数,其中,所述每一个顶点的计算公式为:
Figure FDA0002770752090000011
Figure FDA0002770752090000012
其中,
Figure FDA0002770752090000013
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure FDA0002770752090000014
为该变形矩阵的旋转部分;
Figure FDA0002770752090000015
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分;
对所述能量函数进行求解,并根据求解结果将所述重建模型和所述参数化人体模板与所述三维点云进行对齐,以得到对齐后的模型;
根据所述求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系,其中,所述转换关系的计算公式为:
Figure FDA0002770752090000016
Figure FDA0002770752090000017
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure FDA0002770752090000018
Figure FDA0002770752090000019
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure FDA00027707520900000110
Figure FDA00027707520900000111
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure FDA00027707520900000112
的插值权重;以及
通过深度图更新和补全所述对齐后的模型,并进一步优化参数化人体模板的人体形态参数。
2.根据权利要求1所述的基于参数化人体模板和惯性测量的三维人体重建方法,其特征在于,所述能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数。
3.根据权利要求1所述的基于参数化人体模板和惯性测量的三维人体重建方法,其特征在于,所述将所述单张深度图像变换为三维点云,进一步包括:
获取深度相机的内参矩阵;
根据所述内参矩阵将所述单张深度图像投影到三维空间中,以生成所述三维点云。
4.根据权利要求3所述的基于参数化人体模板和惯性测量的三维人体重建方法,其特征在于,所述单张深度图像的投影公式为:
Figure FDA0002770752090000021
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure FDA0002770752090000022
为深度相机内参矩阵。
5.根据权利要求1所述的基于参数化人体模板和惯性测量的三维人体重建方法,其特征在于,求解骨骼的变形矩阵的公式为:
Figure FDA0002770752090000023
其中,K是骨骼
Figure FDA0002770752090000024
的父骨骼集合,θk表示骨骼k求解出来的转角,
Figure FDA0002770752090000025
则为已知的转轴,
Figure FDA0002770752090000026
为将求解出来的骨骼转动转换为变形矩阵。
6.一种基于参数化人体模板和惯性测量的三维人体重建系统,其特征在于,包括:
深度相机,用于对人体进行深度图拍摄以得到单张深度图像;
多个惯性测量单元,所述多个惯性测量单元绑定于人体四肢上,以读取转向测量读数;
匹配模块,用于将所述单张深度图像变换为三维点云,以获取所述三维点云和重建模型顶点及参数化人体模型顶点之间的匹配点对;
运动结算模块,用于根据所述匹配点对和所述转向测量读数得到能量函数,以共同求解所述重建模型上每一个顶点的非刚性运动位置变换参数和参数化人体模板的人体姿态参数,其中,所述每一个顶点的计算公式为:
Figure FDA0002770752090000027
Figure FDA0002770752090000028
其中,
Figure FDA0002770752090000029
为作用于顶点vi的变形矩阵,包括旋转和平移两部分,
Figure FDA00027707520900000210
为该变形矩阵的旋转部分,
Figure FDA00027707520900000211
为对顶点vi有驱动作用的骨骼的集合,αi,j为第j个骨骼对第i个模型顶点的驱动作用的权重,Tbj为第j个骨骼自身的运动变形矩阵,rot(Tbj)为该变形矩阵的旋转部分;
坐标转换关系优化模块,用于对所述能量函数进行求解,并根据求解结果将所述重建模型和所述参数化人体模板与所述三维点云进行对齐,以得到对齐后的模型,并根据所述求解结果优化相机坐标与大地坐标系之间的转换关系及惯性测量单元坐标系与对应骨骼坐标系之间的转换关系,其中,所述转换关系的计算公式为:
Figure FDA0002770752090000031
Figure FDA0002770752090000032
其中,RI2C为当前的相机坐标系到大地坐标系的转换矩阵,
Figure FDA0002770752090000033
Figure FDA0002770752090000034
分别为当前帧和第一帧得到的传感器在大地坐标系下转向的测量值,
Figure FDA0002770752090000035
Figure FDA0002770752090000036
分别为当前帧和第一帧求解出来的骨骼的旋转,SLERP{·}为球面插值函数,w和ω分别为RI2C
Figure FDA0002770752090000037
的插值权重;以及
模型更新模块,用于通过深度图更新和补全所述对齐后的模型,并进一步优化参数化人体模板的人体形态参数。
7.根据权利要求6所述的基于参数化人体模板和惯性测量的三维人体重建系统,其特征在于,所述能量函数为:
Emot=λIMUEIMUdepthEdepthbindEbindregEregpriorEprior
其中,Emot为运动求解总能量项,EIMU为IMU数据对应的数据项,Edepth为深度图数据对应的数据项,包含非刚性运动跟踪的数据项和参数化人体模板人体姿态求解的数据项,Ebind为参数化人体姿态和非刚性运动一致性约束项,Ereg为局部刚性运动约束项,Epri为人体姿态运动的正则项,λIMU、λdepth、λreg、λbind和λpri分别为对应各个约束项的权重系数;
所述单张深度图像的投影公式为:
Figure FDA0002770752090000038
其中,u,v为像素坐标,d(u,v)为深度图像上像素(u,v)位置处的深度值,
Figure FDA0002770752090000039
为深度相机内参矩阵;
求解骨骼的变形矩阵的公式为:
Figure FDA00027707520900000310
其中,K是骨骼
Figure FDA00027707520900000311
的父骨骼集合,θk表示骨骼k求解出来的转角,
Figure FDA00027707520900000312
则为已知的转轴,
Figure FDA00027707520900000313
为将求解出来的骨骼转动转换为变形矩阵。
8.根据权利要求6所述的基于参数化人体模板和惯性测量的三维人体重建系统,其特征在于,所述匹配模块还用于获取深度相机的内参矩阵,并根据所述内参矩阵将所述单张深度图像投影到三维空间中,以生成所述三维点云。
CN201810317074.6A 2018-04-10 2018-04-10 基于参数化人体模板和惯性测量的三维人体重建方法及系统 Expired - Fee Related CN108629831B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810317074.6A CN108629831B (zh) 2018-04-10 2018-04-10 基于参数化人体模板和惯性测量的三维人体重建方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810317074.6A CN108629831B (zh) 2018-04-10 2018-04-10 基于参数化人体模板和惯性测量的三维人体重建方法及系统

Publications (2)

Publication Number Publication Date
CN108629831A CN108629831A (zh) 2018-10-09
CN108629831B true CN108629831B (zh) 2021-03-12

Family

ID=63705039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810317074.6A Expired - Fee Related CN108629831B (zh) 2018-04-10 2018-04-10 基于参数化人体模板和惯性测量的三维人体重建方法及系统

Country Status (1)

Country Link
CN (1) CN108629831B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110097639B (zh) * 2019-03-18 2023-04-18 北京工业大学 一种三维人体姿态估计方法
CN110148209B (zh) * 2019-04-30 2023-09-15 深圳市重投华讯太赫兹科技有限公司 人体模型生成方法、图像处理设备及具有存储功能的装置
CN110307788B (zh) * 2019-07-12 2021-05-11 清华大学 多深度相机人体参数测量方法及装置
CN110827342B (zh) * 2019-10-21 2023-06-02 中国科学院自动化研究所 三维人体模型重建方法及存储设备、控制设备
CN112883757B (zh) * 2019-11-29 2023-03-24 北京航空航天大学 生成跟踪姿态结果的方法
CN111179418B (zh) * 2019-12-31 2024-02-27 中国人民解放军陆军工程大学 一种无需用户裸体的三维人体测量方法和装置
CN111612887B (zh) * 2020-04-30 2021-11-09 北京的卢深视科技有限公司 一种人体测量方法和装置
CN112037319A (zh) * 2020-08-19 2020-12-04 上海佑久健康科技有限公司 人体测量方法、系统及计算机可读存储介质
CN112132971B (zh) * 2020-09-08 2024-04-23 合肥的卢深视科技有限公司 三维人体建模方法、装置、电子设备和存储介质
CN112233223A (zh) * 2020-09-29 2021-01-15 深圳市易尚展示股份有限公司 基于三维点云的人体参数化模型自动变形方法和装置
CN112330813A (zh) * 2020-11-12 2021-02-05 首都师范大学 一种基于单目深度摄像头的着装下三维人体模型重建方法
CN112365589B (zh) * 2020-12-01 2024-04-26 东方梦幻虚拟现实科技有限公司 一种虚拟三维场景展示方法、装置及系统
CN112837406B (zh) * 2021-01-11 2023-03-14 聚好看科技股份有限公司 一种三维重建方法、装置及系统
CN113689539B (zh) * 2021-07-06 2024-04-19 清华大学 基于隐式光流场的动态场景实时三维重建方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102800103A (zh) * 2012-06-18 2012-11-28 清华大学 基于多视角的深度相机的无标记动作捕捉方法及装置
CN102842148A (zh) * 2012-07-10 2012-12-26 清华大学 一种无标记运动捕捉及场景重建方法及装置
CN103198523A (zh) * 2013-04-26 2013-07-10 清华大学 一种基于多深度图的非刚体三维重建方法及系统
KR101471761B1 (ko) * 2013-05-30 2014-12-10 가톨릭대학교 산학협력단 3d 콘텐츠를 이용한 안구 능력 강화 장치

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103712604B (zh) * 2013-12-20 2016-04-06 清华大学深圳研究生院 一种光学多目标三维空间定位方法及系统
CN107490364A (zh) * 2017-09-01 2017-12-19 中国科学院长春光学精密机械与物理研究所 一种大角度倾斜成像航空相机对地目标定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102800103A (zh) * 2012-06-18 2012-11-28 清华大学 基于多视角的深度相机的无标记动作捕捉方法及装置
CN102842148A (zh) * 2012-07-10 2012-12-26 清华大学 一种无标记运动捕捉及场景重建方法及装置
CN103198523A (zh) * 2013-04-26 2013-07-10 清华大学 一种基于多深度图的非刚体三维重建方法及系统
KR101471761B1 (ko) * 2013-05-30 2014-12-10 가톨릭대학교 산학협력단 3d 콘텐츠를 이용한 안구 능력 강화 장치

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维人体模型的生成技术研究;陈寅;《中国博士学位论文全文数据库 信息科技辑》;20170215(第2期);I138-198,正文正文第27页、第41-76页 *
基于惯性传感器的人体运动姿态三维重构;蹇彪;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140415(第4期);I138-1180,摘要、正文第58页 *

Also Published As

Publication number Publication date
CN108629831A (zh) 2018-10-09

Similar Documents

Publication Publication Date Title
CN108629831B (zh) 基于参数化人体模板和惯性测量的三维人体重建方法及系统
CN108665537B (zh) 联合优化人体体态与外观模型的三维重建方法及系统
CN110296691B (zh) 融合imu标定的双目立体视觉测量方法与系统
CN108711185B (zh) 联合刚性运动和非刚性形变的三维重建方法及装置
US20190371003A1 (en) Monocular vision tracking method, apparatus and non-volatile computer-readable storage medium
AU2011312140B2 (en) Rapid 3D modeling
CN112894832A (zh) 三维建模方法、装置、电子设备和存储介质
CN108053437B (zh) 基于体态的三维模型获取方法及装置
CN110617814A (zh) 单目视觉和惯性传感器融合的远距离测距系统及方法
CN110969665B (zh) 一种外参标定方法、装置、系统及机器人
JP2013539147A5 (zh)
JP2017182695A (ja) 情報処理プログラム、情報処理方法および情報処理装置
CN111739137A (zh) 一种三维姿态估计数据集的生成方法
CN113205603A (zh) 一种基于旋转台的三维点云拼接重建方法
US20230351625A1 (en) A method for measuring the topography of an environment
CN111060006A (zh) 一种基于三维模型的视点规划方法
CN110517209A (zh) 数据处理方法、装置、系统以及计算机可读存储介质
CN111932637B (zh) 一种车身相机外参自适应标定方法和装置
CN103778610A (zh) 一种星载线阵传感器垂轨摆扫影像的几何预处理方法
CN116597080A (zh) 一种多源空间数据的完整场景3d精细模型构建系统及方法
WO2023133125A1 (en) System and method for mapping the skin
CN110307788B (zh) 多深度相机人体参数测量方法及装置
JP2005063012A (ja) 全方位カメラ運動と3次元情報の復元方法とその装置及びプログラム並びにこれを記録した記録媒体
CN107705267B (zh) 基于控制矢量的光学卫星影像几何校正方法
JP2004108836A (ja) 撮像装置の方位角計算方法及び装置、撮像装置の姿勢検出装置、撮像装置の傾きセンサ、コンピュータプログラム、並びに三次元モデル構成装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210312

CF01 Termination of patent right due to non-payment of annual fee