CN113052138B - 一种舞蹈与运动动作的智能对比矫正的方法 - Google Patents

一种舞蹈与运动动作的智能对比矫正的方法 Download PDF

Info

Publication number
CN113052138B
CN113052138B CN202110448369.9A CN202110448369A CN113052138B CN 113052138 B CN113052138 B CN 113052138B CN 202110448369 A CN202110448369 A CN 202110448369A CN 113052138 B CN113052138 B CN 113052138B
Authority
CN
China
Prior art keywords
dance
teacher
video
student
frame
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
CN202110448369.9A
Other languages
English (en)
Other versions
CN113052138A (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.)
Guanghai Art Technology Innovation Shenzhen Co ltd
Original Assignee
Guanghai Art Technology Innovation Shenzhen Co ltd
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 Guanghai Art Technology Innovation Shenzhen Co ltd filed Critical Guanghai Art Technology Innovation Shenzhen Co ltd
Priority to CN202110448369.9A priority Critical patent/CN113052138B/zh
Publication of CN113052138A publication Critical patent/CN113052138A/zh
Application granted granted Critical
Publication of CN113052138B publication Critical patent/CN113052138B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/20Movements or behaviour, e.g. gesture recognition
    • G06V40/23Recognition of whole body movements, e.g. for sport training
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/40Scenes; Scene-specific elements in video content
    • G06V20/46Extracting features or characteristics from the video content, e.g. video fingerprints, representative shots or key frames
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/40Scenes; Scene-specific elements in video content
    • G06V20/48Matching video sequences
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/103Static body considered as a whole, e.g. static pedestrian or occupant recognition

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Human Computer Interaction (AREA)
  • Social Psychology (AREA)
  • Health & Medical Sciences (AREA)
  • Psychiatry (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种舞蹈与运动动作的智能对比矫正的方法,包括以下步骤:采用人体姿态估计模型AlphaPose提取学生舞蹈视频与老师舞蹈视频中的人体关节点;将学生舞蹈视频与老师舞蹈视频进行动作对齐,分为有音乐情形和无音乐情形的动作对齐,并提取学生舞蹈视频与老师舞蹈视频的关键帧;根据提取的人体关节点、学生舞蹈视频的关键帧与老师舞蹈视频的关键帧,对动作对齐后的学生舞蹈视频与老师舞蹈视频进行动作相似对比;根据动作相似对比结果,找出差异区间,再对学生动作进行可视化矫正,得到智能对比矫正结果。

Description

一种舞蹈与运动动作的智能对比矫正的方法
技术领域
本发明属于计算机视频分析和人体动作识别领域,具体涉及一种舞蹈与运动动作的智能对比矫正的方法。
背景技术
近年来,随着科学技术的快速发展,人们对生活水平要求的提高,深度学习得到快速的发展。研究人员希望研发出更加高端的技术以提高机器对人类的理解能力。正是这样的研究热潮下,运动估计与识别算法得到快速的发展。
运动估计与行为技术,当前处于理论研究并逐步应用于现实生活的阶段。它一方面用于行为预判功能的实现,通过对当前的行为分析,推测目标在下一时间轴的若干行为;另一方面用于实现高质量的人机交互体验,如用于评价和指导人的动作姿态是否准确,这对于体育、健身、舞蹈等运动领域的学习和练习有着重大的意义。
目前的人体运动识别系统中,都仅限于对人体的连续性动作进行检测和相关行为的一些识别,还未有过对两组或多组图像或视频进行姿态动作智能对比评价和智能矫正的高精准度的系统性解决方案。
发明内容
针对现有技术中的上述不足,本发明提供的一种舞蹈与运动动作的智能对比矫正的方法解决了现有技术中存在的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种舞蹈与运动动作的智能对比矫正的方法,包括以下步骤:
S1、采用人体姿态估计模型AlphaPose提取学生舞蹈视频与老师舞蹈视频中的人体关节点;
S2、将学生舞蹈视频与老师舞蹈视频进行动作对齐,并提取学生舞蹈视频与老师舞蹈视频的关键帧;
S3、根据提取的人体关节点、学生舞蹈视频的关键帧与老师舞蹈视频的关键帧,对动作对齐后的学生舞蹈视频与老师舞蹈视频进行动作相似对比;
S4、根据动作相似对比结果,找出动作差异区间,对学生动作进行可视化矫正,得到智能对比矫正结果。
进一步地,所述步骤S2中将学生舞蹈视频与老师舞蹈视频进行动作对齐包括有音乐情况的对齐和无音乐情况的对齐。
进一步地,所述有音乐情况的对齐具体为:
S2.1.1、提取学生舞蹈视频与老师舞蹈视频的音频,并提取音频中前N帧没有音乐的时间,并将前N帧时间中的噪音强度取平均值,以平均值构建估计噪声频谱;
S2.1.2、根据包含噪声的音频频谱Ps(w)以及估计噪声频谱Pn(w),获取频谱差值D(w)为:
D(w)=Ps(w)-αPn(w)
其中,α表示过减因子参数;
S2.1.3、将频谱差值D(w)中的负值进行修改,得到纯净音频Ps'(w)为:
其中,β表示频谱下限阈值参数;
S2.1.4、设定伴乐音频和纯净音频采样率分别为Fs1和Fs2,并对伴乐音频进行采样,获取伴乐音频波动的时间点为s1
S2.1.5、获取伴乐音频和纯净音频之间的延迟值d,并根据延迟值d、采样率Fs1、采样率Fs2以及时间点s1,获取纯净音频波动的时间点s2为:
S2.1.6、根据纯净音频波动的时间点s2的获取方法,分别获取学生舞蹈视频与老师舞蹈视频中伴乐的起点,并提取起点后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
进一步地,所述步骤S2.1.2中过减因子参数α具体为:
其中,SNR表示视频帧的信噪比。
进一步地,所述无音乐情况的对齐具体:
S2.2.1、设定计数器T=1,并设定动作相似距离阈值为K,并将两帧之间的Fréchet距离作为两视频帧之间的舞蹈动作相似距离;
S2.2.2、将视频第T帧设定为关键帧;
S2.2.3、判断视频第T帧与第T+1帧之间的舞蹈动作相似距离是否小阈值K,若是,则令T的计数值加一,并重复步骤S2.2.3,否则令T的计数值加一,并返回步骤S2.2.2;
S2.2.4、根据步骤S2.2.2-S2.2.3所述的方法,遍历视频的每一帧,提取视频的关键帧;
S2.2.5、根据步骤S2.2.1-S2.2.4所述的方法,分别提取学生舞蹈视频和老师舞蹈视频关键帧;
S2.2.6、获取学生舞蹈视频关键帧与老师舞蹈视频中起始动作帧的舞蹈动作相似距离,并将舞蹈动作相似距离最小的关键帧作为学生舞蹈视频的起始动作帧;
S2.2.7、分别提取学生舞蹈视频与老师舞蹈视频中起始动作帧后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
进一步地,所述步骤S3具体为:
S3.1、分别提取学生和老师对应关键帧中的人体关节点;
S3.2、骨架大小归一化;
S3.3、将学生舞蹈视频关键帧和老师舞蹈视频关键帧逐帧进行骨架匹配;
S3.4、判断学生与老师骨架匹配是否正确,若是,则进入步骤S3.5,否则结束动作相似对比流程;
S3.5、分别提取学生与老师的手部关节点,并进行手部关节点匹配,完成动作相似对比。
进一步地,所述步骤S3.2具体为:
S3.2.1、分别读取学生和老师对应关键帧中人体关节点的位置坐标,所述人体关节点包括颈部节点、头部中心节点、两个肩部节点以及两个腰部节点;
S3.2.2、根据人体关节点对应的位置坐标,分别获取学生对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为Xs1、Xs2、Xs3、Xs4和Xs5
S3.2.3、根据人体关节点对应的位置坐标,分别获取老师对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为Xt1、Xt2、Xt3、Xt4和Xs5
S3.2.4、使用距离Xsi除以距离Xti,得到距离比例Ri,i=1,2,3,4,5;
S3.2.5、根据距离比例Ri,求取比例平均值为R;
S3.2.6、根据比例平均值为R,将学生骨架缩放至老师骨架大小,完成骨架归一化。
进一步地,所述步骤S3.3中骨架匹配的具体方法为:
S3.3.1、根据学生和老师颈部节点的位置坐标,获取水平距离dis_x和垂直距离dis_y;
S3.3.2、通过比例平均值为R和老师颈部节点的位置坐标获取老师关节点的位置坐标调节值为:
V_Xi=(R-1)*(Xi-tch_x)
V_Yi=(R-1)*(Yi-tch_y)
其中,V_Xi表示第i个关节点的x位置坐标调节值,V_Yi表示第i个关节点的y位置坐标调节值,Xi表示第i个关节点调整前的x坐标,Yi表示第i个关节点调整前的y坐标,tch_x表示老师颈部节点的x坐标,tch_y表示老师颈部节点的y坐标;
S3.3.3、将老师每个关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成骨架匹配。
进一步地,所述步骤S3.5具体为:
S3.5.1、分别读取学生和老师对应关键帧中手部关节点的位置坐标;
S3.5.2、根据比例平均值为R,将学生手部关节点之间的距离进行缩放,使之与老师手部关节点之间的距离相同;
S3.5.3、将老师每个手部关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的手部关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成手部关节点匹配。
进一步地,所述步骤S4中对学生动作进行可视化矫正包括骨架迁移可视化矫正和视频拼接可视化矫正;
所述骨架迁移可视化矫正具体为:根据动作相似对比结果,以学生不匹配的关键帧为中心,向前后扩展一定时长的视频帧,将老师和学生的颈部关节点对齐,以颈部关节点为基准,将老师骨架迁移至学生舞蹈视频中,得到学生舞姿与老师骨架的对比结果,根据对比结果对学生动作进行矫正;
所述视频拼接可视化矫正具体为:通过视频对齐后,将老师和学生在同一律动下的视频帧一一拼接,得到的新帧集合拼接成一个视频,根据找出的动作差异区间,对照老师的标准动作进行矫正。
本发明的有益效果为:
(1)本发明提供了一种舞蹈与运动动作的智能对比矫正的方法,通过在有音乐和无音乐两种情形下的视频对齐技术,找出学习者和教学者的舞姿与动作的差异,进而对学习者的舞姿与动作进行智能可视化矫正。
(2)本发明提取人体运动动作的各个关节点及其运动的轨迹数据,对人体运动姿态的各个关节点的数据进行求方差,在程序中设定一个方差阈值,当所求得的人体姿态的各个关节点数据方差大于设定的方差阈值的时候,则评价该动作是有差异的,在找出所有差异动作后,再根据教学视频的总时长形成时间轴,区别显示出这段视频中学习者的所有动作差异的视频区间。
(3)本发明提供了对于舞蹈与运动动作的两种可视化智能矫正方式,对应学习者的差异视频区间,一是将老师的所有骨架关节点归一化后平移至对应的学生舞蹈视频的关键帧中,观察动作之间的差异;二是根据需要回放有差异部分的或全部的教学者与学习者的拼接视频,使学习者对比不同,矫正自已的不正确动作。
附图说明
图1为本发明提出的一种舞蹈与运动动作的智能对比矫正的方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
下面结合附图详细说明本发明的实施例。
如图1所示,一种舞蹈与运动动作的智能对比矫正的方法,包括以下步骤:
S1、采用人体姿态估计模型AlphaPose提取学生舞蹈视频与老师舞蹈视频中的人体关节点;
S2、将学生舞蹈视频与老师舞蹈视频进行动作对齐,并提取学生舞蹈视频与老师舞蹈视频的关键帧;
S3、根据提取的人体关节点、学生舞蹈视频的关键帧与老师舞蹈视频的关键帧,对动作对齐后的学生舞蹈视频与老师舞蹈视频进行动作相似对比;
S4、根据动作相似对比结果,找出动作差异区间,对学生动作进行可视化矫正,得到智能对比矫正结果。
所述步骤S2中将学生舞蹈视频与老师舞蹈视频进行动作对齐包括有音乐情况的对齐和无音乐情况的对齐。
人一般是伴随着音乐起舞的,所以为了找到视频中起舞阶段的起点,可以将舞蹈视频中所使用的伴乐纯净音频作为基准去匹配对齐视频中的音频,从而找到视频中伴乐的起点,因此也就是舞蹈的起点。同时因为视频录制时一般包含环境噪声,所以为了更好地保证对齐效果,需要对视频中的音频进行去噪处理。
所述有音乐情况的对齐具体为:
S2.1.1、提取学生舞蹈视频与老师舞蹈视频的音频,并提取音频中前N帧没有音乐的时间,并将前N帧时间中的噪音强度取平均值,以平均值构建估计噪声频谱;
S2.1.2、根据包含噪声的音频频谱Ps(w)以及估计噪声频谱Pn(w),获取频谱差值D(w)为:
D(w)=Ps(w)-αPn(w)
其中,α表示过减因子参数;
S2.1.3、将频谱差值D(w)中的负值进行修改,得到纯净音频Ps'(w)为:
其中,β表示频谱下限阈值参数,0<β<<1;
S2.1.4、设定伴乐音频和纯净音频采样率分别为Fs1和Fs2,并对伴乐音频进行采样,获取伴乐音频波动的时间点为s1
S2.1.5、获取伴乐音频和纯净音频之间的延迟值d,并根据延迟值d、采样率Fs1、采样率Fs2以及时间点s1,获取纯净音频波动的时间点s2为:
S2.1.6、根据纯净音频波动的时间点s2的获取方法,分别获取学生舞蹈视频与老师舞蹈视频中伴乐的起点,并提取起点后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
在本实施例中,根据伴乐纯净音频的时间起点与采样率,可以得到伴乐纯净音频采样上的对应起点。之后与延迟值d作和就可以计算出视频音频中伴乐开始的样本点位置,然后除以视频音频的采样频率就可以计算出视频音频开始波动的时间点s2
为了求得s2,需要估计伴乐纯净音频开始波动的样本点s1·Fs1,一种方法可以人工测量直接给出s1,乘以采样率就可以得到样本点。但由于人工测量本身耗费人力,并且往往只能精确到秒,无法进一步精确到更微小的单位,所以采用另一种方法。该方法基于这样一个前提:因为伴乐原曲是在录音棚条件下录制并经过后期处理的,所以其噪声基本为0。则可以设定一个阈值,找出第一个大于该阈值的样本的位置,该位置即为s1·Fs1。这种方法准度与阈值的设定有关,阈值过大就会导致偏差过大,经过大量试验对比,阈值设定为0.01效果较好,准确度较高。
为了求得s2,需要求出两段信号之间的延迟值d,这可以通过计算两段信号之间的互相关性得出。互相关运算计算的是两段信号之间的相关程度,将伴乐纯净音频作为X,去噪后的视频音频作为Y,将其中样本数据较短的信号补零实现长度的对等。假设该长度为N,则计算[1-N,N-1]范围内X、Y的偏移乘积,结果中最大值对应的下标即为延迟值d。从此计算过程中能够看出,X,Y不需要是彼此精确的延迟副本,这意味着对于不同采样率的音频也可以实现对齐。
同时由于视频录制时,基于录制设备的不同,其音频通道可能为单声道或双声道,而伴乐纯净音频一般为双声道。对此情况在计算延迟值d时,当两者都为双声道时,则分别使左声道与左声道对齐,右声道与右声道对齐;当视频音频为单声道时,则用视频的单声道分别与音频的双声道对齐。无论哪种情况都会产生可能的dl、dr,进而计算出s2时有两个结果,而视频时间点都是非负的,所以对非负的值做算术平均作为最终的舞蹈的时间起点。如果值全负,则对齐失败。大量试验分析后发现,在少部分失败的例子中大多数是视频本身录制质量较差,其中的伴乐本身不全,或者伴乐忽强忽弱,少数是伴乐纯净音频录制质量较差。
所述步骤S2.1.2中过减因子参数α具体为:
其中,SNR表示视频帧的信噪比。
当一段运动动作视频中不存在音乐时,显然需要思考另一种解决方案。我们用教师视频的起始动作去匹配学生用户视频中的起始动作,以此找出学生视频的起点。为此,我们需要先提取用户视频中的关键帧,然后将关键帧与教师的起始动作进行对比,找出最相似的一帧,将此帧对应的时间视作时间起点。
不等长无音乐的视频,以较短视频的最后一个关键动作去较长视频中匹配,匹配成功之后,删除较长视频的后续内容,从而实现对齐效果。
所述无音乐情况的对齐具体:
S2.2.1、设定计数器T=1,并设定动作相似距离阈值为K,并将两帧之间的Fréchet距离作为两视频帧之间的舞蹈动作相似距离;
S2.2.2、将视频第T帧设定为关键帧;
S2.2.3、判断视频第T帧与第T+1帧之间的舞蹈动作相似距离是否小阈值K,若是,则令T的计数值加一,并重复步骤S2.2.3,否则令T的计数值加一,并返回步骤S2.2.2;
S2.2.4、根据步骤S2.2.2-S2.2.3所述的方法,遍历视频的每一帧,提取视频的关键帧;
S2.2.5、根据步骤S2.2.1-S2.2.4所述的方法,分别提取学生舞蹈视频和老师舞蹈视频关键帧;
S2.2.6、获取学生舞蹈视频关键帧与老师舞蹈视频中起始动作帧的舞蹈动作相似距离,并将舞蹈动作相似距离最小的关键帧作为学生舞蹈视频的起始动作帧;
S2.2.7、分别提取学生舞蹈视频与老师舞蹈视频中起始动作帧后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
在本实施例中,提取关键帧的目的是为了去除重复的动作,通过去除掉相似的帧以减少计算量,从而提高效率,例如:一个动作保持了t秒钟,则这t秒内所有帧的动作都是相似的,则将这t秒视频的第1帧作为关键帧。
对于视频关键帧的提取,首先循环读取视频的每一帧,然后判断当前帧和后一帧的相似性,如果两帧的相似度在设定的相似性阈值之内,则继续比较下一帧视频,直到两帧的相似性大于设定的阈值为止,然后将其设为新的关键帧。
例如,设定的阈值相似度为t,首先将第1帧设为关键帧,第1帧和第2帧的相似度为d2,小于设定的阈值相似度t,则第2帧不是关键帧;第1帧再和第3帧进行比较,第1帧和第3帧的相似度为d3,小于设定的阈值相似度t,则第3帧不是关键帧;第1帧再和第4帧进行比较,第1帧和第4帧的相似度为d4,大于设定的阈值相似度t,则第4帧是关键帧;然后第4帧再和第5、6、7、8……帧进行比较。若第4帧和第5、6、7、8帧的相似度小于t,而第4帧和第9帧的相似度大于t,则将第9帧设为关键帧。第9帧之后的比较方法类似,直到所有关键帧提取完毕。所以我们认为第1、2、3帧是相似的,是同一个动作;第4、5、6、7、8帧是相似的,是同一个动作,所以取有代表性的第1帧和第4帧为关键帧,第1帧代表了,第1、2、3帧;第4帧代表了,第4、5、6、7、8帧。
对于每一个教师和学生的视频,都将视频转换为关键帧序列。
将提取出来的学生视频的关键帧与教师的起始动作进行对比,计算其相似度,然后找出其中相似度最高的一帧,将此帧作为学生视频的舞蹈起点。
首先计算两帧之间关节点的Fréchet距离,通过Fréchet距离进一步判断视频帧的相似性,Fréchet距离越小,说明这些关节点的形状越相似,然后设定一个阈值参数,如果视频帧的关节点数据的Fréchet距离大于或者等于设定的阈值参数则认为这两帧动作相似,如果视频帧的关节点数据的Fréchet距离小于设定的阈值参数则认为这两帧动作不相似,例如:设定阈值参数为t,而视频帧的关节点数据的Fréchet距离为d,若d大于t,则认为这两帧动作不相似;若d小于或等于t,则认为这两帧动作相似。
Fréchet距离的定义是:两条有方向的曲线,不能回溯,这两条曲线之间最短的最大距离。两条曲线分别用函数f:[a,b]→V和函数g:[a',b']→V来表示。Fréchet距离计算公式:
δF(f,g)=infmax d(f(α(t)),g(β(t)))
α(t),β(t)∈[0,1]
对此数学表达式的理解是,对于每一对可能的描述函数α(t)和β(t),我们总能找到整个运动过程中曲线最短的最长距离,所以这个最小的距离即为Fréchet距离。
可以采用递归的方法来计算离散型Fréchet距离。我们用coupling distance来求Fréchet距离的上界,也就是说两条连续曲线,找到其中的关键点作为endpoint,然后coupling distance是端点间四个点的距离。Ca(i,j)为i,j点以及这两点之前的Fréchet距离。d(ui,vj)为i,j两点之间的欧氏距离。其中:
ca(i,j)=max{min(c(i-1,j),c(i-1,j-1),c(i,j-1)),d(ui,vj)}
离散曲线可以从上面三种情况到达i,j点。所以我们需要找到这三个位置最小的那一个距离与i,j点的距离进行比较,然后选择其中最大的距离作为i,j点的Fréchet距离。
所述步骤S3具体为:
S3.1、分别提取学生和老师对应关键帧中的人体关节点;
在本实施例中,提取关键帧对应的人体关节点具体方法为:
(1)Load input video:读取输入视频。
(2)Load detection loader:加载yolo目标检测模型,检测图像中待检测的人体。
(3)Load pose model:选择Mscoco作为评估数据集,使用MSCOCO keypoints的训练集和验证集,来fine-tune SPPE。
(4)Pose Estimation:经过上述的初始化之后,然后对视频各帧中的人体进行姿态估计,识别出人体的骨架关节点,最后所有帧的识别结果将会被保存进json文件中。结果中包含提取出的视频各帧所有关节点的位置信息和置信度。人体姿态检测过程如下:
a)首先使用yolo做目标检测,得到人体所在的区域框。
b)然后将检测到的人体区域框作为输入,使用STN+SPPE模块完成人体姿态的自动检测。STN空间变换网络,对于不准确的输入,得到准确的人的框。输入候选区域,用于获取高质量的候选区域。SPPE得到估计的姿态。SDTN空间逆变换网络,将估计的姿态映射回原始的图像坐标。
c)再通过Parametrics Pose NMS进行调整改善,消除额外的估计到的姿态。
d)在训练的同时,通过Parallel SPPE防止局部最优并进而提升SSTN的效果。
不准确的输入图像经过SSTN进行姿态估计,并将估计的结果映射在原图上,以此来调整原本的框,使框变的精准。
S3.2、骨架大小归一化;
S3.3、将学生舞蹈视频关键帧和老师舞蹈视频关键帧逐帧进行骨架匹配;
S3.4、判断学生与老师骨架匹配是否正确,若是,则进入步骤S3.5,否则结束动作相似对比流程;
S3.5、分别提取学生与老师的手部关节点,并进行手部关节点匹配,完成动作相似对比。
所述步骤S3.2具体为:
S3.2.1、分别读取学生和老师对应关键帧中人体关节点的位置坐标,所述人体关节点包括颈部节点、头部中心节点、两个肩部节点以及两个腰部节点;
S3.2.2、根据人体关节点对应的位置坐标,分别获取学生对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为Xs1、Xs2、Xs3、Xs4和Xs5
S3.2.3、根据人体关节点对应的位置坐标,分别获取老师对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为Xt1、Xt2、Xt3、Xt4和Xs5
S3.2.4、使用距离Xsi除以距离Xti,得到距离比例Ri,i=1,2,3,4,5;
S3.2.5、根据距离比例Ri,求取比例平均值为R;
S3.2.6、根据比例平均值为R,将学生骨架缩放至老师骨架大小,完成骨架归一化。
所述步骤S3.3中骨架匹配的具体方法为:
S3.3.1、根据学生和老师颈部节点的位置坐标,获取水平距离dis_x和垂直距离dis_y;
S3.3.2、通过比例平均值为R和老师颈部节点的位置坐标获取老师关节点的位置坐标调节值为:
V_Xi=(R-1)*(Xi-tch_x)
V_Yi=(R-1)*(Yi-tch_y)
其中,V_Xi表示第i个关节点的x位置坐标调节值,V_Yi表示第i个关节点的y位置坐标调节值,Xi表示第i个关节点调整前的x坐标,Yi表示第i个关节点调整前的y坐标,tch_x表示老师颈部节点的x坐标,tch_y表示老师颈部节点的y坐标;
S3.3.3、将老师每个关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成骨架匹配。
所述步骤S3.5具体为:
S3.5.1、分别读取学生和老师对应关键帧中手部关节点的位置坐标;
S3.5.2、根据比例平均值为R,将学生手部关节点之间的距离进行缩放,使之与老师手部关节点之间的距离相同;
S3.5.3、将老师每个手部关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的手部关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成手部关节点匹配。
所述步骤S4中对学生动作进行可视化矫正包括骨架迁移可视化矫正和视频拼接可视化矫正;
所述骨架迁移可视化矫正具体为:根据动作相似对比结果,以学生不匹配的关键帧为中心,向前后扩展一定时长的视频帧,将老师和学生的颈部关节点对齐,以颈部关节点为基准,将老师骨架迁移至学生舞蹈视频中,得到学生舞姿与老师骨架的对比结果,根据对比结果对学生动作进行矫正;
所述视频拼接可视化矫正具体为:通过视频对齐后,将老师和学生在同一律动下的视频帧一一拼接,得到的新帧集合拼接成一个视频,根据找出的动作差异区间,对照老师的标准动作进行矫正。
将输入的学生视频提取为若干关键帧,并对学生舞蹈视频的关键帧进行姿态识别。为了将老师的标准舞蹈动作的姿态骨架与学生的舞蹈动作进行对比,首先将老师的骨架进行归一化确定需要进行放缩的比例,然后以学生姿态识别结果的中心关节点为基准把老师骨架的中心节点定位匹配到学生舞蹈视频的关键帧中,最后将老师舞蹈姿态的所有骨架关节点平移至对应的学生舞蹈视频的关键帧中。
至此,老师舞蹈姿态的骨架以完整迁移至学生的舞蹈视频关键帧中,可在学生的舞蹈视频关键帧中看到学生的舞蹈动作与老师的舞蹈动作之间的存在差异。通过对学生和老师舞蹈姿态关节位置信息的计算,系统还可指出学生舞蹈动作与老师舞蹈动作差异较大的关键帧,使得学生可以较直观的看出动作哪里不标准,从而进行改正。
在本实施例中,通过对集中测试12组视频(每组视频包括舞蹈老师的标准视频和优秀、良好、较差三组学生的4个视频)进行测试,每个视频的时长在30秒至1分钟左右不等,总时长约在半小时;本发明检测的准确度可达到90%的正确率。
拼接视频对比矫正:通过S2视频对齐后,得到老师和学生在同一律动下的视频帧。将两者的视频帧一一拼接,得到新的一个对比视频,具体操作如下:
创建一张新的空白图片,将老师的视频帧在左边写入,学生的视频帧在右边写入,该图片的宽度即为老师和学生视频帧的宽度之和,得到的新的帧的集合拼接成一个视频。最后在拼接的视频中加入经过对齐截取后的标准伴奏音乐。
在该视频播放过程中,老师的舞蹈视频位于左侧,学生自己的舞蹈视频位于右侧,便于学生对照老师标准的动作、节奏观察学习,直接清晰地看到自身舞蹈的不足之处。

Claims (8)

1.一种舞蹈与运动动作的智能对比矫正的方法,其特征在于,包括以下步骤:
S1、采用人体姿态估计模型AlphaPose提取学生舞蹈视频与老师舞蹈视频中的人体关节点;
S2、将学生舞蹈视频与老师舞蹈视频进行动作对齐,并提取学生舞蹈视频与老师舞蹈视频的关键帧;
S3、根据提取的人体关节点、学生舞蹈视频的关键帧与老师舞蹈视频的关键帧,对动作对齐后的学生舞蹈视频与老师舞蹈视频进行动作相似对比;
S4、根据动作相似对比结果,找出动作差异区间,对学生动作进行可视化矫正,得到智能对比矫正结果;
所述步骤S2中将学生舞蹈视频与老师舞蹈视频进行动作对齐包括有音乐情况的对齐和无音乐情况的对齐;
所述有音乐情况的对齐具体为:
S2.1.1、提取学生舞蹈视频与老师舞蹈视频的音频,并提取音频中前N帧没有音乐的时间,并将前N帧时间中的噪音强度取平均值,以平均值构建估计噪声频谱;
S2.1.2、根据包含噪声的音频频谱以及估计噪声频谱/>,获取频谱差值为:
其中,表示过减因子参数;
S2.1.3、将频谱差值中的负值进行修改,得到纯净音频/>为:
其中,表示频谱下限阈值参数;
S2.1.4、设定伴乐音频和纯净音频采样率分别为和/>,并对伴乐音频进行采样,获取伴乐音频波动的时间点为/>
S2.1.5、获取伴乐音频和纯净音频之间的延迟值d,并根据延迟值d、采样率、采样率以及时间点/>,获取纯净音频波动的时间点/>为:
S2.1.6、根据纯净音频波动的时间点的获取方法,分别获取学生舞蹈视频与老师舞蹈视频中伴乐的起点,并提取起点后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
2.根据权利要求1所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S2.1.2中过减因子参数具体为:
其中,表示视频帧的信噪比。
3.根据权利要求2所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述无音乐情况的对齐具体:
S2.2.1、设定计数器T=1,并设定动作相似距离阈值为K,并将两帧之间的Fréchet距离作为两视频帧之间的舞蹈动作相似距离;
S2.2.2、将视频第T帧设定为关键帧;
S2.2.3、判断视频第T帧与第T+1帧之间的舞蹈动作相似距离是否小阈值K,若是,则令T的计数值加一,并重复步骤S2.2.3,否则令T的计数值加一,并返回步骤S2.2.2;
S2.2.4、根据步骤S2.2.2-S2.2.3所述的方法,遍历视频的每一帧,提取视频的关键帧;
S2.2.5、根据步骤S2.2.1-S2.2.4所述的方法,分别提取学生舞蹈视频和老师舞蹈视频关键帧;
S2.2.6、获取学生舞蹈视频关键帧与老师舞蹈视频中起始动作帧的舞蹈动作相似距离,并将舞蹈动作相似距离最小的关键帧作为学生舞蹈视频的起始动作帧;
S2.2.7、分别提取学生舞蹈视频与老师舞蹈视频中起始动作帧后的视频,完成学生舞蹈视频与老师舞蹈视频的动作对齐。
4.根据权利要求3所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S3具体为:
S3.1、分别提取学生和老师对应关键帧中的人体关节点;
S3.2、骨架大小归一化;
S3.3、将学生舞蹈视频关键帧和老师舞蹈视频关键帧逐帧进行骨架匹配;
S3.4、判断学生与老师骨架匹配是否正确,若是,则进入步骤S3.5,否则结束动作相似对比流程;
S3.5、分别提取学生与老师的手部关节点,并进行手部关节点匹配,完成动作相似对比。
5.根据权利要求4所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S3.2具体为:
S3.2.1、分别读取学生和老师对应关键帧中人体关节点的位置坐标,所述人体关节点包括颈部节点、头部中心节点、两个肩部节点以及两个腰部节点;
S3.2.2、根据人体关节点对应的位置坐标,分别获取学生对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为、/>、/>、/>和/>
S3.2.3、根据人体关节点对应的位置坐标,分别获取老师对应人体关节点中颈部节点到头部中心节点、第一肩部节点、第二肩部节点、第一腰部节点以及第二腰部节点的距离为、/>、/>、/>和/>
S3.2.4、使用距离除以距离/>,得到距离比例/>i=1,2,3,4,5;
S3.2.5、根据距离比例,求取比例平均值为R
S3.2.6、根据比例平均值为R,将学生骨架缩放至老师骨架大小,完成骨架归一化。
6.根据权利要求5所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S3.3中骨架匹配的具体方法为:
S3.3.1、根据学生和老师颈部节点的位置坐标,获取水平距离dis_x和垂直距离dis_y;
S3.3.2、通过比例平均值为R和老师颈部节点的位置坐标获取老师关节点的位置坐标调节值为:
其中,表示第i个关节点的x位置坐标调节值, />表示第i个关节点的y位置坐标调节值,/>表示第i个关节点调整前的x坐标,/>表示第i个关节点调整前的y坐标,表示老师颈部节点的x坐标,/>表示老师颈部节点的y坐标;
S3.3.3、将老师每个关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成骨架匹配。
7.根据权利要求6所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S3.5具体为:
S3.5.1、分别读取学生和老师对应关键帧中手部关节点的位置坐标;
S3.5.2、根据比例平均值为R,将学生手部关节点之间的距离进行缩放,使之与老师手部关节点之间的距离相同;
S3.5.3、将老师每个手部关节点坐标加上位置坐标调节值进行调节,以颈部节点为粘合节点,将老师所有调节后的手部关节点坐标分别减去水平距离dis_x和垂直距离dis_y,完成手部关节点匹配。
8.根据权利要求7所述的舞蹈与运动动作的智能对比矫正的方法,其特征在于,所述步骤S4中对学生动作进行可视化矫正包括骨架迁移可视化矫正和视频拼接可视化矫正;
所述骨架迁移可视化矫正具体为:根据动作相似对比结果,以学生不匹配的关键帧为中心,向前后扩展一定时长的视频帧,将老师和学生的颈部关节点对齐,以颈部关节点为基准,将老师骨架迁移至学生舞蹈视频中,得到学生舞姿与老师骨架的对比结果,根据对比结果对学生动作进行矫正;
所述视频拼接可视化矫正具体为:通过视频对齐后,将老师和学生在同一律动下的视频帧一一拼接,得到的新帧集合拼接成一个视频,根据找出的动作差异区间,对照老师的标准动作进行矫正。
CN202110448369.9A 2021-04-25 2021-04-25 一种舞蹈与运动动作的智能对比矫正的方法 Active CN113052138B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110448369.9A CN113052138B (zh) 2021-04-25 2021-04-25 一种舞蹈与运动动作的智能对比矫正的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110448369.9A CN113052138B (zh) 2021-04-25 2021-04-25 一种舞蹈与运动动作的智能对比矫正的方法

Publications (2)

Publication Number Publication Date
CN113052138A CN113052138A (zh) 2021-06-29
CN113052138B true CN113052138B (zh) 2024-03-15

Family

ID=76520452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110448369.9A Active CN113052138B (zh) 2021-04-25 2021-04-25 一种舞蹈与运动动作的智能对比矫正的方法

Country Status (1)

Country Link
CN (1) CN113052138B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114241595A (zh) * 2021-11-03 2022-03-25 橙狮体育(北京)有限公司 数据处理方法、装置、电子设备及计算机存储介质
CN114648810B (zh) * 2022-03-11 2022-10-14 北京中科睿医信息科技有限公司 一种交互式步态识别方法、装置及电子设备
CN114742889A (zh) * 2022-03-16 2022-07-12 北京工业大学 一种基于九轴姿态传感器和机器视觉的人体舞蹈动作检测与纠正方法
CN114582029B (zh) * 2022-05-06 2022-08-02 山东大学 一种非专业舞蹈运动序列增强方法及系统

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2461067C1 (ru) * 2011-05-04 2012-09-10 Государственное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ГОУ ВПО ТГТУ Способ обнаружения движущихся объектов
JP2014023745A (ja) * 2012-07-27 2014-02-06 Seek Japan Corp ダンス教習装置
CN104992712A (zh) * 2015-07-06 2015-10-21 成都云创新科技有限公司 能识别音乐自动成谱的方法
CN105989853A (zh) * 2015-02-28 2016-10-05 科大讯飞股份有限公司 一种音频质量评测方法及系统
CN108734104A (zh) * 2018-04-20 2018-11-02 杭州易舞科技有限公司 基于深度学习图像识别的健身动作纠错方法及系统
CN108777081A (zh) * 2018-05-31 2018-11-09 华中师范大学 一种虚拟舞蹈教学方法及系统
CN109508656A (zh) * 2018-10-29 2019-03-22 重庆中科云丛科技有限公司 一种舞蹈评级自动判别方法、系统及计算机可读存储介质
CN109522850A (zh) * 2018-11-22 2019-03-26 中山大学 一种基于小样本学习的动作相似度评估方法
CN109545191A (zh) * 2018-11-15 2019-03-29 电子科技大学 一种歌曲中人声起始位置的实时检测方法
CN109589563A (zh) * 2018-12-29 2019-04-09 南京华捷艾米软件科技有限公司 一种基于3d体感摄像头的舞蹈姿态教辅方法及系统
CN109833608A (zh) * 2018-12-29 2019-06-04 南京华捷艾米软件科技有限公司 一种基于3d体感摄像头的舞蹈动作教辅方法及系统
CN110782908A (zh) * 2019-11-05 2020-02-11 广州欢聊网络科技有限公司 一种音频信号处理的方法及装置
CN110841287A (zh) * 2019-11-22 2020-02-28 腾讯科技(深圳)有限公司 视频处理方法、装置、计算机可读存储介质和计算机设备
CN111144217A (zh) * 2019-11-28 2020-05-12 重庆邮电大学 一种基于人体三维关节点检测的动作评价方法
CN111554256A (zh) * 2020-04-21 2020-08-18 华南理工大学 一种基于强弱标准的钢琴视奏能力评价系统
CN111563487A (zh) * 2020-07-14 2020-08-21 平安国际智慧城市科技股份有限公司 基于姿态识别模型的舞蹈评分方法及相关设备
CN112309432A (zh) * 2020-10-27 2021-02-02 暨南大学 一种基于数据驱动的音符起始点检测方法
CN112560665A (zh) * 2020-12-13 2021-03-26 同济大学 基于深度迁移学习实现人体姿态检测的专业舞蹈评价方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140118522A1 (en) * 2012-11-01 2014-05-01 Josh Heath Zuniga Dance learning system using a computer
CN112530450A (zh) * 2019-09-17 2021-03-19 杜比实验室特许公司 频域中的样本精度延迟识别

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2461067C1 (ru) * 2011-05-04 2012-09-10 Государственное образовательное учреждение высшего профессионального образования "Тамбовский государственный технический университет" ГОУ ВПО ТГТУ Способ обнаружения движущихся объектов
JP2014023745A (ja) * 2012-07-27 2014-02-06 Seek Japan Corp ダンス教習装置
CN105989853A (zh) * 2015-02-28 2016-10-05 科大讯飞股份有限公司 一种音频质量评测方法及系统
CN104992712A (zh) * 2015-07-06 2015-10-21 成都云创新科技有限公司 能识别音乐自动成谱的方法
CN108734104A (zh) * 2018-04-20 2018-11-02 杭州易舞科技有限公司 基于深度学习图像识别的健身动作纠错方法及系统
CN108777081A (zh) * 2018-05-31 2018-11-09 华中师范大学 一种虚拟舞蹈教学方法及系统
CN109508656A (zh) * 2018-10-29 2019-03-22 重庆中科云丛科技有限公司 一种舞蹈评级自动判别方法、系统及计算机可读存储介质
CN109545191A (zh) * 2018-11-15 2019-03-29 电子科技大学 一种歌曲中人声起始位置的实时检测方法
CN109522850A (zh) * 2018-11-22 2019-03-26 中山大学 一种基于小样本学习的动作相似度评估方法
CN109589563A (zh) * 2018-12-29 2019-04-09 南京华捷艾米软件科技有限公司 一种基于3d体感摄像头的舞蹈姿态教辅方法及系统
CN109833608A (zh) * 2018-12-29 2019-06-04 南京华捷艾米软件科技有限公司 一种基于3d体感摄像头的舞蹈动作教辅方法及系统
CN110782908A (zh) * 2019-11-05 2020-02-11 广州欢聊网络科技有限公司 一种音频信号处理的方法及装置
CN110841287A (zh) * 2019-11-22 2020-02-28 腾讯科技(深圳)有限公司 视频处理方法、装置、计算机可读存储介质和计算机设备
CN111144217A (zh) * 2019-11-28 2020-05-12 重庆邮电大学 一种基于人体三维关节点检测的动作评价方法
CN111554256A (zh) * 2020-04-21 2020-08-18 华南理工大学 一种基于强弱标准的钢琴视奏能力评价系统
CN111563487A (zh) * 2020-07-14 2020-08-21 平安国际智慧城市科技股份有限公司 基于姿态识别模型的舞蹈评分方法及相关设备
CN112309432A (zh) * 2020-10-27 2021-02-02 暨南大学 一种基于数据驱动的音符起始点检测方法
CN112560665A (zh) * 2020-12-13 2021-03-26 同济大学 基于深度迁移学习实现人体姿态检测的专业舞蹈评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Lighten Up! Postural Instructions Affect Static and Dynamic Balance in Healthy Older Adults;Cohen, RG等;《NNOVATION IN AGING》;20201231(第2期);1-5 *
运动动作图像正确姿态轮廓提取方法对比研究;邹文君;《现代电子技术》;20181231(第15期);65-68+72 *

Also Published As

Publication number Publication date
CN113052138A (zh) 2021-06-29

Similar Documents

Publication Publication Date Title
CN113052138B (zh) 一种舞蹈与运动动作的智能对比矫正的方法
CN113762133B (zh) 基于人体姿态识别的自重健身辅助教练系统、方法、终端
CN108764120B (zh) 一种人体规范动作评价方法
Folgado et al. Time alignment measurement for time series
CN105426827B (zh) 活体验证方法、装置和系统
Papandreou et al. Adaptive multimodal fusion by uncertainty compensation with application to audiovisual speech recognition
CN109815795A (zh) 基于人脸监测的课堂学生状态分析方法及装置
CN110448870B (zh) 一种人体姿态训练方法
CN109508661B (zh) 一种基于物体检测和姿态估计的举手者检测方法
Chiang et al. Kinect-based in-home exercise system for lymphatic health and lymphedema intervention
Auephanwiriyakul et al. Thai sign language translation using scale invariant feature transform and hidden markov models
CN112541870A (zh) 一种视频处理的方法、装置、可读存储介质和电子设备
CN113158914B (zh) 一种舞蹈动作姿态、节奏和表情的智能评估方法
Wei et al. Performance monitoring and evaluation in dance teaching with mobile sensing technology
CN112101315A (zh) 一种基于深度学习的运动评判指导方法和系统
CN112200138B (zh) 基于计算机视觉的课堂学情分析方法
CN117542121B (zh) 基于计算机视觉的智能化训练考核系统及方法
CN116844084A (zh) 一种融合区块链的体育运动动作分析与纠正方法及系统
CN116343335A (zh) 一种基于动作识别的运动姿态矫正方法
Guo et al. PhyCoVIS: A visual analytic tool of physical coordination for cheer and dance training
CN114333063A (zh) 一种基于人体姿态估计的武术动作纠正方法及装置
CN117137435B (zh) 一种基于多模态信息融合的康复动作识别方法和系统
CN103778439B (zh) 基于动态时空信息挖掘的人体轮廓重构方法
KR101621304B1 (ko) 마우스맵을 이용한 능동형태모델 기반 입술 형태 추정 방법 및 시스템
CN116704603A (zh) 一种基于肢体关键点分析的动作评估纠正方法及系统

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